def geometric_matrix_multiplication(matrix1, matrix2):
    new_matrix = np.ones((len(matrix1), len(matrix2[0])))
jlen = len(matrix2[0])
klen = len(matrix2)
for i in range(len(matrix1)):
    for j in range(jlen):
        for k in range(klen):
            new_matrix[i, j] *= matrix1[i, k] * matrix2[k, j]
return new_matrix
Geometric Matrix Multiplication