def matrix_multiply_elementwise(A, B):
    if A.shape != B.shape:
        raise ShapeError()
    shape = A.shape
    return classof(A, B)._new(shape[0], shape[1],
        lambda i, j: A[i, j]*B[i, j])
Return the Hadamard product (elementwise product) of A and B
dataset/ETHPy150Open sympy/sympy/sympy/matrices/dense.py/matrix_multiply_elementwise