Hi <br><br>I think so. At the beginning I wrote a function in deal.ii to "manually" delete the repeated points. Later I found that I can use a command in Gmsh to do the same.<br>The procedure I followed was:<br><br>
1 - write the geometry script into .geo file<br>2 - import .geo file into Gmsh. Generate mesh (2D or 3D). Save mesh into .mesh file.<br>3 - write another .geo file with the commands<br><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote">
Merge "filename.msh";<br>Coherence Mesh;<br></blockquote>4 - import last .geo file into gmsh. In the message console you should read the amount of repeated points deleted.<br>5 - Save last mesh into another .mesh file.<br>
<br>I tried this with a much more complex 3D geometry and it worked quite well. Although I didn't tried with the original example I posted.<br><br>By the way, if you're using Gmsh with Deal.ii don't forget to check the last tutorial example posted "step-49". It deals with meshes in general and specifically it offers some hints about the interaction between Gmsh and Deal.ii.<br>
<br>Hope this helps.<br><br>Regards<br>Javier<br><br><br><div class="gmail_quote">2013/7/10 <span dir="ltr"><<a href="mailto:ragab.fayez@gmail.com" target="_blank">ragab.fayez@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi <span><span style="color:rgb(34,34,34)">Javier <br>Have you got a solution for this problem because i faced the same problem <br><br>Thanks <br></span></span><div><div class="h5"><br>On Saturday, 29 October 2011 00:55:33 UTC+2, Javier Muñoz wrote:</div>
</div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5">Hi<p>I'm trying to extrude two lines:</p><p>lc = 0.01;<br>p1 = newp; Point (p1) = {0,0,0,lc};<br>
p2 = newp; Point (p2) = {0.5,0,0,lc};<br>p3 = newp; Point (p3) = {1.0,0,0,lc};</p><p>l1 = newl; Line (l1) = {p1,p2};<br>l2 = newl; Line (l2) = {p2,p3};</p><p>p1 -------------- p2 ---------------p3</p><p><br>Extrusion height is 1 for line 1 (with 1 layer). Line two has a height<br>
of z = 0.5 (with 1 layer) and then I want to extrude the 'top' face<br>generated to a height of z = 1 (also with 1 layer).</p><p><br>out1[] = Extrude {0,1,0} {<br> Line {1};<br> Layers {{1},{1}};<br>
Recombine;<br> };<br>out2[] = Extrude {0,0.5,0} {<br> Line {2};<br> Layers {{1},{1}};<br> Recombine;<br> };<br>out3[] = Extrude {0,0.5,0} {<br>
Line {out2[0]};<br> Layers {{1},{1}};<br> Recombine;<br> };<br>Coherence;</p><p>|---------------||------------<u></u>--|<br>| || <u></u> |<br>
| || <u></u> |<br>| || <u></u> |<br>| || <u></u> |<br>| || <u></u> |<br>| || <u></u> |<br>
| ||-----------<u></u>---|<br>| || <u></u> |<br>| || <u></u> |<br>| || <u></u> |<br>| || <u></u> |<br>
| || <u></u> |<br>| || <u></u> |<br>p1 ------------ p2 ------------ p3</p><p>The problem is: I have lines overlapping in the interface of surface<br>out1. The line on the right side of out1 overlaps with the lines on the<br>
left side of out2 and out3. Is there anyway to avoid this? As you can<br>see I already try to use Coherence (maybe in the wrong way?).</p><p>I want to import the mesh generated into deal.ii but as long as I have<br>overlapping lines, deal.ii appears to recognize the interfaces as<br>
boundaries.</p><p>I appreciate your help.</p><p>Regards<br>Javier Munoz</p></div></div><p><br>______________________________<u></u>_________________<br>gmsh mailing list<br><a>gm...@geuz.org</a><br><a href="http://www.geuz.org/mailman/listinfo/gmsh" target="_blank">http://www.geuz.org/mailman/<u></u>listinfo/gmsh</a><br>
</p><p></p><p></p><p></p><p></p><p></p><p></p><p></p><p></p><p></p><p></p><p></p></blockquote><div></div></blockquote></div><br><br clear="all"><br>-- <br>José Javier Muñoz Criollo<br>PhD Student | Cardiff University School of Engineering<br>
E-mail: <a href="mailto:josejmcriollo@gmail.com" target="_blank">josejmcriollo@gmail.com</a><br>