On Wed, May 4, 2011 at 15:29, <domenico.borzacchiello at univ-st-etienne.fr>wrote:
> -fieldsplit_pressure_pc_type lsc \ > -fieldsplit_pressure_lsc_ksp_type fgmres \ > -fieldsplit_pressure_lsc_ksp_max_it 100 \ > -fieldsplit_pressure_lsc_ksp_constant_null_space \ > -fieldsplit_pressure_lsc_pc_type lu \ > -fieldsplit_pressure_lsc_pc_factor_mat_solver_package mumps > Perhaps this direct solve is the problem? That operator is expected to be singular due to the constant null space so a direct solve is expected to produce nonsense. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20110504/0381a492/attachment.html>
