function res = blkmulti(A0,b)
res = [];
n = size(A0,2);
for ii = 1:(size(b,1)/n)
    res = [res;A0*b(((ii-1)*n+1):ii*n)];
end
