>
>
>
>       if (i == 0) {
>         for (f = 0; f < numFields; ++f) {
>           ierr = PCFieldSplitSetIS(pc, fieldNames[f],
> fields[f]);CHKERRQ(ierr);
>           ierr = PetscFree(fieldNames[f]);CHKERRQ(ierr);
>           ierr = ISDestroy(&fields[f]);CHKERRQ(ierr);
>         }
>
> numFields should be 2.  Compare to src/snes/examples/tutorials/ex28.c
> (adapting runex28_3),
>

The matrix/solver just has one field. The error comes from giving it field
split parameters by mistake.

Reply via email to