How to minimize time complexity when executing FEM on a highly refined mesh?

[dokken]

Did you try:

 solve(a==L, u, bc, solver_parameters={"linear_solver": "mumps"})

as in:

Out of memory error

Many things can be done, like using “mumps” ( solve(a==l, deformation, bcs, solver_parameters={"linear_solver": "mumps"}) as the solver or running in parallel using MPI (mpirun -n 2 python3 script.py)