Further change to make sure it works on linux.

- sheesh There are problems with the tolerances and internal convergence
     requirements within the code.
This commit is contained in:
Harry Moffat 2009-07-14 00:17:52 +00:00
parent 3088a0b365
commit d1a4a8e7c1

View file

@ -119,7 +119,7 @@ for n in range(NReactors):
sim = ReactorNet([upstream, r, downstream])
# set relative and absolute tolerances on the simulation
sim.setTolerances(rtol = 1.0e-3, atol = 1.0e-11)
sim.setTolerances(rtol = 1.0e-4, atol = 1.0e-11)
time = 0
while 1 > 0:
@ -142,7 +142,7 @@ for n in range(NReactors):
for ks in range(nsurf):
ratio = sdot[ks]/(cdot[ks] + ddot[ks])
if ratio < 0.0: ratio = -ratio
if ratio > 1.0e-11 or time < 10*dt:
if ratio > 1.0e-9 or time < 10*dt:
alldone = 0
if alldone: break