29 Nov
2016
29 Nov
'16
9:36 p.m.
On 29 Nov 2016, at 21:25, Francis Poulin <fpoulin@uwaterloo.ca> wrote:
Thanks Lawrence.
That did seem to help. Now I get another errror but now I get that it fails to converge after 0 iterations.
It suggests that I add some options. Should I try doing this in firedrake? I am not sure how to use these options but if you think it would help and I can figure out how I am happy to give it a try.
Oh, I should have said. The default LU in petsc does no pivoting and you have a whacking great zero block in your system. Try additionally with: "pc_factor_mat_solver_package": "mumps" The error you saw comes because petsc encountered a zero pivot in the factorisation stage. Lawrence