> I pasted the gmsh commands into a file and ran gmsh on it. It threw a
> weird error
>
> ** On entry to DGESVD parameter number 6 had an illegal value
>
> What error did you get? Did you adapt the gmsh commands from
> somewhere? Maybe we could work from that and check that those commands
> work correctly. Googling the error indicates that it has something to
> do with the tolerance for recombination. Try rescaling the problem so
> the edge lengths are close to 1. Also cut and paste the gmsh commands
> into a geo file and run that against gmsh rather than using fipy
> making it easier to debug.
Hi Daniel,
I have set the edge lengths to 1 and placed the gmsh commands in a .geo
file as you suggested. I've attached the .geo file to this message. I also
removed the 4 lines and the last line of the gmsh commands as I think they
may not be required.
This is the error I get when I run the resulting code:
"Traceback (most recent call last):
File "C:\Users\SL\Documents\Python27\Gmsh Stuff\3DCylinder.py", line 8,
in <module>
mesh = GmshImporter3D("pipe.geo")
File "C:\Program
Files\Python27\lib\site-packages\fipy-2.1.3-py2.7.egg\fipy\meshes\numMesh\gmshImport.py",
line 454, in __init__
mesh.Mesh.__init__(self, **_DataGetter(mshfile.getFilename(),
dimensions=3).getData())
File "C:\Program
Files\Python27\lib\site-packages\fipy-2.1.3-py2.7.egg\fipy\meshes\numMesh\gmshImport.py",
line 196, in getData
vertexCoords = self._calcVertexCoords(self.coordDimensions)
File "C:\Program
Files\Python27\lib\site-packages\fipy-2.1.3-py2.7.egg\fipy\meshes\numMesh\gmshImport.py",
line 253, in _calcVertexCoords
maxNode = max(nodeToVertexIDdict.keys())
ValueError: max() arg is an empty sequence"
Point(1) = {0, 0, 0, cl};
Point(2) = {R, 0, 0, cl};
Point(3) = {0, R, 0, cl};
Point(4) = {-R, 0, 0, cl};
Point(5) = {0, -R, 0, cl};
Point(6) = {r, 0, 0, cl};
Point(7) = {0, r, 0, cl};
Point(8) = {-r, 0, 0, cl};
Point(9) = {0, -r, 0, cl};
Circle(11) = {3,1,4};
Circle(12) = {4,1,5};
Circle(13) = {5,1,2};
Circle(14) = {2,1,3};
Circle(15) = {7,1,8};
Circle(16) = {8,1,9};
Circle(17) = {9,1,6};
Circle(18) = {6,1,7};
Line(19) = {7,3};
Line(20) = {9,5};
Line Loop(21) = {19,11,12,-20,-16,-15};
Plane Surface(22) = {21};
Line Loop(23) = {20,13,14,-19,-18,-7};
Plane Surface(24) = {23};
Transfinite Line{11:18} = 6/2;
Transfinite Line{19,20} = 6;
Transfinite Surface{22,24} = {23,24,22,21};
Recombine Surface '*';
Extrude{0,0,length}{Surface{22,24}; Layers{4,1}; Recombine};
_______________________________________________
fipy mailing list
[email protected]
http://www.ctcms.nist.gov/fipy
[ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
