Hi Krishna,
According to this, fixed point iteration only gets you linear
convergence in general,
http://www.math-cs.gordon.edu/courses/ma342/handouts/rate.pdf
I think you need Newton for quadratic convergence.
Another issue to be aware of is that the error and residual are two
separate things.
Thanks a lot, now it perfectly works ! :-)
Damien
Message d'origine De : "Guyer, Jonathan E. Dr.
(Fed)" Date :05/10/2016 17:21
(GMT+01:00) À : FIPY Objet : Re: problem
with my initial conditions
I suspect the problem is that by doing `from math import *`, you are
importin
Good suggestion, but `u.value = ___` and `u.setValue(___)` *should* do the same
thing in this instance.
> On Oct 5, 2016, at 11:16 AM, Raymond Smith wrote:
>
> Hi, Damien.
>
> Try using u.setValue(0.5 + 0.5*sin(pi*x)*sin(pi*y)) and similarly for v.
> These lines could naturally be placed righ
I see the standard python math library imported a lot in examples provided by
fipy users.
Please don't do this.
The math library supplies nothing useful for fipy (anything useful in math is
available from numpy, which we provide as numerix). As seen in the example of
M. Chudeau, math can ove
I suspect the problem is that by doing `from math import *`, you are importing
a version of the `sin` function that only operates on single numbers, not on
fields. As a result, u and v get assigned uniform values.
`from numpy import *` would be better.
`from fipy.tools.numerix import *` would be
Hi, Damien.
Try using u.setValue(0.5 + 0.5*sin(pi*x)*sin(pi*y)) and similarly for v.
These lines could naturally be placed right after the definition of u and v
(certainly before the time stepping begins).
See
http://www.ctcms.nist.gov/fipy/examples/diffusion/generated/examples.diffusion.mesh1D.h
Hello everybody,
I am a French student and I am trying to model the emergence of patterns
on the fur of animals using the Turing equations (involving the
concentrations u and v) and to code it thanks to the FiPy library. A
pattern appearing only if the system is unstable, it is necessary that
