<div dir="ltr">Can you please give more details on what 'does not work' is.<div>More detail on how you judge what works would also be useful.</div><div>Mark</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Dec 7, 2021 at 11:19 AM Marco Cisternino <<a href="mailto:marco.cisternino@optimad.it">marco.cisternino@optimad.it</a>> wrote:<br></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="IT" style="overflow-wrap: break-word;">
<div class="gmail-m_-4790067142172085774WordSection1">
<p class="MsoNormal">Good morning,<u></u><u></u></p>
<p class="MsoNormal"><span lang="EN-GB">I’m still struggling with the Poisson equation with Neumann BCs.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">I discretize the equation by finite volume method and I divide every line of the linear system by the volume of the cell. I could avoid this division, but I’m trying to understand.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">My mesh is not uniform, i.e. cells have different volumes (it is an octree mesh).<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">Moreover, in my computational domain there are 2 separated sub-domains.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">I build the null space and then I use MatNullSpaceTest to check it.<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">If I do this:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">MatNullSpaceCreate(getCommunicator(), PETSC_TRUE, 0, nullptr, &nullspace);<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">It works<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">If I do this:<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="FR">Vec nsp;<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="FR">VecDuplicate(m_rhs, &nsp);<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="FR">VecSet(nsp,1.0);<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="FR">VecNormalize(nsp, nullptr);<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">MatNullSpaceCreate(getCommunicator(), PETSC_FALSE, 1, &nsp, &nullspace);<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">It does not work<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">Probably, I have wrong expectations, but should not it be the same?<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">Thanks<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
<p class="MsoNormal"><span>Marco Cisternino, PhD<br>
<a href="mailto:marco.cisternino@optimad.it" target="_blank"><span style="color:blue">marco.cisternino@optimad.it</span></a><u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">______________________<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">Optimad Engineering Srl<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB">Via Bligny 5,
</span><span lang="EN-GB" style="font-family:"Segoe UI",sans-serif;color:rgb(32,31,30);background:white">Torino, Italia.</span><span lang="EN-GB" style="font-family:"Segoe UI",sans-serif;color:rgb(32,31,30)"><br>
<span style="background:white">+3901119719782</span><br>
</span><span><a href="http://www.optimad.it/" target="_blank"><span lang="EN-GB" style="font-family:"Segoe UI",sans-serif;color:blue;border:1pt none windowtext;padding:0cm;background:white">www.optimad.it</span></a></span><span lang="EN-GB"><u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-GB"><u></u> <u></u></span></p>
</div>
</div>

</blockquote></div>