On Feb 20, 2015, at 12:08 PM, Daniel Wheeler wrote:
> On Fri, Feb 20, 2015 at 6:39 AM, Ronghai Wu 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 nee
On Feb 20, 2015, at 6:39 AM, Ronghai Wu wrote:
> (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)?
See the discussion on linearization of the source in
http://www.ctcms.nist.gov/fipy/exampl
On Fri, Feb 20, 2015 at 6:39 AM, Ronghai Wu 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 hav
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)?
(2) when solving the equation, sweep() controls the residual but solve()
does not. I am won
