<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
<br>
      if (i == 0) {<br>
        for (f = 0; f < numFields; ++f) {<br>
          ierr = PCFieldSplitSetIS(pc, fieldNames[f], fields[f]);CHKERRQ(ierr);<br>
          ierr = PetscFree(fieldNames[f]);CHKERRQ(ierr);<br>
          ierr = ISDestroy(&fields[f]);CHKERRQ(ierr);<br>
        }<br>
<br>
numFields should be 2.  Compare to src/snes/examples/tutorials/ex28.c<br>
(adapting runex28_3),<br></blockquote><div><br></div><div>The matrix/solver just has one field. The error comes from giving it field split parameters by mistake.</div></div></div></div>