fn couple(
local: Vec<DMatrix<Complex>>,
coupling: DMatrix<Complex>,
expansion: ParticleExpansion,
) -> Result<ParticleClusterResidual>Expand description
Solve the interaction of the local matrices through the singular expansion
coupling. The solve owns the only coupling matrix; the expansion pullback needs
just the geometry.