<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
Hello,
<div class=""><br class="">
</div>
<div class=""><span style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" class="">So I am trying to make a simple 5 pt stencil finite difference laplacian operator (</span><a href="https://akmitra.public.iastate.edu/aero361/design_web/Laplace.pdf" class="">https://akmitra.public.iastate.edu/aero361/design_web/Laplace.pdf</a><span style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0);" class="">)
work on a dmplex mesh, but all of the examples that I have found in ksp/examples/tutorials all use DMDA. The only way that I have </span><font color="#000000" class="">managed to get the KSP to even run is by creating the operator matrix myself by hand. The
workflow for the working <span style="caret-color: rgb(0, 0, 0);" class="">implementation is as follows:</span></font></div>
<div class=""><font color="#000000" class=""><span style="caret-color: rgb(0, 0, 0);" class=""><br class="">
</span></font></div>
<div class=""><font color="#000000" class=""><span style="caret-color: rgb(0, 0, 0);" class="">DMPlexMakeBoxMesh -> KSPCreate -> KSPsetDM -> KSPDeactiveDM -> Make the matrix by hand -> KSPSetOperator(matrix) -> KSPSetComputeRHS -> Solve</span></font></div>
<div class=""><font color="#000000" class=""><span style="caret-color: rgb(0, 0, 0);" class=""><br class="">
</span></font></div>
<div class=""><font color="#000000" class="">Given that I tether KSP to the DM by KSPSetDM, but then immediately deactivate it (because then it complains about mallocs), this doesn’t work when reading in an actual unstructured mesh. Does KSP just not support
dmplex?</font></div>
<div class=""><font color="#000000" class=""><br class="">
</font></div>
<div class="">
<div class="">
<div>Best,</div>
<div><br class="">
</div>
<div>Jacob</div>
</div>
<br class="">
</div>
</body>
</html>