[Numpy-discussion] Kronecker product

[Charles G Waldman]

Nils Wagner writes:
 > 
 > How can I rearrange the result of outerproduct to the result of
 > Kronecker product with numpy ?
 > 

def kron(a,b):
    o = outerproduct(a,b)
    o.shape = a.shape + b.shape
    return concatenate(concatenate(o, axis=1), axis=1)