About src/ksp/pc/impls/fieldsplit/fieldsplit.c

Ryan Yan vyan2000 at gmail.com
Fri Jul 17 17:18:08 CDT 2009

Hi All,
I have some difficulty of understanding the struct PC_FieldSplit. From the
definition I can see that it has a data structure like a "list" and each
member in the list is representing an object Field(or PC_Field).

Suppose that my matrix has a blocksize of 5.   I want to Set 3 Fields,  i.e.
field_0 as {0,1},  field_1 as{2,3},and field_3 as {4}. Then can anyone
please help me to fill in the following parameter.

PC_FieldSplit *jac = (PC_FieldSplit*)pc->data;

what is jac->nsplit?     3, is my guess, since we have split the matrix into
3 split, namely field_0, field_1, field_2.

PC_FieldSplitLink ilink = jac->head;

what is ilink->nfields? 2, is my guess, since the fields_0 has 2 fields

Then ilink=ilink->next; ilink=ilin->next; ilink->nfileds should be 1 right?

Thank you very much in advance.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090717/40c6726b/attachment.htm>

More information about the petsc-users mailing list