1 function onvec = orthonormalize(vec,A,epsilon,method)
2 %
function onvec = orthonormalize(vec[,A,epsilon,method])
4 % orthonormalization of vector a wrt
'A'-scalarproduct. Epsilon can
5 % be set and as methods ''gram-schmidt' (default) and 'qr' are supported.
7 % Bernard Haasdonk 13.6.2002
10 method = 'gram-schmidt';
13 if nargin < 3 || isempty(epsilon)
14 epsilon = 1e-7; % => for nonlinear evolution this value is required.
17 if nargin<2 || isempty(A)
23 warning('gram-schmidt used in orthonormalize might be inaccurate!');
24 onvec = orthonormalize_gram_schmidt(vec,A,epsilon);
26 onvec = orthonormalize_qr(vec,A,epsilon);
28 error('orthonormalization method unknown');