pub(crate) fn interaction_blocks(
blocks: Vec<DMatrix<Complex>>,
coupling: DMatrix<Complex>,
) -> Result<InteractionResidual>Expand description
interaction for a block-diagonal T given by its particle blocks, which are
never expanded to a dense matrix for the factorization.