pub fn reflect<A>(w: &ArrayRef<A, Ix1>, a: &mut ArrayRef<A, Ix1>)where A: Scalar + Lapack,
Take a reflection P = I - 2ww^T
P = I - 2ww^T
w
a