Re: [firedrake] problem solving a linear system

On 21/06/17 17:00, Francis Poulin wrote:

Saving file at /QG2L-Stommel-psi-u-q.ipynb Linear firedrake_0_ solve did not converge due to DIVERGED_PCSETUP_FAILED iterations 0 PCSETUP_FAILED due to FACTOR_NUMERIC_ZEROPIVOT

This is your problem. You have a 0 on the diagonal of your operator. If you are sure it is non-singular, try adding: "pc_factor_mat_solver_package": "mumps" to your solver options. Lawrence