`i ddlZdZn #e$rdZYnwxYwddlZddlmZddlmZddlmZddlm Z Gdd ej Z d Z dS) NTF)_core)_csc)_data)_utilc\eZdZdZdZddZddZddZd Zdd Z dd Z dd Z dd Z ddZ dS) dia_matrixaSparse matrix with DIAgonal storage. Now it has only one initializer format below: ``dia_matrix((data, offsets))`` Args: arg1: Arguments for the initializer. shape (tuple): Shape of a matrix. Its length must be two. dtype: Data type. It must be an argument of :class:`numpy.dtype`. copy (bool): If ``True``, copies of given arrays are always used. .. seealso:: :class:`scipy.sparse.dia_matrix` diaNFc,trRtj|r3|}|j}|j}|j}|j}d}n;t|tr|\}}|tdntdtj |||}tj|}tj |d|}tj|}|jdkrtd|jdkrtd |jd t#|kr,td |jd t#|fztj|}|dd |ddkrtd ||_||_t)j|stdt-|d t-|df|_dS)NFzexpected a shape argumentz,unrecognized form for dia_matrix constructor)dtypecopyizoffsets array must have rank 1zdata array must have rank 2rzBnumber of diagonals (%d) does not match the number of offsets (%d)z&offset array contains duplicate valuesz(invalid shape (must be a 2-tuple of int))_scipy_availablescipysparseissparsetodiadataoffsetsshaper isinstancetuple ValueErrorcupyarray atleast_2d atleast_1dndimlensortanyrisshapeint_shape) selfarg1rr r xrrsorted_offsetss k/home/jaya/work/projects/VOICE-AGENT/VIET/agent-env/lib/python3.11/site-packages/cupyx/scipy/sparse/_dia.py__init__zdia_matrix.__init__%s   @ 5 5d ; ; @ A6DiGGEGEDD e $ $ @ MD'} !<===>@@ @z$e$777t$$*WCd;;;/'** <1  =>> > 9>>:;; ; :a=CLL ( (:a=#g,,/011 1 7++ 3B3 >!""#5 5 : : < < GEFF F  }U## IGHH H%(mmSq]]2 Tc|r/t||jf|jSt||jf|jS)zReturns a matrix with the same sparsity structure as self, but with different data. By default the structure arrays are copied. r)rrr r)r'rr s r+ _with_datazdia_matrix._with_dataQsQ  FtT\%6%6%8%89LLL LtT\2$*EEE Er-ctstd|j|}|j|}t j||f|jS)a:Returns a copy of the array on host memory. Args: stream (cupy.cuda.Stream): CUDA stream object. If it is given, the copy runs asynchronously. Otherwise, the copy is synchronous. Returns: scipy.sparse.dia_matrix: Copy of the array on host memory. zscipy is not availabler/) r RuntimeErrorrgetrrrrr&)r'streamrrs r+r3zdia_matrix.getZsd  9788 8y}}V$$,""6**|&&gdk&JJJr-c|jS)zcReturns the shape of the matrix. Returns: tuple: Shape of the matrix. )r&)r's r+ get_shapezdia_matrix.get_shapeks {r-c |td|j\}}tjddddddd |j||}t |S) zReturns the number of stored values, including explicit zeros. Args: axis: Not supported yet. Returns: int: The number of stored values. Nz5getnnz over an axis is not implemented for DIA formatzint32 offsets, int32 m, int32 nz int32 nnzz7offsets > 0 ? min(m, n - offsets) : min(m + offsets, n)za + bznnz = a0dia_nnz)NotImplementedErrorrrReductionKernelrr%)r'axismnnnzs r+getnnzzdia_matrix.getnnzssz  %GII Iz10e# -{ E YY0015 aDD3xxr-cT|||S)z3Returns a dense matrix representing the same value.)orderout)tocsctoarray)r'rBrCs r+rEzdia_matrix.toarrays#zz||##%S#999r-c|jjdkr tj|j|jS|j\}}|jj\}}t jdddd||jdddf|||j\}}tj |dzd }tj | d |d|dz<||||dzd<|j |j d d } |jj |j } tj| | |f|j|j S)a{Converts the matrix to Compressed Sparse Column format. Args: copy (bool): If ``False``, it shares data arrays as much as possible. Actually this option is ignored because all arrays in a matrix cannot be shared in dia to csc conversion. Returns: cupyx.scipy.sparse.csc_matrix: Converted matrix. rr zGint32 offset_len, int32 offsets, int32 num_rows, int32 num_cols, T datazint32 row, bool maskz int offset_inds = i % offset_len; row = offset_inds - offsets; mask = (row >= 0 && row < num_rows && offset_inds < num_cols && data != T(0)); cupyx_scipy_sparse_dia_tocscNrr )r<F)r )rr )rsizer csc_matrixrr rElementwiseKernelrrzeroscumsumsumTastype) r'r num_rowsnum_cols num_offsets offset_lenrowmaskindptrindicesrs r+rDzdia_matrix.tocscsX 9>Q  ?4:TZ@@@ @!Z("&)/ Z ,E+ % "  + , ,-7 QQQW8M,4h  K K THqL444$(Ka0@0@$A$Aq*q. !"("4zA~%-&&s&77y{46" 7F #4:TZIII Ir-cN|S)axConverts the matrix to Compressed Sparse Row format. Args: copy (bool): If ``False``, it shares data arrays as much as possible. Actually this option is ignored because all arrays in a matrix cannot be shared in dia to csr conversion. Returns: cupyx.scipy.sparse.csc_matrix: Converted matrix. )rDtocsr)r'r s r+rZzdia_matrix.tocsrszz||!!###r-rc~|j\}}|| ks||kr tjd|jjStj|j|k\}td|t||z|}}|j dkr#tj ||z |jjS|j|d||fS)aReturns the k-th diagonal of the matrix. Args: k (int, optional): Which diagonal to get, corresponding to elements a[i, i+k]. Default: 0 (the main diagonal). Returns: cupy.ndarray : The k-th diagonal. rrG) rremptyrr nonzerormaxminrIrL)r'krowscolsidx first_collast_cols r+diagonalzdia_matrix.diagonalsZ d ::d:aty777 7|DLA-..!!QiiTAXt)<)<8 8q==:h2$)/JJJ JyQ8!3344r-)NNF)T)N)NN)F)r)__name__ __module__ __qualname____doc__formatr,r0r3r6r@rErDrZrfr-r+rrs"F*3*3*3*3XFFFFKKKK"*::::$I$I$I$IL $ $ $ $555555r-rc,t|tS)zChecks if a given matrix is of DIA format. Returns: bool: Returns if ``x`` is :class:`cupyx.scipy.sparse.dia_matrix`. )rr)r)s r+isspmatrix_diarns a $ $$r-) scipy.sparserr ImportErrorrrcupyx.scipy.sparserrr _data_matrixrrnrlr-r+rss ######$$$$$$$$$$$$ A5A5A5A5A5#A5A5A5H%%%%%s