<div dir="ltr">Hello everybody,<div><br></div><div>I am using Gmsh inside my own code, through the API, to mesh a lot of very small geometries that are independent one from the other. My goal is to produce for every geometry the coarsest mesh possible.</div><div><br></div><div>In order to speedup the process, I've tried to create a single GModel for every geometry and to process them in parallel (using openmp). However, due to certain runtime errors I am getting, it seems that either I am doing something wrong in my code or it is not thread-safe to create multiple instances of the GModel class and to work with them in parallel.</div><div><br></div><div>On the other hand, I've seen that Gmsh offers multithreading parallelization.</div><div>My question is: can I benefit from that?</div><div>If different regions (corresponding to different geometries) are meshed in parallel, my code would speedup; if the parallelism is perform inside the mesh process of every region/geometry, probably the code would remain mostly serial.</div><div><br></div><div>After meshing all the geometries, I check them, looking for possible negative jacobians. For the ones that present negative jacobians I either refined them or use the high-order optimization algorithm. But this process is only executed for the geometries presenting negative jacobians.</div><div><br></div><div>So, my second question is: after meshing the first GModel (hopefully in parallel), can I remove or deactivate some of the regions / elements (present in the first GModel) that are already valid (their jacobians are positive everywhere)? Then, I would refine/optimize only the ones that present negative jacobians.</div><div><br></div><div><br></div><div><div>Thanks for the support.</div><div>Best regards,</div><div><br></div><div><br></div><div>-- </div><div><br></div><div><div>Pablo Antolin</div><div>Researcher</div><div><span style="font-size:12.8px">Chair of Numerical Modelling and Simulation </span><br><span style="font-size:12.8px">Institute of Mathematics</span></div><div><span style="font-size:12.8px">Ecole polytechnique fédérale de Lausanne</span><br></div></div></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px"><br></span></div></div>