GreenPolyProducts/cantera
5
0
Fork 0
mirror of https://github.com/Cantera/cantera.git synced 2025-02-25 18:55:29 -06:00
Code Issues Packages Projects Releases Wiki Activity

[Reactor] Fix getting sensitivities after calls to advance

Browse Source 
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.
This commit is contained in:
Ray Speth
2023-08-04 13:51:07 -04:00
committed by Ingmar Schoegl
parent 222b0ebcf6
commit e0a2f112f7
2 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
src/numerics/CVodesIntegrator.cpp
View File

@@ -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;
}

2
src/numerics/IdasIntegrator.cpp
View File

@@ -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;
}