function [diff, tdiff] = read_lambdas(H1, H2,fDim)
  
diff = [1:length(fDim)]*0;
for i=1:length(fDim),
    diff(i)= sum(abs(H1(i,:)-H2(i,:)))/fDim(i);
    diff(i) = diff(i)/sum(abs(H1(i,:)));
end;
tdiff = sum(sum(abs(H1-H2)));
return;
