On Fri, Aug 17, 2012 at 3:10 AM, Thomas Witkowski <span dir="ltr"><<a href="mailto:thomas.witkowski@tu-dresden.de" target="_blank">thomas.witkowski@tu-dresden.de</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">
I want to solve some (weakly) coupled system of equations of the following form:<br>
<br>
A  B     u<br>
        .     =   .....<br>
0  C     v<br>
<br>
<br>
so, C is the discrete Laplacian and A and B are some more complicated operators (I make use of linear finite elements). All boundary conditions are periodic, so the unknown v is determined only up to a constant. A and B contain both the identity operator, so u is fixed. Now I want to solve the system on the whole (there are reasons to do it in this way!) and I must provide information about the nullspace to the solver. When I am right, to provide the correct nullspace I must solve one equation with A. Is there any way in PETSc to circumvent the problem?</blockquote>
<div><br></div><div>If I understand you correctly, your null space vector is (0 I). I use the same null space for SNES ex62.</div><div><br></div><div>Note that you should probably use FieldSplit multiplicative here since you have a block triangular system to start with.</div>
<div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="HOEnZb"><font color="#888888"><br>
Thomas</font></span></blockquote></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>