pub fn operator(
basis: impl Into<MultipoleBasis>,
points: Vec<[f64; 3]>,
ks: [Complex; 2],
helicity: bool,
radial: Radial,
) -> Result<(DMatrix<Complex>, OperatorResidual)>Expand description
Matrix mapping multipole amplitudes to Cartesian samples. Rows pack (sample, component). The residual keeps the basis, wavenumbers and points, without the output or any derivative matrix.
Upstream: treams.efield, which returns the matrix with shape (..., 3, modes).