If I change the element order of a geometry-file (which works with 1st order elements) to 2nd order ones, Elmer crashes with a segfault, so there seems to be a bug hidden somewhere.
Here's the output snippet:
Code: Select all
MAIN: -------------------------------------
MAIN: Steady state iteration: 1
MAIN: -------------------------------------
MAIN:
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
CRS_IncompleteLU: ILU(0) (Real), Starting Factorization:
CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 919722
CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100
CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.04
62 0.8251E-11
ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.73219875 2.0000000 ) :: mgdynamics
Segmentation fault