in tutorials/LinearAlgebra/testing.py [0:0]
def tensor_product_ref(a, b):
n = len(a)
m = len(a[0])
k = len(b)
l = len(b[0])
ans = []
for i in range(n * k):
row = [0] * m * l
ans.append(row)
for i in range(n):
for j in range(m):
for o in range(k):
for p in range(l):
ans[i * k + o][j * l + p] = a[i][j] * b[o][p]
return ans