[petsc-users] Fixing a non Delaunay mesh

[Dharmendar Reddy]

I should have worded the problem differently. I intended to say, fix
an element to have its orthocenter inside the element to avoid errors
in approximating the flux along an edge.

Thanks
Reddy

On Fri, Jan 17, 2014 at 5:01 AM, Karl Rupp <rupp at mcs.anl.gov> wrote:
> 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
>