<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:arial, helvetica, sans-serif;font-size:12pt"><div><div><div><span class="Apple-style-span" style="font-size: x-large;"><i>Here is the output with opt 1) &nbsp;VecSet(b,h):</i></span></div><div><br></div><div>&nbsp;0 KSP preconditioned resid norm 5.656854249492e+00 true resid norm 5.656854249492e+00 ||Ae||/||Ax|| 1.000000000000e+00</div><div>&nbsp;&nbsp; &nbsp;1 KSP preconditioned resid norm 1.256073966947e-15 true resid norm 8.881784197001e-16 ||Ae||/||Ax|| 1.570092458684e-16</div><div>Linear solve converged due to CONVERGED_RTOL iterations&nbsp;</div><div>Process [0]</div><div>0.666667</div><div>0.666667</div><div>Process [1]</div><div>0.666667</div><div>0.666667</div><div>Process [2]</div><div>0.666667</div><div>0.666667</div><div>Process [3]</div><div>0.666667</div><div>0.666667</div><div>&nbsp;&nbsp; &nbsp;0 KSP preconditioned resid norm
 5.656854249492e+00 true resid norm 5.656854249492e+00 ||Ae||/||Ax|| 1.000000000000e+00</div><div>&nbsp;&nbsp; &nbsp;1 KSP preconditioned resid norm 1.256073966947e-15 true resid norm 8.881784197001e-16 ||Ae||/||Ax|| 1.570092458684e-16</div><div>Linear solve converged due to CONVERGED_RTOL iterations 1</div><div>&nbsp;&nbsp; &nbsp;0 KSP preconditioned resid norm 5.656854249492e+00 true resid norm 5.656854249492e+00 ||Ae||/||Ax|| 1.000000000000e+00</div><div>&nbsp;&nbsp; &nbsp;1 KSP preconditioned resid norm 1.256073966947e-15 true resid norm 8.881784197001e-16 ||Ae||/||Ax|| 1.570092458684e-16</div><div>Linear solve converged due to CONVERGED_RTOL iterations 1</div><div><br></div></div><div><br></div><div><span class="Apple-style-span" style="font-size: x-large; font-style: italic; ">Here is the output with opt 2) Copied from another function:</span></div><div><br></div><div><br></div><div>Linear solve converged due to CONVERGED_ATOL iterations
 0</div><div>Process [0]</div><div>0</div><div>0</div><div>Process [1]</div><div>0</div><div>0</div><div>Process [2]</div><div>0</div><div>0</div><div>Process [3]</div><div>0</div><div>0</div><div>&nbsp;&nbsp; &nbsp;0 KSP preconditioned resid norm 7.515606252360e+01 true resid norm 7.515606252360e+01 ||Ae||/||Ax|| 1.000000000000e+00</div><div>&nbsp;&nbsp; &nbsp;1 KSP preconditioned resid norm 2.014114995652e+01 true resid norm 2.014114995652e+01 ||Ae||/||Ax|| 2.679910213523e-01</div><div>&nbsp;&nbsp; &nbsp;2 KSP preconditioned resid norm 6.102231836870e+00 true resid norm 6.102231836870e+00 ||Ae||/||Ax|| 8.119413966044e-02</div><div>&nbsp;&nbsp; &nbsp;3 KSP preconditioned resid norm 3.582243925876e-01 true resid norm 3.582243925876e-01 ||Ae||/||Ax|| 4.766407134156e-03</div><div>&nbsp;&nbsp; &nbsp;4 KSP preconditioned resid norm 9.891345584245e-02 true resid norm 9.891345584245e-02 ||Ae||/||Ax|| 1.316107477176e-03</div><div>&nbsp;&nbsp; &nbsp;5 KSP
 preconditioned resid norm 3.103414625576e-14 true resid norm 2.787743265007e-14 ||Ae||/||Ax|| 3.709272640689e-16</div><div>Linear solve converged due to CONVERGED_RTOL iterations 5</div><div>&nbsp;&nbsp; &nbsp;0 KSP preconditioned resid norm 7.567218624719e+01 true resid norm 7.567218624719e+01 ||Ae||/||Ax|| 1.000000000000e+00</div><div>&nbsp;&nbsp; &nbsp;1 KSP preconditioned resid norm 2.027946648035e+01 true resid norm 2.027946648035e+01 ||Ae||/||Ax|| 2.679910213523e-01</div><div>&nbsp;&nbsp; &nbsp;2 KSP preconditioned resid norm 6.144138058565e+00 true resid norm 6.144138058565e+00 ||Ae||/||Ax|| 8.119413966044e-02</div><div>&nbsp;&nbsp; &nbsp;3 KSP preconditioned resid norm 3.606844483858e-01 true resid norm 3.606844483858e-01 ||Ae||/||Ax|| 4.766407134156e-03</div><div>&nbsp;&nbsp; &nbsp;4 KSP preconditioned resid norm 9.959273013421e-02 true resid norm 9.959273013421e-02 ||Ae||/||Ax|| 1.316107477176e-03</div><div>&nbsp;&nbsp; &nbsp;5 KSP
 preconditioned resid norm 2.724152110884e-14 true resid norm 3.587687035382e-14 ||Ae||/||Ax|| 4.741090767039e-16</div><div>Linear solve converged due to CONVERGED_RTOL iterations 5</div><div><br></div><div><br></div><div>Here the printed zeros are the print of DMMGGetx()</div><div><br></div><div>Thanks.</div></div><div style="position:fixed"></div>


</div><br>

      </body></html>