Column-wise AND operation in all column combinations of two matrices

Given matrix a of size MxN and matrix b of size MxL I want to get a matrix C of size M x(N*L) which contains the column-wise logical operation AND of all combinations of columns of matrix a and matrix b.

a = [1     0     0
     1     0     0
     1     0     0
     1     0     0
     0     1     0
     0     1     0
     0     1     0
     0     1     0
     0     0     1
     0     0     1
     0     0     1
     0     0     1];

 b= [0     0     1
     0     1     0
     1     0     0
     0     0     1
     0     1     0
     1     0     0
     0     0     1
     0     1     0
     1     0     0
     0     0     1
     0     1     0
     1     0     0];

C=[];
for ii=1:size(a, 2)
    for jj=1:size(b, 2)        
        C = [C b(:,jj) & a(:,ii)];
    end
end

I am searching for an alternative syntax of the above loops. Thank you!

[Jos (10584)]

A matrix of size MxN, means M rows and N columns, rather than the other way around ...

Answer 1

[Iain]

This should do it:

 C = bsxfun(@and, a, reshape(b,[12,1,3]);
 theresult = C(down,column_a,column_b);

Answer 2

[Azzi Abdelmalek]

   out=reshape(repmat(a,size(b,2),1),size(a,1),[]) & repmat(b,1,size(a,2))

Answer 3

[Jos (10584)]

% Small example: two logical arrays
A = rand(3,3) > .25
B = rand(2,3) > .25

% general engine
[ib,ia] = ndgrid(1:size(B,1), 1:size(A,1))
C = A(ia(:),:) & B(ib(:),:)

Answer 4

[Giorgos Papakonstantinou]

Thank you all for your replies!!