usnistgov / fipy

FiPy is a Finite Volume PDE solver written in Python
http://pages.nist.gov/fipy/en/latest
Other
489 stars 148 forks source link

examples.diffusion.mesh1D No-flux - steady-state doesn't always give zero #1010

Closed agoncalvescoel closed 5 months ago

agoncalvescoel commented 5 months ago

When I run mesh1D.py , the other parts of the script run fine, but the part on the figure below gives me different results image

The blue line does not appear at all image

And when it uns the validation of numerix, gives me False image

guyer commented 5 months ago

So, the code is working properly, but the explanation in the text is not completely correct. Those lines are only tested with PySparse, which does return zero everywhere. Other solvers can return infinity or NaN. Regardless, the solution isn't right, which is the point of the illustration. I'll work on improving the explanation.