The edge length computation in Mmg3d contains some inconsistencies and needs to be tested and cleaned:
- for isotropic metric, it seems that we always compute the "straight" edge length;
- for anisotropic metric:
- the "straight edge length is computed along internal edges
- along boundary edges:
- we may compute the "curve" length along edges marked by the
MG_BDY tag
- or the straight length:
- because the tag is missing (we are not able to ensure the update of this tag everywhere)
- because we are on a boundary edge belonging the a tetra without any boundary faces and without an xtetra (I am sure that such edges exists but not sure that we attempt to compute their lengths however)
- because we are inside an operator when tags are not yet updated
Moreover, the readibility of this part the code needs to be improved (function names not very meaningful and complex function pointer digram)
Steps
The edge length computation in Mmg3d contains some inconsistencies and needs to be tested and cleaned:
MG_BDYtagMoreover, the readibility of this part the code needs to be improved (function names not very meaningful and complex function pointer digram)
Steps