Portable efficient solvers for adaptive finite element simulations of elastostatics in two and three dimensions

Adaptive finite element methods (FEM), generate linear equation systems that require dynamic and irregular patterns of data storage, access and computation, making their parallelization very difficult. Moreover, constantly evolving computer architectures often require new algorithms altogether. We describe here several solvers for solving such systems efficiently in two and three dimensions on multiple parallel architectures.