From 39a6819bf38d2c55d8233e4a749ad688ca01f306 Mon Sep 17 00:00:00 2001 From: Ray Speth Date: Fri, 11 Dec 2015 11:36:00 -0500 Subject: [PATCH] [1D] Only reset Jacobian age after successful *steady-state* solve --- src/oneD/MultiNewton.cpp | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/src/oneD/MultiNewton.cpp b/src/oneD/MultiNewton.cpp index 4f2101c53..00efb9491 100644 --- 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