Hi 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. The procedure I followed was:
1 - write the geometry script into .geo file 2 - import .geo file into Gmsh. Generate mesh (2D or 3D). Save mesh into .mesh file. 3 - write another .geo file with the commands > Merge "filename.msh"; > Coherence Mesh; > 4 - import last .geo file into gmsh. In the message console you should read the amount of repeated points deleted. 5 - Save last mesh into another .mesh file. 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. 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. Hope this helps. Regards Javier 2013/7/10 <[email protected]> > Hi Javier > Have you got a solution for this problem because i faced the same problem > > Thanks > > On Saturday, 29 October 2011 00:55:33 UTC+2, Javier Muñoz wrote: > >> Hi >> >> I'm trying to extrude two lines: >> >> lc = 0.01; >> p1 = newp; Point (p1) = {0,0,0,lc}; >> p2 = newp; Point (p2) = {0.5,0,0,lc}; >> p3 = newp; Point (p3) = {1.0,0,0,lc}; >> >> l1 = newl; Line (l1) = {p1,p2}; >> l2 = newl; Line (l2) = {p2,p3}; >> >> p1 -------------- p2 ---------------p3 >> >> >> Extrusion height is 1 for line 1 (with 1 layer). Line two has a height >> of z = 0.5 (with 1 layer) and then I want to extrude the 'top' face >> generated to a height of z = 1 (also with 1 layer). >> >> >> out1[] = Extrude {0,1,0} { >> Line {1}; >> Layers {{1},{1}}; >> Recombine; >> }; >> out2[] = Extrude {0,0.5,0} { >> Line {2}; >> Layers {{1},{1}}; >> Recombine; >> }; >> out3[] = Extrude {0,0.5,0} { >> Line {out2[0]}; >> Layers {{1},{1}}; >> Recombine; >> }; >> Coherence; >> >> |---------------||------------**--| >> | || ** | >> | || ** | >> | || ** | >> | || ** | >> | || ** | >> | || ** | >> | ||-----------**---| >> | || ** | >> | || ** | >> | || ** | >> | || ** | >> | || ** | >> | || ** | >> p1 ------------ p2 ------------ p3 >> >> The problem is: I have lines overlapping in the interface of surface >> out1. The line on the right side of out1 overlaps with the lines on the >> left side of out2 and out3. Is there anyway to avoid this? As you can >> see I already try to use Coherence (maybe in the wrong way?). >> >> I want to import the mesh generated into deal.ii but as long as I have >> overlapping lines, deal.ii appears to recognize the interfaces as >> boundaries. >> >> I appreciate your help. >> >> Regards >> Javier Munoz >> >> >> ______________________________**_________________ >> gmsh mailing list >> [email protected] >> http://www.geuz.org/mailman/**listinfo/gmsh<http://www.geuz.org/mailman/listinfo/gmsh> >> >> -- José Javier Muñoz Criollo PhD Student | Cardiff University School of Engineering E-mail: [email protected]
_______________________________________________ gmsh mailing list [email protected] http://www.geuz.org/mailman/listinfo/gmsh
