<div dir="ltr"><div>Hi Aaron, <br></div><div><br></div><div>I might be mistaken, but it looks like your in-file Line Loop stuff is good. <br></div><div>That might mean the Surface Loop command is giving the warning because the two Line Loops don't match up (i.e. form a closed surface)?</div><div><br></div><div>Sincerely, <br></div><div>Max<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Oct 11, 2018 at 1:25 PM, Jensen, Aaron J. <span dir="ltr"><<a href="mailto:Aaron.J.Jensen@leidos.com" target="_blank">Aaron.J.Jensen@leidos.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Christophe,<br>
<br>
Thank you for the fast response.  The latest snapshot works great and fixes the crash.<br>
<br>
I am still setting up the problem incorrectly.  Here is what I have now:<br>
<br>
<br>
Merge "square";<br>
Line Loop(1) = {1};<br>
<br>
squareWidth = 1.5;<br>
halfSquareWidth = squareWidth/2.0;<br>
<br>
lc = 0.3;<br>
Point(1001) = {-halfSquareWidth,-<wbr>halfSquareWidth,0.0,lc};<br>
Point(1002) = {-halfSquareWidth, halfSquareWidth,0.0,lc};<br>
Point(1003) = { halfSquareWidth, halfSquareWidth,0.0,lc};<br>
Point(1004) = { halfSquareWidth,-<wbr>halfSquareWidth,0.0,lc};<br>
Line(1001) = {1004,1003};<br>
Line(1002) = {1003,1002};<br>
Line(1003) = {1002,1001};<br>
Line(1004) = {1001,1004};<br>
<br>
Line Loop(1005) = {1002,1003,1004,1001};<br>
<br>
Surface Loop(1006) = {1005, 1};<br>
<br>
<br>
I think it is close to working but I get the warning message: "1D mesh not forming a closed loop".<br>
<br>
Any help would be greatly appreciated.<br>
<br>
Thanks,<br>
Aaron<br>
<br>
-----Original Message-----<br>
From: Christophe Geuzaine [mailto:<a href="mailto:cgeuzaine@uliege.be">cgeuzaine@uliege.be</a>] <br>
Sent: Thursday, October 11, 2018 8:11 AM<br>
To: Jensen, Aaron J. [US-US]<br>
Cc: <a href="mailto:gmsh@onelab.info">gmsh@onelab.info</a><br>
Subject: EXTERNAL: Re: [Gmsh] Embedded Mesh<br>
<br>
<br>
<br>
> On 10 Oct 2018, at 21:12, Jensen, Aaron J. <<a href="mailto:Aaron.J.Jensen@leidos.com">Aaron.J.Jensen@leidos.com</a>> wrote:<br>
> <br>
> Hello,<br>
>  <br>
> I would like to merge an existing mesh into a geometry and then mesh the unmeshed region using gmsh (something like a boundary mesh around a pre-existing mesh generated externally by another program).  Ideally, this would result in a mesh that was partially created in gmsh and partially created externally to gmsh.  Ultimately, I would like to do this in both 2D and 3D but am starting with the following in 2D:<br>
>  <br>
>  <br>
>  <br>
> Merge "square";<br>
>  <br>
> SetFactory("OpenCASCADE");<br>
>  <br>
> Surface Loop(1) = {1};<br>
> //sExternal = news;<br>
> //Rectangle(sExternal) = {0, 0, 0.0, 0.5, 0.5};<br>
>  <br>
> s1 = news;<br>
> Rectangle(s1) = {-0.25, -0.25, 0.0, 1.0, 1.0};<br>
>  <br>
> //BooleanDifference (2000) = { Surface{s1}; Delete; }{ Surface{sExternal}; };<br>
> BooleanDifference (2000) = { Surface{s1}; Delete; }{ Surface{1}; };<br>
>  <br>
<br>
This cannot work: boolean operations can only be applied to CAD entities. It should not crash though - and it does not on my machine. Can you try with the latest snapshot?<br>
<br>
You can still produce a mesh like the one you want. With the built-in engine, simply <br>
<br>
- load the mesh, which should contain the bounding surface mesh<br>
<br>
- define the exterior "Surface Loop" where the interior surfaces are the bounding surfaces (discrete)<br>
<br>
Christophe<br>
<br>
<br>
>  <br>
>  <br>
> OpenCascade doesn’t recognize surface 1 and BooleanDifference causes gmsh to crash.  I assume I am doing something wrong?  Can gmsh embed an existing mesh into a gmsh geometry or CAD model?<br>
>  <br>
> Thanks,<br>
> Aaron<br>
> <square>______________________<wbr>_________________________<br>
> gmsh mailing list<br>
> <a href="mailto:gmsh@onelab.info">gmsh@onelab.info</a><br>
> <a href="http://onelab.info/mailman/listinfo/gmsh" rel="noreferrer" target="_blank">http://onelab.info/mailman/<wbr>listinfo/gmsh</a><br>
<br>
— <br>
Prof. Christophe Geuzaine<br>
University of Liege, Electrical Engineering and Computer Science <br>
<a href="http://www.montefiore.ulg.ac.be/~geuzaine" rel="noreferrer" target="_blank">http://www.montefiore.ulg.ac.<wbr>be/~geuzaine</a><br>
<br>
Free software: <a href="http://gmsh.info" rel="noreferrer" target="_blank">http://gmsh.info</a> | <a href="http://getdp.info" rel="noreferrer" target="_blank">http://getdp.info</a> | <a href="http://onelab.info" rel="noreferrer" target="_blank">http://onelab.info</a><br>
<br>
______________________________<wbr>_________________<br>
gmsh mailing list<br>
<a href="mailto:gmsh@onelab.info">gmsh@onelab.info</a><br>
<a href="http://onelab.info/mailman/listinfo/gmsh" rel="noreferrer" target="_blank">http://onelab.info/mailman/<wbr>listinfo/gmsh</a><br>
</blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>Max Orok</div><div>Contractor<br></div></div><div dir="ltr"><a href="http://www.mevex.com" target="_blank">www.mevex.com</a><br><br><img src="https://docs.google.com/uc?export=download&id=1fHTIiW4OMUjQr1iOkspQ7wiEsxunoOs0&revid=0B6x5w-5zVaEjSkpwbm5oY29jbG1XMzJoYldXTmJpNGFtb3dVPQ" width="164" height="42"><br></div></div></div></div></div></div>
</div>