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)
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)