From 02ff1a6071c43778918f69074859705b7f88cacf Mon Sep 17 00:00:00 2001 From: Harry Moffat Date: Tue, 4 Dec 2007 17:37:41 +0000 Subject: [PATCH] Decreased the tolerance in the setState_UV() routine. This tolerance is needed for the accurate calculation of Jacobians. --- Cantera/src/zeroD/Reactor.cpp | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) diff --git a/Cantera/src/zeroD/Reactor.cpp b/Cantera/src/zeroD/Reactor.cpp index c1e56e26c..d8eadc4da 100644 --- a/Cantera/src/zeroD/Reactor.cpp +++ b/Cantera/src/zeroD/Reactor.cpp @@ -148,8 +148,12 @@ namespace CanteraZeroD { mix.setTemperature(temp); if (m_energy) { - m_thermo->setState_UV(u/mass,m_vol/mass,1.0e-4); - temp = mix.temperature(); //mix.setTemperature(temp); + // Decreased the tolerance on delta_T to 1.0E-7 so that T is + // accurate to 9 sig digits, because this is + // used in the numerical jacobian routines where relative values + // of 1.0E-7 are used in the deltas. + m_thermo->setState_UV(u/mass,m_vol/mass, 1.0e-7); + temp = mix.temperature(); //mix.setTemperature(temp); } //m_state[0] = temp;