<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sat, May 30, 2015 at 2:37 AM, Justin Chang <span dir="ltr"><<a href="mailto:jychang48@gmail.com" target="_blank">jychang48@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Hello everyone,<br><br></div>In SNES ex62, it seems to solve the Stokes equation using TH elements (i.e., using P2-P1 elements). However, this formulation is not locally conservative, and as a result many people use the P2-P1/P0 elements where the pressure space is enriched with a piece-wise constant to ensure mass conservation at the element level. <br><br></div>Is such an element combination possible with the current implementation within DMPlex?<br></div></div></div></blockquote><div><br></div><div>I have implemented P1_disc:</div><div><br></div><div>  <a href="https://bitbucket.org/petsc/petsc/src/779682a0c231963f14b692ce9dcbb6520a47bf75/config/builder.py?at=master#cl-497">https://bitbucket.org/petsc/petsc/src/779682a0c231963f14b692ce9dcbb6520a47bf75/config/builder.py?at=master#cl-497</a></div><div><br></div><div>I believe its working correctly based on my tests for function interpolation, but I have not tried anything hard with it.</div><div><br></div><div> Thanks,</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr"><div><div></div>Thanks,<br></div>Justin<br></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>