On Tue, Oct 25, 2016 at 12:34 PM, Zhekai Deng <[email protected]> wrote: > Thank you. I implemented as you suggested and it works very good now. This > approach seems also applicable for the mesh that are not aligned in the > interface.
Exactly. One thing though, it is probably a good idea to cache the mapping from one grid to the other as this is expensive to calculate (especially of you do this at every time step or sweep). We provided this functionality in FiPY, I'm not sure if Scipy lets you do that. -- Daniel Wheeler _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
