Mesh generation with MMG2D ignoring edge table


#21

Ok the error comes from a wrong path given for loading the mesh. So the right log is now this one:

– READING DATA FILE (null)
**WARNING NO GIVEN TRIANGLE
MAXIMUM MEMORY AUTHORIZED (MB) 2000
_MMG2D_NPMAX 50000
_MMG2D_NTMAX 100000
NUMBER OF VERTICES 551 CORNERS 0
NUMBER OF TRIANGLES 0
NUMBER OF EDGES 551
mmg_mesh_reload_ptr->npmax: 50000
mmg_mesh_reload_ptr->np: 551
mmg_mesh_ptr->npmax: 50000
mmg_mesh_ptr->np: 551
– MMG2D, Release 5.3.8 (Apr. 10, 2017)
Copyright © IMB-LJLL, 2004-
Jun 11 2018 11:55:38

– MMG2DMESH: INPUT DATA
– INPUT DATA COMPLETED. 0.000s

&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE MMG2D-IMB/LJLL : 5.3.8 (Apr. 10, 2017) ISO
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MAXIMUM NUMBER OF POINTS (NPMAX) : 50000
MAXIMUM NUMBER OF TRIANGLES (NTMAX) : 100000

– PHASE 1 : MESH GENERATION
449 vertex inserted 102 not inserted
34 vertex inserted 68 not inserted
2 vertex inserted 66 not inserted
0 vertex inserted 66 not inserted
Error: MMG2_insertpointdelone: some vertex are not inserted 3.

Warning: MMG2_pack: unexpected edge table… Ignored data.
NUMBER OF VERTICES 555 CORNERS 0
NUMBER OF TRIANGLES 1098
NUMBER OF EDGES 4

Sorry for the mistake
Best
Pierre


#22

Have you rebuild and reinstall the Mmg2d library since the last debugging? Maybe you need to remove the library by hand and to rebuild it.

I think that you run an old version of Mmg2d.

BR


#23

You are right: I have rebuild the lib but forgot to reinstall.

OK log are not the same. Thank you for for help it works now on the mesh I gave you previously. Hurray!
Here is the corresponding log: without_line_remeshing_OK.log (8.5 KB)

But I would like to generate a mesh with constant mesh size of 1 meter. So I resampled the given mesh:
here are the new mesh and sol files : before.sol (7.7 KB)
before.mesh (148.7 KB)

Using this meshes there are issues: see attached log here: with_line_remeshing_KO.log (1.9 KB)

Thank for all the help :wink:
I feel that we are very close to solve the problem.

Best
Pierre


#24

Hi,

I have founded another wrong epsilon in the computation of the cavity of Delaunay (the modification is pushed on the develop branch). Now it works on this mesh.

Regards,

Algiane


#25

Hi Algiane,

I have just tried with the updated develop branch and it works perfectly. :grin:
Thanks a lot for your help.

Best regards,
Pierre


#26

Hi Algiane,

thank again for all your help.
I have a subsidiary question about bugs you have fixed. Do mmg use exact predicates for geometrical operations?

Regards
Pierre


#27

Hi Pierre,

You’re welcome.

Mmg doesn’t use exact geometry reconstruction. We build a boundary approximation using fourth order bezier curves computed from the normals at nodes.

Details are available section 3.2 of this paper.

Regards,

Algiane