Recent advances have shown that block recursive approximate factorizations provide among the best preconditioners for solving elasticity equations on large 3D meshes such as arising in geomechanical applications. However, to get full efficiency of these methods, a wide range of implementation...