On Fri, Feb 20, 2015 at 6:39 AM, Ronghai Wu <[email protected]> wrote: > Thanks, I worked it out. But still have three questions: > (1) The 4th-order is split by "psi = d2fdphi2(phi-phiold) + dfdphi - > epsilon**2*laplace phi". I do not understand why we need > d2fdphi2(phi-phiold)?
Where do you see this? Do you have a link to this? > (2) when solving the equation, sweep() controls the residual but solve() > does not. I am wondering what this "residual" exactly is? it is the > residual of what? It's just "Ax-b", the residual from the linearized system of equations (before the solution to the linearized system).l > (3) I noticed the recent issue "coupled 1D PDE on a periodic grid". I > want to implement 2D periodic boundary condition. The question is, in > addition to use PeriodicGrid2D(), should I do anything else to ensure > the periodic boundary condition? Probably not. Unless you are doing something customized then things should work right. -- Daniel Wheeler _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
