Add support for face elements in interpolator class
Clean up + add error check in SumOfMatrices