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

CylindricalGrid1D mesh volumes

2020-01-15 Thread Pavel Aleynikov
Hi, How are "mesh.cellVolumes" defined in fp.CylindricalGrid1D case? The surface of a circle grid is not equal pi*r^2. import fipy as fp L= 1.; nx = 1000 mesh = fp.CylindricalGrid1D(dx=L/nx, nx=nx) print(mesh.cellVolumes.sum()) >> 0.5 Why not pi? How should I integrate Variables on

Re: CylindricalGrid1D mesh volumes

2020-01-17 Thread Pavel Aleynikov
t;theta = 1" is not a problem, as long as it's clear from the documentation. -- Pavel > On 15 Jan 2020, at 17:34, Daniel Wheeler wrote: > > On Wed, Jan 15, 2020 at 4:21 AM Pavel Aleynikov > wrote: >> >> Hi, >> >> How are "mesh.cellVolumes" de