Sorry to bother you, just to understand if I got it.
If your variables are interlaced like > > u0,v0,w0,u1,v1,w1,... > > Then you can (a) set a block size for the matrix with MatSetBlockSize, > (b) use BAIJ in which case the block size is built in, or (c) use > -pc_fieldsplit_block_size. Then to define the splits, use > > -pc_fieldsplit_0_fields 0,3 -pc_fieldsplit_1_fields 1 > -pc_fieldsplit_2_fields 2,4,5 > > to split six fields into three splits of sizes 2,1,3. > > > If my variables are u v w p (3D stokes equation) and I have an equal order approximation Can I just do -pc_fieldsplit_0_fields 0,1,2 -pc_fieldsplit_1_fields 3 -pc_fieldsplit_type schur -pc_filedsplit_block_size <n_dofs_per_variable> Thank you very much Lorenzo ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
