30 Mar
2017
30 Mar
'17
9:24 a.m.
Hi Francis,
On 30 Mar 2017, at 02:50, Francis Poulin <fpoulin@uwaterloo.ca> wrote:
$ python QG1L_Stommel_psiuq.py Linear firedrake_1_ solve did not converge due to DIVERGED_PCSETUP_FAILED iterations 0 PCSETUP_FAILED due to FACTOR_NUMERIC_ZEROPIVOT
The PETSc LU by default doesn't do any pivoting, and you may have a zero on the diagonal somewhere? Try using mumps instead, by adding 'pc_factor_mat_solver_package': 'mumps' to your nonlinear solver options. Cheers, Lawrence