view general/numerical/orthogonalise.m @ 4:e44f49929e56

Adding reorganised general toolbox, now in several subdirectories.

function B1=orthogonalise(B),
% orthogonalise - Orthogonalise a basis matrix
%
% orthogonalise :: [[N,M]] -> [[N,M]].
%
% Works using SVD.

[U,S,V] = svd(B,0);
B1 = U*V';

% alternative method, seems to be slower
% B1 = B*real((B'*B)^(-0.5));