On Wednesday, October 5, 2016 at 3:36:06 PM UTC-4, John Gibson wrote: > > I would do this via QR decomposition. Produce an orthogonal basis for span > v1, ..., vn via QR, then see if v changes when you project onto the space > spanned by the basis. E.g. > > V = [v1 v2 ... vn] > (Q,R) = qr(V) >
Even better would be to use F = qr(V) Q = F[:Q] which uses an implicit representation of Q (e.g. a compact WY representation of the Householder reflectors).
