On Thu, May 30, 2013 at 11:48 AM, Michelangelo Formisano < [email protected]> wrote:
> Hi Daniel, > the output has this form: > > Info : Running 'gmsh tmp.geo -format msh -2 -nopopup tmp.msh' [1 > node(s), max. 4 thread(s)] > Info : Started on Thu May 30 17:40:41 2013 > > Info : Done meshing 1D (0 s) > Info : Meshing 2D... > Info : Meshing surface 6 (Plane, Delaunay) > Info : Done meshing 2D (0.072004 s) > Info : 4370 vertices 8744 elements > That's weird. You seem to have a completely different number of vertices and elements. At least it didn't give zero vertices and elements. I have noticed differences in meshes between versions of Gmsh, but such a large difference seems strange. I'm using version 2.6.2 by the way. Anyway, try running "examples/circle.py" ( https://matforge.org/fipy/browser/fipy/examples/diffusion/circle.py) and see what happens. Does it show something? Try also just running the following at the python command prompt >>> 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 Just to make sure that the mesh is being built correctly.
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
