I could not use streamplot, so updated the firedrake yesterday to fix the problem. It still doesn't work and additionally, I encounter a strange nonlinear solver divergence error for the same problem that it was able to solve before the package was updated. I don't know if the problem is due to petsc configuration or something else. Would you please help me resolve the issue?
Two figures are attached to show the problem could be solved without error before the update. The text file attached shows the SNES monitor output. It is worth to mention the problem can still be solved when the gradient of the function at the top and bottom boundaries of the domain is not large. Thank you very much in advance.