Fix temperature delta for IdealGasConstPressureMoleReactor

Cantera · Aug 30, 2022 · 357902c · 357902c
1 parent 3fc279e
commit 357902c
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/zeroD/IdealGasConstPressureMoleReactor.cpp b/src/zeroD/IdealGasConstPressureMoleReactor.cpp
@@ -167,7 +167,7 @@ Eigen::SparseMatrix<double> IdealGasConstPressureMoleReactor::jacobian()
     if (m_energy) {
         // getting perturbed state for finite difference
         double deltaTemp = m_thermo->temperature()
-            * std::numeric_limits<double>::epsilon();
+            * std::sqrt(std::numeric_limits<double>::epsilon());
         // finite difference temperature derivatives
         vector_fp ydotNext(m_nv);
         vector_fp yCurrent(m_nv);