test_orthonormalize.m

Go to the documentation of this file.

function OK = test_orthonormalize
% function performing a test of the orthonormalization routine by
% generating a time-sequence of data and orthogonalizing this.
% returns 1, if test is OK, 0 otherwise

% Bernard Haasdonk

OK = 1;

load('test_rb_lin_evol_data');
model_data = gen_model_data(model);
sim_data = detailed_simulation(model,model_data);
A = model_data.W;
UON = orthonormalize_qr(sim_data.U,A);
B = UON' * A * UON;

E = eye(size(B));
i_zero = find(E==0);
i_nonzero = find(E~=0);

if (find(B>0.9999)~=i_nonzero)
  disp('gram-matrix not the identity!');
  OK = 0;
end;

if (find(B<0.0001)~=i_zero)
  disp('gram-matrix not the identity!');
  OK = 0;
end;

%pcolor(B);

%| \docupdate