I was able to solve successfully some problems in elasticity in 3D.
Mesh is extremely simple (a revolved square, so forming a ring).
However, as soon as I increase the number of nodes and elements, I get the error :
“Error umfpack_di_solve status -3”
I noticed that other users experienced the same issue but I was able to find no answers.
It cannot be a memory problem since I have RAM 32 Gb and Windows 10 64bit and the memory occupied by the solver is just a tiny fraction before it gives the error.
To understand the problem, I tried to repeat the Example “beam-3d.edp” that I find in the examples/3d subdirectory.
I replaced line 4 :
with the following:
int[int] Nxyz=[20 * factor, 5 * factor, 5 * factor];
with factor up to 3.5 there is no crash, and very little memory occupation. However, it crashes with factor from 4 upwards, even if the number of elements and nodes is still very low.
Do you have any solution?
Many thanks in advance