Changes

Numpy linalg

218 bytes added, 12:18, 19 June 2014

a2_e = np.einsum('...ij,...jk', a2, a2)

print a2_e

# Expand one axis in start, and tile up 2 times.

a3 = np.tile(a2[None,:], (2, 1, 1, 1))

print "a3 shape", a3.shape

print "einsum dot product over higher dimensions"

a3_e = np.einsum('...ij,...jk', a3, a3)

print a3_e

</source>

1,660

edits