<div class="gmail_quote">On Mon, Mar 12, 2012 at 02:56, Klaij, Christiaan <span dir="ltr"><<a href="mailto:C.Klaij@marin.nl">C.Klaij@marin.nl</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">>Declare this as KSP *subksp;<br></div></blockquote><div><br></div><div>Compare this ^^^^^^^^^^^</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">
><br>
> Matt<br>
><br>
<br>
</div>Thanks Matt, but when I do that I cannot use KSPSetNullSpace:<br>
<div class="im"><br>
KSP *subksp[2];<br></div></blockquote><div><br></div><div>To this ^^^^^^^^^^</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">
PetscInt n=2;<br>
ierr = PCFieldSplitGetSubKSP(pc,&n,subksp); CHKERRQ(ierr);<br>
ierr = KSPSetNullSpace(subksp[1],subnullsp); CHKERRQ(ierr);<br>
<br>
</div>This gives the error:<br>
<div class="im"><br>
error: argument of type "KSP *" is incompatible with parameter of type<br>
"KSP"<br>
ierr = KSPSetNullSpace(subksp[1],subnullsp); CHKERRQ(ierr);</div></blockquote></div><br>