pub fn qmat(
destination: Vec<Mode>,
source: Vec<Mode>,
surface: Surface,
ks: [[Complex; 2]; 2],
zs: [Complex; 2],
radial: Radial,
radial_area_factor: bool,
) -> Result<(DMatrix<Complex>, QmatResidual)>Expand description
Integrate Q between regular inner test waves and regular or singular outer waves.
The surface element of node i is (r rhat - r' theta-hat) r sin(theta) w_i, where
w_i is the quadrature weight. With radial_area_factor false the factor r is
omitted. Each entry sums its nodes in order, so the result does not depend on the
thread count.
Upstream: treams.ebcm.qmat. Differences: treams omits the radial area factor r,
and radial_area_factor = false reproduces treams; treams integrates each entry
with adaptive quadrature, while this sums over the fixed nodes of surface.