pub fn field(
vectors: Vec<[Complex; 3]>,
polarizations: Vec<u8>,
points: Vec<[f64; 3]>,
coefficients: Option<Vec<Complex>>,
helicity: bool,
fixed_vectors: bool,
) -> Result<(DMatrix<Complex>, FieldResidual)>Expand description
Evaluate weighted plane fields, or the full operator when coefficients are absent.
The rows hold the three Cartesian components of each point in turn. With
fixed_vectors, the pullback holds the wavevectors fixed and returns zero wavevector
gradients.
Upstream: sums of treams.special.vpw_M, vpw_N or vpw_A over the plane modes.