On Thu, May 30, 2013 at 12:22 PM, Michelangelo Formisano <
[email protected]> wrote:
> The circle.py example works. Firstly it show me a green uniform circle
> and than diffusion starts (last image attached).
>
That seems to be working.
> Running the python code you suggested me, I find:
>
> Info : 0 vertices 0 elements
> Info : Writing '/tmp/tmpTBjIcy.msh'...
> Info : Done writing '/tmp/tmpTBjIcy.msh'
> Info : Stopped on Thu May 30 18:19:53 2013
>
I don't know what is going on here. Maybe the string being passed to the
GmshMesh is messed up. This is so strange. What happens when you put
cellSize = 0.05
radius = 1.
import fipy as fp
mesh = fp.Gmsh2D('''
cellSize = %(cellSize)g;
radius = %(radius)g;
Point(1) = {0, 0, 0, cellSize};
Point(2) = {-radius, 0, 0, cellSize};
Point(3) = {0, radius, 0, cellSize};
Point(4) = {radius, 0, 0, cellSize};
Point(5) = {0, -radius, 0, cellSize};
Circle(6) = {2, 1, 3};
Circle(7) = {3, 1, 4};
Circle(8) = {4, 1, 5};
Circle(9) = {5, 1, 2};
Line Loop(10) = {6, 7, 8, 9};
Plane Surface(11) = {10};
''' % locals())
print mesh.cellCenters
in a file rather than running at the command line?
Can you also run the full test suite "python setup.py test"?
--
Daniel Wheeler
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]