[petsc-users] Fixing a non Delaunay mesh

[Karl Rupp]

Hi Reddy,

 >           I there a way to fix a non Delauny mesh via call to any
> DMPLex routine ?
> I am looking at this scenario.
>
> I calculate the orthocenter of each element. If i find that for some
> elements, the othocenter is outside the element, is there any way i
> can fix it via a call to DMPLEx routine or do i have to go back to the
> mesher and regenerate the mesh ?

keep in mind that the Delaunay property only defines the duality of the 
simplex mesh and its dual (Voronoi) grid. As far as I remember, it does 
*not* require that each orthocenter is in the interior of the 
corresponding element.

Best regards,
Karli