Re: [Gmsh] Issues with gmsh GUI after updating to 3.0.7

2018-02-12 Thread Christophe Geuzaine 


> On 12 Feb 2018, at 11:47, James Gross  wrote:
> 
> Hi,
> 
> It had been suggested to me update my gmsh to the latest snapshot in order to 
> fix a problem I was having with boundary layers.
> 
> Fortunately, after updating the problem was fixed. Unfortunately, after I 
> updated I was unable to use the gmsh GUI. Now every time I call gmsh at 
> command line, I am presented with all the command line options. Furthermore, 
> if I attempt to use command line to open a geo file, I am presented with the 
> standard info:
> 

You probably compiled Gmsh without GUI, i.e. without FLTK...


> Info: Running 'gmsh ../Gmsh/Cases/FreeCAD2D/Impeller.geo' [Gmsh 3.0.7, 1 
> node, max. 1 thread]
> Info: Started on Mon Feb 12 00:45:50 2018
> Info: Reading '../Gmsh/Cases/FreeCAD2D/Impeller.geo'...
> Info: Reading '../Gmsh/Cases/FreeCAD2D/Impeller.brep'...
> Info: Done reading '../Gmsh/Cases/FreeCAD2D/Impeller.brep'
> Info: Done reading '../Gmsh/Cases/FreeCAD2D/Impeller.geo'
> 
> However, I am also presented with two extra lines:
> Info: Stopped on Mon Feb 12 00:45:50 2018
> Info: No current model available: creating one
>  
> As a note, I enabled Python wrappers in the new build, however I do not think 
> this would cause an issue with the GUI. 
> 
> I am using Ubuntu 16.04. After a little research, I saw that some Mac users 
> have had issues with using the GUI in the past, but I couldn't find anything 
> similar for Linux users.
> 
> I still have the older version (3.0.6) in my /usr/bin/ directory so I was 
> able to make a temporary workaround by creating an alias that allowed me to 
> call the older version to use the GUI. However it would be greatly preferred 
> to be able to call the GUI with the same command as with other operations.
> 
> Any ideas as to what the issue could be?
> 
> Thanks,
> James
> ___
> 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] refine existing mesh

2018-02-12 Thread Christophe Geuzaine 


> On 12 Feb 2018, at 15:42, Juan E. Sanchez  wrote:
> 
> Hello,
> 
> Unfortunately, I do not have the geometrical description.  Is there a way to 
> convert the mesh to a geometrical description?
> 

No ; we are working on a new implementation of our “reparametrization” 
algorithm, which will allow you to do some of this, but it’s not ready yet...


> Regards,
> 
> Juan
> 
> On 2/12/18 1:21 AM, Christophe Geuzaine wrote:
>>> On 10 Feb 2018, at 22:06, Juan E. Sanchez >> > wrote:
>>> 
>>> Hello,
>>> 
>>> I have an existing msh file, and a pos file containing a background mesh.  
>>> Is there a script example showing how to use these files to create a 
>>> refined mesh using the background field?
>>> 
>> Start from the geometrical description (not the mesh): this way Gmsh will be 
>> able to place the new vertices correctly, i.e. on the CAD
>> gmsh -3 test.geo -bgm test.pos
>> CG
>>> Using
>>> 
>>> gmsh -3 test.msh -bgm test.pos
>>> 
>>> results in a mesh where there is the same number of nodes and elements, but 
>>> it appears that the physical entity number is lost:
>>> 
>>> last element in original file:
>>> 
>>> 36210 4 3 15 15 0 5786 5791 5792 5800
>>> 
>>> last element in new file:
>>> 
>>> 36210 4 2 0 15 2603 5822 2604 2605
>>> 
>>> Now all of the mesh elements belong to the same physical entity, 0.
>>> 
>>> If I perform a "refine by splitting" in the gui, the physical entity 
>>> information is maintained:
>>> 
>>> 266256 4 2 15 15 42825 10009 41452 9710
>>> 
>>> The file is a 3d tetrahedral mesh, and it also contains triangular surfaces 
>>> to impose boundary conditions.
>>> 
>>> Regards,
>>> 
>>> Juan
>>> 
>>> 
>>> ___
>>> 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] embedded surface not meshed ?

2018-02-12 Thread Christophe Geuzaine 


> On 12 Feb 2018, at 21:01, Winfried Boxleitner <winfried.boxleit...@ait.ac.at> 
> wrote:
> 
> dear gmsh developers,
> when running the tutorial file: t15.geo, the embedded surface defined
> with "Surface{s} In Volume{1};" is not included in the 3D mesh; only
> the 2D mesh is determined correctly.
> 
> I used this feature some years ago, and at that time it worked as
> expected.
> 
> what is missing now ?
> 

Make sure you use the 3D Delaunay algorithm (not 3D Frontal, which does not 
handle embedded entities or complex size fields).


> I'm ruinning an up2date version (git pull at: Mon Feb 12 10:10:16 CET 2018)
> Version      : 3.0.7
> Build OS : Linux64
> Build date   : 20180212
> Build host   : gandalf
> Build options: 64Bit Ann Bamg Bfgs Blas(ATLAS) Blossom Cairo DIntegration 
> Dlopen Fltk GMP Jpeg Kbipack Lapack(ATLAS) LinuxJoystick MPI MUMPS MathEx 
> Mesh Metis Mmg3d Mpeg NativeFileChooser Netgen ONELAB ONELABMetamodel 
> OpenCASCADE OpenGL OptHom PETSc Parser Plugins Png Poppler Post SLEPc Solver 
> TetGen/BR Voro++ Zlib
> FLTK version : 1.3.4
> PETSc version: 3.7.5 (real arithmtic)
> OCC version  : 7.1.0
> 
> regards,
> winfried
> 
>  __
> /
> | WINFRIED BOXLEITNER
> | Scientist
> | Security & Communication Technologies
> | Center for Digital Safety & Security
> |
> | AIT Austrian Institute of Technology GmbH
> | Donau-City-Straße 1 | 1220 Vienna | Austria
> | T +43 50550-4164 | M +43 664 8251061 | F +43 50550-4190
> | winfried.boxleit...@ait.ac.at | www.ait.ac.at
> |
> | FN: 115980 i HG Wien | UID: ATU14703506
> | www.ait.ac.at/Email-Disclaimer
> \_
> 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] error: Cannot bind multiple faces to single tag

2018-02-12 Thread Christophe Geuzaine 


> On 12 Feb 2018, at 20:16, Sathyanarayan Rao  
> wrote:
> 
> Dear Gmsh team,
> 
> I am following the same syntax as 
> https://gitlab.onelab.info/gmsh/gmsh/blob/master/demos/boolean/boolean.geo 
> 
> 
> for generating boolean unioun of surfaces. But however it gives me error as 
> in the subject line. Can you please
>  help me to correctly combine surfaces 1:3 into one surface by a boolean 
> unioun ?
> 

You can only use the declarative syntax

  BooleanUnion(7) = { Surface {1}; Delete; }{ Surface {2:3}; Delete; };

If the operation results in the creation of a single shape. If not, use

 BooleanUnion { Surface {1}; Delete; }{ Surface {2:3}; Delete; }

or

 out() = BooleanUnion { Surface {1}; Delete; }{ Surface {2:3}; Delete; };

If you need the result.


> SetFactory("OpenCASCADE");
> 
> // mesh size description
> 
> cl_1   =  5;
> cl_2   =  0.01;
> 
> Geometry.NumSubEdges = 100; // nicer display of curve
> 
> var1 = 0.02; 
> var2 = 0.1;
> 
> 
> // Define Segment coordinates
> 
> X1 
> ={0,0.000811,0.001097,0.001233,0.00129,0.001325,0.001335,0.001343,0.001346,0.001347,0.001348,0.001348,0.001349,0.001349,0.001349,0.001349,0.001349,0.001349,0.001349,0.0013487};
> Y1 
> ={0,-2.4,-4,-5.6,-7.2,-8.8,-10.4,-12,-13.6,-15.2,-16.8,-18.4,-20,-21.6,-23.2,-24.8,-26.4,-28,-28.96,-29.118};
> 
> Point(newp)  ={X1[1], Y1[1], 0, cl_2};
> 
> // Line 1
> Extrude {var1, var1, 0} {
>   Point{newp-1}; 
> }
> 
> // Line 2
> Extrude {var2, var2, 0} {
>   Point{newp-1}; 
> }
> 
> // Line 3
> Extrude {var1, var1, 0} {
>   Point{newp-1}; 
> }
> 
> 
> // Define 1 segment
> 
> l0 = newl;
> l  = l0;
> nR = #X1[];
> p0  =  newp;
> p   =  p0;
> For i In {0:nR-1}
> Point(newp)  ={X1[i], Y1[i], 0, cl_2};
> EndFor
> p2  =  newp-1;
> Spline(l0)   =  {p:p2};
> 
> Wire(l0) = {l0};
> Extrude { Line{l0-3}; } Using Wire {l0}   //left
> Extrude { Line{l0-2}; } Using Wire {l0}   // center
> Extrude { Line{l0-1}; } Using Wire {l0}   // right
> 
> BooleanUnion(7) = { Surface {1}; Delete; }{ Surface {2:3}; Delete; };
> 
> 
> 
> 
> 
> Best Regards,
> 
> Sathyanarayan Rao, PhD student
> Earth and Life Institute/Environmental Sciences (ELI-e)
> Université catholique de Louvain
> c.037, Croix du Sud 2, Louvain-la-Neuve, Belgium
> Phone: 010473827 ( intercom 73827)
> ___
> 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

[Gmsh] error: Cannot bind multiple faces to single tag

2018-02-12 Thread Sathyanarayan Rao 
Dear Gmsh team,


I am following the same syntax as 
https://gitlab.onelab.info/gmsh/gmsh/blob/master/demos/boolean/boolean.geo


for generating boolean unioun of surfaces. But however it gives me error as in 
the subject line. Can you please

 help me to correctly combine surfaces 1:3 into one surface by a boolean unioun 
?


SetFactory("OpenCASCADE");

// mesh size description

cl_1   =  5;
cl_2   =  0.01;

Geometry.NumSubEdges = 100; // nicer display of curve

var1 = 0.02;
var2 = 0.1;


// Define Segment coordinates

X1 
={0,0.000811,0.001097,0.001233,0.00129,0.001325,0.001335,0.001343,0.001346,0.001347,0.001348,0.001348,0.001349,0.001349,0.001349,0.001349,0.001349,0.001349,0.001349,0.0013487};
Y1 
={0,-2.4,-4,-5.6,-7.2,-8.8,-10.4,-12,-13.6,-15.2,-16.8,-18.4,-20,-21.6,-23.2,-24.8,-26.4,-28,-28.96,-29.118};

Point(newp)  ={X1[1], Y1[1], 0, cl_2};

// Line 1
Extrude {var1, var1, 0} {
  Point{newp-1};
}

// Line 2
Extrude {var2, var2, 0} {
  Point{newp-1};
}

// Line 3
Extrude {var1, var1, 0} {
  Point{newp-1};
}


// Define 1 segment

l0 = newl;
l  = l0;
nR = #X1[];
p0  =  newp;
p   =  p0;
For i In {0:nR-1}
Point(newp)  ={X1[i], Y1[i], 0, cl_2};
EndFor
p2  =  newp-1;
Spline(l0)   =  {p:p2};

Wire(l0) = {l0};
Extrude { Line{l0-3}; } Using Wire {l0}   //left
Extrude { Line{l0-2}; } Using Wire {l0}   // center
Extrude { Line{l0-1}; } Using Wire {l0}   // right

BooleanUnion(7) = { Surface {1}; Delete; }{ Surface {2:3}; Delete; };






Best Regards,

Sathyanarayan Rao, PhD student
Earth and Life Institute/Environmental Sciences (ELI-e)
Université catholique de Louvain
c.037, Croix du Sud 2, Louvain-la-Neuve, Belgium
Phone: 010473827 ( intercom 73827)


___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

[Gmsh] embedded surface not meshed ?

2018-02-12 Thread Winfried Boxleitner 

dear gmsh developers,
when running the tutorial file: t15.geo, the embedded surface defined
with "Surface{s} In Volume{1};" is not included in the 3D mesh; only
the 2D mesh is determined correctly.

I used this feature some years ago, and at that time it worked as
expected.

what is missing now ?

I'm ruinning an up2date version (git pull at: Mon Feb 12 10:10:16 CET 2018)
Version  : 3.0.7
Build OS : Linux64
Build date   : 20180212
Build host   : gandalf
Build options: 64Bit Ann Bamg Bfgs Blas(ATLAS) Blossom Cairo DIntegration 
Dlopen Fltk GMP Jpeg Kbipack Lapack(ATLAS) LinuxJoystick MPI MUMPS MathEx Mesh 
Metis Mmg3d Mpeg NativeFileChooser Netgen ONELAB ONELABMetamodel OpenCASCADE 
OpenGL OptHom PETSc Parser Plugins Png Poppler Post SLEPc Solver TetGen/BR 
Voro++ Zlib
FLTK version : 1.3.4
PETSc version: 3.7.5 (real arithmtic)
OCC version  : 7.1.0

regards,
winfried

  __
 /
| WINFRIED BOXLEITNER
| Scientist
| Security & Communication Technologies
| Center for Digital Safety & Security
|
| AIT Austrian Institute of Technology GmbH
| Donau-City-Straße 1 | 1220 Vienna | Austria
| T +43 50550-4164 | M +43 664 8251061 | F +43 50550-4190
| winfried.boxleit...@ait.ac.at | www.ait.ac.at
|
| FN: 115980 i HG Wien | UID: ATU14703506
| www.ait.ac.at/Email-Disclaimer
 \_
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

[Gmsh] Merging two interesting surfaces as one surface

2018-02-12 Thread Sathyanarayan Rao 
Dear Gmsh team,


I am trying to transform a set of spline function into thin surfaces that can 
be given one physical group.
I am using extrude function to make these splines to have specific thickness so 
that it becomes a surface.
However the resulting surface from extrude cannot be fused together. My problem 
in a simplified code:
Is there a easy way in gmsh to convert set of splines to become surface with 
specific thickness ?  Kindly
let me know ?

SetFactory("OpenCASCADE");

// mesh size description

cl_1   =  5;
cl_2   =  3.00e-01;
r   =  0.5;


// Define Segment coordinates

X1 
={0,0.000811,0.001097,0.001233,0.00129,0.001325,0.001335,0.001343,0.001346,0.001347,0.001348,0.001348,0.001349,0.001349,0.001349,0.001349,0.001349,0.001349,0.001349,0.0013487};
Y1 
={0,-2.4,-4,-5.6,-7.2,-8.8,-10.4,-12,-13.6,-15.2,-16.8,-18.4,-20,-21.6,-23.2,-24.8,-26.4,-28,-28.96,-29.118};


X2 ={0.0002703,1.276,1.564,1.433,1.302,1.171,1.04,0.9087,0.7122,0.71171};
Y2 ={-0.8,-0.8009,-1.366,-2.351,-3.335,-4.32,-5.305,-6.29,-7.767,-7.7709};


// Define 1 segment

LN = 90;
nR = #X1[];
p0  =  newp;
p   =  p0;
For i In {0:nR-1}
Point(newp)  ={X1[i], Y1[i], 0, cl_2};
EndFor
p2  =  newp-1;
Spline(LN)   =  {p:p2};
Extrude {r, r, 0} { Line{LN};}

// Define 2 segment

LN = 95;
nR = #X2[];
p0  =  newp;
p   =  p0;
For i In {0:nR-1}
Point(newp)  ={X2[i], Y2[i], 0, cl_2};
EndFor
p2  =  newp-1;
Spline(LN)   =  {p:p2};
Extrude {r, r, 0} { Line{LN};}

BooleanUnion(100) = { Surface {1}; Delete; }{ Surface {2}; Delete; };
Physical Surface(1) = {100};


Best Regards,

Sathyanarayan Rao, PhD student
Earth and Life Institute/Environmental Sciences (ELI-e)
Université catholique de Louvain
c.037, Croix du Sud 2, Louvain-la-Neuve, Belgium
Phone: 010473827 ( intercom 73827)


___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh

Re: [Gmsh] Merging two interesting surfaces as one surface

2018-02-12 Thread Sathyanarayan Rao 
I also tried using wire function which looks better than extrusion but I am 
unable to fuse two surfaces together to form a single surface


Wire code:


SetFactory("OpenCASCADE");

// mesh size description

cl_1   =  5;
cl_2   =  3.00e-01;
r   =  0.5;


// Define Segment coordinates

X1 
={0,0.000811,0.001097,0.001233,0.00129,0.001325,0.001335,0.001343,0.001346,0.001347,0.001348,0.001348,0.001349,0.001349,0.001349,0.001349,0.001349,0.001349,0.001349,0.0013487};
Y1 
={0,-2.4,-4,-5.6,-7.2,-8.8,-10.4,-12,-13.6,-15.2,-16.8,-18.4,-20,-21.6,-23.2,-24.8,-26.4,-28,-28.96,-29.118};


X2 ={0.0002703,1.276,1.564,1.433,1.302,1.171,1.04,0.9087,0.7122,0.71171};
Y2 ={-0.8,-0.8009,-1.366,-2.351,-3.335,-4.32,-5.305,-6.29,-7.767,-7.7709};


// Define 1 segment

LN = 90;
nR = #X1[];
p0  =  newp;
p   =  p0;
For i In {0:nR-1}
Point(newp)  ={X1[i], Y1[i], 0, cl_2};
EndFor
p2  =  newp-1;
Spline(LN)   =  {p:p2};

// Define 2 segment

LN = 95;
nR = #X2[];
p0  =  newp;
p   =  p0;
For i In {0:nR-1}
Point(newp)  ={X2[i], Y2[i], 0, cl_2};
EndFor
p2  =  newp-1;
Spline(LN)   =  {p:p2};

//+
Translate {0, -0.1, 0} {
  Duplicata { Point{21}; }
}
//+
Line(96) = {21, 31};
Wire(1) = {95};
Extrude { Line{96}; } Using Wire {1}

Translate {0.1, 0.1, 0} {
  Duplicata { Point{1}; }
}

Line(101) = {1, 36};
Wire(10) = {90};
Extrude { Line{101}; } Using Wire {10}

Mesh.CharacteristicLengthMax = 0.1;
BooleanUnion(1000) = { Surface {1}; Delete; }{ Surface {2}; Delete; };
Physical Surface(1) = {1000};




Best Regards,

Sathyanarayan Rao, PhD student
Earth and Life Institute/Environmental Sciences (ELI-e)
Université catholique de Louvain
c.037, Croix du Sud 2, Louvain-la-Neuve, Belgium
Phone: 010473827 ( intercom 73827)




From: Sathyanarayan Rao
Sent: 12 February 2018 12:07:49
To: gmsh@onelab.info
Subject: Merging two interesting surfaces as one surface


Dear Gmsh team,


I am trying to transform a set of spline function into thin surfaces that can 
be given one physical group.
I am using extrude function to make these splines to have specific thickness so 
that it becomes a surface.
However the resulting surface from extrude cannot be fused together. My problem 
in a simplified code:
Is there a easy way in gmsh to convert set of splines to become surface with 
specific thickness ?  Kindly
let me know ?

SetFactory("OpenCASCADE");

// mesh size description

cl_1   =  5;
cl_2   =  3.00e-01;
r   =  0.5;


// Define Segment coordinates

X1 
={0,0.000811,0.001097,0.001233,0.00129,0.001325,0.001335,0.001343,0.001346,0.001347,0.001348,0.001348,0.001349,0.001349,0.001349,0.001349,0.001349,0.001349,0.001349,0.0013487};
Y1 
={0,-2.4,-4,-5.6,-7.2,-8.8,-10.4,-12,-13.6,-15.2,-16.8,-18.4,-20,-21.6,-23.2,-24.8,-26.4,-28,-28.96,-29.118};


X2 ={0.0002703,1.276,1.564,1.433,1.302,1.171,1.04,0.9087,0.7122,0.71171};
Y2 ={-0.8,-0.8009,-1.366,-2.351,-3.335,-4.32,-5.305,-6.29,-7.767,-7.7709};


// Define 1 segment

LN = 90;
nR = #X1[];
p0  =  newp;
p   =  p0;
For i In {0:nR-1}
Point(newp)  ={X1[i], Y1[i], 0, cl_2};
EndFor
p2  =  newp-1;
Spline(LN)   =  {p:p2};
Extrude {r, r, 0} { Line{LN};}

// Define 2 segment

LN = 95;
nR = #X2[];
p0  =  newp;
p   =  p0;
For i In {0:nR-1}
Point(newp)  ={X2[i], Y2[i], 0, cl_2};
EndFor
p2  =  newp-1;
Spline(LN)   =  {p:p2};
Extrude {r, r, 0} { Line{LN};}

BooleanUnion(100) = { Surface {1}; Delete; }{ Surface {2}; Delete; };
Physical Surface(1) = {100};


Best Regards,

Sathyanarayan Rao, PhD student
Earth and Life Institute/Environmental Sciences (ELI-e)
Université catholique de Louvain
c.037, Croix du Sud 2, Louvain-la-Neuve, Belgium
Phone: 010473827 ( intercom 73827)


___
gmsh mailing list
gmsh@onelab.info
http://onelab.info/mailman/listinfo/gmsh