On Fri, Oct 12, 2012 at 9:03 AM, Jonathan Guyer <[email protected]> wrote: > The script below illustrates the way I would do this. > > Note that it crashes once h == hs because C is then zero and the > DiffusionTerm gets an infinite coefficient. > Multiplying the equation > through by C does not work for some reason; Wheeler might know why.
Thanks! That certainly is simpler. I had seen something like that approach in the archives of the mailing list, but without the conditions incorporated. I was confused by the rules or convention that FiPy uses to do lazy evaluation of expressions - allowing functional dependence. I guess I assumed it also handled functions lazily. I had started off with C multiplied through, and I was having problems with the convection term. I believe the remaining problems are related to my expressions for the parameters, and not their implementation. Kris _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
