function R = radave(A,P)
global SIGMA
[n,m] = size(A);
%R = P'*max(abs(SIGMA * A),[],2)/n;
R = P'*max((SIGMA * A),[],2)/n; % w/o abs it's zero!!!
return

R = 0;
for xi = 1:2^n
  x = i2c(xi,n);
  R = R + max(abs(x*A))*P(xi);
end
R=R/n;
return


function c=i2c(i,n)
m=2;
dims = repmat(m,[1 n]);
c = ind2coord(i,dims)-1;
c = 2*c-1;
return