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