<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Jan 16, 2014 at 2:41 PM, Dharmendar Reddy <span dir="ltr"><<a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello,<br>
         I there a way to fix a non Delauny mesh via call to any<br>
DMPLex routine ?<br>
I am looking at this scenario.<br>
<br>
I calculate the orthocenter of each element. If i find that for some<br>
elements, the othocenter is outside the element, is there any way i<br>
can fix it via a call to DMPLEx routine or do i have to go back to the<br>
mesher and regenerate the mesh ?<br></blockquote><div><br></div><div>"Flipping" to enforce the Delaunay condition is useful, but can be costly and involve</div><div>large changes to the structure. Thus you want something like the heap structures used</div>
<div>in Triangle. Thus I do not think it belongs in PETSc.</div><div><br></div><div>However, you can use DMRefine() with a large cell tolerance to just do flipping.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

Thanks<br>
<span class="HOEnZb"><font color="#888888">Reddy<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>