<div dir="ltr"><div class="gmail_extra">On Thu, Jan 10, 2013 at 2:35 PM, David Scott <span dir="ltr"><<a href="mailto:d.scott@ed.ac.uk" target="_blank">d.scott@ed.ac.uk</a>></span> wrote:<br><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div id=":4e6">Thanks for the replies. I obviously need to think about the problem more. I was aware that the solution was not fully determined and I had created a MatNullSpace to deal with the arbitrary constant.<br>

<br>
I have contacted the original author of the code (I am just modifying it to use PETSc) to find out more about the formulation of the problem. I should point out that the the actual problem I want to deal with does not have a zero RHS.<br>

<br>
This systems of equations is only part of the complete system and only the gradient of the solution is required elsewhere so the arbitrary constant has no physical significance.</div></blockquote></div><br>You are in 2D and the dimension of your null space is the number of grid points on that "non-boundary". It can change your solution in a myriad of ways.</div>
</div>