Hi Anna, On 10 Nov 2014, at 19:39, Anna Kalogirou <a.kalogirou@leeds.ac.uk> wrote:
Hi Lawrence,
Does it matter that I am solving a linear problem with the nonlinear solver? I am just trying to reproduce the linear one that works first, before introducing the nonlinear terms. I tried also adding the parameter {'snes_type': 'ksponly'} and the code works, but the result is not as good.
Do you have any idea what might be happening?
P.S. The linear code can be found attached.
Thanks, I think I have found the problem, it is not in your code. The problem appears to be a subtle issue in which we interact with PETSc, which means that it is not computing the initial residual correctly. I will prepare a fix tomorrow and let you know when things are fixed. Cheers, Lawrence