loosen solver tolerances for some tests

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)