Re: [Numpy-discussion] Kronecker product
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)
Charles G Waldman schrieb:
What is the difference between kron(z,eye) and kron(yt,eye) ? Nils from Numeric import * def kron(a,b): o = outerproduct(a,b) o.shape = a.shape + b.shape return concatenate(concatenate(o, axis=1), axis=1) y = array(([5,6,7],[8,9,10])) print 'y' print print y print z = array(([5,8],[6,9],[7,10])) print 'z' print print z print yt = transpose(y) print 'yt' print print yt print eye = identity(2) a = kron(y,eye) print print 'kron(y,eye)' print print a print b = kron(z,eye) print print 'kron(z,eye)' print print b print c = kron(yt,eye) print print c print
Charles G Waldman schrieb:
What is the difference between kron(z,eye) and kron(yt,eye) ? Nils from Numeric import * def kron(a,b): o = outerproduct(a,b) o.shape = a.shape + b.shape return concatenate(concatenate(o, axis=1), axis=1) y = array(([5,6,7],[8,9,10])) print 'y' print print y print z = array(([5,8],[6,9],[7,10])) print 'z' print print z print yt = transpose(y) print 'yt' print print yt print eye = identity(2) a = kron(y,eye) print print 'kron(y,eye)' print print a print b = kron(z,eye) print print 'kron(z,eye)' print print b print c = kron(yt,eye) print print c print
participants (2)
-
Charles G Waldman -
Nils Wagner