fn coupling_matrix(
particles: usize,
modes_per_particle: usize,
pair: impl Fn(usize, usize) -> Result<Vec<Complex>> + Sync,
) -> Result<DMatrix<Complex>>Expand description
Assemble complete particle columns in parallel, for either coupling values or their directional derivative. Self-couplings are exactly zero in both cases.