Re: [Gmsh] Mesh Statistics
> On 11 Jan 2016, at 11:56, Paul Comino wrote:
>
> Hi there,
>
> Simple question: Is it possible to export the mesh statistics (File -> Save
> As -> Post-Processing - Mesh Statistics (.pos)) through the command line or
> in the .geo file after meshing?
>
Sure; for example:
lc = 1e-2;
Point(1) = {0, 0, 0, lc};
Point(2) = {.1, 0, 0, lc} ;
Point(3) = {.1, .3, 0, lc} ;
Point(4) = {0, .3, 0, lc} ;
Line(1) = {1,2} ;
Line(2) = {3,2} ;
Line(3) = {3,4} ;
Line(4) = {4,1} ;
Line Loop(5) = {4,1,-2,3} ;
Plane Surface(6) = {5} ;
// mesh
Mesh 2;
// what to print:
Print.PostElementary = 1;
Print.PostElement = 0;
Print.PostGamma = 0;
Print.PostEta = 0;
Print.PostRho = 1;
Print.PostDisto = 0;
// save (using automatic format detection using the extension)
Save "aa.pos";
> If it is possible how can it be done? If not is it possible to print Gamma,
> Rho, Eta to a file using printf?
>
> Thank you.
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Gmsh exits immediately after run, "Requesting more near neighbors than data points"
(0 , 1 , 0 ) > Debug : t2 : (1 , 0 , 0 ) > Debug : pt : (0 , 0 , 15.24 ) > Debug : Surface: 303 > Debug : SVD: 243.84,15.24,0 (min=2) > Debug : Plane : (1 x + 0 y + 0 z = 0) > Debug : Normal : (1 , 0 , 0 ) > Debug : t1 : (0 , 0 , 1 ) > Debug : t2 : (0 , 1 , 0 ) > Debug : pt : (0 , 0 , 0 ) > Debug : Surface: 304 > Debug : SVD: 45.72,15.24,0 (min=2) > Debug : Plane : (-0 x + -1 y + -0 z = -243.84) > Debug : Normal : (-0 , -1 , -0 ) > Debug : t1 : (0 , 0 , 1 ) > Debug : t2 : (1 , 0 , 0 ) > Debug : pt : (0 , 243.84 , 0 ) > Debug : Surface: 305 > Debug : SVD: 243.84,15.24,-0 (min=2) > Debug : Plane : (1 x + 0 y + 0 z = 45.72) > Debug : Normal : (1 , 0 , 0 ) > Debug : t1 : (0 , 0 , 1 ) > Debug : t2 : (0 , 1 , 0 ) > Debug : pt : (45.72 , 0 , 0 ) > Debug : Surface: 306 > Debug : SVD: 45.72,15.24,0 (min=2) > Debug : Plane : (-0 x + -1 y + -0 z = -0) > Debug : Normal : (-0 , -1 , -0 ) > Debug : t1 : (0 , 0 , 1 ) > Debug : t2 : (1 , 0 , 0 ) > Debug : pt : (0 , 0 , 0 ) > Debug : Surface: 307 > Debug : SVD: 15.24,2.54,0 (min=2) > Debug : Plane : (1 x + 0 y + 0 z = 21.59) > Debug : Normal : (1 , 0 , 0 ) > Debug : t1 : (0 , 0 , 1 ) > Debug : t2 : (0 , 1 , 0 ) > Debug : pt : (21.59 , 0 , 0 ) > Debug : Surface: 308 > Debug : SVD: 15.24,2.54,0 (min=2) > Debug : Plane : (-0 x + -1 y + -0 z = -123.19) > Debug : Normal : (-0 , -1 , -0 ) > Debug : t1 : (0 , 0 , 1 ) > Debug : t2 : (1 , 0 , 0 ) > Debug : pt : (0 , 123.19 , 0 ) > Debug : Surface: 309 > Debug : SVD: 15.24,2.54,-0 (min=2) > Debug : Plane : (1 x + 0 y + 0 z = 24.13) > Debug : Normal : (1 , 0 , 0 ) > Debug : t1 : (0 , 0 , 1 ) > Debug : t2 : (0 , 1 , 0 ) > Debug : pt : (24.13 , 0 , 0 ) > Debug : Surface: 310 > Debug : SVD: 15.24,2.54,0 (min=2) > Debug : Plane : (-0 x + -1 y + -0 z = -120.65) > Debug : Normal : (-0 , -1 , -0 ) > Debug : t1 : (0 , 0 , 1 ) > Debug : t2 : (1 , 0 , 0 ) > Debug : pt : (0 , 120.65 , 0 ) > Debug : Gmsh model (GModel) imported: > Debug : 24 Vertices > Debug : 28 Edges > Debug : 10 Faces > Debug : 1 Regions > Info: Done reading > 'C:\Programs\woodemc\lap_water-master\models\oak_hole.geo' > Info: Meshing 1D... > Info: Meshing curve 101 (Line) > ANN: ERROR--->Requesting more near neighbors than data > points<-ERROR > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Editing imported file in GUI
> On 03 Feb 2016, at 03:28, Aaron Kiefer wrote: > > GMSH Users, > > I’ve suffered through a few hours of removing volumetric parts of my.geo file > so it can be meshed in 2D. From working through the menus with GMSH, it > appears that the .geo file cannot be displayed with solid surfaces as a .step > file can? If anyone has ideas, please let me know – I’m going cross eyed > staring at the little crossed surface regions trying to determine which > surface belongs where. Indeed... As a workaround, just mesh the surfaces - you can then select them by clicking the mesh. > > Thanks, > > > Aaron J. Kiefer MSME, PE > Accident Research Specialists, PLLC > 1631 NW Maynard Road, Suite 101 > Cary, NC 27513 > Phone: (919) 467-8134 > Fax:(919) 678-1261 > Cell:(919) 810-8429 > Email:akie...@accident-research.com > Web: www.accident-research.com > Email Confidentiality Notice: The information contained in this transmission > is confidential, proprietary, and/or privileged. The message is intended for > the sole use of the individual or entity to whom it is addressed. If you are > not the intended recipient, please be advised that any use, distribution, or > copying of the message is strictly prohibited and may be subject to > penalties, If you received this transmission in error, please contact the > sender immediately and delete this material from any computer. > > > > From: Aaron Kiefer > Sent: Tuesday, February 2, 2016 9:06 AM > To: 'gmsh@onelab.info' > Subject: Editing imported file in GUI > > Gents, > > I’m attempting to edit a .geo file that I imported as a .step from Solidworks > using the GUI in GMSH 2.11. I’d like to display the surfaces as shaded > solids so I know which ones I’m selecting. The ‘solid’ surface display is > selected in the geometry/visibility tab but GMSH appears to be showing the > ‘wireframe’ or maybe the ‘cross’ surface display? > > Thanks in advance, > > Aaron J. Kiefer MSME, PE > Accident Research Specialists, PLLC > 1631 NW Maynard Road, Suite 101 > Cary, NC 27513 > Phone: (919) 467-8134 > Fax:(919) 678-1261 > Cell:(919) 810-8429 > Email:akie...@accident-research.com > Web: www.accident-research.com > Email Confidentiality Notice: The information contained in this transmission > is confidential, proprietary, and/or privileged. The message is intended for > the sole use of the individual or entity to whom it is addressed. If you are > not the intended recipient, please be advised that any use, distribution, or > copying of the message is strictly prohibited and may be subject to > penalties, If you received this transmission in error, please contact the > sender immediately and delete this material from any computer. > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] periodic mesh
> On 02 Feb 2016, at 20:59, Josefin Ahlkrona wrote:
>
> Hi,
>
> I have a cuboid (x,y,z) domain for which i have created tetrahedral,
> structured mesh. Now, I want to apply periodic boundary conditions in x, and
> y direction and in order to do that I need the mesh to be identical for the
> x=0 and x=xmax surfaces, as well as for the y=0 and y=ymax surfaces. It is
> not sufficient for only the nodes to match but I also want the edges to match!
>
> I have tried to use the "Periodic Surface" command without success. I have
> seen several posts on using the "Periodic Surface" commands without finding a
> solution for my problem.
>
The Periodic meshing command only apply to unstructured grids. In your case the
grid is set structured (meshed by extrusion with the "Layers" option).
> Could anyone help me?
>
> My geofile looks like this:
>
> cl=1.0;
> //cuboid which is 8000,8000,1000 big
>
> Point(1) = {0,0,0,cl};
> Extrude {8000,0,0} {
> Point{1}; Layers{10};
> }
> Extrude {0,8000,0} {
> Line{1}; Layers{10};
> }
> Extrude {0,0,1000} {
> Surface{5}; Layers{5};
> }
>
>
> boundMaster[] = Boundary{26};
> boundSlave[] = Boundary{-8}; //tried plus and minus sign here..
> boundMaster2[] = Boundary{22};
> boundSlave2[] = Boundary{-14}; //tried plus and minus sign here..
>
> Periodic Surface 18 {boundSlave[]} = 26 {boundMaster[]};
> Periodic Surface 14 {boundSlave2[]} = 22 {boundMaster2[]};
>
> Physical Surface(28) = {14}; //y=y0
> Physical Surface(29) = {18}; //x=xmax
> Physical Surface(30) = {22}; //y=ymax
> Physical Surface(31) = {26}; //x=x0
> Physical Surface(32) = {5};//top
> Physical Surface(33) = {27}; //bottom
>
> Physical Volume(34) = {1};
>
> Best,
>
> Josefin
>
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Bug: "print messages on terminal" option not perstent.
> On 02 Feb 2016, at 20:04, Jason C wrote: > > gmsh 2.11.0, windows. > > Minor issue: The value of the "print messages on terminal" option is not > saved with the other options. It is always initially unchecked on program > start. > Indeed - this is now fixed in SVN. Thanks for the report! > By contrast, options such as "message verbosity" save with no issue. > > Jason > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] periodic mesh
> On 03 Feb 2016, at 18:06, Josefin Ahlkrona wrote:
>
> Hi,
>
> Thank you for the answer. Is there a way to make structured meshes periodic?
> Since it's structured of course the nodes match, but I would like the edges
> to match eachother too!
>
It's possible, albeit tedious: you need to use the "Transfinite" algorithm for
that, and specify the orientions by hand (the automatic mode will not ensure
edge matching). Here's an example:
lc = 0.3;
// example of a purely hexahedral mesh using only transfinite
// mesh constraints
Point(1) = {-2,0,0,lc};
Point(2) = {-1,0,0,lc};
Point(3) = {-1,1,0,lc};
Point(4) = {-2,1,0,lc};
Point(5) = {-2,0,1,lc};
Point(6) = {-1,0,1,lc};
Point(7) = {-1,1,1,lc};
Point(8) = {-2,1,1,lc};
Line(1) = {4,3};
Line(2) = {3,2};
Line(3) = {2,1};
Line(4) = {1,4};
Line(6) = {5,6};
Line(7) = {6,7};
Line(8) = {7,8};
Line(9) = {8,5};
Line(10) = {1,5};
Line(11) = {4,8};
Line(12) = {2,6};
Line(13) = {3,7};
Line Loop(14) = {3,4,1,2};
Plane Surface(15) = {14};
Line Loop(16) = {6,7,8,9};
Plane Surface(17) = {16};
Line Loop(18) = {10,-9,-11,-4};
Plane Surface(19) = {18};
Line Loop(20) = {8,-11,1,13};
Plane Surface(21) = {20};
Line Loop(22) = {12,7,-13,2};
Plane Surface(23) = {22};
Line Loop(24) = {6,-12,3,10};
Plane Surface(25) = {24};
Surface Loop(1) = {17,-25,-23,-21,19,15};
Volume(1) = {1};
Transfinite Line{1:4,6:13} = 5;
// need to specify orientation (through vertex list) by hand to have correctly
// matching volume/surface edges
Transfinite Surface {15} = {1,2,3,4};
Transfinite Surface {17} = {5,6,7,8};
Transfinite Surface {19} = {1,5,8,4};
Transfinite Surface {21} = {4,8,7,3};
Transfinite Surface {23} = {2,6,7,3};
Transfinite Surface {25} = {1,5,6,2};
Transfinite Volume{1} = {1,2,3,4,5,6,7,8};
> Best Regards,
>
> Josefin
>
> On 02/03/2016 06:02 PM, Christophe Geuzaine wrote:
>>> On 02 Feb 2016, at 20:59, Josefin Ahlkrona
>>> wrote:
>>>
>>> Hi,
>>>
>>> I have a cuboid (x,y,z) domain for which i have created tetrahedral,
>>> structured mesh. Now, I want to apply periodic boundary conditions in x,
>>> and y direction and in order to do that I need the mesh to be identical for
>>> the x=0 and x=xmax surfaces, as well as for the y=0 and y=ymax surfaces. It
>>> is not sufficient for only the nodes to match but I also want the edges to
>>> match!
>>>
>>> I have tried to use the "Periodic Surface" command without success. I have
>>> seen several posts on using the "Periodic Surface" commands without finding
>>> a solution for my problem.
>>>
>> The Periodic meshing command only apply to unstructured grids. In your case
>> the grid is set structured (meshed by extrusion with the "Layers" option).
>>
>>
>>> Could anyone help me?
>>>
>>> My geofile looks like this:
>>>
>>> cl=1.0;
>>> //cuboid which is 8000,8000,1000 big
>>>
>>> Point(1) = {0,0,0,cl};
>>> Extrude {8000,0,0} {
>>> Point{1}; Layers{10};
>>> }
>>> Extrude {0,8000,0} {
>>> Line{1}; Layers{10};
>>> }
>>> Extrude {0,0,1000} {
>>> Surface{5}; Layers{5};
>>> }
>>>
>>>
>>> boundMaster[] = Boundary{26};
>>> boundSlave[] = Boundary{-8}; //tried plus and minus sign here..
>>> boundMaster2[] = Boundary{22};
>>> boundSlave2[] = Boundary{-14}; //tried plus and minus sign here..
>>>
>>> Periodic Surface 18 {boundSlave[]} = 26 {boundMaster[]};
>>> Periodic Surface 14 {boundSlave2[]} = 22 {boundMaster2[]};
>>>
>>> Physical Surface(28) = {14}; //y=y0
>>> Physical Surface(29) = {18}; //x=xmax
>>> Physical Surface(30) = {22}; //y=ymax
>>> Physical Surface(31) = {26}; //x=x0
>>> Physical Surface(32) = {5};//top
>>> Physical Surface(33) = {27}; //bottom
>>>
>>> Physical Volume(34) = {1};
>>>
>>> Best,
>>>
>>> Josefin
>>>
>>>
>>> ___
>>> gmsh mailing list
>>> gmsh@onelab.info
>>> http://onelab.info/mailman/listinfo/gmsh
>
> --
> Josefin Ahlkrona
> PhD Student
> Division of Scientific Computing
> Department of Information Technology
> Uppsala University
> Box 337
> SE-751 05 Uppsala
> Sweden
>
> Email: josefin.ahlkr...@it.uu.se
> Webpage: http://www.it.uu.se/katalog/josah992
> Ph: +46-18-4712980
> Fax: +46 18 523049, +46 18 511925
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Unexpected behavior (bug?)
> On 09 Feb 2016, at 21:21, Kópházi József wrote:
>
> Dear list members,
>
> I have encountered some unusual behavior while trying to generate a very
> simple structured tetrahedral mesh for two boxes. The attached images show
> the surface (triangle) and volumetric (tetrahedron) elements of the same
> mesh. However, some of the triangles do not "sit" properly on the face of the
> corresponding tetrahedron but they are "laid across" on two tetrahedra. This
> is marked with red arrows on the images. I have also attached the script
> resulting this mesh.
>
Yes : the "Transfinite Volume" algorithm in automatic mode (when the corner
vertices are not explicitly specified) will not try to match the surface mesh
(because in some cases it's impossible, depending on the surface mesh). Without
"Recombine" you should thus specify the orientions by hand:
lc = 0.3;
// example of a purely hexahedral mesh using only transfinite
// mesh constraints
Point(1) = {-2,0,0,lc};
Point(2) = {-1,0,0,lc};
Point(3) = {-1,1,0,lc};
Point(4) = {-2,1,0,lc};
Point(5) = {-2,0,1,lc};
Point(6) = {-1,0,1,lc};
Point(7) = {-1,1,1,lc};
Point(8) = {-2,1,1,lc};
Line(1) = {4,3};
Line(2) = {3,2};
Line(3) = {2,1};
Line(4) = {1,4};
Line(6) = {5,6};
Line(7) = {6,7};
Line(8) = {7,8};
Line(9) = {8,5};
Line(10) = {1,5};
Line(11) = {4,8};
Line(12) = {2,6};
Line(13) = {3,7};
Line Loop(14) = {3,4,1,2};
Plane Surface(15) = {14};
Line Loop(16) = {6,7,8,9};
Plane Surface(17) = {16};
Line Loop(18) = {10,-9,-11,-4};
Plane Surface(19) = {18};
Line Loop(20) = {8,-11,1,13};
Plane Surface(21) = {20};
Line Loop(22) = {12,7,-13,2};
Plane Surface(23) = {22};
Line Loop(24) = {6,-12,3,10};
Plane Surface(25) = {24};
Surface Loop(1) = {17,-25,-23,-21,19,15};
Volume(1) = {1};
Transfinite Line{1:4,6:13} = 5;
// need to specify orientation (through vertex list) by hand to have correctly
// matching volume/surface edges
Transfinite Surface {15} = {1,2,3,4};
Transfinite Surface {17} = {5,6,7,8};
Transfinite Surface {19} = {1,5,8,4};
Transfinite Surface {21} = {4,8,7,3};
Transfinite Surface {23} = {2,6,7,3};
Transfinite Surface {25} = {1,5,6,2};
Transfinite Volume{1} = {1,2,3,4,5,6,7,8};
Note that it's often less tedious to use "Extrude" with "Layers" instead of
Transfinite.
Christophe
> I generated these pictures using gmsh version 2.9.1 but some other versions
> also produce the same problem.
>
> Thank you for your help,
>
> Jozsef Kophazi
>
>
> marked.png>_______
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Issue with if statements and planes
> On 09 Feb 2016, at 08:08, Jimmie Lahti wrote:
>
> Hello,
>
> I have been trying to solve an issue where I can decide if the 2d plane shall
> have a hollow structure or not. Please note that I have tried changing
> numbering and if I comment any half of the if statement it works. So the
> issue seems to be the double definition of the Plane.
>
> The error message I get is:
> “””
> Info: Running 'gmsh -2 hollowbowtie.geo' [Gmsh 2.10.1, 1 node, max. 1
> thread]
The "Else" command was only added in Gmsh 2.11...
> Info: Started on Mon Feb 8 14:47:23 2016
> Info: Reading 'hollowbowtie.geo'...
> Error : 'hollowbowtie.geo', line 118 : syntax error (Plane)
> Info: Done reading 'hollowbowtie.geo'
> “””
>
> The code:
> “””
> 106 If (hollow_length > 0)
> 107 Circle(9) = {9, 21, 10};
> 108 Circle(10) = {11, 19, 12};
> 109 Circle(11) = {12, 18, 13};
> 110 Circle(12) = {13, 20, 14};
> 111 Line(13) = {9, 14};
> 112 Line(14) = {11, 10};
> 113 Line Loop(16) = {13, -12, -11, -10, 14, -9};
> 114
> 115 Plane Surface(17) = {15, 16};
> 116 Physical Surface(18) = {17};
> 117 Else
> 118 Plane Surface(30) = {15};
> 119 Physical Surface(31) = {30};
> 120 EndIf
> “””
>
>
> Thanks for your help,
> Jimmie Lahti
> Electronic engineer
>
> Mobile No number
> Direct +46 953 34 565
> jimmie.la...@guidelinegeo.com
>
> Guideline Geo AB
> http://www.guidelinegeo.com
> http://www.abem.se
> http://www.malags.se
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] visualization problem ?
> On 07 Feb 2016, at 19:25, walter steffe wrote: > > The attached file is a mesh generated by Gmsh library linked to EmCAD. > The mesh associated with Elementary Entity Surface number 4 (and also e > few other faces) is incomplete. > This surface mesh should copied from an other surface. EmCAD makes use > of GFace::setMeshMaster to enforce this condition. The affine > transformation relating the two faces is a rotation about the cylinder > axis. > > It is quite strange that the original face is properly meshed while the > copy is not. It is even more strange that other 18 copies of the > same original face (with different rotation angles) are properly meshed > while only 3 faces are incomplete. > > I do not understand well the mesh file format but I think that the slave > face meshes should just refer to the master face and the related > transformation. > > So I do not understand why there is such big difference between master > and slave meshes. Is it a visualization problem ?. > I don't think so: the actual mesh is indeed incomplete. The problem is most probably in the generation of the slave mesh. Quite a bit of work has been going on on periodic meshes in recent releases: I would try with the latest stable release to see if the problem persists. > > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] physical group labels
> On 28 Jan 2016, at 09:13, Владимир Кузнецов wrote:
>
> Hi, i have problem with physical group labels.
> I open gmsh.exe and create physical surface. In geo file i see something like
> this:
>
> Physical Surface(35) = {18, 20, 22, 24};
>
> But, i need label for this surface like this:
>
> Physical Surface("label", 35) = {18, 20, 22, 24};
>
> Can i label this surface without edit geo file in text editor?
>
We have added this is recent nightly snapshots. Give it a try...
> -
>
> Second question:
> Ok, i have geo file with physical surfaces without labels:
> Physical Surface(35) = {18, 20, 22, 24};
>
> can i get id for this physical surface?
> If physical surface has label, i just use this function:
>
> Gmodel *m;
> m->getPhysicalName(dimension, id);
> m->getPhysicalNumber(dimension, name);
>
> How can i get id for surfaces without labels?
>
> thanks in advice
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] visualization problem ?
> On 11 Feb 2016, at 15:20, walter.ste...@alice.it wrote: > > > > I al already using the last stable version (gmsh 2.11.0) > > I have investigated a little bit more and I have seen that the msh file > contains a few lines as the followings > 2 4 259 > Affine -0.9994 1.110223024625157e-16 0 0 -1.110223024625157e-16 > -0.9994 0 0 0 0 1 0 0 0 0 1 > 64 > 14624 23986 > -1 23987 > -1 23990 > > The problem is related with the -1 value associated with the tags of some > slave vertices. > I have used GDB with a breackpoint inside of the function > writeMSHPeriodicNodes (GModelIO_MSH.cpp line 574) > I have also used a break point condition (v1->getIndex() < 0) in order to > spot the problematic data. > So I have discovered that the slave vertex is good (for what concerns the > geometrical position) but has an index= -1. > > The problem now is that I do not know where the wrong index was generated. > The slave vertex was created inside of the copyMesh(...) function: > MVertex *vt =new MFaceVertex ... > But this function seems not to be the responsible for the index value of -1 > (the MFaceVertex constructor sets that value to 0). > The indexing is done in GModel::indexMeshVertices(). It looks like if the vertices keep the -1 index, it's because they do not belong to elements (triangles) that need to get saved !? Let us know if you find the issue... > And why most of the slave indices were properly renumbered with the exeption > of those associated with face 4 (having master face=259) and a few others ? > > > >I don't think so: the actual mesh is indeed incomplete. The problem is most > >probably in the generation of the slave mesh. Quite a bit of work has been > >going on on periodic meshes in recent releases: I would try with the latest > >stable release to see if the problem >persists. > > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Question about plugin crack
> On 06 Feb 2016, at 00:18, Babak Hassas Irani wrote: > > Hi, > > I want to insert edge cracks [notches] with *different lengths and angles* in > a 2D mesh. Is plugin crack capable to do it? > > I have read the documentation, but I'm not sure yet. > A priori, yes. Give it a try... > Best, > Babak > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] fatal error: 'MTrihedron.h' file not found
> On 15 Feb 2016, at 07:42, Felicity Graham wrote: > > Hi, > > I'm trying to install gmshpy on Mac 10.10 and get the error > src/gmshGeoPYTHON_wrap.cxx:3731:12: fatal error: 'MTrihedron.h' file not found > > #include "MTrihedron.h" > > I've downloaded http://gmsh.info/bin/MacOSX/gmsh-svn-MacOSX-dynamic.tgz, > untarred it, and followed the install directions at > http://qssi.cs.umt.edu/wiki/index.php/Setup. > > I can find "MTriangle.h", etc., in > gmsh-2.11.0-dynamic-svn-MacOSX/include/gmsh, however, there is no file called > "MTrihedron.h" in the bundle at all. > Indeed - this will be fixed in the next snapshot. Thanks for the report. > Any help would be appreciated. > > Cheers, > Felicity > > _______ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Bug in Gmsh 22699
> On 19 Feb 2016, at 15:26, Theler German Guillermo > wrote: > > Hello > > I've found that in Gmsh compiled from SVN version 22699 when merging a > BREP, General.FileName gets replaced with the location of the merged > file so the addition of physical groups fail because they are added to > the brep instead of being added to the geo. > Indeed - this should now be fixed. There's a lot of stuff going on in the SVN in order to enable fully interactive model construction (i.e. constructing and solving finite element problems with Gmsh without editing any files, by defining materials, boundary conditions, etc., interactively) that might affect previous behavior. Let us know if you see any other regressions. Christophe > The stable Gmsh version published in the web page does not have this > problem. > > > > > -- > Germán Theler :: CTO Ingeniería & TICs > > CITES – Centro de Innovación Tecnológica Empresarial y Social S.A. > Dirección General Sancor Seguros > Grupo Sancor Seguros > tel +54 3493 –428 500 – Int.: 3374 > gthe...@cites-gss.com > www.cites-gss.com - www.gruposancorseguros.com > > > > > Imprima este mensaje sólo si es absolutamente necesario. > Para imprimir, en lo posible utilice el papel de ambos lados. > El Grupo Sancor Seguros se compromete con el cuidado del medioambiente. > > > > AVISO DE CONFIDENCIALIDAD > > El Grupo Sancor Seguros comunica que: > > Este mensaje y todos los archivos adjuntos a el son para uso exclusivo del > destinatario y pueden contener información confidencial o propietaria, cuya > divulgación es sancionada por ley. Si usted recibió este mensaje > erróneamente, por favor notifíquenos respondiendo al remitente, borre el > mensaje original y destruya las copias (impresas o grabadas en cualquier > medio magnético) que pueda haber realizado del mismo. Todas las opiniones > contenidas en este mail son propias del autor del mensaje. La publicación, > uso, copia o impresión total o parcial de este mensaje o documentos adjuntos > queda prohibida. > > Disposición DNDP 10-2008. El titular de los datos personales tiene la > facultad de ejercer el derecho de acceso a los mismos en forma gratuita a > intervalos no inferiores a seis meses, salvo que acredite un interés legítimo > al efecto conforme lo establecido en el artículo 14, inciso 3 de la Ley > 25.326. La DIRECCIÓN NACIONAL DE PROTECCIÓN DE DATOS PERSONALES, Organo de > Control de la Ley 25.326, tiene la atribución de atender las denuncias y > reclamos que se interpongan con relación al incumplimiento de las normas > sobre la protección de datos personales. > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] trouble when setting Rotation center > Select
Dear Jean-Pierre, Indeed, the change in rotation center will introduce a change in the viewport unless no rotation was set before the change... Not sure how to fix this. Christophe > On 27 Jan 2016, at 08:55, jean pierre aubry wrote: > > hello > > when i do Rotation center > Select, after selecting the point > Gmsh changes the point of view and jumps to a "fit all" zoom view > this is a little annoying when one is fiddling for details in a very > close (almost microscopic) view > > is it possible to correct that so the point of view does not change > after the Rotation center > Select? > > jean pierre aubry -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] unexpected results
> On 11 Mar 2016, at 19:28, Theler German Guillermo
> wrote:
>
> Dear all
>
> I wrote the following .geo file but got the unexpected results shown in
> the attached images. I am using verions 2.12.0 as compiled from svn.
> Any suggestion?
It's one of the annoying quirks of Gmsh's internal CAD engine, for backward
compatibility. Setting "Geometry.ExactExtrusion = 0;" will do the trick.
>
>
> e = 3;
> theta = 1;
> R = e*e;
> h = 2*e*Pi;
>
> lc = 0.4;
>
>
> Point(1) = {0, 0, 0, lc};
> Point(2) = {R, 0, 0, lc};
> Point(3) = {0, -R, 0, lc};
> Point(4) = {0, R, 0, lc};
> Point(5) = {-R, 0, 0, lc};
>
> Circle(1) = {3, 1, 2};
> Circle(2) = {2, 1, 4};
> Circle(3) = {4, 1, 5};
> Circle(4) = {5, 1, 3};
> Line Loop(5) = {3, 4, 1, 2};
> Plane Surface(6) = {5};
>
> Extrude {0, 0, h} {
> Surface{6};
> }
>
> Rotate {{0.1, 0.2, 0.3}, {0, 0, 0}, theta} {
> Volume{1};
> }
>
>
>
> --
> Germán Theler :: CTO Ingeniería & TICs
>
> CITES – Centro de Innovación Tecnológica Empresarial y Social S.A.
> Dirección General Sancor Seguros
> Grupo Sancor Seguros
> tel +54 3493 –428 500 – Int.: 3374
> gthe...@cites-gss.com
> www.cites-gss.com - www.gruposancorseguros.com
>
>
>
>
> Imprima este mensaje sólo si es absolutamente necesario.
> Para imprimir, en lo posible utilice el papel de ambos lados.
> El Grupo Sancor Seguros se compromete con el cuidado del medioambiente.
>
>
>
> AVISO DE CONFIDENCIALIDAD
>
> El Grupo Sancor Seguros comunica que:
>
> Este mensaje y todos los archivos adjuntos a el son para uso exclusivo del
> destinatario y pueden contener información confidencial o propietaria, cuya
> divulgación es sancionada por ley. Si usted recibió este mensaje
> erróneamente, por favor notifíquenos respondiendo al remitente, borre el
> mensaje original y destruya las copias (impresas o grabadas en cualquier
> medio magnético) que pueda haber realizado del mismo. Todas las opiniones
> contenidas en este mail son propias del autor del mensaje. La publicación,
> uso, copia o impresión total o parcial de este mensaje o documentos adjuntos
> queda prohibida.
>
> Disposición DNDP 10-2008. El titular de los datos personales tiene la
> facultad de ejercer el derecho de acceso a los mismos en forma gratuita a
> intervalos no inferiores a seis meses, salvo que acredite un interés legítimo
> al efecto conforme lo establecido en el artículo 14, inciso 3 de la Ley
> 25.326. La DIRECCIÓN NACIONAL DE PROTECCIÓN DE DATOS PERSONALES, Organo de
> Control de la Ley 25.326, tiene la atribución de atender las denuncias y
> reclamos que se interpongan con relación al incumplimiento de las normas
> sobre la protección de datos personales.
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] periodic mesh problem
Surface{22};
> Line Loop(5018)={5015,5016,5017};
> Plane Surface(5019)={5018};
> Physical Surface(12)={5019};
> Surface{5019} In Volume{1};
> ...
> best regards,
> Xiaoxin
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Options for Packing Quads and Bunin Opt
> On 07 Mar 2016, at 20:47, Chaman Singh Verma wrote: > > Hello, > > > > In the geo file, how should we use Packing Parallellogram option, and Bunin > Optimization ? > These options have actually been removed in 2.12 > Thanks, > csv > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Patch: added export menu
Hi Rui, Thanks for the proposed patch. I thing the "File->Save As" menu should actually be replaced completely with "File->Export", as we currently never save the model (.geo) file with Gmsh (since it's created on the fly), and we only save (export) representations of the model (meshes, images, etc.) which can not be read back and produce the current state. Maybe file a bug with the idea so we won't forget? Christophe > On 04 Mar 2016, at 12:16, Rui Maciel wrote: > > Hi, everyone. > > Exporting data to a file, whether it is a mesh definition or image, is > one of the main use cases of a mesh generator. I also noticed that > the question of how to export data is posted often in Gmsh's mailing > list. Right now, Gmsh provides this functionality through the "Save > As..." menu, and although this works well I've found that sometimes it > would be preferable if this feature could be used in a tad more > straight-forward way. > > To make this use case a bit more straight-forward, I tweaked Gmsh's > file menu to add an "Export" menu, which includes options to export > the mesh, postprocessing info, image, or a movie. This menu offers > users a clear way to do what they want to do, and the ability to > unequivocally discover where they need to go in the GUI to export > specific bits of information. > > I've attached the patch that implements this feature. The patch was > created with regards the current stable version (v2.11.0) as was made > available through Gmsh's site. > > The patch adds the "Export" menu, the respective submenu entries, and > a callback for each export option. The export callbacks were based on > the file_save_as_cb() callback function. > > The "Save as" menu was left untouched due to potential implications > regarding backward compatibility. > > If this patch is of any interest, feel free to use it. > > > Best regards, > Rui Maciel > <0001-Added-export-menu-and-menu-items.patch>___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Msh-file element type sequence
> On 29 Feb 2016, at 16:31, Benjamin Isbarn wrote: > > Hi, > > I'm wondering if it is safe to assume a constant sequence regarding > the element types in the MSH2 output files?! > > Judging by the source code in the file 'GModelIO_MSH2.cpp' the > sequence is as follows (with Mesh.SaveAll=1, or no physical entities): > > $Elements > > points > lines > triangles > quads > polygons > tets > hexas > prisms > ... > > I just wanted to verify that i got this right and if this sequence is > somehow future proof, so that newer versions won't break this, since > I'd like to rely on this sequence for my parser. > Gmsh will continue to save elements in that order; but the format does not enforce it... > Best regards, > > Benjamin Isbarn > > _______ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Question sur GMSH (PERRIN Pierre - CEREMA/DTerEst/Infra/DOA)
>>> *Pierre PERRIN *
>>> *Ing?nieur d'?tudes Ouvrages d'Art
>>> D?partement Conception et Exploitation des Infrastructures
>>> Division Ouvrages d'Art *
>>> *T?l.: : +33(0)3 87 20 46 39
>>> * */Mobilis? pour sauver le site de Metz/*
>>>
>>>
>>> Centre d??tudes et d?expertise sur les risques, l?environnement, la
>>> mobilit? et l?am?nagement www.cerema.fr <http://www.cerema.fr>
>>> Direction territoriale Est
>>> 1, boulevard Solidarit? BP 85230 57076 METZ CEDEX 3 - T?l : +33 (0)3 87
>>> 20 43 00
>>> Si?ge social : Cit? des Mobilit?s - 25, avenue Fran?ois Mitterrand - CS
>>> 92 803 - F-69674 Bron Cedex - T?l : +33 (0)4 72 14 30 30
>>>
>>>
>>>
>>> -- next part --
>>> An HTML attachment was scrubbed...
>>> URL:
>>> <http://onelab.info/pipermail/gmsh/attachments/20160318/1e8e28a6/attachment-0001.html>
>>>
>>> -- next part --
>>> //PARAMETRES
>>> h=0.1;
>>> cote_x=20.0 ;
>>> cote_y=10.0 ;
>>> deb_impact_x=5.0;
>>> deb_impact_y=5.0;
>>> impact_x=0.6 ;
>>> impact_y=0.4 ;
>>>
>>>
>>> //POINTS
>>> Point(1)={0,0,0,h};
>>> Point(2)={cote_x,0,0,h};
>>> Point(3)={cote_x,cote_y,0,h};
>>> Point(4)={0,cote_y,0,h};
>>> Point(5)={deb_impact_x,deb_impact_y,0,h};
>>> Point(6)={deb_impact_x+impact_y,deb_impact_y,0,h};
>>> Point(7)={deb_impact_x+impact_y,deb_impact_y+impact_y,0,h};
>>> Point(8)={deb_impact_x,deb_impact_y+impact_y,0,h};
>>>
>>> //LINES
>>> Line(1)={1,2};
>>> Line(2)={2,3};
>>> Line(3)={3,4};
>>> Line(4)={4,1};
>>> Line(5)={5,6};
>>> Line(6)={6,7};
>>> Line(7)={7,8};
>>> Line(8)={8,5};
>>>
>>> //SURFACES
>>> Line Loop(1)={1,2,3,4};
>>> Plane Surface(1)={1};
>>> Transfinite Surface{1};
>>> Recombine Surface{1};
>>>
>>> // incorporation des lignes dans les surfaces :
>>> // En GMSH ancien (v1.60)
>>> //??
>>> // En GMSH moderne (v2.12)
>>> //Line {5,6,7,8} In Surface {1};
>>>
>>> Line {5} In Surface {1};
>>> Line {6} In Surface {1};
>>> Line {7} In Surface {1};
>>> Line {8} In Surface {1};
>>>
>>>
>>> Line Loop(2)={5,6,7,8};
>>> Plane Surface(2)={2};
>>> Transfinite Surface{2};
>>> Recombine Surface{2};
>>>
>>>
>>> // Essais inutile au final :
>>> //Compound Surface(3)={1,2};
>>> //Transfinite Surface{3};
>>> //Mesh.RemeshAlgorithm=1;
>>> Coherence;
>>> Mesh.Format=1;
>>>
>>> //PHYSICAL
>>> Physical Line(1)={1, 3}; //section d'encastrement
>>> Physical Line(2)={2, 4}; //section d'encastrement
>>> Physical Surface(1)={1};
>>> Physical Surface(2)={2};
>>> Physical Surface(3)={1,2};
>>>
>>> //Color Yellow { Surface {1} ; }
>>> //Color Red { Line {2} ; }
>>>
>>>
>>>
>>>
>>>
>>> --
>>>
>>> Subject: Digest Footer
>>>
>>> ___
>>> gmsh mailing list
>>> gmsh@onelab.info
>>> http://onelab.info/mailman/listinfo/gmsh
>>>
>>>
>>> --
>>>
>>> End of gmsh Digest, Vol 158, Issue 15
>>> *
>>
>>
>>
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] detect negative jacobians
> On 24 Mar 2016, at 19:06, Christophe Geuzaine wrote: > >> >> On 24 Mar 2016, at 12:27, Bernd Hahnebach wrote: >> >> >> What is the recommended way to avoid such negative jacobians? Attached the >> geometry of the mesh above. >> And to avoid such negative jacobians, you can use the "High order tools" (in the Mesh menu). This is still experimental, but works quite well. The algorithms are described (in part) in T. Toulorge, C. Geuzaine, J.-F. Remacle, J. Lambrechts. Robust untangling of curvilinear meshes. Journal of Computational Physics 254, pp. 8-26, 2013. http://gmsh.info/doc/preprints/gmsh_untangling_preprint.pdf Christophe -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] which med version is supported with gmsh .2.11.0 and above
Hi Christophe, I don't think we ever tried compiling with MED >= 3.0.6. Where can we get the latest MED release? (Any big changes in this release?) Christophe > On 25 Mar 2016, at 09:35, Christophe Trophime > wrote: > > Hi, > I used to compile gmsh 2.11.0 with MED support. > It works for until I move to MED 3.1.0. > > With MED 3.1.0 when I save a MED mesh into a GMSH mesh > the elements (ie tetra) are not saved in .msh file!!! > > That's very annoying. > Did someone experience that issue? > Is it fix in recent gmsh version? > > Best > C > > > > Christophe TROPHIME > Research Engineer > > CNRS - LNCMI > 25, rue des Martyrs > BP 166 > 38042 GRENOBLE Cedex 9 > FRANCE > > Tel : +33 (0)4 76 88 90 02 > Fax : +33 (0) 4 76 88 10 01 > Office U 19 > M@il : christophe.troph...@lncmi.cnrs.fr > > _______ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] which med version is supported with gmsh .2.11.0 and above
> On 30 Mar 2016, at 09:15, Christophe Trophime > wrote: > > > - Original Message - >> From: "Christophe Geuzaine" >> To: "Christophe Trophime" >> Cc: g...@geuz.org >> Sent: Tuesday, March 29, 2016 10:33:57 PM >> Subject: Re: [Gmsh] which med version is supported with gmsh .2.11.0 and >> above >> >> >> Hi Christophe, >> >> I don't think we ever tried compiling with MED >= 3.0.6. Where can we get the >> latest MED release? (Any big changes in this release?) > > Hi Christophe, > I've tried with MED 3.8.0 and it works like a charm > The changes happen with 3.10 version. > You could retreive latest version from Salome > http://www.salome-platform.org/downloads/current-version > Weird, I recompiled MED 3.1.0 and everything seems to work ok here. Can you send a small test case with instructions on how to reproduce the issue? Thanks, Christophe >> >> Christophe >> >>> On 25 Mar 2016, at 09:35, Christophe Trophime >>> wrote: >>> >>> Hi, >>> I used to compile gmsh 2.11.0 with MED support. >>> It works for until I move to MED 3.1.0. >>> >>> With MED 3.1.0 when I save a MED mesh into a GMSH mesh >>> the elements (ie tetra) are not saved in .msh file!!! >>> >>> That's very annoying. >>> Did someone experience that issue? >>> Is it fix in recent gmsh version? >>> >>> Best >>> C >>> >>> >>> >>> Christophe TROPHIME >>> Research Engineer >>> >>>CNRS - LNCMI >>> 25, rue des Martyrs >>> BP 166 >>> 38042 GRENOBLE Cedex 9 >>> FRANCE >>> >>> Tel : +33 (0)4 76 88 90 02 >>> Fax : +33 (0) 4 76 88 10 01 >>> Office U 19 >>> M@il : christophe.troph...@lncmi.cnrs.fr >>> >>> ___ >>> gmsh mailing list >>> gmsh@onelab.info >>> http://onelab.info/mailman/listinfo/gmsh >> >> -- >> Prof. Christophe Geuzaine >> University of Liege, Electrical Engineering and Computer Science >> http://www.montefiore.ulg.ac.be/~geuzaine >> >> Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be >> Free software: http://gmsh.info | http://getdp.info | http://onelab.info >> >> -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] OneLab (socket?) remote controllable?
> On 26 Mar 2016, at 20:35, Todd Pierce wrote: > > Hi Experts, > > I have developed a natural language processor which is supposed to be the > front end to a simulation environment. OneLab would be perfect, if I knew > how to redirect all of the inputs and outputs. I've seen people control > GNUplot seamlessly from other software and it would be nice to do this with > OneLab. > > Here's what I would like to do: > > I'd like to feed the code to Gmsh (pipe? socket?) and see the results in the > CAD window as it is being illustrated. > I would also like to also be able to select and configure the solvers > remotely. > When that is done, perform the rendering and have it visible on the screen. > Have a look at the examples in the source distribution: gmsh/utils/solvers/c++ for simple C++ clients gmsh/utils/solvers/python for simple Python clients Christophe > Any ideas? > > Thanks, > > -Todd > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] How to add huge number of facets to a physical surface group quickly?
> On 25 Mar 2016, at 17:03, Yang, Jun wrote: > > Hi, > I'm trying to build a nano-structure model using Gmsh. The model contains > hundreds of cubes inside. Is there a way that I can add all the facets of > these cubes to a physical surface group quickly? It is painful to choose the > facets one by one with a mouse. Ctrl+left click to start a rectangular selection; left click to add; shift+left click to remove. > Thank you! > > Jun > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] detect negative jacobians
> On 30 Mar 2016, at 07:30, Bernd Hahnebach wrote: > > Thanks Christophe! > > I got the problem of the high order meshes by reading a big part of the > paper. :-) :-) Is there any user documentation or an example how to use the > HighOrderTools around? > > Just to get it the right way done. How does it work in the gui? > 1) Load the geometry 2) Mesh (e.g. press "3" to do a 3D mesh) 3) Open "High order tools" - Choose the polynomial order and press "Generate" (this will create high order elements) - Press "Regularize" (this will optimize the high order elements) > - Load the geometry (in my case a BRep) > - make the adjustments on Tools --> Options --> Geometry and Mesh --> one of > them is 2nd order mesh > - In GMSH Tree press Mesh 1D, 2D, 3D > - on 3D I get negative Jacobians if the mesh order in options is set to two > - do not use any optimization, do not use the order buttons (mesh is 2nd > order already) > - open HighOrderTools > - Generate 2nd order vertices > - Regularize with given values > > - Or do I have to start with a 1st order mesh and use Generate 2nd order > vertices on a 1st order mesh? > - Should any of the optimizations be done before the HighOrderTools? > > kind regards bernd > > > Zitat von Christophe Geuzaine : > >> >>> On 24 Mar 2016, at 19:06, Christophe Geuzaine wrote: >>> >>>> >>>> On 24 Mar 2016, at 12:27, Bernd Hahnebach wrote: >>>> >>>> >>>> What is the recommended way to avoid such negative jacobians? Attached the >>>> geometry of the mesh above. >>>> >> >> And to avoid such negative jacobians, you can use the "High order tools" (in >> the Mesh menu). This is still experimental, but works quite well. The >> algorithms are described (in part) in >> >> T. Toulorge, C. Geuzaine, J.-F. Remacle, J. Lambrechts. Robust untangling of >> curvilinear meshes. Journal of Computational Physics 254, pp. 8-26, 2013. >> http://gmsh.info/doc/preprints/gmsh_untangling_preprint.pdf >> >> Christophe >> >> >> -- >> Prof. Christophe Geuzaine >> University of Liege, Electrical Engineering and Computer Science >> http://www.montefiore.ulg.ac.be/~geuzaine >> >> Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be >> Free software: http://gmsh.info | http://getdp.info | http://onelab.info >> >> >> > > > > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] bug conversion maillage Gmsh vers Abaqus
Hello Julien, Indeed, HEX27 export is not coded yet for Abaqus. Could you have a look at the MHexahedron27 class in https://onelab.info/svn/gmsh/trunk/Geo/MHexahedron.h (login and passwd: gmsh) and provide the little "virtual MVertex *getVertexINP(int num)" function? Thanks, Christophe > On 04 Apr 2016, at 16:02, troufflard wrote: > > Bonjour, > > je pense avoir trouvé un bug (Gmsh 2.11.0 sous Mac OS X). Je n'ai pas pu > accédé à la création d'un ticket via le lien : > https://onelab.info/trac/gmsh/newticket > > Je poste donc ici le problème rencontré à propos de la conversion Gmsh vers > Abaqus d'un maillage constitué d'un élément type 12 (27-node hexaedron). > ci-joint un exemple avec un Readme pour expliquer. > > cordialement, > Julien Troufflard > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] compile with med support
> On 04 Apr 2016, at 17:49, Bernd Hahnebach wrote:
>
> According CMakeLists.txt is MED support activated by default
>
> opt(MED "Enable MED mesh and post file formats" ${DEFAULT})
>
>
> But at my cmake output it is missing, but cmake does not complain some
> library is missing ?!? How do I set the path to libmedc1 for cmake gmsh?
>
Bernd - Gmsh looks for a library named libmed (.a, .so, .dylib, etc.): not sure
why yours is named libmedc1?
You can modify the line
find_library(MED_LIB med)
into
find_library(MED_LIB medc1)
in the CMakeLists.txt file to see if this fixes it.
> cheers bernd
>
>
> cmake gmsh
>
> -- Gmsh 2.12.1 has been configured for Linux
> --
> -- * Build options: Ann Bamg Bfgs Blas(Generic) Blossom Cairo Chaco
> DIntegration Dlopen Fltk GMP Gmm Jpeg Kbipack Lapack(Generic) LinuxJoystick
> MathEx Mesh Metis Mmg3d Mpeg NativeFileChooser Netgen ONELAB ONELABMetamodel
> OpenCascade OpenGL OptHom Parser Plugins Png Post Salome Solver Taucs
> TetGen/BR Tetgen1.5 Voro3D Zlib
>
>
>
>
> Zitat von Bernd Hahnebach :
>
>> Hi gmsh folks,
>>
>> just compiled gmsh on debian jessie. It worked like a charm. I have libmedc1
>> installed but gmsh compiled without med support. How is support for med
>> export activated in cmake?
>>
>> cheers bernd
>>
>>
>>
>> ___
>> gmsh mailing list
>> gmsh@onelab.info
>> http://onelab.info/mailman/listinfo/gmsh
>>
>
>
>
>
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Scaling problem
> On 04 Apr 2016, at 22:30, David Colignon wrote: > > Hi, > > Can you try with: > > Geometry.ScalingFactor > >Global geometry scaling factor >Default value: 1 >Saved in: General.OptionsFileName > You can also select multiple points in the interface by pressing "Ctrl+left mouse button", or simply do a For loop in the .geo file and iterate over the points. > Regards, > > Dave > > -- > David Colignon, Ph.D. > 1er logisticien de recherche > Université de Liège > ACE - Applied & Computational Electromagnetics > Sart-Tilman B28 > 10, Grande Traverse > 4000 Liège - BELGIQUE > Tél: +32 (0)4 366 37 32 > http://www.ulg.ac.be/nic4 > > On 04/04/16 21:03, nbouv...@student.ulg.ac.be wrote: >> Hello, >> >> I was wondering if it's possible to change the scaling of every point >> defined in my code with a factor 100 immediately in the .geo file? >> >> I am actually looking for an iterative methode as the scaling command of the >> Gmsh interface requires (if i'm right) to select each point one by one. The >> underlying problem of my situation is that i've already created one thousand >> of points... >> >> Kind regards. >> >> Nicolas Bouvier. >> >> ___ >> gmsh mailing list >> gmsh@onelab.info >> http://onelab.info/mailman/listinfo/gmsh >> > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] detect negative jacobians
> On 31 Mar 2016, at 19:26, Bernd Hahnebach wrote: > > Thanks Christophe, got again the point! Works quit well for me now. > > Is it possible to script the hight order tools in a geo file? I could not > find anything in this regard (only the window position ?!?) > You can use: Mesh.ElementOrder = 2; Mesh.HighOrderOptimize = 1; > bernd > > > Zitat von Christophe Geuzaine : > >> >>> On 30 Mar 2016, at 07:30, Bernd Hahnebach wrote: >>> >>> Thanks Christophe! >>> >>> I got the problem of the high order meshes by reading a big part of the >>> paper. :-) :-) Is there any user documentation or an example how to use >>> the HighOrderTools around? >>> >>> Just to get it the right way done. How does it work in the gui? >>> >> >> 1) Load the geometry >> 2) Mesh (e.g. press "3" to do a 3D mesh) >> 3) Open "High order tools" >> - Choose the polynomial order and press "Generate" (this will create high >> order elements) >> - Press "Regularize" (this will optimize the high order elements) >> >> >>> - Load the geometry (in my case a BRep) >>> - make the adjustments on Tools --> Options --> Geometry and Mesh --> one >>> of them is 2nd order mesh >>> - In GMSH Tree press Mesh 1D, 2D, 3D >>> - on 3D I get negative Jacobians if the mesh order in options is set to two >>> - do not use any optimization, do not use the order buttons (mesh is 2nd >>> order already) >>> - open HighOrderTools >>> - Generate 2nd order vertices >>> - Regularize with given values >>> >>> - Or do I have to start with a 1st order mesh and use Generate 2nd order >>> vertices on a 1st order mesh? >>> - Should any of the optimizations be done before the HighOrderTools? >>> >>> kind regards bernd >>> >>> >>> Zitat von Christophe Geuzaine : >>> >>>> >>>>> On 24 Mar 2016, at 19:06, Christophe Geuzaine wrote: >>>>> >>>>>> >>>>>> On 24 Mar 2016, at 12:27, Bernd Hahnebach wrote: >>>>>> >>>>>> >>>>>> What is the recommended way to avoid such negative jacobians? Attached >>>>>> the geometry of the mesh above. >>>>>> >>>> >>>> And to avoid such negative jacobians, you can use the "High order tools" >>>> (in the Mesh menu). This is still experimental, but works quite well. The >>>> algorithms are described (in part) in >>>> >>>> T. Toulorge, C. Geuzaine, J.-F. Remacle, J. Lambrechts. Robust untangling >>>> of curvilinear meshes. Journal of Computational Physics 254, pp. 8-26, >>>> 2013. http://gmsh.info/doc/preprints/gmsh_untangling_preprint.pdf >>>> >>>> Christophe >>>> >>>> >>>> -- >>>> Prof. Christophe Geuzaine >>>> University of Liege, Electrical Engineering and Computer Science >>>> http://www.montefiore.ulg.ac.be/~geuzaine >>>> >>>> Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be >>>> Free software: http://gmsh.info | http://getdp.info | http://onelab.info >>>> >>>> >>>> >>> >>> >>> >>> >>> >>> ___ >>> gmsh mailing list >>> gmsh@onelab.info >>> http://onelab.info/mailman/listinfo/gmsh >> >> -- >> Prof. Christophe Geuzaine >> University of Liege, Electrical Engineering and Computer Science >> http://www.montefiore.ulg.ac.be/~geuzaine >> >> Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be >> Free software: http://gmsh.info | http://getdp.info | http://onelab.info >> >> >> > > > > -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] compile with med support
> On 05 Apr 2016, at 19:19, Bernd Hahnebach wrote:
>
> Found it. It seams a known problem on Debian Jessie (I do not know about
> other debian). See https://bugs.debian.org/cgi-bin/bugreport.cgi?bug=756472
>
> Makeing this changes to the CMakeLists.txt did the trick. HDF5 where found.
>
Thanks - I've merge the Debian patch into our CMakeLists.txt.
> cheers bernd
>
>
> Zitat von Bernd Hahnebach :
>
>> For med support on Debian Jessie package libmedc-dev has to be installed.
>> Watchout for the c !
>>
>> But the culprit is somewhere else. Cmake only searches for libmed if HDF5
>> was found. OK run cmake ../ again.
>> YEAH
>> -- HDF5_LIB-NOTFOUND
>> -- HDF5 not found
>>
>> ok I searched for them, and they seam installed. OK
>>
>> cmake \
>> -DHDF5_LIB=/usr/lib/x86_64-linux-gnu/libhdf5_cpp.so.8.0.2 \
>> ../
>>
>> BUT
>>
>> [ 1%] Building CXX object CMakeFiles/gmsh.dir/Common/CommandLine.cpp.o
>> In file included from
>> /home/hugo/Documents/dev/gmsh/gmsh-dev/Common/CommandLine.cpp:41:0:
>> /usr/include/med.h:22:18: fatal error: hdf5.h: Datei oder Verzeichnis nicht
>> gefunden
>> #include
>> ^
>> compilation terminated.
>>
>>
>> BUT
>>
>> $ locate hdf5.h
>> /usr/include/hdf5/serial/hdf5.h
>> /usr/include/vtk-6.1/vtk_hdf5.h
>>
>>
>> Which results in how do I set the include Dir for hdf5 ?
>>
>> Bernd
>>
>>
>> Zitat von Christophe Geuzaine :
>>
>>>
>>>> On 04 Apr 2016, at 17:49, Bernd Hahnebach wrote:
>>>>
>>>> According CMakeLists.txt is MED support activated by default
>>>>
>>>> opt(MED "Enable MED mesh and post file formats" ${DEFAULT})
>>>>
>>>>
>>>> But at my cmake output it is missing, but cmake does not complain some
>>>> library is missing ?!? How do I set the path to libmedc1 for cmake gmsh?
>>>>
>>>
>>> Bernd - Gmsh looks for a library named libmed (.a, .so, .dylib, etc.): not
>>> sure why yours is named libmedc1?
>>>
>>> You can modify the line
>>>
>>> find_library(MED_LIB med)
>>>
>>> into
>>>
>>> find_library(MED_LIB medc1)
>>>
>>> in the CMakeLists.txt file to see if this fixes it.
>>>
>>>> cheers bernd
>>>>
>>>>
>>>> cmake gmsh
>>>>
>>>> -- Gmsh 2.12.1 has been configured for Linux
>>>> --
>>>> -- * Build options: Ann Bamg Bfgs Blas(Generic) Blossom Cairo Chaco
>>>> DIntegration Dlopen Fltk GMP Gmm Jpeg Kbipack Lapack(Generic)
>>>> LinuxJoystick MathEx Mesh Metis Mmg3d Mpeg NativeFileChooser Netgen ONELAB
>>>> ONELABMetamodel OpenCascade OpenGL OptHom Parser Plugins Png Post Salome
>>>> Solver Taucs TetGen/BR Tetgen1.5 Voro3D Zlib
>>>>
>>>>
>>>>
>>>>
>>>> Zitat von Bernd Hahnebach :
>>>>
>>>>> Hi gmsh folks,
>>>>>
>>>>> just compiled gmsh on debian jessie. It worked like a charm. I have
>>>>> libmedc1 installed but gmsh compiled without med support. How is support
>>>>> for med export activated in cmake?
>>>>>
>>>>> cheers bernd
>>>>>
>>>>>
>>>>>
>>>>> ___
>>>>> gmsh mailing list
>>>>> gmsh@onelab.info
>>>>> http://onelab.info/mailman/listinfo/gmsh
>>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> ___
>>>> gmsh mailing list
>>>> gmsh@onelab.info
>>>> http://onelab.info/mailman/listinfo/gmsh
>>>
>>> --
>>> Prof. Christophe Geuzaine
>>> University of Liege, Electrical Engineering and Computer Science
>>> http://www.montefiore.ulg.ac.be/~geuzaine
>>>
>>> Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
>>> Free software: http://gmsh.info | http://getdp.info | http://onelab.info
>>>
>>>
>>>
>>
>>
>>
>>
>>
>> ___
>> gmsh mailing list
>> gmsh@onelab.info
>> http://onelab.info/mailman/listinfo/gmsh
>>
>
>
>
>
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] constant mesh despite background field
> On 05 Apr 2016, at 17:34, Andrea Cimatoribus > wrote: > > Hi everyone, > my first post here, I hope it is not too silly a question. > > I am trying to generate a gis-based 2D mesh for some large-scale geophysical > fluid dynamics model. My .geo file is attached, and has been generated by > this QGIS plugin: > https://github.com/ccorail/qgis-gmsh/wiki > > The essential part is this: > NF = newf; > Field[NF] = Structured; > Field[NF].TextFormat = 0; > Field[NF].FileName = "/home/cimatori/Data/Mesh/mesh.dat"; > Background Field = NF; > > As far as I understand, this should mean that the mesh follows the scales > given in the mesh.dat file. However, the mesh I get is almost entirely fixed > at the smallest scale (apart from some small region close to complex > geometry). Maybe the size prescribed on the boundary is smaller? In this case set Mesh.CharacteristicLengthExtendFromBoundary = 0; > I did check the content of this file, and it contains the numbers I expect to > find (range 30-300m). The data in mesh.dat seems to be completely ignored. > Could it be related to the fact that on all the boundaries the scale is > expected to be the same (30m)? But how to fix it? > > Thanks for your help. > > -- > > Andrea Cimatoribus > post-doctoral researcher > EPFL ENAC IIE ECOL > people.epfl.ch/andrea.cimatoribus > _______ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] CFD mesh with prismatic boundary layer based on STL
Hi Johannes - Here's the corrected file. wing_gmsh_cg.geo Description: Binary data > On 05 Apr 2016, at 09:50, Johannes Dillinger wrote: > > Dear all, > > still didn't have any luck in meshing the STL with prismatic sublayer.. > Either the tetrahedra mesh directly connects to the wing surface rather than > to the prisms, or the tetrahedra are not generated at all due to intersection > issues. > Maybe someone has another hint on how to modify Ruths file to make it work.. > Thanks a lot in advance, kind regards, > > Johannes > > > > > 2016-04-02 2:35 GMT+02:00 Johannes Dillinger : > Dear Ruth, > > intermdeiate result: > I tried changing line 54 from > > slbndwing=newsl; Surface Loop(newsl) = 1; > > to > > slbndwing=newsl; Surface Loop(newsl) = num[0]; > > but then a lot of warnings occur and no tetrahedra mesh outside the boundary > layer is generated.. As I understand, num[0] is the outer surface of the > boundary layer.. > > Maybe you have another idea..:) > Kind regards, > > Johannes > > > > > > > > 2016-04-01 19:15 GMT+02:00 Johannes Dillinger : > Dear Ruth, > thanks a lot for the modification! There is is just one problem left, the > tetrahedra mesh does not connect to the outer surface of the prism sublayer, > but directly to the wing surface; the meshes intersect (see picture attached). > I got a similar result with Delauney-3D-meshing. > Do you have a hint on how to solve this? > Thank you so much again!! > Kind regards, > Johannes > > > > > 2016-04-01 11:59 GMT+02:00 Ruth Vazquez Sabariego > : > Hope the attached patched version of your file, does what you want. > There was a small problem with the surface loops. > > Regards, > Ruth > > > > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] compile with med support
> On 06 Apr 2016, at 07:54, Bernd Hahnebach wrote:
>
> Zitat von Christophe Geuzaine :
>
>>
>>> On 05 Apr 2016, at 19:19, Bernd Hahnebach wrote:
>>>
>>> Found it. It seams a known problem on Debian Jessie (I do not know about
>>> other debian). See https://bugs.debian.org/cgi-bin/bugreport.cgi?bug=756472
>>>
>>> Makeing this changes to the CMakeLists.txt did the trick. HDF5 where found.
>>>
>>
>> Thanks - I've merge the Debian patch into our CMakeLists.txt.
>
> Where? https://onelab.info/svn/gmsh/trunk/CMakeLists.txt has not been changed
> !
>
Cf. https://onelab.info/trac/gmsh/timeline: I merged the patch, then reverted
it as it broke our nightly builds... (maybe this requires a too recent version
of cmake)
Doesn't the old version work if you specify the path where the hdf5 stuff is
installed using e.g.
cmake -DCMAKE_PREFIX_PATH=/path/to/hdf5/stuff ..
> Bernd
>
>
>
>>
>>
>>> cheers bernd
>>>
>>>
>>> Zitat von Bernd Hahnebach :
>>>
>>>> For med support on Debian Jessie package libmedc-dev has to be installed.
>>>> Watchout for the c !
>>>>
>>>> But the culprit is somewhere else. Cmake only searches for libmed if HDF5
>>>> was found. OK run cmake ../ again.
>>>> YEAH
>>>> -- HDF5_LIB-NOTFOUND
>>>> -- HDF5 not found
>>>>
>>>> ok I searched for them, and they seam installed. OK
>>>>
>>>> cmake \
>>>> -DHDF5_LIB=/usr/lib/x86_64-linux-gnu/libhdf5_cpp.so.8.0.2 \
>>>> ../
>>>>
>>>> BUT
>>>>
>>>> [ 1%] Building CXX object CMakeFiles/gmsh.dir/Common/CommandLine.cpp.o
>>>> In file included from
>>>> /home/hugo/Documents/dev/gmsh/gmsh-dev/Common/CommandLine.cpp:41:0:
>>>> /usr/include/med.h:22:18: fatal error: hdf5.h: Datei oder Verzeichnis
>>>> nicht gefunden
>>>> #include
>>>> ^
>>>> compilation terminated.
>>>>
>>>>
>>>> BUT
>>>>
>>>> $ locate hdf5.h
>>>> /usr/include/hdf5/serial/hdf5.h
>>>> /usr/include/vtk-6.1/vtk_hdf5.h
>>>>
>>>>
>>>> Which results in how do I set the include Dir for hdf5 ?
>>>>
>>>> Bernd
>>>>
>>>>
>>>> Zitat von Christophe Geuzaine :
>>>>
>>>>>
>>>>>> On 04 Apr 2016, at 17:49, Bernd Hahnebach wrote:
>>>>>>
>>>>>> According CMakeLists.txt is MED support activated by default
>>>>>>
>>>>>> opt(MED "Enable MED mesh and post file formats" ${DEFAULT})
>>>>>>
>>>>>>
>>>>>> But at my cmake output it is missing, but cmake does not complain some
>>>>>> library is missing ?!? How do I set the path to libmedc1 for cmake gmsh?
>>>>>>
>>>>>
>>>>> Bernd - Gmsh looks for a library named libmed (.a, .so, .dylib, etc.):
>>>>> not sure why yours is named libmedc1?
>>>>>
>>>>> You can modify the line
>>>>>
>>>>> find_library(MED_LIB med)
>>>>>
>>>>> into
>>>>>
>>>>> find_library(MED_LIB medc1)
>>>>>
>>>>> in the CMakeLists.txt file to see if this fixes it.
>>>>>
>>>>>> cheers bernd
>>>>>>
>>>>>>
>>>>>> cmake gmsh
>>>>>>
>>>>>> -- Gmsh 2.12.1 has been configured for Linux
>>>>>> --
>>>>>> -- * Build options: Ann Bamg Bfgs Blas(Generic) Blossom Cairo Chaco
>>>>>> DIntegration Dlopen Fltk GMP Gmm Jpeg Kbipack Lapack(Generic)
>>>>>> LinuxJoystick MathEx Mesh Metis Mmg3d Mpeg NativeFileChooser Netgen
>>>>>> ONELAB ONELABMetamodel OpenCascade OpenGL OptHom Parser Plugins Png Post
>>>>>> Salome Solver Taucs TetGen/BR Tetgen1.5 Voro3D Zlib
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> Zitat von Bernd Hahnebach :
>>>>>>
>>>>>>> Hi gmsh folks,
>>>>>>>
>>>>>>> just compiled gmsh on debian jessie. It worked like a charm. I have
>>>>>>> libmedc1 installed but gmsh com
Re: [Gmsh] OneLab VirtualBox or individual components?
Hi Todd - On Linux I recommend installing the ONELAB bundle (Gmsh+GetDP) from http://onelab.info and installing Elmer through your package manager. Christophe > On 06 Apr 2016, at 16:53, Todd Pierce wrote: > > Hi All, > > I like this idea that gmsh, elmer and getdp are all getting bundled together > into OneLab. That's really helpful for people like me who are armchair > scientists and barely know what we're doing. > > So, naturally, at some point I want to install OneLab. My question is, > should I really go for this VirtualBox installation where it's all one huge > thing or should I just install the components and get them working together > myself? > > I'm pretty good at Linux system administration so I don't think it would blow > my mind to configure them individually, but if, say, these software > components are constantly being upgraded or if OneLab is held together with > Scotch tape and dental floss, maybe the VirtualBox installation is better. > > I ask since I'd rather not install VirtualBox, but if I have to I'll do it. > > Any ideas? > > Thanks, > > -Todd > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] EssaiGmsh_v1.java:775: error: method add_physical in class WrapGmsh cannot be applied to given types;
> On 28 Apr 2016, at 19:36, Nico Schlömer wrote: > > Hi everyone, > > When compiling the nightly sources from Gmsh, I'm getting > ``` > «PKGBUILDDIR»/wrappers/java/WrappingJava/src/main/java/com/artenum/sample/EssaiGmsh_v1.java:775: > error: method add_physical in class WrapGmsh cannot be applied to given > types; > WrapGmsh.add_physical(str, lst1, m.getFileName()); > ^ > > ``` > and more such errors. Full build log at [1]. > I don't think anybody is actually maintaining or using the Java bindings. I propose to simply remove them: please react if you think this is a bad idea? > Cheers, > Nico > > PS: Is there a bug tracker for gmsh? > Yes: https://onelab.info/trac/gmsh/report/1 (login and passwd: gmsh) > > [1] > https://launchpadlibrarian.net/256925558/buildlog_ubuntu-trusty-amd64.gmsh_2.12.1~20160428182441-trusty1_BUILDING.txt.gz > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Importing .geo file to .stl file
> On 28 Apr 2016, at 12:03, Swathi Kanuku wrote: > > Hello sir, > I have created geometry using gmsh. Now I want to convert .geo file to > .stl file to simulate it in paraview. Can you please tell me how to convert > it. gmsh file.geo -2 -o file.stl > Thanks in advance > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] gmsh 2.12.0 : gmshpy : installation directory
> On 27 Apr 2016, at 11:02, Fabien Rozar wrote: > > Hello, > > I am a postdoc at LMGC in Montpellier with Frédéric Dubois as advisor. For > some user case of our software, LMGC90, we use the python module 'gmshpy'. > > In order to use it on a new computer, I need to install it. I downloaded the > source of gmsh 2.12.0 from the web page : > http://gmsh.info/ > > Before the compilation, I configured the generation of 'gmshpy' with : > $ cd build > $ cmake -DENABLE_WRAP_PYTHON=ON .../gmsh-2.12.0-source > $ make > > The module 'gmshpy' is well build under 'build/wrappers/gmshpy'. > > The point is that, after the installation step : > $ sudo make install > > the directory 'build/wrappers/gmshpy' is copied in '/usr/local' whereas > this directory is not part of python default search directories : > $ python -c "import sys; print '\n'.join(sys.path)" > /usr/lib/python2.7 > /usr/lib/python2.7/plat-x86_64-linux-gnu > /usr/lib/python2.7/lib-tk > /usr/lib/python2.7/lib-old > /usr/lib/python2.7/lib-dynload > /usr/local/lib/python2.7/dist-packages > /usr/lib/python2.7/dist-packages > /usr/lib/pymodules/python2.7 > > Question : Could you provide a release where the directory > 'build/wrappers/gmshpy' is installed (copied) > in one of the python default search directories? > Why not just cmake -DCMAKE_INSTALL_PREFIX=/where/you/want/to/install/stuff ? > (the directory '/usr/local/lib/python2.7/dist-packages' seems to be the more > suitable) > > Best regards, > ROZAR Fabien > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] 4 pole wound synchronous machine example
Dear Alexandre, It's indeed possible (the examples from the paper you are referring to were solved with GetDP), but there are currently no examples for multiharmonic simulations available online yet. It's on our TODO list... Christophe > On 28 Apr 2016, at 06:10, Alexander Shendi wrote: > > Hi folks, > > I have compiled gmsh and getdp for OpenBSD 4.9 > (amd64 and i386 platforms). Only small problems > there (I could resolve them). Now I have some > questions regarding using the gmsh/getdp combo > for the analysis of rotating electric machines. > > I am particularily interested in the 4 pole > salient pole synchronous machine in the example > "wfsm_4p.pro". In particular I would like to > calculate the rotor iron losses in the frequency > domain using the Harmonic Balance Finite Element > Method. The method is outlined in the paper > "Steady State Finite Element Analysis of a > Salient-Pole Synchronous Machine in the Frequency > Domain" paper by Gyselinck et al. , which is > referenced on the onelab.info web site. > > Now to my questions: > * Is it possible to use the above method with > gmsh/getdp? > * I would appreciate some guidance how to > write a *.pro file for the problem. Any > pointers welcome! > > Many thanks in advance for your help. > > Best Regards, > > Alexander > > _______ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Remeshing with GMesh in 2D
> On 27 Apr 2016, at 01:49, Dr. Stephan Schmidt > wrote: > > Dear all, > > I am trying to use gmsh to remesh a 2D geometry. > > So far, following the 3D examples (that is having a 2D stl surface as input > to generate a 3D mesh) work fabulously just following the examples on > compound surfaces, such as the one found here: > http://people.sc.fsu.edu/~jburkardt/examples/gmsh/t13.geo > > But for my problem, I would need to remesh a 2D geometry (that is having a > non-overlapping planar graph in a .msh file as input) seeking a 2D triangle > mesh as output with reparameterized boundary graphs. > > Is there any way of doing this as elegantly as the above link does for > surfaces? > Sure, you can use the exact same procedure... > Regards and many thanks, > Stephan Schmidt > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Interactive use, bad performance
> On 19 Apr 2016, at 22:29, R. Vestergaard wrote: > > Greetings! > > I am trying to visualize a large-ish mesh (~500MB), and the interface is > really slow. > When rotating the mesh, I seems gmsh maxes out one CPU core, while the > rest remain idle. Adding a Nvidia video card did not help speed things > up. > > Any general tips on optimizing performance for interactive use? > Remove the display of invisible things (e.g. volume edges/faces)? Gmsh does not simplify the mesh when displaying it (e.g. like Paraview): everything that's in the mesh is put in the OpenGL vertex arrays... even if it's not visible. > Regards, > Ruben Vestergaard > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Gmsh parallel delaunay algorithm (3d)
Dear Andrew, The results from the paper were obtained using a standalone prototype. This new parallel 3D algorithm is currently being integrated in Gmsh, but the work is not complete yet. We plan to make this available in Gmsh 3.0, hopefully this summer. Christophe > On 19 Apr 2016, at 12:24, Андрей Демченко wrote: > > Hello! > I am student from Belarus and I am interested in gmsh implementation of > parallel delaunay algorithm (using openmp) for my course project. I have read > this article (http://www.imr.sandia.gov/papers/imr24/01_IMR24_Remacle.pdf) > and wanted to test speedup of gmsh implementation comparing to tetgen 1.5 and > to my own implementation but faced some problems. > I have downloaded gmsh sources, extracted files related to 3d delaunay > triangulation (such as Mesh/delaunay3d.cpp, Mesh/delaunay3d_private.h etc.), > compiled it but recognized that there are some problems and checks (for > exampe under FIXME in delaunay3d_private.h) that doesn't allow to run that > code without modification for more than one thread using openmp. Am i right > that gmsh implementation is not ready for parallel processing yet? > Sorry for my bad English and thank you. > Best regards, Andrew Demchenko. > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] gmsh Digest, Vol 159, Issue 32
> On 29 Apr 2016, at 17:19, benjamin JEANTY-RUARD wrote: > > Dear Christophe, dear Nico, > > Artenum, my company, proposes the java wrapping in Gmsh as a > contribution about 5 years ago. Currently we use this java wrapping on > old Gmsh versions and we do not try it on the new Gmsh versions. > > I will try to test it during the month of May and if the problem can be > solved easily I will discuss about it with my boss to do it and I will > send you the correction. > > Tell me if it is ok on your side, Great! > > Regards, > > Benjamin Jeanty Ruard > > On Thu, 2016-04-28 at 21:51 +0200, > gmsh-requ...@ace20.montefiore.ulg.ac.be wrote: >> Message: 4 >> Date: Thu, 28 Apr 2016 21:48:43 +0200 >> From: Christophe Geuzaine >> To: Nico Schl?mer >> Cc: "g...@geuz.org" >> Subject: Re: [Gmsh] EssaiGmsh_v1.java:775: error: method add_physical >>in class WrapGmsh cannot be applied to given types; >> Message-ID: >> Content-Type: text/plain; charset=utf-8 >> >> >>> On 28 Apr 2016, at 19:36, Nico Schl?mer >> wrote: >>> >>> Hi everyone, >>> >>> When compiling the nightly sources from Gmsh, I'm getting >>> ``` >>> ?PKGBUILDDIR?/wrappers/java/WrappingJava/src/main/java/com/artenum/sample/EssaiGmsh_v1.java:775: >>> error: method add_physical in class WrapGmsh cannot be applied to given >>> types; >>>WrapGmsh.add_physical(str, lst1, m.getFileName()); >>>^ >>> >>> ``` >>> and more such errors. Full build log at [1]. >>> >> >> I don't think anybody is actually maintaining or using the Java >> bindings. I propose to simply remove them: please react if you think >> this is a bad idea? >> >>> Cheers, >>> Nico >>> >>> PS: Is there a bug tracker for gmsh? >>> >> >> Yes: https://onelab.info/trac/gmsh/report/1 (login and passwd: gmsh) >> >>> >>> [1] >> https://launchpadlibrarian.net/256925558/buildlog_ubuntu-trusty-amd64.gmsh_2.12.1~20160428182441-trusty1_BUILDING.txt.gz >>> ___ >>> gmsh mailing list >>> gmsh@onelab.info >>> http://onelab.info/mailman/listinfo/gmsh >> > > -- > > > Benjamin JEANTY-RUARD > > > Artenum Toulouse - Science & Groupware > http://www.artenum.com > > Bâtiment Calfocenter > 10, rue Marguerite Long > 31320 Castanet-Tolosan > France > Phone: +33 (0)5 82 95 51 97 > > > > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Remeshing with GMesh in 2D
> On 29 Apr 2016, at 17:57, Dr. Stephan Schmidt > wrote: > > Dear Prof. Geuzaine, > > many thanks for the feedback! I have tried to use exactly the same procedure, > but gmsh simply crashes without an error message when I try to e.g. refine > the remeshed mesh. > > Attached is a very minimal working example: „channel_orig.msh“ is the > original mesh created by gmsh itself. It appears to be valid in all regards, > my flow solver works flawlessly on this. „Remesh_channel.geo“ is the file I > actually open for remeshing and I manage to remesh the geometry successfully. > > But whenever I try to refine the rather coarse remeshed cannel using the gui, > gmsh crashes without throwing any error messages and I am stuck on how to > make a finer mesh. I am using the following build and version of gmsh: > Can you try with a recent nightly build? It works fine on my machine. > Info: --- > Info: Gmsh version : 2.11.0 > Info: Build OS : MacOSX > Info: Build options : Ann Bamg Bfgs Blas(Custom) Blossom Chaco > DIntegration Dlopen Fltk Gmm Jpeg(Fltk) Kbipack Lapack(Custom) MathEx Med > Mesh Metis Mmg3d Mpeg NativeFileChooser Netgen ONELAB ONELABMetamodel > OpenCascade OpenGL OptHom PETSc Parser Plugins Png(Fltk) Post SLEPc Salome > Solver Tetgen(1.5) Voro3D Zlib > Info: Build date : 20151107 > Info: Build host : Christophes-Mac.local > Info: Packager : geuzaine > Info: Executable : /Applications/Gmsh.app/Contents/MacOS/gmsh > Info: Home directory : /Users/schmidt/ > Info: Launch date: Fri Apr 29 17:47:49 2016 > Info: Command line : /Applications/Gmsh.app/Contents/MacOS/gmsh > Info: --- > > I would be very thankful if you could kindly provide feedback with respect to > this: > > - How can I generate a finer mesh? Ideally, I would love to carry over the > Delauney-weights associated with the original points into the remeshing > information. You can use any of the standard methods to define element sizes: global factor, fields (e.g. defined on a mesh), etc. > - Does the remeshed mesh also contain the same physical boundaries as the > original mesh or do I have to re-identify the different boundaries afterwards > manually? I know from the .stl 3D example how to do this when each physical > boundary is contained within its own Currently the "compounds" create new geometrical entities - you thus need to define the Physical entities using the new tags. However, we are working on a new implementation of the remeshing tools for Gmsh 3.0, which will change the workflow: the "Compound" command will become a meshing constraint command, and will not create new entities. Christophe > .stl file, but I would much appreciate hints on how to match a compound > quantity with a physical surface and consequently a boundary marker in the > case here, where everything stems from a single .msh file. > > Regards and many thanks, > Stephan Schmidt > > > > >> Am 28.04.2016 um 21:53 schrieb Christophe Geuzaine : >> >>> >>> On 27 Apr 2016, at 01:49, Dr. Stephan Schmidt >>> wrote: >>> >>> Dear all, >>> >>> I am trying to use gmsh to remesh a 2D geometry. >>> >>> So far, following the 3D examples (that is having a 2D stl surface as input >>> to generate a 3D mesh) work fabulously just following the examples on >>> compound surfaces, such as the one found here: >>> http://people.sc.fsu.edu/~jburkardt/examples/gmsh/t13.geo >>> >>> But for my problem, I would need to remesh a 2D geometry (that is having a >>> non-overlapping planar graph in a .msh file as input) seeking a 2D triangle >>> mesh as output with reparameterized boundary graphs. >>> >>> Is there any way of doing this as elegantly as the above link does for >>> surfaces? >>> >> >> Sure, you can use the exact same procedure... >> >> >>> Regards and many thanks, >>> Stephan Schmidt >>> ___ >>> gmsh mailing list >>> gmsh@onelab.info >>> http://onelab.info/mailman/listinfo/gmsh >> >> -- >> Prof. Christophe Geuzaine >> University of Liege, Electrical Engineering and Computer Science >> http://www.montefiore.ulg.ac.be/~geuzaine >> >> Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be >> Free software: http://gmsh.info | http://getdp.info | http://onelab.info > >
Re: [Gmsh] make a hexa second order mesh out of a hexa first order mesh
> On 13 May 2016, at 07:16, Bernd Hahnebach wrote: > > Found the culprit :-) but it does not solve my problem :-( > > If either "Mesh --> set Order 2" or "Mesh -> Hight order tools" is used on a > hexahedron8 it does not make a hexahedron20 but a hexahedron27 instead. Means > a node is not only made on any edge center point but on any face center point > too. > Use "Mesh.SecondOrderIncomplete = 1;" in your script (or "Tools->Options->Mesh->General : Use incomplete elements" in the GUI). > Neither unv nor inp export of GMSH support hexahedron27. A hexahedron20 can > be exported to unv and inp and these files are fine. > > This means back to the first question. Is it possible to make a hexahedron20 > element out of an hexahedron8 in gmsh? > > kind regards bernd > > attached a hexa20.unv and hexa8.unv to test with. > > > Zitat von Bernd Hahnebach : > >> More informations in this regard ... >> >> - open file b5__hexa8_coarse.msh >> - in GMSH Tree click on Mesh --> Set order 2 >> - Menue Tools --> Options --> Mesh --> Visibility --> activate Nodes >> - It seams we have got a second order mesh ?! >> >> - Menue Tools --> Visibility --> Numeric --> click on hide elements (put in >> * before) --> click on show elements (put in 50 before) >> - screen is attached --> it really seams to be a hexa20 mesh >> >> - save as msh --> reload --> still a hexa20 mesh >> - save as unv --> reload --> no hexa20 mesh, it is a hexa8 mesh >> >> Could it be hexa20 ist not supported in GMSH unv export ?? Could someone >> else give more informations in this regard? >> >> GMSH 2.12.0 >> >> kind regards bernd >> >> >> >> Zitat von Bernd Hahnebach : >> >>> Dear gmsh folks, >>> >>> attache a coarse hexa first order mesh. I'd like to get an second orde >>> hexa mesh out of this. I found the hight order tools can be used to >>> create the appropriate nodes for a second order hexa mesh. But how do I >>> create elements out of the nodes. >>> >>> kind regards bernd >>> >>> >> >> >> >> > > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Help
PS : from the log it seems that you are using a very old version of Gmsh (that uses ./configure, and not cmake!) > On 15 May 2016, at 09:39, David Colignon wrote: > > Hi, > > If you really need to compile Gmsh yourself (instead of just using the > available precompiled libraries), you will have to install Fltk , > http://www.fltk.org/index.php > > Regards, > > Dave > > -- > David Colignon, Ph.D. > 1er logisticien de recherche > Université de Liège > ACE - Applied & Computational Electromagnetics > Sart-Tilman B28 > 10, Grande Traverse > 4000 Liège - BELGIQUE > Tél: +32 (0)4 366 37 32 > http://www.ulg.ac.be/nic4 > > On 14/05/16 19:26, krish...@msu.edu wrote: >> When I compile gmsh, I got a following error >> checking for gcc... /usr/bin/gcc >> checking for C compiler default output... a.out >> checking whether the C compiler works... yes >> checking whether we are cross compiling... no >> checking for suffix of executables... >> checking for suffix of object files... o >> checking whether we are using the GNU C compiler... yes >> checking whether /usr/bin/gcc accepts -g... yes >> checking whether we are using the GNU C++ compiler... yes >> checking whether /usr/bin/g++ accepts -g... yes >> checking how to run the C preprocessor... /usr/bin/gcc -E >> checking for ranlib... ranlib >> checking for ar... /usr/bin/ar >> checking for main in -lm... yes >> checking for fltk-config... no >> configure: error: Could not find fltk-config. Try --with-fltk-prefix? >> please help me to fix the error >> thanks >> Research Scholar >> >> >> ___ >> gmsh mailing list >> gmsh@onelab.info >> http://onelab.info/mailman/listinfo/gmsh > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Segfault at partitionning / 2.12 freshly compiled
(Delaunay) > Delaunizing vertices… > Delaunay seconds: 0.527774 > Creating surface mesh … > Surface mesh seconds: 0.00 > Recovering boundaries… > Boundary recovery seconds: 0.004269 > Removing exterior tetrahedra … > Exterior tets removal seconds: 0.000612 > Recovering Delaunayness… > Delaunay recovery seconds: 0.146909 > Optimizing mesh… > Optimization seconds: 0.00224 > > Writing nodes. > Writing elements. > Writing faces. > Writing edges. > > Output seconds: 0.002485 > Total running seconds: 0.68894 > > Statistics: > > Input points: 1603 > Input facets: 3202 > Input segments: 4803 > Input holes: 0 > Input regions: 0 > > Mesh points: 1603 > Mesh tetrahedra: 4448 > Mesh faces: 10497 > Mesh edges: 7651 > Mesh faces on exterior boundary: 3202 > Mesh edges on exterior boundary: 4803 > Mesh faces on input facets: 3202 > Mesh edges on input segments: 4803 > > Info : 1603 points 7651 edges and 3202 faces in the initial mesh > Info : Found region 30 > Info : 0 points created - Worst tet radius is 10.7622 (PTS removed 0 0) > Info : 3D point insertion terminated (4111 points created): > Info : - 0 Delaunay cavities modified for star shapeness > Info : - 0 points could not be inserted > Info : - 20740 tetrahedra created in 0.213642 sec. (97078 tets/sec.) > Info : Done meshing 3D (0.972973 s) > Info : 4112 vertices 24141 elements > Info : Building graph… > [1] 57253 segmentation fault (core dumped) gmsh -3 ./demos/sphere.geo -part 3 > > -- > Docteur Ingénieur de recherche > CeMoSiS - vincent.hu...@cemosis.fr > Tel: +33 (0)3 68 85 02 06 > IRMA - 7, rue René Descartes > 67 000 Strasbourg > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Variable characteristic length over a surface
> On 26 May 2016, at 21:50, Karl Kästner wrote:
>
> Hello,
>
> I want to vary the element size over a surface. More specifically I want to
> geometrically increase the element size with increasing distance from the
> boundary. I tried to use the boundary layer option but this did not work, as
> the surface shape is nontrivial and the edge length has also to vary along
> the boundary.
>
> Therefore I first generate an initial mesh, then compute the local element
> size at the interior points of the initial mesh, add the interior points of
> the initial mesh to the geometry file and finally generate a second mesh.
>
> Unfortunately the points of the first mesh are also vertices of the second
> mesh, thus the element size of the second mesh cannot be lower than that of
> the first mesh. I tried to add the points to an attractor field, but gmsh
> displayed following warning and crashes:
>
> Error : Wrong mesh element size lc = 0 (lcmin = 0, lcmax = 1e+22)
>
(the attractor returns distances, which can here be equal to zero when you are
on the points)
Replace everything after
Plane Surface(1) = {1};
in your file with this:
// field 1 computes the distance to your boundary
Field[1] = Attractor;
Field[1].NNodesByEdge = 100;
Field[1].EdgesList = {1:32};
// field 2 processes this distance to create a size field from 0.01 to 0.1
Field[2] = Threshold;
Field[2].IField = 1;
Field[2].LcMin = 0.01;
Field[2].LcMax = 0.1;
Field[2].DistMin = 0.02;
Field[2].DistMax = 0.5;
Background Field = 2;
// don't extend the size field from the boundary
Mesh.CharacteristicLengthExtendFromBoundary = 0;
which will do what you want.
> Attached a trivial example reproducing the behaviour. Any help is appreciated.
>
> Kind regards,
>
> --
> Karl Kästner
>
> Phd Student at Wageningen University & Research centre
> Hydrology and Quantitative Water Management Group
> Droevendaalsesteeg 3, Lumen building, 6708 PB Wageningen, The Netherlands
> Phone: +31617876970 (Netherlands) +491627351075 (Germany)
> Mail: kastner.k...@gmail.com
> _______
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] gmsh crash
ng surface 321 (P"..., 48Info: Meshing
> surface 321
> (Plane, Delaunay)
>
> ) = 48
>
> write(1, "Info: Meshing surface 322 (P"..., 48Info: Meshing
> surface 322
> (Plane, Delaunay)
>
> ) = 48
>
> write(1, ""...,
> 47**
> *** NULL
>
> ) = 47
>
> write(2, "terminate called without an acti"..., 45terminate called
> without an ac
> tive exception
>
> ) = 45
>
> rt_sigprocmask(SIG_UNBLOCK, [ABRT], NULL, 8) = 0
>
> tgkill(15476, 15476, SIGABRT) = 0
>
> --- SIGABRT {si_signo=SIGABRT, si_code=SI_TKILL, si_pid=15476,
> si_uid=1000} ---
>
> +++ killed by SIGABRT (core dumped) +++
>
> Aborted (core dumped)
> ---
>
> Don't know if that helps any..
>
> - kv, Karen
>
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] #288: Invalid cells for a fine mesh around a sharp leading edge
Hi Xavier, There seem to be different issues with the file: - some curves are duplicated: when refining the curve visualization (Tools->Geometry->Aspect->Curve subdivisions) for example, curves 18 and 36 seem identical. This is likely a problem with the IGES representation: I would suggest using STEP (or even better the native BREP format since the model seems to come from OpenCASCADE). - the Characteristic Length command is only applied to some of these duplicate curves. Could you re-send with the BREP instead of the IGES ? Thanks, Christophe > On 30 May 2016, at 15:57, GARNAUD Xavier (SAFRAN) > wrote: > > Here is the attachement. > > -Message d'origine- > De : GARNAUD Xavier (SAFRAN) > Envoyé : lundi 30 mai 2016 15:57 > À : 'gmsh@onelab.info' > Objet : RE: [Gmsh] #288: Invalid cells for a fine mesh around a sharp leading > edge > > Here is an example that gives a similar error (it is not exactly the same but > I had to modify the geometry). > > In this case, replacing the cylinder by an edge attractor solves the problem. > In my other case it does not. > Best regards, > > Xavier > > -Message d'origine- > De : Gmsh [mailto:g...@ace20.montefiore.ulg.ac.be] > Envoyé : dimanche 29 mai 2016 09:52 > Objet : Re: [Gmsh] #288: Invalid cells for a fine mesh around a sharp leading > edge > > #288: Invalid cells for a fine mesh around a sharp leading edge > ---+--- > Reporter: gmsh| Owner: geuzaine >Type: defect | Status: accepted > Priority: minor | Resolution: > Keywords: | > ---+--- > Changes (by geuzaine): > > * owner: => geuzaine > * status: new => accepted > > > Comment: > > Could you attach a small example? > > -- > Ticket URL: <https://onelab.info/trac/gmsh/ticket/288#comment:1> > Gmsh <http://gmsh.info> > A three-dimensional finite element mesh generator with built-in pre- and > post-processing facilities > # > " Ce courriel et les documents qui lui sont joints peuvent contenir des > informations confidentielles, être soumis aux règlementations relatives au > contrôle des exportations ou ayant un caractère privé. S'ils ne vous sont pas > destinés, nous vous signalons qu'il est strictement interdit de les > divulguer, de les reproduire ou d'en utiliser de quelque manière que ce soit > le contenu. Toute exportation ou réexportation non autorisée est interdite Si > ce message vous a été transmis par erreur, merci d'en informer l'expéditeur > et de supprimer immédiatement de votre système informatique ce courriel ainsi > que tous les documents qui y sont attachés." > ** > " This e-mail and any attached documents may contain confidential or > proprietary information and may be subject to export control laws and > regulations. If you are not the intended recipient, you are notified that any > dissemination, copying of this e-mail and any attachments thereto or use of > their contents by any means whatsoever is strictly prohibited. Unauthorized > export or re-export is prohibited. If you have received this e-mail in error, > please advise the sender immediately and delete this e-mail and all attached > documents from your computer system." > # > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Help
> On 15 Jun 2016, at 22:58, Anusha Vs wrote:
>
> Hello
>
> My name is Anusha, and I am masters student studying in Germany. I have
> been using gmsh for quite a few months now for my academic work. I need some
> help regarding meshing where i have to mesh a surface using quadratic
> elements.,i.e. triangles having 6 nodes(3 vertices and 3 midpoints). It would
> be very helpful if you could guide me how I can generate such mesh.
>
Here we go:
Point(1) = {0, 0, 0};
Point(2) = {1, 0, 0} ;
Point(3) = {1, 1, 0} ;
Point(4) = {0, 1, 0} ;
Line(1) = {1,2} ;
Line(2) = {3,2} ;
Line(3) = {3,4} ;
Line(4) = {4,1} ;
Line Loop(5) = {4,1,-2,3} ;
Plane Surface(6) = {5} ;
Transfinite Line{1:4} = 10;
Transfinite Surface{6} Alternate;
Mesh.ElementOrder = 2;
Christophe
> NOTE::Please find attached, the picture of the mesh i wanted to generate
> using GMSH.
>
> --
> Thanks and Regards
> ANUSHA
> Department of Computer Science & Electrical
> Fachhochschule,
> Kiel,Germany.
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] [12.12.0] compiler error related to Metis
e/gmsh-2.12.0-hw4jsttw4y7lzozl5m4e3kz32fqbf6sw' > '-DCMAKE_BUILD_TYPE=RelWithDebInfo' '-DCMAKE_FIND_FRAMEWORK=LAST' > '-DCMAKE_INSTALL_RPATH_USE_LINK_PATH=FALSE' > '-DCMAKE_INSTALL_RPATH=/Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/gmsh-2.12.0-hw4jsttw4y7lzozl5m4e3kz32fqbf6sw/lib:/Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/gmsh-2.12.0-hw4jsttw4y7lzozl5m4e3kz32fqbf6sw/lib64:/Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/zlib-1.2.8-olur32o4bwaagk27ppelhlggws4pxknx/lib:/Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/petsc-3.6.4-24g5zsrebs3qba5uqt6l25chhnvwa5j7/lib:/Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/openmpi-1.10.3-5nfj7fkjw6wvypo5weuvso3evcw4kttg/lib:/Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/gmp-6.1.0-5zolcd74gboeemtykebegzqdsoipv5l2/lib:/Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/openblas-0.2.18-ligmia5aqpjb4hlunalo5xwfn5igmmx6/lib' > '-DENABLE_OS_SPECIFIC_INSTALL=OFF' > '-DGMSH_BIN=/Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/gmsh-2.12.0-hw4jsttw4y7lzozl5m4e3kz32fqbf6sw/bin' > > '-DGMSH_LIB=/Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/gmsh-2.12.0-hw4jsttw4y7lzozl5m4e3kz32fqbf6sw/lib' > > '-DGMSH_DOC=/Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/gmsh-2.12.0-hw4jsttw4y7lzozl5m4e3kz32fqbf6sw/share/doc' > > '-DGMSH_MAN=/Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/gmsh-2.12.0-hw4jsttw4y7lzozl5m4e3kz32fqbf6sw/man' > > '-DCMAKE_INSTALL_NAME_DIR:PATH=/Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/gmsh-2.12.0-hw4jsttw4y7lzozl5m4e3kz32fqbf6sw/lib' > '-DENABLE_OS_SPECIFIC_INSTALL=OFF' '-DENABLE_OCC=OFF' '-DENABLE_PETSC=ON' > '-DENABLE_SLEPC=OFF' '-DENABLE_BUILD_SHARED:BOOL=ON' > '-DENABLE_BUILD_DYNAMIC:BOOL=ON' '-DENABLE_MPI:BOOL=ON' > '-DENABLE_COMPRESSED_IO:BOOL=ON' > -- The CXX compiler identification is AppleClang 7.3.0.7030031 > -- The C compiler identification is AppleClang 7.3.0.7030031 > -- Check for working CXX compiler: > /Users/davydden/spack/lib/spack/env/clang/clang++ > -- Check for working CXX compiler: > /Users/davydden/spack/lib/spack/env/clang/clang++ -- works > -- Detecting CXX compiler ABI info > -- Detecting CXX compiler ABI info - done > -- Detecting CXX compile features > -- Detecting CXX compile features - done > -- Check for working C compiler: > /Users/davydden/spack/lib/spack/env/clang/clang > -- Check for working C compiler: > /Users/davydden/spack/lib/spack/env/clang/clang -- works > -- Detecting C compiler ABI info > -- Detecting C compiler ABI info - done > -- Detecting C compile features > -- Detecting C compile features - done > -- Looking for sys/types.h > -- Looking for sys/types.h - found > -- Looking for stdint.h > -- Looking for stdint.h - found > -- Looking for stddef.h > -- Looking for stddef.h - found > -- Check size of void* > -- Check size of void* - done > -- Found Blas(VecLib) > -- Found Lapack(VecLib) > -- Found Mesh > -- Found Solver > -- Found Post > -- Found Plugins > -- Found Parser > -- Found OpenGL: /System/Library/Frameworks/OpenGL.framework > -- Looking for XOpenDisplay in > /opt/X11/lib/libX11.dylib;/opt/X11/lib/libXext.dylib > -- Looking for XOpenDisplay in > /opt/X11/lib/libX11.dylib;/opt/X11/lib/libXext.dylib - found > -- Looking for gethostbyname > -- Looking for gethostbyname - found > -- Looking for connect > -- Looking for connect - found > -- Looking for remove > -- Looking for remove - found > -- Looking for shmat > -- Looking for shmat - found > -- Looking for IceConnectionNumber in ICE > -- Looking for IceConnectionNumber in ICE - found > -- Found X11: /opt/X11/lib/libX11.dylib > -- Could NOT find FLTK (missing: FLTK_LIBRARIES) > -- Found NativeFileChooser > -- Found ONELAB > -- Found ONELABMetamodel > -- System ANN not found: using contrib/ANN instead > -- Found Ann > -- Found Bfgs > -- Found DIntegration > -- Found OptHom > -- Found Kbipack > -- Found GMP > -- Found MathEx > -- Found MPI_C: > /Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/openmpi-1.10.3-5nfj7fkjw6wvypo5weuvso3evcw4kttg/lib/libmpi.dylib > > -- Found MPI_CXX: > /Users/davydden/spack/opt/spack/darwin-elcapitan-x86_64/clang-7.3.0-apple/openmpi-1.10.3-5nfj7fkjw6wvypo5weuvso3evcw4kttg/lib/libmpi_cxx.dylib;/Users/davydden/spack/opt/spack/darwin-elcapitan-
Re: [Gmsh] Translating points after extrusion
Hi Jose,
Extrude + Layers will not work if you modify the geometry : the extrusion
direction (translation, rotation, etc.) is used as-is when creating the
structured mesh. You will have to resort to the "Transfinite" command in that
case.
Christophe
> On 05 Jul 2016, at 17:06, Jose A. Abell M. wrote:
>
> Dear gmsh devs,
>
> Suppose I want to create a skewed hexahedral geometry to be meshed using a
> structured method. I want to use extrusions to create a regular cube first
> and then translating the points created by the extrusion to 'skew' it.
>
> Translating the extruded points and lines works, but the surfaces and volumes
> created by extrusion seem attached to their old geometry, generating errors.
>
> Please see example below:
>
> //Create 1x1x1 cube using extrusions
> Point(1) = {0, 0, 0, 1};
>
> lines[] = Extrude {1, 0, 0}{
> Point{1};
> Layers{10};
> Recombine;
> };
>
> surfaces[] = Extrude {0, 1, 0}{
> Line{lines[]};
> Layers{10};
> Recombine;
> };
>
> volumes[] = Extrude {0, 0, 1}{
> Surface{surfaces[]};
> Layers{10};
> Recombine;
> };
>
> // Distort the geometry using a translation for point on a face
>
> Translate {1, 0, 0} {
> Point{2, 1, 4, 3};
> }
>
>
>
>
> Saludos/Regards,
>
> --
> José A. Abell M., PhD
> Ingeniero Civil
> Profesor Investigador
> Facultad de Ingeniería y Ciencias Aplicadas
> Universidad de los Andes
> (+56)2-2412-9321
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Gmsh crash
Hi Martin, Give 2.13.1 a try: the issue is that compounds cannot currently be meshed high-order (this will be supported with the new compound implementation in Gmsh 3.0); also, the error message was misleading. Christophe > On 15 Jul 2016, at 20:45, Martin Kraska > wrote: > > Hi, > > when I try to create a compound surface in gmsh 2.13.0-svn-linux 64, I get > this error: > > Info: --- > Info: Gmsh version : 2.13.0-svn > Info: Build OS : Linux > Info: Build options : Ann Bamg Bfgs Blas(Generic) Blossom Chaco > DIntegration Dlopen Fltk Gmm Jpeg(Fltk) Kbipack Lapack(Generic) LinuxJoystick > MathEx Med Mesh Metis Mmg3d Mpeg Netgen ONELAB ONELABMetamodel OpenCascade > OpenGL OptHom PETSc Parser Plugins Png(Fltk) Post SLEPc Salome Solver Taucs > TetGen/BR Tetgen1.5 Voro3D Zlib > Info: Build date : 20160715 > Info: Build host : debian > Info: Packager : geuzaine > Info: Executable : > /home/user/A_Inst_gmsh/gmsh-2.13.0-svn-Linux/bin/gmsh > Info: Home directory : /home/user/ > Info: Launch date: Fri Jul 15 20:29:15 2016 > Info: Command line : > /home/user/A_Inst_gmsh/gmsh-2.13.0-svn-Linux/bin/gmsh part.geo > Info: --- > Info: Reading 'part.geo'... > Info: Reading 'part.step'... > Info: Done reading 'part.step' > Error : Unknown surface 23 > Warning : Could not make generatrices list for compound surface 30 > > When I try to mesh this, gmsh crashes. Without the compound surface I can > mesh the part without problems, see > > https://github.com/mkraska/CalculiX-Examples/tree/master/CAD/OnshapeTutorial > > What can I do to properly create a compound surface? > > Best regards, Martin Kraska > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Compound surfaces and lines
> On 17 Jul 2016, at 14:35, Martin Kraska
> wrote:
>
> Hi Christophe,
>
> Thanks for the hint. I can indeed create a mesh with compound surfaces in
> 2.13.1. I observed some restrictions:
>
> - the midside nodes of the generated elements don't follow the compound
> surface (except on the edges of the surface).
> - Successful meshing relies on introduction of an upper element size limit
> - If I also combine the adjacent lines in order to get rid of the short edges
> (they still attract small elements), then 3D meshing fails.
>
With the current compound implementation you will also need to define the
surfaces that touch a compound line as compound. Here's the corrected script:
Merge "part.step";
// Mesh control
Mesh.ElementOrder=2;
Mesh.CharacteristicLengthMax=10;
Mesh.Optimize=1;
// Display control
Mesh.SurfaceEdges = 1;
Mesh.SurfaceFaces = 0;
Mesh.VolumeEdges = 1;
Mesh.VolumeFaces = 0;
Compound Line(1000)={63,64};
Compound Line(1001)={65,66};
Compound Surface(30) = {23,25};
Compound Surface(31) = {9,24,26};
Compound Surface(32) = {21};
Compound Surface(33) = {22};
Mesh 3;
Physical Surface("support") = {5};
Physical Surface("load") = {17};
Physical Volume("part") = {1};
Mesh.SaveGroupsOfNodes = 1;
Save "gmsh.inp";
> The console complains about intersecting surfaces.
>
> Also, switching on all geometry healing options and setting a high geometry
> tolerance does not remove the small edges.
>
> Martin
>
>>
>> 1. Re: Gmsh crash (Christophe Geuzaine)
>>
>>
>> --
>>
>> Message: 1
>> Date: Sat, 16 Jul 2016 22:44:23 +0200
>> From: Christophe Geuzaine
>> To: gmsh@onelab.info
>> Subject: Re: [Gmsh] Gmsh crash
>> Message-ID: <42cdb7be-5537-4658-8de2-b90a2d3c4...@ulg.ac.be>
>> Content-Type: text/plain; charset=us-ascii
>>
>>
>> Hi Martin,
>>
>> Give 2.13.1 a try: the issue is that compounds cannot currently be meshed
>> high-order (this will be supported with the new compound implementation in
>> Gmsh
>> 3.0); also, the error message was misleading.
>>
>> Christophe
>>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Compound surfaces and lines
> On 17 Jul 2016, at 19:28, Martin Kraska
> wrote:
>
> Hi Christophe,
>
> Thanks again. The script even works without setting the max element size.
> y
> However, the problem of some of the midside nodes not obeying the geometric
> surface remains. Is there a way to get this right?
>
Not with the current implementation: currently compounds cannot be meshed with
curved elements. This will be fixed in the new implementation... hopefully by
the end of the summer.
> Martin
>
>>>
>>
>> With the current compound implementation you will also need to define the
>> surfaces that touch a compound line as compound. Here's the corrected script:
>>
>>
>> Merge "part.step";
>>
>> // Mesh control
>> Mesh.ElementOrder=2;
>> Mesh.CharacteristicLengthMax=10;
>> Mesh.Optimize=1;
>> // Display control
>> Mesh.SurfaceEdges = 1;
>> Mesh.SurfaceFaces = 0;
>> Mesh.VolumeEdges = 1;
>> Mesh.VolumeFaces = 0;
>>
>> Compound Line(1000)={63,64};
>> Compound Line(1001)={65,66};
>> Compound Surface(30) = {23,25};
>> Compound Surface(31) = {9,24,26};
>> Compound Surface(32) = {21};
>> Compound Surface(33) = {22};
>>
>> Mesh 3;
>> Physical Surface("support") = {5};
>> Physical Surface("load") = {17};
>> Physical Volume("part") = {1};
>> Mesh.SaveGroupsOfNodes = 1;
>> Save "gmsh.inp";
>>
>>
>>
>>> The console complains about intersecting surfaces.
>>>
>>> Also, switching on all geometry healing options and setting a high geometry
>>> tolerance does not remove the small edges.
>>>
>>> Martin
>>>
>>>>
>>>> 1. Re: Gmsh crash (Christophe Geuzaine)
>>>>
>>>>
>>>> --
>>>>
>>>> Message: 1
>>>> Date: Sat, 16 Jul 2016 22:44:23 +0200
>>>> From: Christophe Geuzaine
>>>> To: gmsh@onelab.info
>>>> Subject: Re: [Gmsh] Gmsh crash
>>>> Message-ID: <42cdb7be-5537-4658-8de2-b90a2d3c4...@ulg.ac.be>
>>>> Content-Type: text/plain; charset=us-ascii
>>>>
>>>>
>>>> Hi Martin,
>>>>
>>>> Give 2.13.1 a try: the issue is that compounds cannot currently be meshed
>>>> high-order (this will be supported with the new compound implementation in
>>>> Gmsh
>>>> 3.0); also, the error message was misleading.
>>>>
>>>> Christophe
>>>>
>>> ___
>>> gmsh mailing list
>>> gmsh@onelab.info
>>> http://onelab.info/mailman/listinfo/gmsh
>>
>> --
>> Prof. Christophe Geuzaine
>> University of Liege, Electrical Engineering and Computer Science
>> http://www.montefiore.ulg.ac.be/~geuzaine
>>
>> Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
>> Free software: http://gmsh.info | http://getdp.info | http://onelab.info
>
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] problem building Gmsh 2.13.2
Dear William - Thanks for the report; this should be fixed in SVN r23446. Christophe > On 19 Aug 2016, at 17:59, William F Mitchell > wrote: > > I am trying to update Gmsh to 2.13.2, but get a compiler error. My build > command is > > make -DCMAKE_PREFIX_PATH=/home/mitchell/apps/fltk > -DCMAKE_INSTALL_PREFIX=/home/mitchell/apps/gmsh .. > > I get the following error: > > [ 32%] Building CXX object CMakeFiles/gmsh.dir/Fltk/FlGui.cpp.o > In file included from > /home/mitchell/src/gmsh/gmsh-2.13.2-source/Fltk/FlGui.cpp:15:0: > /home/mitchell/src/gmsh/gmsh-2.13.2-source/Fltk/drawContextFltk.h: In member > function ‘virtual int drawContextFltk::getFontEnum(int)’: > /home/mitchell/src/gmsh/gmsh-2.13.2-source/Fltk/drawContextFltk.h:66:15: > error: ‘intptr_t’ was not declared in this scope > return (intptr_t)menu_font_names[index].user_data(); > ^ > /home/mitchell/src/gmsh/gmsh-2.13.2-source/Fltk/drawContextFltk.h:66:24: > error: expected ‘;’ before ‘menu_font_names’ > return (intptr_t)menu_font_names[index].user_data(); >^ > make[2]: *** [CMakeFiles/gmsh.dir/Fltk/FlGui.cpp.o] Error 1 > make[1]: *** [CMakeFiles/gmsh.dir/all] Error 2 > make: *** [all] Error 2 > > Thanks, > Bill > > -- > William F. Mitchell > Applied and Computational Mathematics Division > National Institute of Standards and Technology > william.mitch...@nist.gov http://math.nist.gov/~WMitchell > > _______ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Erroneous command line to call GetDP from Gmsh
Dear Paul - This is indeed quite strange. Could you try the latest (precompiled) ONELAB bundle on Linux and tell us if you can reproduce the problem? Christophe > On 17 Aug 2016, at 22:18, Paul Zander wrote: > > I am having the exact problem. I am running Debian. > > getdp and gmsh were installed from the packages for "testing". > gmsh version 2.12.0 > Build OS: Linux > Build date: 20160531 > Build host: binet > > One a different PC running Windows8, the models downloaded from onelab all > run well. > With Linux, I can't get beyond building the mesh. > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Gmsh error for high order element generation
> On 16 Aug 2016, at 16:34, acss vt wrote: > > Dear Gmsh developers, > I'm wondering if anyone has ever experienced program crash when generating > element order > 9. When I tried to generate 10th order hex element in Gmsh > using high-order tool, the program crashed. Based on user manual and other > Gmsh materials, it is supposed to support any order element. In theory yes :-) but indeed currently the Gmsh file format only defines elements up to order 9. > > Thanks for sharing any advice on that! > > Best, > Lu > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Can arcs be exactly pi radians?
> On 14 Aug 2016, at 16:49, Jim Monte wrote:
>
> Hi All,
>
> I am very new here. I need to make many circles. The documentation says in
> section 5.1.2
> "Circle ( expression ) = { expression, expression, expression };
> Creates a circle arc (strictly) smaller than Pi..."
>
> This would mean at least three perimeter points would be needed to make each
> circle. I tried with two points making arcs exactly pi radians each. They
> display properly, but I am wondering what problems I would have making all of
> the circles with two points (and a center point) instead of more. Thanks for
> any information.
The little internal Gmsh CAD engine is indeed limited to circle arcs < Pi. When
= Pi the specification of the plane in which the circle lives cannot be
computed in a univoque manner.
(It's not documented but you can actually specify the plane by its normal
vector explicitly as follows: Circle(num) = {p1,p2,p3} Plane {nx,ny,nz};)
Christophe
> Jim
>
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Proper way to make a volume from a line-loop
> On 10 Aug 2016, at 23:24, Brian Drawert wrote:
>
> Hi all,
> I am having trouble with something that is probably quite simple. I hope
> this community can help. I have a 1D shape, and I need to rotate it to make
> a 3D volume. I start by defining the points along the arc, connect them with
> their mirror image across the y/z=0 axis, then extrude by rotating. This
> gives me the proper shape, but not a volume mesh. Can anyone help?
>
The volume generated this way is "twisted" - extrude e.g. with an angle < Pi to
see what I mean.
You should build half the surface, and extrude that. Beware that you should
extrude by angles < Pi. This gives e.g.
p1 = newp;
Point(p1) = {0, 0, 0, 0.05};
p2 = newp;
Point(p2) = {0.25, 0.5, 0, 0.05};
p3 = newp;
Point(p3) = {0.5, 0.75, 0, 0.05};
p4 = newp;
Point(p4) = {1, 1, 0, 0.05};
p5 = newp;
Point(p5) = {2, 1, 0, 0.05};
p6 = newp;
Point(p6) = {2, 0, 0, 0.05};
l1 = newl;
Line(l1) = {p1, p2};
l2 = newl;
Line(l2) = {p2, p3};
l3 = newl;
Line(l3) = {p3, p4};
l4 = newl;
Line(l4) = {p4, p5};
l5 = newl;
Line(l5) = {p5, p6};
l6 = newl;
Line(l6) = {p1, p6};
ll7 = newll;
Line Loop(ll7) = {l2, l3, l4, l5, -l6, l1};
s1 = news;
Plane Surface(s1) = {ll7};
ex1[] = Extrude{{1,0,0}, {0,0,0}, 2*Pi/3}{Surface{s1};};
ex2[] = Extrude{{1,0,0}, {0,0,0}, 2*Pi/3}{Surface{ex1[0]};};
ex3[] = Extrude{{1,0,0}, {0,0,0}, 2*Pi/3}{Surface{ex2[0]};};
> Thanks in advance!
> Brian
>
>
> // Gmsh code here:
> p1 = newp;
> Point(p1) = {0, 0, 0, 0.05};
> p2 = newp;
> Point(p2) = {0.25, 0.5, 0, 0.05};
> p3 = newp;
> Point(p3) = {0.5, 0.75, 0, 0.05};
> p4 = newp;
> Point(p4) = {1, 1, 0, 0.05};
> p5 = newp;
> Point(p5) = {2, 1, 0, 0.05};
> p6 = newp;
> Point(p6) = {2, 0, 0, 0.05};
> p7 = newp;
> Point(p7) = {2, -1, 0, 0.05};
> p8 = newp;
> Point(p8) = {1, -1, 0, 0.05};
> p9 = newp;
> Point(p9) = {0.5, -0.75, 0, 0.05};
> p10 = newp;
> Point(p10) = {0.25, -0.5, 0, 0.05};
> l1 = newl;
> Line(l1) = {p1, p2};
> l2 = newl;
> Line(l2) = {p2, p3};
> l3 = newl;
> Line(l3) = {p3, p4};
> l4 = newl;
> Line(l4) = {p4, p5};
> l5 = newl;
> Line(l5) = {p5, p6};
> l6 = newl;
> Line(l6) = {p6, p7};
> l7 = newl;
> Line(l7) = {p7, p8};
> l8 = newl;
> Line(l8) = {p8, p9};
> l9 = newl;
> Line(l9) = {p9, p10};
> l10 = newl;
> Line(l10) = {p10, p1};
> ll1 = newll;
> Line Loop(ll1) = {l1,l2,l3,l4,l5,l6,l7,l8,l9,l10};
> surf1 = news;
> Plane Surface(surf1) = {ll1};
> ex1[] = Extrude{{0,0,0}, {1,0,0}, {0,0,0}, Pi}{Surface{surf1};};
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] too slow coherence function
> On 19 Aug 2016, at 21:34, andrew wrote:
>
> Hi,
>
> I try to create a 3d mesh from an extruded 2d mesh. The mesh is a simple box
> but with cuts in it of zero width. The geometry is created ok and the meshing
> even when I have 2000-3000 cuts is done in less than a few minutes. However
> when I convert the mesh for use with openfoam I was getting errors. By adding
> the command 'Coherence mesh;' after the creation of the mesh, the mesh is
> converted correctly for use with openfoam by using the ultility 'gmshToFoam'.
Indeed, it's a limitation of our extrusion algorithm with embedded curves
("Line In Surface"): it currently creates duplicate vertices. The real fix is
thus for us to enhance extrusion of embedded curves, but it's not trivial due
to the way extruded meshes are generated.
> The problem is in the use of the 'Coherence Mesh'. It takes almost two hours
> to execute while the meshing takes a few minutes. Is there a way to improve
> the Coherence command or to make the same geometry without the duplicate
> nodes that need so much time to be removed?
>
> I attached a file that needs 2.5 minutes on a quad core for meshing and more
> than ten minutes to remove 1600 duplicates.
> I have to create thousands of these meshes so a speedup would be welcomed.
>
That seems a bit slow: on my laptop meshing and duplicate removal take about
the same time (2 minutes each). Which version of Gmsh do you use?
> kind regards
>
> Andrew tsiantis
>
> _______
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] windows 7 text editor
> On 23 Aug 2016, at 14:44, Mark Smith wrote: > > Hi, > I set the text editor in gmsh_options to > General.TextEditor = """C:\Program Files (x86)\Notepad++\notepad++.exe" %s"; > // System command to launch a text editor > > Then when I launch gmsh it shows in tools>general>advancednotepad.exe %s > in the text editor command, > it completely ignores the setting in the gmsh_options file? > the location of gmsh_options is C:\Users\\AppData\Roaming/gmsh-options as > reported from "show file path" > What is wrong? Not sure; but the double quotes will clearly upset the parser. Could you try changing the setting in the GUI (Tools->Options->Advanced) with 'C:\Program Files (x86)\Notepad++\notepad++.exe' %s (note the single quotes around the executable), then File->Save options as default? > > If you hit restore all options to default settings where are the settings > obtained from is there a default_options file somewhere & if so where? Restore all options actually deletes that file. Christophe > > Thanks in advance > mark > > > _______ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] too slow coherence function
Andrew,
The problem actually seems to come from the way you define the surface: can you
regenerate your geometry script after
- removing all the zero-area line loops (e.g. "Line Loop(6) = {5, -5};")
- removing these loops from the definition of the main surface (i.e. simply
have "Plane Surface(1) = {205};")
With this you should not have duplicate vertices anymore, which should make
"Coherence Mesh" unnecessary.
Let us know,
Christophe
> On 21 Aug 2016, at 17:33, andrew wrote:
>
> hi,
>
> Thank you for the response. The file I submitted was a small one and the
> problem was not so obvious. I am attaching a 'real' file that I use on my
> simulations. It has 3000 cuts in it. This takes 2 hours on a quad core and 1+
> hour on a i7 @4ghz. I am using the version 2.13.1.
>
> Kind regards
>
> Andrew Tsiantis
>
>
> 2016-08-21 17:50 GMT+03:00 Christophe Geuzaine :
>
> > On 19 Aug 2016, at 21:34, andrew wrote:
> >
> > Hi,
> >
> > I try to create a 3d mesh from an extruded 2d mesh. The mesh is a simple
> > box but with cuts in it of zero width. The geometry is created ok and the
> > meshing even when I have 2000-3000 cuts is done in less than a few minutes.
> > However when I convert the mesh for use with openfoam I was getting errors.
> > By adding the command 'Coherence mesh;' after the creation of the mesh, the
> > mesh is converted correctly for use with openfoam by using the ultility
> > 'gmshToFoam'.
>
> Indeed, it's a limitation of our extrusion algorithm with embedded curves
> ("Line In Surface"): it currently creates duplicate vertices. The real fix is
> thus for us to enhance extrusion of embedded curves, but it's not trivial due
> to the way extruded meshes are generated.
>
> > The problem is in the use of the 'Coherence Mesh'. It takes almost two
> > hours to execute while the meshing takes a few minutes. Is there a way to
> > improve the Coherence command or to make the same geometry without the
> > duplicate nodes that need so much time to be removed?
> >
> > I attached a file that needs 2.5 minutes on a quad core for meshing and
> > more than ten minutes to remove 1600 duplicates.
> > I have to create thousands of these meshes so a speedup would be welcomed.
> >
>
> That seems a bit slow: on my laptop meshing and duplicate removal take about
> the same time (2 minutes each). Which version of Gmsh do you use?
>
> > kind regards
> >
> > Andrew tsiantis
> >
> > ___
> > gmsh mailing list
> > gmsh@onelab.info
> > http://onelab.info/mailman/listinfo/gmsh
>
> --
> Prof. Christophe Geuzaine
> University of Liege, Electrical Engineering and Computer Science
> http://www.montefiore.ulg.ac.be/~geuzaine
>
> Tetrahedron V, July 4-5 2016: http://tetrahedron.montefiore.ulg.ac.be
> Free software: http://gmsh.info | http://getdp.info | http://onelab.info
>
>
> <3k_cuts.zip>___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Close current mesh and open new
> On 26 Aug 2016, at 10:01, Praveen Kumar wrote: > > Hi, > > I am looking for a command in .geo file, to clear or close current mesh and > open new .msh file (using merge) > Use either "NewModel;" (keep the old model, add a new empty one), "Delete Model;" (destroys the current model) or "Delete All;" (destroys everything). > Any help on this would be great. > > Best Regards > -- > Praveen Kumar R | Technical Manager > Centre for Computational Technologies - CCTech > 1 Akshay Residancy | 50 Anand Park | Aundh | Pune | 411007 | India > Ph: +91 20 40098381/82 | M: +91-9923406816 | www.cctech.co.in > CCTech – "Transforming human life by democratization of technology" > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] too slow coherence function
> On 28 Aug 2016, at 13:57, andrew wrote:
>
> Hi,
>
> Thank you for your response.
>
> I have tried the suggestion you made but this way I don't get the internal
> surfaces which I need. The surface mesh is constructed correctly but as you
> can see from the first image I don't have the internal regions. I do get the
> points but openfoam does not understand them no matter what I do. In the
> second image is the mesh that I want with the internal surfaces.
>
> In the third image is a mesh I get with the loops but without the 'Line In
> Surface' set of commands where you can see that I get a bad mesh near the
> narrow areas where the mesh crosses the boundaries of the loop lines. I would
> expect that since I construct the surface with the loops, even if they were
> of zero width, they would be embedded in the surface since its creation
> making the need for the 'Line In Surface' commands obsolete.
>
I see; indeed this is a valid approach.
I have recoded the "Coherence Mesh;" function to make it much faster. Could you
download the next build and give it a try?
Thanks for the detailed feedback,
Christophe
> Anyway I found a workaround but it involves the Salome mesh generator.
>
> In the first step I create the mesh in Gmsh without the coherence function. I
> save it in unv format. I reload it in Salome and I remove the duplicate nodes
> there with 'merge nodes'. I delete the 'vol' since I don't need it and I save
> the mesh with only the surfaces I need again in unv format.
>
> It sounds cybersome but all these are done by scripts and Salome needs a few
> minutes to load a 3+ million cells mesh and save it back and less than a
> minute to remove the duplicates. Gmsh needs a few minutes for the mesh and
> 2+hrs for the removal so the procedure through Salome is way faster.
>
> Salome is written in python which is no faster than C++ so I guess the
> coherence routine in gmsh could take a refinement to a faster way.
>
> Kind regards
>
> Andrew Tsiantis
>
>
> 2016-08-28 12:10 GMT+03:00 Christophe Geuzaine :
>
> Andrew,
>
> The problem actually seems to come from the way you define the surface: can
> you regenerate your geometry script after
>
> - removing all the zero-area line loops (e.g. "Line Loop(6) = {5, -5};")
> - removing these loops from the definition of the main surface (i.e. simply
> have "Plane Surface(1) = {205};")
>
> With this you should not have duplicate vertices anymore, which should make
> "Coherence Mesh" unnecessary.
>
> Let us know,
>
> Christophe
>
>
> > On 21 Aug 2016, at 17:33, andrew wrote:
> >
> > hi,
> >
> > Thank you for the response. The file I submitted was a small one and the
> > problem was not so obvious. I am attaching a 'real' file that I use on my
> > simulations. It has 3000 cuts in it. This takes 2 hours on a quad core and
> > 1+ hour on a i7 @4ghz. I am using the version 2.13.1.
> >
> > Kind regards
> >
> > Andrew Tsiantis
> >
> >
> > 2016-08-21 17:50 GMT+03:00 Christophe Geuzaine :
> >
> > > On 19 Aug 2016, at 21:34, andrew wrote:
> > >
> > > Hi,
> > >
> > > I try to create a 3d mesh from an extruded 2d mesh. The mesh is a simple
> > > box but with cuts in it of zero width. The geometry is created ok and the
> > > meshing even when I have 2000-3000 cuts is done in less than a few
> > > minutes. However when I convert the mesh for use with openfoam I was
> > > getting errors. By adding the command 'Coherence mesh;' after the
> > > creation of the mesh, the mesh is converted correctly for use with
> > > openfoam by using the ultility 'gmshToFoam'.
> >
> > Indeed, it's a limitation of our extrusion algorithm with embedded curves
> > ("Line In Surface"): it currently creates duplicate vertices. The real fix
> > is thus for us to enhance extrusion of embedded curves, but it's not
> > trivial due to the way extruded meshes are generated.
> >
> > > The problem is in the use of the 'Coherence Mesh'. It takes almost two
> > > hours to execute while the meshing takes a few minutes. Is there a way to
> > > improve the Coherence command or to make the same geometry without the
> > > duplicate nodes that need so much time to be removed?
> > >
> > > I attached a file that needs 2.5 minutes on a quad core for meshing and
> > > more than ten minutes to remove 1600 duplicates.
> > > I have t
Re: [Gmsh] Gmsh may mistakenly treats an interface surface between two regions as a boundary surface
Hi Jun, It seems that your last volume definition is wrong. Attached is a file where the volume is defined interactively - does it work? Christophe pipe2.geo Description: Binary data > On 29 Aug 2016, at 21:22, Junhui Liu wrote: > > Hi All, > > I am generating structured tetra meshes (without the final recombine) inside > and outside a pipe. Please see the attache plot. Using structured tetra > meshes outside the pipe is to control the time step size. There are three > regions generated. Two are structured tetra meshes and another one is the > regular unstructured tetra mesh. Somehow our internal code complains that > "Boundary face construction failed with: xxx faces remaining”. if I use a > version that does not care about boundary construction, It appears that the > results of the cylindrical surface pointed by blue arrows behave very > strangely and look like those of a boundary face. If I take region III away, > the problem disappears. I have also attached the pipe.geo file. > > I have to say GMSH has lots of functionalities to control mesh and I really > like it. I hope this problem can be resolved soon. Thanks a lot for your help. > > Jun > > > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] set different order to different region of structure
Dear Nan, Normally this is not supported (Gmsh produces meshes of a given order globally). But in your specific case you can make it work by using a combination of "Mesh.MeshOnlyVisible" and "SetOrder": see attached file. Christophe mesh_hybrid_order.geo Description: Binary data > On 29 Aug 2016, at 19:35, Nan Li wrote: > > Hello everyone, > > I met a problem seems pretty simple but have bothered me for two days. > I have a structure constructed with beam and plate. So I want to use 1d > element (for the beam) and 2d element (for of plate) at the same time in > gmsh. As I have already programed in matlab an element 1d (order 1) with 2 > nodes and an element 2d with 8 nodes (order 2), I want to use them directly. > Then I find it difficult in gmsh to set different order to different region > of my structure. If I use Mesh.ElementOrder =2; Mesh.SecondOrderIncomplete = > 1, then the 1d element will have 3 nodes. If I delete these lines, the 2d > element will have 4 nodes. > > Does anyone have any idea on this question? > > Thank you very much. > > By the way, I use gmsh 2.8.5 and the file .geo is in the attachment. This is > only a simplified exemple of my structure. There is only three cells in this > exemple and I plan to have 50 cells in x direction and 30 in y direction. > > -- > Li Nan > Doctorante > Laboratoire Roberval > Université de Technologie de Compiègne, France > Email: nan...@utc.fr > _______ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] patching CMakeLists.txt for TAUCS
Dear Joseph,
The Taucs timing code uses ftime, which indeed is obsolete. I've changed our
Taucs config to disable timings - we don't use them anyway. The fix is in SVN.
Christophe
> On 01 Sep 2016, at 16:09, Joseph Mingrone wrote:
>
> Hello,
>
> There is a linking issue with the FreeBSD gmsh package when TAUCS is enabled.
> Here is a build snippet.
>
>
> CMakeFiles/gmsh.dir/contrib/taucs/G/taucs_timer.c.o: In function
> `taucs_wtime':
> /wrkdirs/usr/ports/cad/gmsh/work/gmsh-2.13.2-source/contrib/taucs/G/taucs_timer.c:(.text+0xd):
> undefined reference to `ftime'
>
>
> I believe adding '-lcompat' when linking gmsh solves the problem. How can
> cmake be instructed to do this?
>
> I tried patching CMakeLists.txt as below, but without success.
>
> Regards,
>
> Joseph
>
> --- CMakeLists.txt.orig 2016-08-31 20:17:30 UTC
> +++ CMakeLists.txt
> @@ -887,6 +887,10 @@ if(HAVE_SOLVER)
> else(WIN32)
> add_definitions(-DOSTYPE_linux)
> endif(WIN32 OR CYGWIN)
> + find_library(COMPAT_LIB compat)
> + if(COMPAT_LIB)
> +list(APPEND LINK_LIBRARIES ${COMPAT_LIB})
> + endif(COMPAT_LIB)
> file(GLOB_RECURSE TAUCS_D
> ${CMAKE_CURRENT_BINARY_DIR}/contrib/taucs/D/*.c)
> file(GLOB_RECURSE TAUCS_S
> ${CMAKE_CURRENT_BINARY_DIR}/contrib/taucs/S/*.c)
> file(GLOB_RECURSE TAUCS_Z
> ${CMAKE_CURRENT_BINARY_DIR}/contrib/taucs/Z/*.c)
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] How to find if a variable exists
> On 07 Sep 2016, at 22:30, Jim Monte wrote: > > Hi All, > > I have written macros returning arrays of various items, such as point IDs. > Depending on the arguments to the function, some arrays can be of different > sizes. To ensure that there is no "junk" at the end of the arrays left from a > previous call, I would like to Delete the arrays before building them. > However, the first time the macro is called, the array does not exist, and > that causes an error. One workaround is to force the array to exist by adding > something like x[0] = 0; before the Delete command, but I was wondering if > there is a more direct way to query if a specific variable exists. Thanks for > any information. > Exists(name) will do the trick > Jim Monte > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Question about early post increment
Jim - Indeed, good catch: our post-increment var++ was actually behaving like
++var. This is now fixed in SVN.
> On 10 Sep 2016, at 00:39, Jim Monte wrote:
>
> Hi All,
>
> I was trying to build arrays of point IDs and related things and used the
> post increment essentially as shown in the simple example below. In this
> example, the array element pid[0] is unassigned while pid[1] has a value,
> although I would have expected pid[0] to be the only assigned value. Am I
> misunderstanding how the increment operator is supposed to work in Gmsh?
>
> Jim Monte
>
> /*** start of example */
> i = 0; /* index to current point ID */
> j = newp; /* get a new point ID */
> pid[i++] = j; /* save point ID */
> Point(j) = {1, 2, 3}; /* assign value to point */
> xyz[] = Point{pid[0]}; /* retrieve coordinates of point */
> Printf("x=%g; y=%g; z=%g", xyz[0], xyz[1], xyz[2]); /* output */
> xyz[] = Point{pid[1]}; /* retrieve coordinates */
> Printf("x=%g; y=%g; z=%g", xyz[0], xyz[1], xyz[2]); /* output */
> /*** end of example */
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Questions related to Extrusions
> On 14 Sep 2016, at 20:44, Jim Monte wrote:
>
> Hi,
>
> I would like to know how to find the order of the "side" elements (lines or
> surfaces) that are generated during an extrusion. Thanks for any information.
When you extrude a line, the surface is created from a line loop defined by the
original line, the line extruded from the end point, the (reversed) "top" line
and (reversed) lines extruded from the starting point.
When you extrude a surface, the surface loop is created with the original
surface, the "top" surface, then all the surfaces extruded from the sides of
the original surface (in the order they appear in the surface definitions).
Note that these resulting entities can be changed if they are found to
duplicate existing entities, through the explicit or implicit call to
"Coherence;".
>
> Also, while doing some related testing, I saw some results that I could not
> explain from the script below.
>
> /*** start ***/
> Printf("next line ID = %g; next surface ID = %g", newl, news);
> Point(1) = {0, 0, 0};
> Point(2) = {0, 1, 0};
> Line(1) = {1, 2};
> Printf("next line ID = %g; next surface ID = %g", newl, news);
>
> out[] = Extrude {0, 0, 1} {Line{1}; };
>
> n = #out[];
> Printf("Extrude has returned %g elements", n);
> n -= 1;
> For i In {0 : n}
> Printf("Extrusion value[%g] = %g.", i, out[i]);
> EndFor
>
> Printf("next line ID = %g; next surface ID = %g", newl, news);
> /*** end ***/
>
> The output is
>
> 1 next line ID = 1; next surface ID = 1
> 2 next line ID = 2; next surface ID = 2
> 3 Extrude has returned 4 elements
> 4 Extrusion value[0] = 2.
> 5 Extrusion value[1] = 5.
> 6 Extrusion value[2] = 4.
> 7 Extrusion value[3] = -3.
> 8 next line ID = 6; next surface ID = 6
>
> Here are the things that did not seem right:
> • The value returned by news in line 2 is too high. No surface was
> created yet, so it should remain as 1.
> • The next free surface has a value of 1, so line 5 should give the
> extruded surface a value of 1 instead of 5.
This is for backward compatibility with old Gmsh versions: if you prefer a more
natural numbering, you can set "Geometry.OldNewReg=0;" at te beginning of your
script.
> • I thought the IDs should all be positive, so line 7 seems wrong.
Curves actually exist with both orientations in Gmsh; use "FAbs()" to get the
ID.
Christophe
> I would appreciate some explanations of what Gmsh is doing in these instances.
>
> Jim Monte
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] boundary layer disconnects at the corner of geometry(stl)
Dear Praveen,
The Boundary Layer field was an experiment: the 3D version was not functional
and has been removed in recent versions. The 2D version is fairly usable. New
implementation of boundary layers are forthcoming. In 3D, I currently advise
simply using the Extrude command: it will not create fans, but for simple
geometries like yours it should be sufficient.
Here's an example:
// create a square
lc = 0.3;
Point(1) = {0.0,0.0,0.0,lc};
Point(2) = {1,0.0,0.0,lc};
Point(3) = {1,1,0.0,lc};
Point(4) = {0,1,0.0,lc};
Line(1) = {4,3};
Line(2) = {3,2};
Line(3) = {2,1};
Line(4) = {1,4};
Line Loop(1) = {2,3,4,1};
Plane Surface(1) = {1};
// create a cube
Extrude {0,0.0,1} { Surface{1}; }
// remove the volume
Delete { Volume{1}; }
// extrude the boundary of the cube inwards by 0.05, with 5 layers of elements
Extrude { Surface{26, -1, 13, 17, 21, 25}; Layers{5, 0.05}; }
// create the inner volume
Surface Loop(159) = {92, 158, 70, 48, 136, 114};
Volume(160) = {159};
General.ExpertMode = 1; // disable warning about mixing Delaunay and extrusion
> On 14 Sep 2016, at 10:52, Praveen Kumar wrote:
>
> Hi,
>
> I am trying to generate boundary layers on all the surfaces of a simple
> rectangular duct geometry (imported stl). The prism layers are not getting
> connected on the edges/corners between each surfaces.
>
>
>
>
> Can someone help me understand what is missing here? Please Note, I am using
> 2.11.0 version of gmsh as this is not working at all in the latest version
>
> --
> Merge "duct.stl";
> CreateTopology;
>
> ll[] = Line "*";
> l = newl;
> For j In {0 : #ll[]-1}
> Compound Line(l+j) = ll[j];
> EndFor
>
> ss[] = Surface "*";
> s = news;
> For i In {0 : #ss[]-1}
> Compound Surface(s+i) = ss[i];
> EndFor
>
> bl = newl;
> bs = news;
>
> Field[1] = MathEval;
> Field[1].F = "0.1";
>
> Background Field = 1;
>
> Mesh.CharacteristicLengthExtendFromBoundary = 0;
>
> Physical Surface("Side-01") = {s};
> Physical Surface("Top") = {s+1};
> Physical Surface("Side-02") = {s+2};
> Physical Surface("Outlet") = {s+3};
> Physical Surface("Bottom") = {s+4};
> Physical Surface("Inlet") = {s+5};
>
> sl = newsl; Surface Loop(sl) = {s:s+5};
> vo = newv; Volume(vo) = {sl};
>
> //Physical Volume("V_Fluid") = {vo};
>
> pp[] = Point "*";
>
> Field[2] = BoundaryLayer;
> //Field[2].AnisoMax = 60;
> Field[2].FacesList = {s:bs-1};
> Field[2].EdgesList = {l:bl-1};
> Field[2].NodesList = {pp[0]:pp[7]};
> //Field[2].FanNodesList = {pp[0]:pp[7]};
> //Field[2].IntersectMetrics = 0;
> //Field[2].Quads = 0;
> //Field[2].hfar = 0.1;
> Field[2].hwall_n = 0.025;
> //Field[2].hwall_t = 0.1;
> Field[2].ratio = 1.1;
> Field[2].thickness = 0.05;
> BoundaryLayer Field = 2;
> //Background Field = 2;
>
> --
>
> --
> Praveen Kumar R | Technical Manager
> Centre for Computational Technologies - CCTech
> 1 Akshay Residancy | 50 Anand Park | Aundh | Pune | 411007 | India
> Ph: +91 20 40098381/82 | M: +91-9923406816 | www.cctech.co.in
> CCTech – "Transforming human life by democratization of technology"
>
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Questions related to Extrusions
> On 15 Sep 2016, at 20:47, Jim Monte wrote:
>
>
> On Thu, Sep 15, 2016 at 3:23 AM, Christophe Geuzaine
> wrote:
>
> > On 14 Sep 2016, at 20:44, Jim Monte wrote:
> >
> > Hi,
> >
> > I would like to know how to find the order of the "side" elements (lines or
> > surfaces) that are generated during an extrusion. Thanks for any
> > information.
>
> When you extrude a line, the surface is created from a line loop defined by
> the original line, the line extruded from the end point, the (reversed) "top"
> line and (reversed) lines extruded from the starting point.
>
> When you extrude a surface, the surface loop is created with the original
> surface, the "top" surface, then all the surfaces extruded from the sides of
> the original surface (in the order they appear in the surface definitions).
>
> Note that these resulting entities can be changed if they are found to
> duplicate existing entities, through the explicit or implicit call to
> "Coherence;".
>
> >
> > Also, while doing some related testing, I saw some results that I could not
> > explain from the script below.
> >
> > /*** start ***/
> > Printf("next line ID = %g; next surface ID = %g", newl, news);
> > Point(1) = {0, 0, 0};
> > Point(2) = {0, 1, 0};
> > Line(1) = {1, 2};
> > Printf("next line ID = %g; next surface ID = %g", newl, news);
> >
> > out[] = Extrude {0, 0, 1} {Line{1}; };
> >
> > n = #out[];
> > Printf("Extrude has returned %g elements", n);
> > n -= 1;
> > For i In {0 : n}
> > Printf("Extrusion value[%g] = %g.", i, out[i]);
> > EndFor
> >
> > Printf("next line ID = %g; next surface ID = %g", newl, news);
> > /*** end ***/
> >
> > The output is
> >
> > 1 next line ID = 1; next surface ID = 1
> > 2 next line ID = 2; next surface ID = 2
> > 3 Extrude has returned 4 elements
> > 4 Extrusion value[0] = 2.
> > 5 Extrusion value[1] = 5.
> > 6 Extrusion value[2] = 4.
> > 7 Extrusion value[3] = -3.
> > 8 next line ID = 6; next surface ID = 6
> >
> > Here are the things that did not seem right:
> > • The value returned by news in line 2 is too high. No surface was
> > created yet, so it should remain as 1.
> > • The next free surface has a value of 1, so line 5 should give the
> > extruded surface a value of 1 instead of 5.
>
> This is for backward compatibility with old Gmsh versions: if you prefer a
> more natural numbering, you can set "Geometry.OldNewReg=0;" at te beginning
> of your script.
>
> > • I thought the IDs should all be positive, so line 7 seems wrong.
>
> Curves actually exist with both orientations in Gmsh; use "FAbs()" to get the
> ID.
>
> Christophe
>
>
> > I would appreciate some explanations of what Gmsh is doing in these
> > instances.
> >
> > Jim Monte
> >
> > ___
> > gmsh mailing list
> > gmsh@onelab.info
> > http://onelab.info/mailman/listinfo/gmsh
>
> Thanks for the details on extrusion. But I am still not seeing how to
> determine which index will receive the ID for a particular surface. Perhaps I
> am trying to do something the hard way. My original interest in this question
> is that I want to know where specific surfaces will be put in the returned
> array so I can associate a physical surface with specific plane surfaces that
> get formed as part of the extrusion. Is there another approach to do this
> programatically (without looking at the GUI for information).
>
> I did more testing with a "brick" shape, made by two extrusions,
> line->square->brick. The Geometry.OldNewReg option did help to make the IDs
> closer to the expected values, but it is still confusing what is being done
> by the program. At a high level, the brick has 8 points, 12 lines, and 6
> plane surfaces. But there are
>
> 49-1-8 = 40 unused points,
> 21-1-12 = 8 unused lines (9, 12, 13, 15, 16, 18, 19, 20), and
> 8-1-6 = 1 unused surface (2)
> generated as part of the extrusion. Is there some way to query what will be
> made by the extrusion?
>
The extrusion works by recursively extruding entities of lower dimension, so
when you extrude a surface, you actually extrude the curves first; each curve
extrusion creates a surface but also new curves and new points. Some of these
are duplicates, and they are merged. Hence all the "unused" points, lines, etc.
You can retrieve the IDs programmatically: you ju
Re: [Gmsh] Instabilities of mesh generation in interactive and batch mode of gmsh
> On 20 Sep 2016, at 17:07, andreas.schus...@dlr.de wrote: > > Hi, > > > I’m using gmsh for creating a finite element model of a full aircraft model > for structural analysis. > For that purpose, I have written a python process which writes the input > geometry of the aircraft as .geo file and pass it to gmsh via a simple batch > process in Python. My problem is that when I run that process in batch mode, > gmsh is running and produces some files called “before.msh” or “after.msh” > again and again and seems not to converge the mesh generation. > If I start the GUI of gmsh with the input geometry and start the mesh > generation inside, gmsh finishes after around 4 seconds with a nice mesh! > The aircraft is mainly defined by regular 4 point areas which I recombine to > get quadrilateral elements. At some places due to the local aircraft > structure definition, there are some areas with 5 lines which produces some > small degenerated elements. > > I already tested different meshing algorithms and options directly written in > the geometry input file to overcome to problem, but the batch process of gmsh > doesn’t finish, though. > > Can anybody help me with this problem? Weird... The "before.msh"/"after.msh" files are created for debugging purposes in the triangle->quad recombination algorithm, when the "-saveall" option is used. (This will be removed in a future version.) What happens if you run Gmsh (without Python), directly from the command line, e.g. gmsh myfile.geo -2 > Thanks in advance! > > Best regards > Andreas > > —— > Deutsches Zentrum für Luft- und Raumfahrt e.V. (DLR) > Institut für Faserverbundleichtbau und Adaptronik | > > > _______ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] redundant nodes when using transfinite
> On 20 Sep 2016, at 14:53, Thomas Unfer wrote: > > Hi there, > > I am new to this list, so first how all greetings to everybody. > I am trying to mesh two adjacent boxes with hexaedra using transfinite (.geo > attached) plus a few surfaces to apply different boundary conditions in my > solver. > The meshing goes alright and the cells are (visually) conforming. > However in the .msh file the physical surfaces and volumes rely on different > point ids which makes it impossible to retrieve proper connectivity. I tryed > Coherence but it did no good. Does any one has any clue how to fix this? > Your additional surfaces (e.g. surface 26) are not topologically connected to the volumes: the quadrangles on the surfaces will thus be completely independent of the hexahedra in the volumes. With recent nightly builds you could use "Coherence Mesh;" to glue the mesh back together after the fact, i.e., end your script with Mesh 3; Coherence Mesh; Save "mymesh.msh"; Not super elegant, but that should do the trick... > Thanks for your help, > > Thomas > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Problem with periodic mesh around cylinder with extrusion
> On 19 Sep 2016, at 18:41, ashish bhole wrote: > > Dear All, > > I am facing problem with 'periodic surface' in gmsh. I have created 3D mesh > around cylinder using extrusion (.geo attached here). > When I use Periodic Surface feature of gmsh it do not write complete periodic > data in .mhs file. > Indeed, with "Layers" the mesh is generated in a structured way, and no periodic information is saved - the periodic mesh generator is completely bypassed. > Am I using this feature correctly, specially when I using along with > extrusion feature? > > I can generate periodic mesh inside cube without any problem. > > -- With Regards > Ashish > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Interpolation of polynomial solutions
> On 19 Sep 2016, at 18:53, GARNAUD Xavier (SAFRAN) > wrote: > > Dear gmsh developpers, > > I’d like to know if it is possible to interpolate a polynomial solution at > given locations through the python interface. Only the adapted (h-refined) piecewise linear approximation I'm afraid, using the plugins. > Best regards, > > Xavier > # > " Ce courriel et les documents qui lui sont joints peuvent contenir des > informations confidentielles, être soumis aux règlementations relatives au > contrôle des exportations ou ayant un caractère privé. S'ils ne vous sont pas > destinés, nous vous signalons qu'il est strictement interdit de les > divulguer, de les reproduire ou d'en utiliser de quelque manière que ce soit > le contenu. Toute exportation ou réexportation non autorisée est interdite.Si > ce message vous a été transmis par erreur, merci d'en informer l'expéditeur > et de supprimer immédiatement de votre système informatique ce courriel ainsi > que tous les documents qui y sont attachés." > ** > " This e-mail and any attached documents may contain confidential or > proprietary information and may be subject to export control laws and > regulations. If you are not the intended recipient, you are notified that any > dissemination, copying of this e-mail and any attachments thereto or use of > their contents by any means whatsoever is strictly prohibited. Unauthorized > export or re-export is prohibited. If you have received this e-mail in error, > please advise the sender immediately and delete this e-mail and all attached > documents from your computer system." > # > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Query
> On 13 Sep 2016, at 18:18, UMESH KUMAR SETH > wrote: > > Hello, > > I am a new user to GMSH. > I downloaded the current stable release (version 2.13.2) 64 bit Linux > version of GMSH, last week, from the gmsh website (http://gmsh.info/). > When I checked the About Gmsh option in Help menu, I found it says Version > 2.8.3. You probably have a path-related issue: if Help->About says "2.8.3" you are clearly runnning the old version. > The file I downloaded from the website is named "gmsh-2.13.2-Linux64.tgz". > Because of this I am not sure if I am using the latest version 2.13.2 . > > The problems I face : > > 1.) As soon as I launch gmsh from the command prompt I find these two lines, > always. > > libGL error: failed to open drm device: Permission denied > libGL error: failed to load driver: i965 > > My system is Linux 3.13.0-77(Ubuntu version 14.04. 4LTS) > > 2) gmsh window hangs too much, when I go for generating the mesh (2D or 3D) > even with simple cube of 3 nodes per line. > 3) Sometimes, When I change the number of nodes on lines in my script and > then Reload and try to mesh with 3D from GUI window of gmsh, the window > CRASHES abruptly. > 4) I see the error messages : > > gmsh: ../../src/XlibInt.c:245: _XPrivSyncFunction: Assertion > `dpy->synchandler == _XPrivSyncFunction' failed. > [h2-efd-seth:03831] *** Process received signal *** > [h2-efd-seth:03831] Signal: Aborted (6) > [h2-efd-seth:03831] Signal code: (-6) > *** Error in `gmsh': double free or corruption (fasttop): 0x7fdefc000960 > *** > Abandon (core dumped) > > > gmsh: ../../src/xcb_io.c:529: _XAllocID: Assertion `ret != inval_id' failed. > [h2-efd-seth:12728] *** Process received signal *** > [h2-efd-seth:12728] Signal: Aborted (6) > [h2-efd-seth:12728] Signal code: (-6) > gmsh: ../../src/xcb_io.c:529: _XAllocID: Assertion `ret != inval_id' failed. > Abandon (core dumped) > > > Could someone help me with this, please ? > > Many Thanks, > Regards, > > > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] reading stl and meshing
> On 19 Sep 2016, at 17:39, Martin Genet wrote: > > Dear all: > > Attaching a gmsh script + stl that used to work fine (loading the stl, > remeshing the surface and then meshing the volume), at least until 2.11 (see > test.2.11.0.log), but not anymore (see test-2.13.2.log). I tried with 2.12.0, > 2.13.0, 2.13.1, 2.13.2, but it just hangs after "Re-orient all 0 face normals > coherently". Any idea on what might be going on? Thanks a lot! > Indeed, that's a regression. We are currently rewriting the STL remesher from scratch - this will soon be ready in "beta" and should be much more robust than the old code. > Martin > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Possible bug when extruding relatively small surfaces
Dear Kate, This is "normal" : by default the tolerance to distinguish separate points/vertices is set to 1e-6 times the size of the model bounding box. Set e.g. Geometry.Tolerance = 1e-12; at the beginning of your script. Christophe > On 21 Sep 2016, at 14:26, look wrote: > > Hello all, > > I think I found a bug related to the extrusion (and recombination) of > relatively small surfaces. > > The problem occurs when you try to extrude very small 2-D cells into the > third dimension while having a huge domain, like it is often necessary for a > resolved boundary layer mesh for a 2-D airfoil simulation. > Gmsh then fails to recombine the extruded mesh into quadrangles. > > You can test this with the attached minimal example. It creates the mesh well > until you uncomment the last line. You can then see in the GUI > statistics-panel, that gmsh created triangles, which is obviously not > intended. > > Please see also the related discussion at cfd-online.com: > http://www.cfd-online.com/Forums/openfoam-solving/164600-gmshtofoam-fatal-error.html > > Best regards, > > Kate > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Interactive gmsh command line?
Hi Todd, Just put all the commands you want in a text file, and run Gmsh on it: either "gmsh file" to show the GUI, or "gmsh file -" to don't show the GUI. If you really want an interactive command line like gnuplot, have a look at gmsh/utils/solvers/c++/interactive.cpp in the Gmsh source code: it's an example on how you can build such an interactive system by communicating with a Gmsh server process. Christophe > On 3 Oct 2016, at 14:51, Todd Pierce wrote: > > Hi Crew, > > I am completely new to gmsh. I may have posted this question before, but > this time I actually have gmsh installed so this time I am prepared. > > I was hoping to use gmsh as a "back end" to a program I am writing, similar > to the way I utilize gnuplot. What I do there is simply start gnuplot and > hurl commands at it through a pipe and the graphic window magically opens up > with results. > > Now, I know gmsh is a totally different animal. However, if I could have a > command line to work with and throw commands at it to define the 3D structure > of something and perform meshing (among other things) this would open up many > doors. Of critical importance would be to see the output in a graphical > window, since my program has no graphical capabilities. > > Does this make sense to anybody? Has anybody done this? > > Thanks, > > -Todd > _______ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] aaaaaaaaaaaaaaargh
Hi Jeremy - I've removed both debug statements from SVN; this is the cost of
living on the bleeding edge ;-)
> On 3 Oct 2016, at 15:37, German Theler wrote:
>
> Hi all
>
> Following up with (still unresolved) hard-coded printfs,
> http://onelab.info/pipermail/gmsh/2016/010837.html
> I am now getting the message 'aaargh':
>
> gtheler@tom:~/codigos/3ros/gmsh$ svn blame Mesh/meshGRegion.cpp | grep
> -C3 aaargh
> 12845 geuzaine
> 12808remacle for(int j = 0; j < 3; j++){
> 22971ajohnen if (out.trifacelist[i * 3 + j] - 1 >=
> initialSize){
> 22971ajohnen printf("aaargh\n");
> 22971ajohnen //if(v[j]->onWhat()->dim() == 3){
> 22971ajohnen v[j]->onWhat()->mesh_vertices.erase
> 22971ajohnen
> (std::find(v[j]->onWhat()->mesh_vertices.begin(),
> gtheler@tom:~/codigos/3ros/gmsh$
>
>
> Any comment or solution is appreciated.
> Thanks!
>
> --
> Jeremy Theler
> www.seamplex.com
>
>
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Adding items to physical entities
Hi Alessandro - you can use e.g. "Physical Volume(1) += {1,2,3};" (or simply
use lists to construct what you need, and affect the list at the end).
Christophe
> On 30 Sep 2016, at 08:48, Alessandro Vicini
> wrote:
>
>
> Hullo everybody,
>
> is there a way to add elementary entities to existing physical groups? Thank
> you.
>
> A.
>
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] How to generate 'plugin crack' for post processing
Hi Xiaoxin,
Plugin(Crack) works on meshes: you thus need to have a mesh before you can run
it. For example, add "Mesh 3;" before calling the plugin.
Christophe
> On 29 Sep 2016, at 17:44, Xiaoxin Lu wrote:
>
> Hello,
>
> I am a PhD student using GMSH to generate crack faces. However, I could not
> find any example for this kind of post processing.
> My .geo file is attached below, I put 'plugin crack' command inside but it
> did not work, I was wondering how to generate a post process file for this
> command.
>
> Thank you for your help.
>
> best regards,
> Xiaoxin
> //start of the file
> lc1 = 3;
> lc2 = 3;
>
> // We proceed by defining some elementary entities describing a cube:
>
> Point(1) = {0,0,0,lc1};
> Point(2) = {40,0,0,lc1};
> Point(3) = {40,40,0,lc1};
> Point(4) = {0,40,0,lc1};
> Point(5) = {0,0,40,lc1};
> Point(6) = {40,0,40,lc1};
> Point(7) = {40,40,40,lc1};
> Point(8) = {0,40,40,lc1};
>
> Line(30)={1,2};
> Line(31)={2,3};
> Line(32)={3,4};
> Line(33)={4,1};
> Line(34)={5,6};
> Line(35)={6,7};
> Line(36)={7,8};
> Line(37)={8,5};
> Line(38)={1,5};
> Line(39)={2,6};
> Line(40)={3,7};
> Line(41)={4,8};
> Line Loop(42)={30,31,32,33};
> Plane Surface(14)={42};
> Line Loop(43)={-37,-36,-35,-34};
> Plane Surface(16)={43};
> Line Loop(44)={38,34,-39,-30};
> Plane Surface(18)={44};
> Line Loop(45)={-32,40,36,-41};
> Plane Surface(20)={45};
> Line Loop(46)={-33,41,37,-38};
> Plane Surface(22)={46};
> Line Loop(47)={39,35,-40,-31};
> Plane Surface(24)={47};
> vol1[]=Extrude{0,0,40}{Surface{14};};
> Delete {Volume{vol1[1]};}
> vol2[]=Extrude{0,40,0}{Surface{18};};
> Delete {Volume{vol2[1]};}
> vol3[]=Extrude{40,0,0}{Surface{22};};
> Delete {Volume{vol3[1]};}
> Physical Surface(1)={14};
> Physical Surface(2)={vol1[0]};
> Physical Surface(3)={18};
> Physical Surface(4)={vol2[0]};
> Physical Surface(5)={22};
> Physical Surface(6)={vol3[0]};
> Surface Loop(1)={14,vol2[0],vol3[0],18,22,vol1[0]};
> Volume(1)={1};
> Physical Volume(1)={1};
>
> //n=1; add a crack surface.
> Point(5000)={1.645000e+01,7.466800e+00,5.890600e+00,lc2};
> Point(5001)={1.728000e+01,2.112300e+01,1.204100e+01,lc2};
> Point(5002)={2.823700e+01,1.636900e+01,2.111600e+01,lc2};
> Point(5003)={2.740700e+01,2.713500e+00,1.496600e+01,lc2};
> Line(5004)={5000,5001};
> Line(5005)={5001,5002};
> Line(5006)={5002,5003};
> Line(5007)={5003,5000};
> Line Loop(5008)={5004,5005,5006,5007};
> Plane Surface(5009)={5008};
> Physical Surface(11)={5009};
> Surface{5009} In Volume{1};
>
> Plugin(Crack).Dimension = 2 ;
> Plugin(Crack).PhysicalGroup = 11 ;
> Plugin(Crack).OpenBoundaryPhysicalGroup = 0 ;
> Plugin(Crack).Run ;
> //end of the file
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Repeated ids with Duplicata
Alexander - Some of your translations are by (0,0,0), hence creating a copy of
the surface at the same location as the original. One of the two copies gets
removed automatically (with Geometry.AutoCoherence=1).
Christophe
> On 13 Oct 2016, at 21:23, Alexander Lindsay wrote:
>
> I'm getting repeated IDs using Duplicata. Here's the relevant excerpt from my
> script:
>
> Plane Surface(11) = {9, 10};
> Plane Surface(12) = {10};
> fuel_surfaces[] = {};
> moder_surfaces[] = {};
> For xtrans In {0:108:12}
> For ytrans In {0:108:12}
> new_f_surface = Translate {xtrans, ytrans, 0} {
> Duplicata { Surface{11}; }
> };
> fuel_surfaces += new_f_surface;
> new_m_surface = Translate {xtrans, ytrans, 0} {
> Duplicata { Surface{12}; }
> };
> moder_surfaces += new_m_surface;
> Printf("New f surface '%g'", new_f_surface);
> Printf("New m surface '%g'", new_m_surface);
> EndFor
> EndFor
>
> And here's the beginning of the output:
>
> New f surface '13'
> New m surface '13'
> New f surface '13'
> New m surface '22'
> New f surface '23'
> New m surface '32'
> New f surface '33'
> New m surface '42'
>
> How do I ensure that everything has a unique ID?
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Problems with version 2.14.0 - Part 2
It would help if you would send the actual files, and maybe describe what you
want to do? Do you want to mix .brep and .geo entities in the same model?
> On 12 Oct 2016, at 11:43, Alessandro Vicini
> wrote:
>
>
> This is another problem I have with version 2.14.0:
>
> 1) I open an existing .geo file and it looks fine, see attached geometry.jpg
> picture;
> 2) I interactively delete a few surfaces, obtaining geometry_2.jpg. If I
> look at the geo file, I have these lines at the end:
>
>Delete {
> Surface{122, 124, 126, 120};
>}
>
> 3) The problem is that, when I issue the “delete” command, I get these
> (apparent) error messages in the console:
>
> Error : Unknown surface 515812064
> Error : Unknown surface 515812632
> Error : Unknown surface 515826792
> Error : Unknown surface 515813200
> Error : Unknown surface 515812064
> Error : Unknown surface 515812632
> Error : Unknown surface 515826792
> Error : Unknown surface 515813200
>
> 4) I see the same kind of (apparent) errors also if I reload the .geo file
> (I say apparent because in the graphics window the geometry looks ok).
>
> Where is the problem? Thank you.
>
> Alessandro
>
>
>
>
> _______
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Command line options
> On 12 Oct 2016, at 15:10, Alessandro Vicini > wrote: > > > Is it possible from command line to load a .geo file and save it in geo > unrolled format? Thanks. > gmsh file.geo -0 > Alessandro > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Recherche information sur l'installation de GMSH (Linux & windows)...
Dear Gregory, There are no dependencies: you can simple copy the binary file (gmsh.exe) on the machines. Christophe > On 11 Oct 2016, at 16:44, GACE, Gregory (SOGETI) > wrote: > > Bonjour, > > Je suis en charge de l’installation des applications d’analyse au sein du > bureau d’études d’Airbus Defence & Space. > Je souhaite installer gmsh sous Linux et Windows et souhaite donc connaitre > les prérequis et la méthode d’installation. > Existe-t-il une documentation notamment pour réaliser des installations > scriptées sous Windows ? > En d’autres terme, est-il possible d’installer gmsh sous windows via des > lignes de commandes de sorte à pouvoir les automatiser ? > > Cordialement, > > Grégory Gace | Product Engineering Support Team | Sogeti High Tech > IM – TOIS – Information Management > Phone +33 (0)5 62 19 59 26 > gregory.gace.exter...@airbus.com > > > > > *** > Ce courriel (incluant ses eventuelles pieces jointes) peut contenir des > informations confidentielles et/ou protegees ou dont la diffusion est > restreinte. Si vous avez recu ce courriel par erreur, vous ne devez ni le > copier, ni l'utiliser, ni en divulguer le contenu a quiconque. Merci d'en > avertir immediatement l'expediteur et d'effacer ce courriel de votre systeme. > Airbus Defence and Space et les sociétés Airbus Group declinent toute > responsabilite en cas de corruption par virus, d'alteration ou de > falsification de ce courriel lors de sa transmission par voie electronique. > This email (including any attachments) may contain confidential and/or > privileged information or information otherwise protected from disclosure. If > you are not the intended recipient, please notify the sender immediately, do > not copy this message or any attachments and do not use it for any purpose or > disclose its content to any person, but delete this message and any > attachments from your system. Airbus Defence and Space and Airbus Group > companies disclaim any and all liability if this email transmission was virus > corrupted, altered or falsified. > - > Airbus Defence and Space SAS (393 341 516 RCS Versailles) - Capital: > 16.587.728 EUR - Siege social: 51-61 Route de Verneuil, 78130 Les Mureaux, > France > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Orient gui camera along axis
> On 14 Oct 2016, at 01:17, Alexander Lindsay wrote: > > Is there a way to orient the gui camera along an arbitrary axis, kind of like > with paraview? E.g. if I've rotated the geometry around, but now my clumsy > mouse skills won't allow me to perfectly orient along the z-axis? > Sure: cf. the "X", "Y", "Z" and "rotate 90 degrees" buttons on the bottom of the GUI. > Alex > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Problems with version 2.14.0 - Part 2
> On 14 Oct 2016, at 08:59, Alessandro Vicini
> wrote:
>
>
> Hullo Christophe,
>
> here are the .geo and step files I am using.
> In this case yes, I am starting from a .geo kind of geometry, and I want to
> merge a step file into it.
>
Here's what's happening:
1) Delete { Surface{122, 124, 126, 120}; } : you cannot delete these surfaces
- they are embedded in the volume. The error message is bad - this is now fixed
in SVN.
2) Disappearing entities when you merge the .step file in the .geo file : this
is "normal" in the current version, you need to explicitly call "SyncModel;"
before the "Merge", so that the .geo model gets imported right away in the
model database. Otherwise the .step gets imported with entity tags that can
clash with those used in the .geo mode. This explicit sync step is indeed a bit
annoying - the SVN version now does it automatically.
Christophe
> Regards,
>
> Alessandro
>
>
> -Messaggio originale-
> Da: Christophe Geuzaine [mailto:cgeuza...@ulg.ac.be]
> Inviato: giovedì 13 ottobre 2016 21.58
> A: Alessandro Vicini
> Cc: gmsh@onelab.info
> Oggetto: Re: [Gmsh] Problems with version 2.14.0 - Part 2
>
>
> It would help if you would send the actual files, and maybe describe what you
> want to do? Do you want to mix .brep and .geo entities in the same model?
>
>> On 12 Oct 2016, at 11:43, Alessandro Vicini
>> wrote:
>>
>>
>> This is another problem I have with version 2.14.0:
>>
>> 1) I open an existing .geo file and it looks fine, see attached
>> geometry.jpg picture;
>> 2) I interactively delete a few surfaces, obtaining geometry_2.jpg. If I
>> look at the geo file, I have these lines at the end:
>>
>> Delete {
>> Surface{122, 124, 126, 120};
>> }
>>
>> 3) The problem is that, when I issue the “delete” command, I get these
>> (apparent) error messages in the console:
>>
>> Error : Unknown surface 515812064
>> Error : Unknown surface 515812632
>> Error : Unknown surface 515826792
>> Error : Unknown surface 515813200
>> Error : Unknown surface 515812064
>> Error : Unknown surface 515812632
>> Error : Unknown surface 515826792
>> Error : Unknown surface 515813200
>>
>> 4) I see the same kind of (apparent) errors also if I reload the .geo file
>> (I say apparent because in the graphics window the geometry looks ok).
>>
>> Where is the problem? Thank you.
>>
>> Alessandro
>>
>>
>>
>>
>> ____
>> ___
>> gmsh mailing list
>> gmsh@onelab.info
>> http://onelab.info/mailman/listinfo/gmsh
>
> --
> Prof. Christophe Geuzaine
> University of Liege, Electrical Engineering and Computer Science
> http://www.montefiore.ulg.ac.be/~geuzaine
>
> Free software: http://gmsh.info | http://getdp.info | http://onelab.info
>
>
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Problems with structured meshing
> On 15 Oct 2016, at 19:02, Kushagra Vidyarthi
> wrote:
>
> Dear Sir/Madam,
>
> I am working on a CFD analysis around an ellipse. For this purpose, I have
> created a circular domain around the geometry, and I am trying to create a
> structured mesh around it, so that I can have some control over the boundary
> layer. My .geo file is below:
>
> // Gmsh project created on Fri Oct 14 10:40:40 2016
>
> //geometry creation
> //ellipse: major axis 0.075, minor axis 0.009375
> // 1st arc
> Point(1) = {-0.0375, 0, 0, 0.01};
> Point(2) = {0, 0, 0, 0.01}; //global center
> Point(3) = {-0.02, 0, 0, 0.1};
> Point(4) = {-0.026516504, 0.003314563, 0, 0.01};
> Ellipse(1) = {1, 2, 3, 4};
>
> //arc 2
> Point(5) = {0, 0.0046875, 0, 0.01};
> Point(6) = {0.02, 0, 0, 0.01};
> Ellipse(2) = {4, 2, 3, 5};
>
> //arc3
> Point(7) = {0.026516504, 0.003314563, 0, 0.01};
> Ellipse(3) = {5, 2, 6, 7};
>
> //arc4
> Point(8) = {0.0375, 0, 0, 0.01};
> Ellipse(4) = {7, 2, 6, 8};
>
> //arc5
> Point(9) = {0.026516504, -0.003314563, 0, 0.01};
> Ellipse(5) = {8, 2, 6, 9};
>
> //arc6
> Point(10) = {0, -0.0046875, 0, 0.01};
> Ellipse(6) = {9, 2, 6, 10};
>
> //arc7
> Point(11) = {-0.026516504, -0.003314563, 0, 0.01};
> Ellipse(7) = {10, 2, 3, 11};
>
> //arc8
> Ellipse(8) = {11, 2, 3, 1};
>
> //circle
> Point(12) = {-1.7677669, -1.7677669, 0, 0.01};
> Point(13) = {-1.7677669, 1.7677669, 0, 0.01};
> Point(14) = {1.7677669, 1.7677669, 0, 0.01};
> Point(15) = {1.7677669, -1.7677669, 0, 0.01};
>
> Circle(9) = {12, 2, 13};
> Circle(10) = {13, 2, 14};
> Circle(11) = {14, 2, 15};
> Circle(12) = {15, 2, 12};
>
> //blocking lines
> Line(13) = {11, 12};
> Line(14) = {4, 13};
> Line(15) = {7, 14};
> Line(16) = {9, 15};
>
> //splitting lines
> Transfinite Line {9} = 100 Using Progression 1;
> Transfinite Line {10} = 100 Using Progression 1;
> Transfinite Line {11} = 100 Using Progression 1;
> Transfinite Line {12} = 100 Using Progression 1;
> Transfinite Line {14} = 100 Using Progression 1;
> Transfinite Line {15} = 100 Using Progression 1;
> Transfinite Line {16} = 100 Using Progression 1;
> Transfinite Line {13} = 100 Using Progression 1;
> Transfinite Line {8} = 100 Using Progression 1;
> Transfinite Line {1} = 100 Using Progression 1;
> Transfinite Line {2} = 100 Using Progression 1;
> Transfinite Line {3} = 100 Using Progression 1;
> Transfinite Line {4} = 100 Using Progression 1;
> Transfinite Line {5} = 100 Using Progression 1;
> Transfinite Line {6} = 100 Using Progression 1;
> Transfinite Line {7} = 100 Using Progression 1;
> Line Loop(17) = {13, 9, -14, -1, -8};
> Plane Surface(1) = {17};
> Line Loop(19) = {15, -10, -14, 2, 3};
> Plane Surface(2) = {19};
> Line Loop(21) = {16, -11, -15, 4, 5};
> Plane Surface(3) = {21};
> Line Loop(23) = {13, -12, -16, 6, 7};
> Plane Surface(4) = {23};
> Transfinite Surface {1} = {4, 13, 12, 11};
> Transfinite Surface {2} = {7, 14, 13, 4};
> Transfinite Surface {3} = {9, 15, 14, 7};
> Transfinite Surface {4} = {11, 12, 15, 9};
>
> Recombine Surface {1};
> Recombine Surface {2};
> Recombine Surface {3};
> Recombine Surface {4};
>
> When I run this file in gmsh, I get errors like: "Error : Surface 4 cannot
> be meshed using the transfinite algo". I also get a warning stating: "Warning
> : Cannot apply Blosson: odd number of triangles (15167) in surface 4". My
> error log and .geo files are attached. Please advise on what I can do to fix
> these errors.
>
Transfinite Surface meshes require matching number of vertices on opposite
sides. Add e.g. this to your file:
N = 10;
Transfinite Line {1:8} = N;
Transfinite Line {9:12} = 2*N-1;
Transfinite Line {13:16} = 50 Using Progression 1.2;
> Thanks and Regards,
>
> Kushagra Vidyarthi
>
> (copy).geo>___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] pos high order interpolation error
> On 18 Oct 2016, at 16:21, Олег Рябков wrote: > > Dear mesh users and developers. > I got very strange error with attached pos file. It contains view with just 2 > elements and when i open it with gmsh and choose 'Adapt visualization grid' > check button gmsh crashes. The strange thing is that this 2 elements could be > visualized separatly (i,e, when i delete from view either one of them it > works). > Indeed, this is current limitation of the high-order visualization code: it cannot handle hybrid meshes (i.e. meshes with more than 1 element type). > Gmsh - 2.14.0 x64, Win 7. > > -- > С уважением, Олег > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] [Gmsh-announce] New Gmsh release: 2.14.1
> On 2 Nov 2016, at 09:24, Florian Blachère > wrote: > > On 30/10/16 09:49, Christophe Geuzaine wrote: >> New in 2.14.1: fixed regression in periodic meshes; small bug fixes and code >> cleanups. >> >> Downloads, mailing lists, etc.: http://gmsh.info >> >> >> ___ >> gmsh-announce mailing list >> gmsh-annou...@onelab.info >> http://onelab.info/mailman/listinfo/gmsh-announce > Hello, > > It seems that the 2.14.1 archive (from > http://gmsh.info/src/gmsh-2.14.1-source.tgz) is updated every day as the > nightly 'gmsh-svn-source.tgz' archive. For instance, today the folder > inside the archived is named 'gmsh-2.14.1-svn-20161102-source'. > > Is it an expected behaviour ? Hi Florian - No, it was a bug due a bad symlink. Thanks for letting us know: it should now be fixed. > > Regards, > Florian Blachère > > _______ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Question about gmsh
> On 22 Oct 2016, at 04:09, Mayank Jog wrote: > > Hello, > I have been using "simnibs" software which uses gmsh to calculate electric > fields over a 3D mesh/geometry (human head). > It outputs a final .msh file, that contains the electric field information, > and can be viewed using gmsh. > > My issue is, I want to convert this file to a list where each line = > [coord location, electric field value ] How do I go about doing it? > > The ultimate goal is to regrid this data onto a regular cartesian grid to do > some funky postprocessing (fourier transforms, followed by funky stuff). Here's an example: Merge "tutorial/view3.pos"; Plugin(CutBox).X0=0; Plugin(CutBox).Y0=0; Plugin(CutBox).Z0=0; Plugin(CutBox).X1=1; Plugin(CutBox).Y1=0; Plugin(CutBox).Z1=0; Plugin(CutBox).X2=0; Plugin(CutBox).Y2=1; Plugin(CutBox).Z2=0; Plugin(CutBox).X3=0; Plugin(CutBox).Y3=0; Plugin(CutBox).Z3=1; Plugin(CutBox).NumPointsU=20; Plugin(CutBox).NumPointsV=20; Plugin(CutBox).NumPointsW=20; Plugin(CutBox).ConnectPoints=0; Plugin(CutBox).Boundary=0; Plugin(CutBox).View=0; Plugin(CutBox).Run; Save View[1] "points.txt"; > I tried saving the .msh file as a .txt ** , and then reading it line by > line.I was thinking of replacing each triangular element with its > centroid. However, I don't know how to proceed with other elementsplease > help me out > > Thank you! > Mayank > > ** "Post processing Generic txt" option in "save as" in gmsh > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Problem about 'Plugin Crack'
> On 25 Oct 2016, at 16:12, Xiaoxin Lu wrote:
>
> Dear GMSH users and developers,
>
> I have a question about 'Plugin crack' command.
> In the manual, it is said that 'The elements touching the crack on the
> “negative” side are modified to use the newly generated vertices'.
> I was wondering how to define the 'negative' side for 2D crack. Is it defined
> by NormalX, Normal Y and Normal Z, which are the components of normal vector
> of the crack surface?
No, Normal{X,Y,Z} is only used to explicitly provide for "1D cracks" (curves)
the normal to the surface in which these cracks live. For 2D cracks the normal
is computed as usual, based on the orientation of the triangles/quadrangles.
> Thank you very much.
>
> Best ragards,
> Xiaoxin
>
> ___
> gmsh mailing list
> gmsh@onelab.info
> http://onelab.info/mailman/listinfo/gmsh
--
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine
Free software: http://gmsh.info | http://getdp.info | http://onelab.info
___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] 2.14.1 tarball change
> On 3 Nov 2016, at 06:09, Joseph Mingrone wrote: > > Hello, > > The FreeBSD package for gmsh is currently failing, because the tarball for > this release has changed. > > fetch: http://gmsh.info/src/gmsh-2.14.1-source.tgz: size mismatch: expected > 10862009, actual 10860112 > *** Error code 1 > > Were there changes other than the renaming of the top level directory? I > don't have an old tarball to do a diff. > The changes were minor (top level directory name change + some one-liners to fix compilation warnings). You should use the tarball currently available on the server as the reference, as the original one was lost. > Thanks, > > Joseph > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Question about gmsh
> On 3 Nov 2016, at 02:01, Mayank Jog wrote: > > Hello Jeremy and Christophe, > Thanks for the help! > > Jeremy: > The data is Current-density magnitude at nodes (attached Fig1.png). I need to > recast this onto a uniformly-spaced 3D cartesian grid (nodes are > non-uniformly spaced in 3D). > > One hack I tried (please correct me :) ), was to export the .msh file as > ascii, then interpolate the data at each element onto the regular cartesian > grid. The coordinates of each element I set as the centroid of the nodes > comprising the element. > I'm new to handling mesh datadoes that sound like a fair approximation? > > Christophe: > I tried what you suggested in the GUI. I didnt understand the first line > (Merge "tutorial/view3.pos";) , but I carried out all the rest steps. I got a > new "view", but it only had a couple of voxels. I'm enclosing Fig1.png which > shows the data, and Fig2.png which shows what I get when I carried out the > step. > You should change the points (X0, Y0, ...) to match your geometry... Interactively, use Tools->Plugins->CutBox. > , > Mayank > > > > > > > > > > > > > > > On Wed, Nov 2, 2016 at 5:24 AM, Christophe Geuzaine > wrote: > > > On 22 Oct 2016, at 04:09, Mayank Jog wrote: > > > > Hello, > > I have been using "simnibs" software which uses gmsh to calculate electric > > fields over a 3D mesh/geometry (human head). > > It outputs a final .msh file, that contains the electric field information, > > and can be viewed using gmsh. > > > > My issue is, I want to convert this file to a list where each line = > > [coord location, electric field value ] How do I go about doing it? > > > > The ultimate goal is to regrid this data onto a regular cartesian grid to > > do some funky postprocessing (fourier transforms, followed by funky stuff). > > Here's an example: > > Merge "tutorial/view3.pos"; > Plugin(CutBox).X0=0; > Plugin(CutBox).Y0=0; > Plugin(CutBox).Z0=0; > Plugin(CutBox).X1=1; > Plugin(CutBox).Y1=0; > Plugin(CutBox).Z1=0; > Plugin(CutBox).X2=0; > Plugin(CutBox).Y2=1; > Plugin(CutBox).Z2=0; > Plugin(CutBox).X3=0; > Plugin(CutBox).Y3=0; > Plugin(CutBox).Z3=1; > Plugin(CutBox).NumPointsU=20; > Plugin(CutBox).NumPointsV=20; > Plugin(CutBox).NumPointsW=20; > Plugin(CutBox).ConnectPoints=0; > Plugin(CutBox).Boundary=0; > Plugin(CutBox).View=0; > Plugin(CutBox).Run; > Save View[1] "points.txt"; > > > > > > I tried saving the .msh file as a .txt ** , and then reading it line by > > line.I was thinking of replacing each triangular element with its > > centroid. However, I don't know how to proceed with other > > elementsplease help me out > > > > Thank you! > > Mayank > > > > ** "Post processing Generic txt" option in "save as" in gmsh > > ___ > > gmsh mailing list > > gmsh@onelab.info > > http://onelab.info/mailman/listinfo/gmsh > > -- > Prof. Christophe Geuzaine > University of Liege, Electrical Engineering and Computer Science > http://www.montefiore.ulg.ac.be/~geuzaine > > Free software: http://gmsh.info | http://getdp.info | http://onelab.info > > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
Re: [Gmsh] Surface In Volume does not work GMSH for Linux
> On 3 Nov 2016, at 14:26, bill wrote: > > Hi, > But it does, strangely enough, for windows! What is wrong? Is this a > bug or do I have a setting wrong? > "Mesh.Algorithm3D = 4;" selects the 3D "Frontal" algorithm, which does not support embedded surfaces. Reset it to the default "Delaunay" (1). > Geometry file and settings dump is attached. > > Cheers, > Bill Slade > > ___ > gmsh mailing list > gmsh@onelab.info > http://onelab.info/mailman/listinfo/gmsh -- Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info ___ gmsh mailing list gmsh@onelab.info http://onelab.info/mailman/listinfo/gmsh
