based on how stamps are performed it seems clear the coefficient matrices are Hermitian and thus we can use a faster Cholesky decomposition for solving