Cannot run Fenics code with fine mesh (memory issue?)

[dokken]

@mekong

UMFPACK V5.7.1 (Oct 10, 2014): ERROR: out of memory

The error is quite clear.
Please search the forum for this message.

UMFPACK error: out of memory despite system having free memory

Hey, this could indeed be a bug of umfpack. As far as I know this is because fenics is using a 32 bit version of this (someone please correct me if I’m wrong). What you can do, however, is to use mumps, which is 64 bit, and can support more memory. For a bilinear form a, linear form L and boundary conditions given in bcs, the solve command for this looks like solve(a==L, u, bcs, solver_parameters={'linear_solver' : 'mumps'}) If you want more freedom in configuring mumps, I’ll refer you to my…

Large Mesh Error PETSc Error 76 - Adjoint Calculation

So this error message is related to the UMFPACK bug i referenced above. However, I notice that you get an actual error when running: solve(A, w.vector(), L, "mumps") Temporary fix would be to change this to: solver_lin = LUSolver("mumps") solver_lin.solve(A, w.vector(), L) Edit: PR added