Dear Dokken,
Thank you a lot for your time and your quick reply.
Sorry, I had no idea it was different; I guess I need to study the documentation a bit better.
In my script substituted:
solve(F == 0, u, form_compiler_parameters={"linear_solver": "mumps"})
with:
solve(F == 0, u, solver_parameters={"newton_solver": {"linear_solver": "mumps"}})
And I stopped getting the out of Memory error message.
My script still didn’t work; I got an error message similar to the one reported in this previous topic, and I solved it as explained in that topic.
Now it looks like it is working, even though it’s taking a lot of time to solve of course.