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>,
Take a reflection P = I - 2ww^T
P = I - 2ww^T
w
a