The problem is, that in line 4893 not the diagonal element of the given matrix is used for this decision, but the first element of the ith row of the diagonal matrix block (A_diag_data[A_diag_I[i]]), ...

The finite element method (FEM) is a numerical method for solving problems of engineering and mathematical physics. Typical problem areas of interest include structural analysis, heat transfer, fluid ...