[Reactor] Fix getting sensitivities after calls to advance

Previously, values would be returned for the sensitivity at the last internal timestep of the integrator, rather than being interpolated to the user-specified integration time.
Cantera · Aug 4, 2023 · 652ed24 · 652ed24
1 parent 59625eb
commit 652ed24
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/numerics/CVodesIntegrator.cpp b/src/numerics/CVodesIntegrator.cpp
@@ -637,7 +637,7 @@ double CVodesIntegrator::sensitivity(size_t k, size_t p)
         return 0.0;
     }
     if (!m_sens_ok && m_np) {
-        int flag = CVodeGetSens(m_cvode_mem, &m_time, m_yS);
+        int flag = CVodeGetSensDky(m_cvode_mem, m_time, 0, m_yS);
         checkError(flag, "sensitivity", "CVodeGetSens");
         m_sens_ok = true;
     }

diff --git a/src/numerics/IdasIntegrator.cpp b/src/numerics/IdasIntegrator.cpp
@@ -512,7 +512,7 @@ double IdasIntegrator::sensitivity(size_t k, size_t p)
         return 0.0;
     }
     if (!m_sens_ok && m_np) {
-        int flag = IDAGetSens(m_ida_mem, &m_time, m_yS);
+        int flag = IDAGetSensDky(m_ida_mem, m_time, 0, m_yS);
         checkError(flag, "sensitivity", "IDAGetSens");
         m_sens_ok = true;
     }