Block SSOR preconditionings for high-order 3D FE systems. II: Incomplete BSSOR preconditionings (Q1176527)

[For part I see BIT 29, 805-823 (1989; Zbl 0689.65016)]. The authors are interested in solving a linear system arising from a conforming hierarchical finite element discretization of the equations for small displacements of an orthotropic elastic material. Orders of the polynomial trial functions range from two to five. The matrix involved is not the original stiffness matrix, but its Schur complement constructed by considering \(2\times 2\times 2\) groupings of elements as superelements. The model problem is defined on the unit cube divided nonuniformly into \(8\times 8\times 8=512\) elements in a manner which is not the Cartesian product of one dimensional meshes. Two implementations of a conjugate gradient iterative solution method preconditioned by the block symmetric successive overrelaxation matrix (BSSOR-CG) are considered, as well as two implementations of an incomplete preconditioner (IBSSOR-CG). Detailed numerical results for each case for different maximum orders of test functions are presented. A theorem bounding the deterioration of the condition number due to the incomplete preconditioner is also presented. The numerical results indicate that the convergence rate of IBSSOR-CG is only slightly worse than that of BSSOR-CG while the incomplete preconditioner can be more efficient, especially for parallel implementations employing higher order finite elements.