Dear Andy, Yes, it works. Thank you very much!
Best Regards, Rui On Sun, Mar 9, 2014 at 4:50 AM, Andrew Reeve <[email protected]> wrote: > Dear Rui, > > Does changing 'L' to a floating point value (40.0 instead of 40) > improve the first script? > > Andy > > > On Fri, Mar 7, 2014 at 7:48 PM, rui tu <[email protected]> wrote: > >> Dear FiPy group, >> >> I am trying to solve a very simple Cahn-Hilliard equation in FiPy, where >> I am using a scaling parameter, L, to adjust the dimensions of the domain. >> When I set the length, L = 40, the equation cannot be solved (test1.py). >> >> However, t if I set L = 0.4 e-5.*Magnification >> >> where >> >> Magnification = 1.0e7 >> >> which numerically is ALSO equal to 40, the equation becomes >> solvable(test2.py). I am baffled by this, and cannot understand why one >> simulation works and the other one doesn't. >> >> By typing, >> >> diff test1.py test2.py >> >> in a shell terminal you can find: >> >> 16c16,17 >> < L = 40 # dimension of simulation domain >> --- >> > Magnification = 1.0e7 # magnification >> > L = 0.4e-5*Magnification # dimension of simulation domain >> >> which proves that the both have the same L value. The only difference is >> how it is being inputed. I am attaching both scripts, test1.py and test2.py >> >> I would appreciate your help. >> >> Thanks, >> Rui >> >> _______________________________________________ >> fipy mailing list >> [email protected] >> http://www.ctcms.nist.gov/fipy >> [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] >> >> > > > -- > -------------------- > Andrew Reeve > School of Earth and Climate Sciences > University of Maine > Orono, ME 04469 > > _______________________________________________ > fipy mailing list > [email protected] > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] > >
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
