[petsc-users] Fixing a non Delaunay mesh

Matthew Knepley knepley at gmail.com
Thu Jan 16 15:22:27 CST 2014

On Thu, Jan 16, 2014 at 2:41 PM, Dharmendar Reddy
<dharmareddy84 at gmail.com>wrote:

> Hello,
>          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 ?

"Flipping" to enforce the Delaunay condition is useful, but can be costly
and involve
large changes to the structure. Thus you want something like the heap
structures used
in Triangle. Thus I do not think it belongs in PETSc.

However, you can use DMRefine() with a large cell tolerance to just do


> Thanks
> Reddy

What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140116/01689f93/attachment.html>

More information about the petsc-users mailing list