U Zʹj @szdZddlmZddlmZddlmZddlmZddl m Z m Z ddl m Z dd d d d d dddddddg Ze dZesGdddZnejZe dddgidZesGdd d ZdZn ejZejZdd Zdd Zdd Zd d Zd!dZd"dZd#dZd$dZd%dZd&d'Zd(d)Z d*d+Z!d,dZ"d-d.Z#d/d0Z$d1d2Z%d3d4Z&d5d6Z'd7dZ(d8d9Z)d:d;Z*dzrrr_scipy_sparse_tensor_products rBcGsFt|dtrt|St|dtr,t|St|dtrBt|SdS)zGCompute the matrix tensor product of sympy/numpy/scipy.sparse matrices.rN)r(rr<r rAr rB)r?rrrrs cCststttj|ddS)znumpy version of complex eye.r#r$)r%r&r)rnrrr _numpy_eyesrEcCststtj||ddS)z$scipy.sparse version of complex eye.r#r$)r r&rrCrrr_scipy_sparse_eyesrFcKsL|dd}|dkrt|S|dkr,t|S|dkrsh