Hi there, While using the AdsorbingSurfactantEquation class, defined in the electroChem example, located in:
fipy/examples/levelSet/electroChem/adsorbingSurfactantEquation.py
I noticed two lines that seem erroneous, both related to the timestep used
for the solution. On line 331, the following code:
self.eq.solve(var,
boundaryConditions=boundaryConditions,
solver = solver,
dt=1.)
should, I think, be:
self.eq.solve(var,
boundaryConditions=boundaryConditions,
solver = solver,
dt=self.dt)
and likewise, line 358, the following code:
return self.eq.sweep(var, solver=solver,
boundaryConditions=boundaryConditions, underRelaxation=underRelaxation,
residualFn=residualFn, dt=1.)
should be:
return self.eq.sweep(var, solver=solver,
boundaryConditions=boundaryConditions, underRelaxation=underRelaxation,
residualFn=residualFn, dt=self.dt)
so that the timestep information set by the user is not ignored. Attached
is a patch (created with git diff --no-prefix) that fixes this. If it
passes the muster, would someone mind checking the change into the git
repository?
Thanks,
Charles
adsorbingSurfactantEquation_dt_fix.patch
Description: Binary data
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
