def matrix_tensor_product(*product):
    if isinstance(product[0], Matrix):
        return _sympy_tensor_product(*product)
    elif isinstance(product[0], numpy_ndarray):
        return _numpy_tensor_product(*product)
    elif isinstance(product[0], scipy_sparse_matrix):
        return _scipy_sparse_tensor_product(*product)
Compute the matrix tensor product of sympy/numpy/scipy.sparse matrices.
dataset/ETHPy150Open sympy/sympy/sympy/physics/quantum/matrixutils.py/matrix_tensor_product