loosen solve tolerances

ReactionMechanismGenerator · Sep 23, 2024 · de4d6ba · de4d6ba
1 parent 7de5980
commit de4d6ba
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/TestReactors.jl b/src/TestReactors.jl
@@ -77,7 +77,7 @@ using SciMLSensitivity
         interfaces = [kLAkHCondensationEvaporationWithReservoir(domain, conds)]
         react = Reactor(domain, y0, (0.0, 140000.01), interfaces; p=p)
 
-        sol1 = solve(react.ode, react.recommendedsolver, abstol=1e-18, reltol=1e-6)
+        sol1 = solve(react.ode, react.recommendedsolver, abstol=1e-16, reltol=1e-6)
 
         phaseDict = readinput("../src/testing/TdependentkLAkH.rms")
         spcs = phaseDict["phase"]["Species"]
@@ -91,7 +91,7 @@ using SciMLSensitivity
         interfaces = [kLAkHCondensationEvaporationWithReservoir(domain, conds)]
         react = Reactor(domain, y0, (0.0, 140000.01), interfaces; p=p) #Create the reactor object
 
-        sol2 = solve(react.ode, react.recommendedsolver, abstol=1e-18, reltol=1e-6)
+        sol2 = solve(react.ode, react.recommendedsolver, abstol=1e-16, reltol=1e-6)
 
         spcnames = getfield.(liq.species, :name)
         octaneind = findfirst(isequal("octane"), spcnames)