[petsc-users] fieldsplit: Unhandled case, must have at least two fields, not 1!
Jed Brown
jedbrown at mcs.anl.gov
Mon Dec 2 08:04:32 CST 2013
Mark Adams <mfadams at lbl.gov> 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.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20131202/fc101ba0/attachment.pgp>
More information about the petsc-users
mailing list