Dear all, I managed to install fipy parallel on supercomputers. The examples/Diffusion/mesh20x20.py runs fine. I use phi.getGlobalValue() to get phi at each time and what I write in ascii "x-y-phi" file. The figures are the same you see on fipy website.
However, if I run the examples/cahnHilliard/mesh2D.py parallel; the script runs, but phi with using phi.getGlobalValue() does not really change in time. The script with pysparse on 1CPU produce the figures with nice spinodal decomposition (see fipy website). Though, if I change the equation in the examples/cahnHilliard/mesh2D.py from eq = (TransientTerm() == DiffusionTerm(coeff=D * a**2 * (1 - 6 * PHI * (1 - PHI))) - DiffusionTerm(coeff=(D, epsilon**2))) to eq = (TransientTerm() == DiffusionTerm(coeff=D * a**2 * (1 - 6 * PHI *(1 - PHI)))) then I see the drastic time evaluation of phi. Do you have some idea how to get then the cahn-hiliard script working? What is wrong in my case? I have already tried different solvers. Thank you for your help! Bests, Ferenc _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
