<div dir="ltr"><div dir="ltr">On Thu, Mar 19, 2020 at 5:05 PM Jeremy L Thompson <<a href="mailto:jeth8984@colorado.edu">jeth8984@colorado.edu</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Good Afternoon,<br>
<br>
I have a working high-order 3D elasticity code built on PETSc + libCEED. <br>
My DM has a three component displacement FE field. I'd like to add a <br>
discontinuous pressure field (P_-1).<br></blockquote><div><br></div><div>By P_{-1} do you mean a discontinuous linear basis? I only ask because FV is normally P_0.</div><div>Anyway, PetscFE has discontinuous bases if you want to use that. You use the flag</div><div><br></div><div>  <a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DUALSPACE/PetscDualSpaceLagrangeSetContinuity.html">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DUALSPACE/PetscDualSpaceLagrangeSetContinuity.html</a></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:1px solid rgb(204,204,204);padding-left:1ex">
Am I missing an obvious way to set up a discontinuous FE field? Or is a <br>
FV field a more natural choice for this sort of thing? Or a different <br>
approach?<br>
<br>
I am managing action of the FDE operator myself via libCEED, so I don't <br>
need to get too into the weeds about anything but the DM. I just need to <br>
find the best way to manage the discontinuous pressure values on each <br>
element in my DM.<br>
<br>
Thanks,<br>
<br>
Jeremy L Thompson<br>
<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>