> > > > Yeah, specifying the fields with options like -pc_fieldsplit_0_field > 0,1,2 requires a collocated discretization where the fields at each node > are interlaced (this gives better performance anyway). It won't work > for mixed discretizations, or where the fields are not interlaced. In > that case, you should set them with PCFieldSplitSetIS(). > > I'm not sure I understand this point. What do you mean with interlaced? Isn't the off-diagonal block realizing the inter-field coupling?
I'm looking at the source code and PCSetUp_FieldSplit() is creating the IS according to the field splits. What kind of situation can I manage with with PCFieldSplitSetIS() ? Thanks Lorenzo ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
