<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear Christophe,</div><div><br></div><div>my example files are scientific data, however originate from processing certain proprietary 3D models that were shared with me under certain restrictions. Therefore, it's difficult to share my original file with you.</div><div>However, I was indeed able to reproduce the issue using only a slight modification of your example file!</div><div><br></div><div>Compare  the following two files:</div><div>test1.geo:</div><div>%%%
%%%


%%%


%%%


%%%


%%%


%%%


%%%


%%%


%%%

</div><div>SetFactory("OpenCASCADE");<br>Box(1) = {0,0,0, 1,1,1};<br>For i In {1:1000}<br>  Point(100+i) = {0.25 + 1e-4*i, 0.1,0.1};<br>  Point{100+i} In Volume{1};<br>EndFor</div><div>
%%%
%%%


%%%


%%%


%%%


%%%


%%%


%%%


%%%


%%% <br></div><div><br></div><div>test2.geo:</div><div>

%%%
%%%


%%%


%%%


%%%


%%%


%%%


%%%


%%%


%%% <br></div><div>SetFactory("OpenCASCADE");<br>Box(1) = {0,0,0, 1,1,1};<br>For i In {1:10000}<br>  Point(100+i) = {0.25 + 1e-5*i, 0.1,0.1};<br>  Point{100+i} In Volume{1};<br>EndFor<br></div><div>

%%%
%%%


%%%


%%%


%%%


%%%


%%%


%%%


%%%


%%% <br></div><div><br></div><div>The first file contains 1000 internal pts, and the second one contains 10 times as much.</div><div>Now run:</div><div>  $ time gmsh -1 test1.geo</div><div>    -> Takes 0.88 sec<br></div><div>
  $ time gmsh -1 test2.geo <br></div><div>
    -> Takes 75 sec <br></div><div>This is a nearly 100 (!!!) times increase in parsing run time, when the model size is increased by only a factor of 10.</div><div>This nicely illustrates what I meant by "quadratic growth of the parsing run time as a function of the number of internal pts."</div><div>As a consequence of this "quadratic growth", the run time for larger models quickly becomes enormous.</div><div><br></div><div>Best regards</div><div><br></div><div>A.S.<br></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Am Do., 31. Jan. 2019 um 14:31 Uhr schrieb Christophe Geuzaine <<a href="mailto:cgeuzaine@uliege.be">cgeuzaine@uliege.be</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
Can you send a test file?<br>
<br>
I tried this, i.e. adding 1000 embedded points in a volume, and it is quite fast (< 2 seconds):<br>
<br>
SetFactory("OpenCASCADE");<br>
Box(1) = {0,0,0, 1,1,1};<br>
For i In {1:1000}<br>
  Point(100+i) = {0.25 + 5e-4*i, 0.1,0.1};<br>
  Point{100+i} In Volume{1};<br>
EndFor<br>
<br>
Maybe you modify the model after each insertion of an embedded point, which would force a rebuild of the topological model?<br>
<br>
Christophe<br>
<br>
<br>
> On 31 Jan 2019, at 13:04, Al Sc <<a href="mailto:al.sc.gmsh@gmail.com" target="_blank">al.sc.gmsh@gmail.com</a>> wrote:<br>
> <br>
> Dear Sir or Madam,<br>
> <br>
> some more details on the previously-reported case of extremely slow parsing of a gmsh file:<br>
> The file contains 34795 points -- of which 23041 points are "Point In Volume" (e.g. embedded_vertices of a GRegion?). Moreover, it contains 4998 "Plane Surface"s and one single "Volume". Almost all plane surfaces are triangles and quads -- except for <10 facets, which each have a very high vertex count (typ. 7000).<br>
> <br>
> I found out that if I remove all "Point In Volume" objects, then the parsing takes only 20 seconds. However, if all the "Point In Volume" objects are not removed, parsing takes around 3000 seconds. This led me to the hypothesis that there is some huge inefficiency with parsing gmsh files with "Point In Volume" (probably quadratic run time in N(PointInVolume)??).<br>
> <br>
> Furthermore, I analyzed the runs using "perf" (linux utility).<br>
> Run 1 (about 3000 seconds):<br>
>   $ time perf record ./gmsh-3.0.6-Linux64/bin/gmsh -1 input_orig.geo<br>
> Run 2 (about 20 seconds):<br>
>   $ time perf record ./gmsh-3.0.6-Linux64/bin/gmsh -1 input_no_point_in_volume.geo<br>
> In Run 1, "perf report" yields the following top-scoring functions:<br>
>   11.57%  gmsh     gmsh                 [.] GModel::getEdgeByTag<br>
>    7.73%  gmsh     gmsh                 [.] GEO_Internals::synchronize<br>
>    6.80%  gmsh     <a href="http://libc-2.12.so" rel="noreferrer" target="_blank">libc-2.12.so</a>         [.] _int_free<br>
>    6.31%  gmsh     gmsh                 [.] GModel::getVertexByTag<br>
>    4.91%  gmsh     <a href="http://libc-2.12.so" rel="noreferrer" target="_blank">libc-2.12.so</a>         [.] malloc<br>
>    4.53%  gmsh     gmsh                 [.] InterpolateCurve<br>
>    4.40%  gmsh     libstdc++.so.6.0.22  [.] std::_Rb_tree_increment<br>
>    3.88%  gmsh     <a href="http://libc-2.12.so" rel="noreferrer" target="_blank">libc-2.12.so</a>         [.] memcpy<br>
>    3.69%  gmsh     gmsh                 [.] List_Nbr<br>
>    3.35%  gmsh     <a href="http://libc-2.12.so" rel="noreferrer" target="_blank">libc-2.12.so</a>         [.] _int_malloc<br>
>    3.24%  gmsh     gmsh                 [.] GEntity::GEntity<br>
>    3.13%  gmsh     gmsh                 [.] avl_lookup<br>
>    2.63%  gmsh     gmsh                 [.] gmshFace::resetNativePtr<br>
>    2.53%  gmsh     gmsh                 [.] GEdge::addFace<br>
>    2.52%  gmsh     gmsh                 [.] CompareVertex<br>
>    1.97%  gmsh     gmsh                 [.] std::_List_base<GEdge*, std::allocator<GEdge*><br>
>    1.88%  gmsh     gmsh                 [.] GModel::getFaceByTag<br>
>    1.76%  gmsh     gmsh                 [.] Tree_Action<br>
>    1.74%  gmsh     gmsh                 [.] gmshEdge::degenerate<br>
>    1.63%  gmsh     gmsh                 [.] CompareCurve<br>
>    1.53%  gmsh     gmsh                 [.] std::_List_base<int, std::allocator<int> >::_M<br>
>    1.38%  gmsh     gmsh                 [.] GModel::deletePhysicalGroups<br>
>    1.23%  gmsh     gmsh                 [.] std::_List_base<GEdgeSigned, std::allocator<GE<br>
>    1.10%  gmsh     gmsh                 [.] GVertex::addEdge<br>
>    0.85%  gmsh     gmsh                 [.] List_Read<br>
>    0.75%  gmsh     gmsh                 [.] GEdge::getBeginVertex<br>
>    0.73%  gmsh     gmsh                 [.] GFace::computeMeanPlane<br>
>    0.73%  gmsh     <a href="http://libc-2.12.so" rel="noreferrer" target="_blank">libc-2.12.so</a>         [.] free<br>
>    0.67%  gmsh     gmsh                 [.] CTX::instance<br>
>    0.65%  gmsh     gmsh                 [.] gmshEdge::resetMeshAttributes<br>
> <br>
> This suggests that maybe GModel::getEdgeByTag is eventually called a quadratic number of times in the number of PointInVolume-objects -- and this causes the drastic slow-down.<br>
> <br>
> Could you please investigate this? Thanks a lot!<br>
> (As already mentioned, this issue also occurs with gmsh-4.x.x)<br>
> <br>
> Best regards<br>
> A. S.<br>
> _______________________________________________<br>
> gmsh mailing list<br>
> <a href="mailto:gmsh@onelab.info" target="_blank">gmsh@onelab.info</a><br>
> <a href="http://onelab.info/mailman/listinfo/gmsh" rel="noreferrer" target="_blank">http://onelab.info/mailman/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.be/~geuzaine</a><br>
<br>
<br>
<br>
</blockquote></div>