<div dir="ltr"><div><div>Hi Bart,<br><br></div>Is it possible for you to send the STL file so that I can play with it ?<br><br></div>Thanks,<br><br>V.<br></div><div class="gmail_extra"><br><div class="gmail_quote">2014-09-10 6:07 GMT-05:00 Bart Vermeulen <span dir="ltr"><<a href="mailto:bart.vermeulen@wur.nl" target="_blank">bart.vermeulen@wur.nl</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
  
    
  
  <div text="#000000" bgcolor="#FFFFFF">
    <div>Hi Vincent,<br>
      <br>
      I get the attached error messages when I remove the CreateTopology
      command.<br>
      <br>
      I think it misses lines and points to build compund line from.<span class="HOEnZb"><font color="#888888"><br>
      <br>
      Bart</font></span><div><div class="h5"><br>
      <br>
      On 09/09/2014 05:27 PM, Vincent Libertiaux wrote:<br>
    </div></div></div><div><div class="h5">
    <blockquote type="cite">
      
      <div dir="ltr">
        <div>Hi Bart,<br>
          <br>
        </div>
        Could you please check what happens if you remove the
        CreateTopology command ?<br>
        <br>
        V.<br>
      </div>
      <div class="gmail_extra"><br>
        <div class="gmail_quote">2014-09-03 5:17 GMT-05:00 Bart
          Vermeulen <span dir="ltr"><<a href="mailto:bart.vermeulen@wur.nl" target="_blank">bart.vermeulen@wur.nl</a>></span>:<br>
          <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear all,<br>
            <br>
            I am trying to create a mesh starting with a STL model. This
            model consists of 6 different surfaces. I use the file
            included below (quite similar to one of the tutorials). When
            I load this into gmsh I can see the geometry with the
            different surfaces. At the corners of each of the surfaces
            there is a point (4 points for each surface). The Compound
            lines connect these points directly instead of following the
            geometry (see attached screenshot). Obviously the meshing
            fails.<br>
            <br>
            Any idea how to properly define the compound lines? (The
            surfaces ar OK)<br>
            <br>
            Thanks in advance,<br>
            Bart Vermeulen<br>
            <br>
            <br>
            <br>
            Merge "mmbendo.stl";<br>
            CreateTopology;<br>
            <br>
            // We can now define a compound line (resp. surface) for
            each discrete<br>
            // line (resp. surface) in the model<br>
            ll[] = Line "*";<br>
            For j In {0 : #ll[]-1}<br>
              Compound Line(newl) = ll[j];<br>
            EndFor<br>
            ss[] = Surface "*";<br>
            s = news;<br>
            For i In {0 : #ss[]-1}<br>
              Compound Surface(s+i) = ss[i];<br>
            EndFor<br>
            <br>
            // And we can create the volume based on the new compound
            entities<br>
            Surface Loop(1) = {s : s + #ss[]-1};<br>
            Volume(1) = {1};<br>
            <br>
            Physical Surface(1) = {s : s + #ss[]-1};<br>
            Physical Volume(1) = 1;<br>
            <br>
            uniform = 1;<br>
            If(uniform)<br>
              // uniform mesh size...<br>
              Mesh.CharacteristicLengthMin = 15;<br>
              Mesh.CharacteristicLengthMax = 15;<br>
            EndIf<br>
            If(!uniform)<br>
              // ... or apply a funny mesh size field, just because we
            can :-)<br>
              Field[1] = MathEval;<br>
              Field[1].F = "2*Sin((x+y)/5) + 3";<br>
              Background Field = 1;<br>
            EndIf<br>
            <br>
            Mesh.RemeshAlgorithm = 0; // (0) no split (1) automatic (2)
            automatic only with metis<br>
            Mesh.RemeshParametrization = 7; // (0) harmonic (1)
            conformal spectral (7) conformal finite element<br>
            Geometry.HideCompounds = 0; // don't hide the compound
            entities<br>
            Mesh.Algorithm = 6; // Frontal<span><font color="#888888"><br>
                <br>
                -- <br>
                Bart Vermeulen<br>
                Wageningen University, Department of Environmental
                Sciences<br>
                Hydrology and Quantitative Water Management Group<br>
                Postal address: Postbus 47, 6700AA Wageningen,The
                Netherlands<br>
                Visiting address: Room 0.C027, Droevendaalsesteeg
                3,6708PB Wageningen, The Netherlands<br>
                Tel. <a href="tel:%2B31-317-482765" value="+31317482765" target="_blank">+31-317-482765</a> Fax <a href="tel:%2B31-317-419000" value="+31317419000" target="_blank">+31-317-419000</a><br>
                email: <a href="mailto:bart.vermeulen@wur.nl" target="_blank">bart.vermeulen@wur.nl</a><br>
                <br>
              </font></span><br>
            _______________________________________________<br>
            gmsh mailing list<br>
            <a href="mailto:gmsh@geuz.org" target="_blank">gmsh@geuz.org</a><br>
            <a href="http://www.geuz.org/mailman/listinfo/gmsh" target="_blank">http://www.geuz.org/mailman/listinfo/gmsh</a><br>
            <br>
          </blockquote>
        </div>
        <br>
      </div>
    </blockquote>
    <br>
    <br>
    <pre cols="72">-- 
Bart Vermeulen
Wageningen University, Department of Environmental Sciences
Hydrology and Quantitative Water Management Group
Postal address: Postbus 47, 6700AA Wageningen,The Netherlands
Visiting address: Room 0.C027, Droevendaalsesteeg 3,6708PB Wageningen, The Netherlands
Tel. <a href="tel:%2B31-317-482765" value="+31317482765" target="_blank">+31-317-482765</a> Fax <a href="tel:%2B31-317-419000" value="+31317419000" target="_blank">+31-317-419000</a>
email: <a href="mailto:bart.vermeulen@wur.nl" target="_blank">bart.vermeulen@wur.nl</a></pre>
  </div></div></div>

</blockquote></div><br></div>