A fast global nodewise mass matrix inversion framework tailored for sparse block-diagonal systems. (March 2022)