%TEST_AGGS constructs random aggregations

p = randperm(100);

AA = zeros(20,100);

for i = 1:20
   
   p1 = 5*(i-1) + 1;
   p2 = p1 + 4;
   
   AA(i,p(p1:p2)) = ones(1,5);
   
end


AA = inv(diag(sum(AA')))*AA;