On Mon, Jul 15, 2013 at 5:34 AM, Edwin Sze Lun Khoo <[email protected]>wrote:
> Hello Daniel and Jonathan, > > Thanks for the suggestions! > > I have looked up Picard sweeps and Newton sweeps. I have implemented them > before in MATLAB but not in FiPy, so I was wondering if you could provide > some toy examples or some references to illustrate how both kinds of sweeps > can be coded in FiPy (I don't recall seeing such examples in the > documentation). > Edwin, The following http://www.ctcms.nist.gov/fipy/examples/reactiveWetting/generated/examples.reactiveWetting.liquidVapor1D.html is an example of Picard iterations at every time step (this uses a coupled solution, but that isn't necessary). Unfortunately, we don't have any examples of Newton steps right now. Sorry about that. We have used FiPy to solve equations using Newton though. It is very similar to Picard in terms of the coding, it's just that you take the variation of the equations with respect to the variables that are being solved for. > Regarding Jonathan's comments, I assume the residual that the sweep method > returns is the absolute, and not relative, residual? > Yes. It is the absolute. You would obviously have to pass in a reference residual in order for the equation to return a relative residual. Good luck and let us know how things proceed. Cheers. -- Daniel Wheeler
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
