RBcheck.m

Go to the documentation of this file.

function res = RBcheck(RB,G,NOS)
%function res = RBcheck(RB)
% => Ergebnis RB' * G * RB = eye;
%    check result



    check=[];
    check=RB'*G*RB;  %should be eye(NOS);
    %trace_check=trace(check);    %should be NOS
    %matrix_check=sum(sum(check)); %should be NOS
   err  = max(max(abs(check - eye(NOS))));
   
   if err < 1e-6
       %disp('RB passed');
       res=1;
   else 
       %disp('RB problem');
       res=0;
   end