[1D] Only reset Jacobian age after successful *steady-state* solve

2015-12-11 11:36:00 -05:00 · 2015-12-11 11:36:00 -05:00 · 39a6819bf3
commit 39a6819bf3
parent f0a1e25cd9
1 changed files with 3 additions and 1 deletions
--- a/src/oneD/MultiNewton.cpp
+++ b/src/oneD/MultiNewton.cpp
@ -349,7 +349,9 @@ int MultiNewton::solve(doublereal* x0, doublereal* x1,
            copy(x1, x1 + m_n, m_x.begin());
        } else if (m == 1) {
            // convergence
-            jac.setAge(0); // for efficient sensitivity analysis
+            if (rdt == 0) {
+                jac.setAge(0); // for efficient sensitivity analysis
+            }
            break;
        } else if (m < 0) {
            // If dampStep fails, first try a new Jacobian if an old one was