pub fn reflect<A, S1, S2>(w: &ArrayBase<S1, Ix1>, a: &mut ArrayBase<S2, Ix1>) where
    A: Scalar + Lapack,
    S1: Data<Elem = A>,
    S2: DataMut<Elem = A>, 
Expand description

Take a reflection P = I - 2ww^T

Panic

  • if the size of w and a mismaches