On Jul 24, 2008, at 12:11 AM, Zhiwen Liang wrote:
Now I am trying to solve a diffusion problem in 3D. And the
diffusivity is anisotropic.
But when I run it, I got the error of "Warning: invalid value
encountered in divide".
At what point? Immediately or after awhile? Your script runs fine for
me, so I can't diagnose. Please post the error trace.
What platform are you running on?
Precisely what version of FiPy?
What supporting libraries (NumPy, PySparse, etc.)?
P.S. I'd like to recommend once again not to start with such a big
mesh when you don't even know if things work yet. You can diagnose a
5**3 mesh a thousand times faster than a 50**3 mesh.
P.P.S. There is no need to "sweep" this problem as there are no non-
linearities
