Mark Adams <[email protected]> writes: >> >> >> >> 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.
So this is a different KSP; not the one you want to use fieldsplit on? Please always use different prefixes for different solvers.
pgpIZz9FlYQ5Y.pgp
Description: PGP signature
