pub fn interaction(
local: DMatrix<Complex>,
coupling: DMatrix<Complex>,
) -> Result<InteractionResidual>Expand description
Solve a multiple-scattering system for the interacting T-matrix
X = (I - T C)⁻¹ T, keeping its LU factorization for the pullback.
local is the block-diagonal T-matrix T of all particles, and coupling the
singular expansion C between them, zero on the diagonal blocks.
Upstream: treams.TMatrix.interaction.solve.