Hi,
I'm wondering how to do with the "dt = " setting in my sweeps when I use exponential time steps.

When I use fixed time step increases, the scope is like:

dt = 200
time = a time variable keeping track of current time
while time() < endTime:
   time.setValue(time() + dt)
   while someResidual < desiredResidual:
       someResidual1 = eq1.sweep(var = ..., dt = dt)
       someResidual2 = eq2.sweep(var = ..., dt = dt)
      ....

What happens if I change time increase according to manual example:

dexp = -5.

while time() < endTime:
   time.setValue(numpy.exp(dexp))
   dexp += .5
   while ....:
      ????


Thanks,
Elizabeth


Reply via email to