On Thu, 29 Apr 2010 13:02:49 +0200, Lorenzo Botti <[email protected]> wrote: > Thank you for the hint, > I had to put also -pc_type fieldsplit, just in case someone else wants to > give it a try... > > Is there a way to monitor the convergence of the solvers?
The most common monitors are -ksp_monitor, -ksp_monitor_true_residual, -ksp_converged_reason These can be set on inner solves with, e.g. -fieldsplit_0_ksp_converged_reason -fieldsplit_1_ksp_monitor_true_residual > Another question, is bs, the block size, always equal to the total > number of fields? 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(). Jed ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
