function avg_top = compute_avg_top( similarity_matrix )

[n1,n2] = size( similarity_matrix); 
if( n1 ~= n2 ) 
  error('Matrix not square!'); 
end 
if( mod(n1,4) ) 
  error('Size not divisible by 4!'); 
end

N = n1/4; 
toprank = 0; 
for i=0:(N-1)
  for j=1:4
[dummy,I] = sort( -similarity_matrix(4*i+j,:));
     dummy = -dummy;
    I = I(1:4); 
t = (I>(4*i)) & (I<=(4*(i+1))) & (dummy(1:4) > dummy(5));
toprank = toprank + sum( t);
  end
end
avg_top = toprank/(4*N);