
<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>

</head>

<body ocsi="0" fpstyle="1">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">I have an elliptic PDE of the form Div[n(x) T Grad(phi)]=f(n,...) subject to periodic BC. (phi is the electrostatic potential in a PIC/fluid model). I expect there to be a constant

 null space but I'm not sure I'm dealing with it properly. <br>

<br>

1) Section 4.6 of the manual suggests MatNullSpaceCreate() followed by KSPSetNullSpace(), but I have also seen examples that use MatSetNullSpace(). Is one preferred over the other?<br>

2) Should I follow XXXSetNullSpace() with MatNullSpaceRemove() on the RHS Vec?<br>

3) Do I need to MatNullSpaceDestroy() for either XXXSetNullSpace() option?<br>

4) If so, should the destruction happen after KSPSolve()?<br>

5) The manual (same section) says "You can run with the additional options -pc_factor_shift_nonzero <dampingfactor> or -pc_factor_shift_nonzero <dampingfactor> to prevent a zero pivot." This seems like a typo.<br>

<br>

I believe I've done my due Google diligence, but I'm happy to be pointed elsewhere if these questions have come up before.<br>

<br>

--A Different Matt<br>

<div><br>

<div style="font-family:Tahoma; font-size:13px">

<div style="font-family:Tahoma; font-size:13px">--------------------------------------------------------------

<div>Matthew Young</div>

<div>Graduate Student</div>

<div>Boston University Dept. of Astronomy</div>

<div>--------------------------------------------------------------</div>

<div><br>

</div>

</div>

</div>

</div>

</div>

</body>

</html>