<div dir="ltr"><div dir="ltr">On Mon, Jul 4, 2022 at 1:29 AM BOYD Andrew via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</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">





<div lang="EN-GB" style="overflow-wrap: break-word;">
<div class="gmail-m_-2113785299130254604WordSection1">
<div>
<p class="MsoNormal"><span lang="EN-US">Dear PETSc users,</span><u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p>
</div>
<p class="MsoNormal"><span lang="EN-US">I am a PhD student looking to solve a Laplacian on a 3D domain which consists of a volume with ellipsoid “holes” in it. I am a newcomer to PETSc, and would like advice on how to begin to write this code. I have studied
 the examples, and example no.12 in SNES seems to share some of the </span>same characteristics that I need.
<u></u><u></u></p>
<p class="MsoNormal">I have also looked at the mail archives, but unfortunately couldn’t find too much on this problem – it may be that I lack the background to understand it properly.
<u></u><u></u></p>
<p class="MsoNormal">Thanks in advance for any advice, it is much appreciated.<u></u><u></u></p>
<p class="MsoNormal">(As an aside, I am also planning to run this in Rust using the PETSc bindings – if anyone has any experience with this it would be interesting to hear from you as well!)</p></div></div></blockquote><div><br></div><div>It seems job 1 for you is forming the matrix with the correct boundary conditions. Ex 12 should be enough for that if you can input a mesh</div><div>with the boundaries marked as you wish. After that I would</div><div><br></div><div>1) Run with a direct solver like MUMPS or SuperLU_dist</div><div><br></div><div>2) Run with algebraic multigrid, like GAMG, or Hypre, or ML</div><div><br></div><div>3) If that is not enough, think about a few steps of geometric or p-refinement.</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"><div lang="EN-GB" style="overflow-wrap: break-word;"><div class="gmail-m_-2113785299130254604WordSection1">
<p class="MsoNormal">Sincerely,<u></u><u></u></p>
<p class="MsoNormal">Andrew Boyd<u></u><u></u></p>
</div>
The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336. Is e buidheann carthannais a th’ ann an Oilthigh Dhùn Èideann, clàraichte an Alba, àireamh clàraidh SC005336.
</div>

</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>