What I get is close to what Michelangelo has got, but still slightly
different:

fausto@fausto-quad:~/Programs/CFD/Gmsh$ gmsh tmp.geo -format msh -2 -nopopup 
tmp.msh
Info    : Running 'gmsh tmp.geo -format msh -2 -nopopup tmp.msh' [1 node(s), 
max. 4 thread(s)]
Info    : Started on Thu May 30 16:24:32 2013
Info    : Reading 'tmp.geo'...
Info    : Done reading 'tmp.geo'
Info    : Meshing 1D...
Info    : Meshing curve 1 (Line)
Info    : Meshing curve 2 (Line)
Info    : Meshing curve 3 (Line)
Info    : Meshing curve 4 (Line)
Info    : Meshing curve 5 (Line)
Info    : Meshing curve 6 (Line)
Info    : Done meshing 1D (0 s)
Info    : Meshing 2D...
Info    : Meshing surface 6 (Plane, Delaunay)
Info    : Done meshing 2D (0.124008 s)
Info    : 4378 vertices 8760 elements
Info    : Writing 'tmp.msh'...
Info    : Done writing 'tmp.msh'
Info    : Stopped on Thu May 30 16:24:33 2013

I'm using version 2.5.1.

Fausto





>________________________________
> From: Daniel Wheeler <[email protected]>
>To: Multiple recipients of list <[email protected]> 
>Sent: Thursday, May 30, 2013 1:11:42 PM
>Subject: Re: Problem with cahnHilliard.sphere example
> 
>
>
>
>
>
>
>
>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 ]
>
>
>
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
  [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]

Reply via email to