[Getfem-users] GMM - orthonormals basis for range of matrix

[Philipp Berger]

To whom it may concern,

sorry for asking something of such simple nature, but I can not figure it out on my own.

Is GMM able to calculate the orthonormal basis for the range of a matrix? Something like "orth()" in Matlab ( http://www.mathworks.de/help/matlab/ref/orth.html ).
I found the gmm_modified_gram_schmidt header, which seems to be what I want but since there is no documentation on the functions I find myself unable to understand how they are meant to be applied.

For example, given a matrix M:
gmm::row_matrix<gmm::wsvector<double>> M1;
M1.resize(4,4);
// Insert some data

And now I would like to get a matrix containg the orthonormal base of M1.

Thank you for your time,

P. Berger

smime.p7s
Description: S/MIME Kryptografische Unterschrift