Re: Integral term

2018-07-13 Thread Daniel Wheeler
On Thu, Jul 12, 2018 at 9:59 AM, Pavel Aleynikov wrote: > How should such equation be arranged? fp.numerix.cumsum seems to return a numpy array so won't automatically update when f1 changes. Probably best to update in the loop like this, import fipy as fp mesh = fp.Grid1D(dx=0.002,

Integral term

2018-07-12 Thread Pavel Aleynikov
Hello, I would like to find a steady-state solution of an integro-differential equation with a spacial variable in the integration limit. So far it does not work. The minimal working example below gives wrong result. I understand that M1 and M2 are not updated in sweeps at all. If I set up M1