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. Traceback (most recent call last): File "linear_stommel_sw.py", line 58, in <module> sw_solver.solve() File "/home/fpoulin/software/firedrake/src/firedrake/firedrake/variational_solver.py", line 199, in solve solving_utils.check_snes_convergence(self.snes) File "/home/fpoulin/software/firedrake/src/firedrake/firedrake/solving_utils.py", line 153, in check_snes_convergence %s""" % (snes.getIterationNumber(), msg)) firedrake.exceptions.ConvergenceError: Nonlinear solve failed to converge after 0 nonlinear iterations. Reason: Inner linear solve failed to converge after 0 iterations with reason: unknown reason (petsc4py enum incomplete?), try with -snes_convered_reason and -ksp_converged_reason Francis ------------------ Francis Poulin Associate Professor Department of Applied Mathematics University of Waterloo email: fpoulin@uwaterloo.ca Web: https://uwaterloo.ca/poulin-research-group/ Telephone: +1 519 888 4567 x32637 ________________________________________ From: firedrake-bounces@imperial.ac.uk [firedrake-bounces@imperial.ac.uk] on behalf of Lawrence Mitchell [lawrence.mitchell@imperial.ac.uk] Sent: Tuesday, November 29, 2016 4:15 PM To: firedrake@imperial.ac.uk Subject: Re: [firedrake] no normal flow Boundary Conditions
On 29 Nov 2016, at 20:56, Francis Poulin <fpoulin@uwaterloo.ca> wrote:
Again, the updated code can be found at the following in case someone wants to look at it.
https://github.com/francispoulin/firedrakeQG/blob/master/linear_stommel_sw.p...
Set "mat_type": "aij" In your solver parameters. Lawrence _______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake