All tests passing locally with both GSL and BOOST.

Author: Dilawar Singh

builtins/Function.cpp

@@ -578,9 +578,9 @@ double Function::getDerivative() const
 
 void Function::setNumVar(const unsigned int num)
 {
+    //cerr << "Deprecated: numVar has no effect. MOOSE can infer number of variables "
+    //     " from the expression. " << endl;
     numVar_ = num;
-    cerr << "Deprecated: numVar has no effect. MOOSE can infer number of variables "
-         " from the expression. " << endl;
 }
 
 unsigned int Function::getNumVar() const

tests/py_rdesigneur/test_72_CICR.py

@@ -22,24 +22,26 @@
 
 sdir_ = os.path.dirname(os.path.realpath(__file__))
 
-E = (np.array([1.090e-07, 7.281e-13, 2.754e-08, 4.094e-01, 5.138e-04, 5.044e-04,
-        5.185e-04, 5.203e-04, 5.203e-04, 5.203e-04, 5.203e-04, 5.203e-04,
-        5.203e-04, 5.203e-04, 5.138e-04, 5.044e-04, 5.185e-04, 5.203e-04,
-        5.203e-04, 5.203e-04, 5.203e-04, 5.203e-04, 5.203e-04, 5.203e-04,
-        4.033e-01, 4.046e-01, 4.038e-01, 4.039e-01, 4.039e-01, 4.039e-01,
-        4.039e-01, 4.039e-01, 4.039e-01, 4.039e-01, 1.081e-06, 1.037e-06,
-        1.046e-06, 1.050e-06, 1.050e-06, 1.050e-06, 1.050e-06, 1.050e-06,
-        1.050e-06, 1.050e-06]),
-    np.array([2.648e-06, 3.539e-12, 1.063e-07, 2.596e-05, 1.508e-03, 1.444e-03,
-        1.465e-03, 1.464e-03, 1.464e-03, 1.464e-03, 1.464e-03, 1.464e-03,
-        1.464e-03, 1.464e-03, 1.508e-03, 1.444e-03, 1.465e-03, 1.464e-03,
-        1.464e-03, 1.464e-03, 1.464e-03, 1.464e-03, 1.464e-03, 1.464e-03,
-        1.268e-02, 1.160e-02, 1.193e-02, 1.201e-02, 1.201e-02, 1.201e-02,
-        1.201e-02, 1.201e-02, 1.201e-02, 1.201e-02, 2.116e-06, 2.063e-06,
-        2.081e-06, 2.086e-06, 2.086e-06, 2.086e-06, 2.086e-06, 2.086e-06,
-        2.086e-06, 2.086e-06]))
-
-
+E = (np.array([
+    1.090e-07, 7.281e-13, 2.754e-08, 4.094e-01, 1.733e-05, 1.733e-05,
+    1.733e-05, 1.733e-05, 1.733e-05, 1.733e-05, 1.733e-05, 1.733e-05,
+    1.733e-05, 1.733e-05, 1.733e-05, 1.733e-05, 1.733e-05, 1.733e-05,
+    1.733e-05, 1.733e-05, 1.733e-05, 1.733e-05, 1.733e-05, 1.733e-05,
+    4.089e-01, 4.089e-01, 4.089e-01, 4.089e-01, 4.089e-01, 4.089e-01,
+    4.089e-01, 4.089e-01, 4.089e-01, 4.089e-01, 0.000e+00, 0.000e+00,
+    0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,
+    0.000e+00, 0.000e+00
+]),
+     np.array([
+         2.648e-06, 3.539e-12, 1.063e-07, 2.596e-05, 2.022e-06, 2.024e-06,
+         2.024e-06, 2.024e-06, 2.024e-06, 2.024e-06, 2.024e-06, 2.024e-06,
+         2.024e-06, 2.024e-06, 2.022e-06, 2.024e-06, 2.024e-06, 2.024e-06,
+         2.024e-06, 2.024e-06, 2.024e-06, 2.024e-06, 2.024e-06, 2.024e-06,
+         2.403e-05, 2.292e-05, 2.291e-05, 2.291e-05, 2.291e-05, 2.291e-05,
+         2.291e-05, 2.291e-05, 2.291e-05, 2.291e-05, 0.000e+00, 0.000e+00,
+         0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00,
+         0.000e+00, 0.000e+00
+     ]))
 
 
 def test():
@@ -76,10 +78,12 @@ def test():
     moose.start(20)
     data = [t.vector for t in moose.wildcardFind('/##[TYPE=Table2]')]
     m, s = np.mean(data, axis=1), np.std(data, axis=1)
+    # print(np.array_repr(m))
+    # print(np.array_repr(s))
     # In multithreaded mode, the numers are not exactly the same as in
     # expected.
-    assert np.allclose(m, E[0], rtol=1e-3), (m,  E[0])
-    assert np.allclose(s, E[1], rtol=1e-3), (s,  E[1])
+    assert np.allclose(m, E[0], rtol=1e-3), (m, E[0])
+    assert np.allclose(s, E[1], rtol=1e-3), (s, E[1])
     print('done')