Now that I've changed it as you mentioned, it says Can't expand MemType 0. What can I do about this? http://pastebin.com/buFS1NRu
On Fri, Dec 20, 2013 at 10:24 AM, Daniel Wheeler <[email protected]>wrote: > On Fri, Dec 20, 2013 at 9:13 AM, Jane Hung <[email protected]> wrote: > > That definitely helps the simulation run, but the behavior is different > from > > what I've seen with other methods. So I'm wondering if you think the way > I'm > > rewriting the equations is alright http://pastebin.com/buFS1NRu > > I didn't check each equation against the maths, but the sweep loop is > definitely wrong. You have > > for sweep in range(sweeps): > res0=res1 = 1e100 > phi.updateOld() > psi.updateOld() > while max(res0, res1)>.1: > res0= eq11.sweep(var=phi, dt=dt) > print(res0) > res1= eq21.sweep(var=psi, dt=dt) > print(res1) > > I'm surprised that it even worked. It should be > > for sweep in range(sweeps): > res0=res1 = 1e100 > for v in all_variables: > v.updateOld() > while res > 0.1: > res= eq.sweep(dt=dt) > print(res) > > > where 'all_variables' is a list of all 6 variables that need to be > solved (regardless of whether you care about their values. > > Hope that helps. > _______________________________________________ > fipy mailing list > [email protected] > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] > -- Jane Hung Graduate Student | MIT Department of Chemical Engineering Hatton Lab 66-325 | Doyle Lab E18-509 [email protected] | 415.952.6325
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
