pub fn matvec(
transa: Transpose,
m: usize,
n: usize,
a: &[f64],
x: &[f64],
y: &mut [f64],
) -> Result<(), SolversError>Expand description
Matrix–vector product y := op(A) · x.
op(A) is m×n; a holds op==No ? m×n : n×m (always m*n elements),
x is length n, y is length m (overwritten). Zero-heap; the dynamic
analogue of super::Matrix4x4::multiply_vector, and the shape the GPU
coop_gemv decode kernel computes per output row.