Numerical experiences of solving elasticity systems by PCG methods

In solving systems of linear equations derived from FEM structural analysis by the conjugate gradient method, an efficient way for preparing effective incomplete factorization preconditioner was published by Saint-Georges [1]. In [2,3], we prepared several variants of the cited method due to Saint-Georges [1], and some of the relevant PCG-solvers were found more efficient than the original ones [1]. Here, we prepare some additional variants for producing good enough preconditioners by applying our spiral ordering [2] and execute some incomplete factorizations (MIC, RIC) with the level of fill-in dropping strategy which was found the best in [3]. In the comparison of the considered solvers, we analyze our computer results when systems from elastic frame structures with 3D beam elements were solved.