From c1c7965a6b81fe93ea52628e21d7fcca3e1c266d Mon Sep 17 00:00:00 2001 From: ssun30 Date: Thu, 7 Nov 2024 22:22:59 -0500 Subject: [PATCH] loosen solver tolerances for some tests --- src/TestReactors.jl | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/TestReactors.jl b/src/TestReactors.jl index d05e97f7..4b05c6d3 100644 --- a/src/TestReactors.jl +++ b/src/TestReactors.jl @@ -76,7 +76,7 @@ using DataFrames 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"] @@ -90,7 +90,7 @@ using DataFrames 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) @@ -136,7 +136,7 @@ using DataFrames domains = (domainliq, domaingas) interfaces = [vl, inletgas, outletgas] react, y0, p = Reactor(domains, (y0liq, y0gas), (0.0, tf), interfaces, (pliq, pgas, pinter)) - sol = solve(react.ode, react.recommendedsolver, abstol=1e-18, reltol=1e-6) + sol = solve(react.ode, react.recommendedsolver, abstol=1e-16, reltol=1e-6) name = "oxygen" ind = findfirst(x -> x == name, liqspcnames) @@ -581,7 +581,7 @@ using DataFrames react, y0, p = Reactor((domaingas, domaincat), (y0gas, y0cat), (0.0, 0.1), (inter,), (pgas, pcat, pinter)) - sol = solve(react.ode, CVODE_BDF(), abstol=1e-20, reltol=1e-6) + sol = solve(react.ode, CVODE_BDF(), abstol=1e-16, reltol=1e-6) ssys = SystemSimulation(sol, (domaingas, domaincat,), (inter,), p)