A= ones(b*c,1)*exp(1i*pi/4); % initial_value
d= 1;
while (1)
    for n=1:length(A)
        A(n)= A(n)*exp(1i*pi/2);
    end
    code=reshape(A(n),b,c);
    r = rank(code);
    if r==c
       codebook(:,:,d)=code
       d=d+1;
    end
    if ( norm(A-A)==0 )
       break
    end
end