Code covered by the BSD License
Andrew Knyazev (view profile)
14 Jan 2000
19 May 2009)
Orthonormalization Relative to matrix A
Watch this File
ORTHA Orthonormalization Relative to matrix A
computes an orthonormal basis Q for the range of X, relative to the scalar product using a positive definite and selfadjoint matrix A. That is, Q'*A*Q = I, the columns of Q span the same space as columns of X, and rank(Q)=rank(X).
This file inspired Subspacea.M.
License change to free software (BSD)