o Th?@sTdZddlZddlmZddlmZmZeej fZ ddZ Gdd d Z d d Z dS) z0Indexing mixin for sparse array/matrix classes. N) isintlike)sparrayissparsecCs0t||\}}|jj|j_|jj|j_||fS)a  Same as np.broadcast_arrays(a, b) but old writeability rules. NumPy >= 1.17.0 transitions broadcast_arrays to return read-only arrays. Set writeability explicitly to avoid warnings. Retain the old writeability rules, as our Cython code assumes the old behavior. )npbroadcast_arraysflags writeable)abxyrS/home/air/segue/gpt/backup/venv/lib/python3.10/site-packages/scipy/sparse/_index.py_broadcast_arrays s  rc@seZdZdZddZddZddZdd Zd d Zd d Z ddZ ddZ ddZ ddZ ddZddZddZddZddZd d!Zd"d#Zd$d%Zd&d'Zd(d)Zd*d+Zd,d-Zd.d/Zd0d1Zd2S)3 IndexMixinzS This class provides common dispatching and validation logic for indexing. cCst||\}}t|dkrq|d}t|tjr |jdkr |}t|tr+||}nt|t r6| |}n| |}t|t sB|S|jdkrl|dkrlt|dkr[|j |g||jdSt|dkrl|j |gg||jdS||S|\}}t|trt|tr|||}nt|t r|||}n|jdkr|||}n|jdkr|||}ntdt|t rt|tr|||}nt|t r|t dkr||kr|}n|||}n|jdkr|||}nstdt|tr|||}nct|t r|||}nV|jdkr.|jddkr.|jdks|jddkr.||dddf|}n+t||\}}|j|jkr@td|jdkrS|j t|j|jd}n| ||}t|t s|dksqt|dkrs|jdkrs|St|dkr~d |}||jkr|S||S|j|kr|j!d krt|dkr|jdkr|j"|g|d S|#|S||S|S) Nrrr)shapedtypezindex results in >2 dimensions'number of row and column indices differr)rlilr)$_validate_indiceslen isinstancerndarrayritem INT_TYPES_get_intslice _get_slice _get_arrayr __class__rreshape _get_intXint_get_intXslicendim_get_intXarray IndexError_get_sliceXintcopy_get_sliceXslice_get_sliceXarray_get_arrayXint_get_arrayXslice_get_columnXarrayravelrsize atleast_2d_get_arrayXarrayformat_coo_containertocoo)selfkeyindex new_shapeidxresrowcolrrr __getitem__s                           $   zIndexMixin.__getitem__cCsp||\}}t|dkr|d}t|r|}ntj||jd}t|tr;|j dkr0t d| ||j ddSt|t rst||jd}t|}|dkrd|j dkrd| |d|j ddStj||jd}|j}n|j}|j|krt||j}|j dkr|||dS|\} } t| trt| trtj||jd}|j dkrt d|| | |j ddSt| t rtj| |jddddf} nt| } t| t rtj| |jddddf} | jdkr| dddf} nt| } t| | \} } | j| jkrtdt|rd|jvr"dS| jdkr0| d} | d} |jddj|jdd}|jddkoJ| jddk} |jddkoY| jddk}| sh|jd| jdkrv|sz|jd| jdkszt d ||| | |dStj||jd}|j| jkrt|| j}|j dkrdS|| j}|| | |dS) Nrrrz&Trying to assign a sequence to an itemrF)r+Tzshape mismatch in assignment) rrrtoarrayrasarrayrrrr2 ValueError_set_intflatr rangeindicesrarangesqueeze broadcast_to _set_array _set_intXint atleast_1dr'rr)r7r$ _shape_as_2dsum_duplicates_set_arrayXarray_sparse_set_arrayXarray)r8r9r r:_r< idx_rangeN idx_shaper>r?ij broadcast_row broadcast_colrrr __setitem__|s          $  "       zIndexMixin.__setitem__cs|turtdf|j|jfSt|ts|g}d}g}d}t|D]Y\}}|tur4|dur1td|}q!|dur>||q!t|tsGt |rQ|||d7}q!t ||j}durf||||j7}q!t |rntd|t ||d7}q!|j|tdg}|r|dur||n|d||||d}g} d} gg} t|D]\}}|dur| dqt|trۈ||j| } tt|| } | | | d7} qt |r|j| }| |kr|ksntd|dt|dkr||n|}|| d7} q|jjdkrK|}| |j}|j| |}|j|kr9td |d |d |j|| t| ||} q|j| }|||}|| | | d7} q| |jkrwtd |jd | dt| dkrtfdd| DtfddddDrdddD}d|}t|tdj| } n t| dkr| d}t|j}| d||| |d} t| }dkrtd|dtt| fS)z8Returns two tuples: (index tuple, requested shape tuple)Nrz(an index can only have a single ellipsisrzoIndexing with sparse matrices is not supported except boolean indexing where matrix and index are equal shapes.zindex (z) out of ranger z bool index z has shape z instead of zinvalid index ndim. Array is zD. Index needs Dc3s|]}|VqdSNr).0rV)r:rr )sz/IndexMixin._validate_indices..c3s |] }dj|jkVqdS)rNrr]ix) idx_arraysrrr^*s css|]}t|jVqdSr\)strrr_rrrr^+szLshape mismatch: indexing arrays could not be broadcast together with shapes rz*Only 1D or 2D arrays allowed. Index makes )Ellipsisr r'rrtuple enumerater)appendr_compatible_boolean_indexrrrBextend_shaperrFrGintrkindnonzero _asindicesranyjoinlist)r8r9 ellps_pos index_1st prelim_ndimrVr<r` ellip_slicesrU index_ndim array_indicesMs len_slicerTtmp_ndim mid_shapeshapesmsg arr_index arr_shaper'r)rar:rrs                           zIndexMixin._validate_indicesc Cszt|}Wntttfy}ztd|d}~ww|jdvr%td|jdkr,|S|}||kr:td|| }|dkrc|| krMtd|||usU|j j sY| }||dk|7<|S)zConvert `idx` to a valid index for an axis with a given length. Subclasses that need special validation can override this method. z invalid indexN)rrzIndex dimension must be 1 or 2rindex (%d) out of range) rrBrC TypeError MemoryErrorr)r'r2maxminrowndatar+)r8r<lengthr emax_indxmin_indxrrrrn9s*      zIndexMixin._asindicescCsP|j\}}t|}|| ks||krtd||dkr ||7}||tdS)zGReturn a copy of row i of the matrix, as a (1 x n) row vector. rrN)rrkr)r&r r8rVMrTrrr_getrowW  zIndexMixin._getrowcCsP|j\}}t|}|| ks||krtd||dkr ||7}|td|S)zMReturn a copy of column i of the matrix, as a (m x 1) column vector. rrN)rrkr)r*r rrrr_getcolbrzIndexMixin._getcolcCtr\NotImplementedErrorr8r<rrrrmzIndexMixin._get_intcCrr\rrrrrr!przIndexMixin._get_slicecCrr\rrrrrr"srzIndexMixin._get_arraycCrr\rr8r>r?rrrr%vrzIndexMixin._get_intXintcCrr\rrrrrr(yrzIndexMixin._get_intXarraycCrr\rrrrrr&|rzIndexMixin._get_intXslicecCrr\rrrrrr*rzIndexMixin._get_sliceXintcCrr\rrrrrr,rzIndexMixin._get_sliceXslicecCrr\rrrrrr-rzIndexMixin._get_sliceXarraycCrr\rrrrrr.rzIndexMixin._get_arrayXintcCrr\rrrrrr/rzIndexMixin._get_arrayXslicecCrr\rrrrrr0rzIndexMixin._get_columnXarraycCrr\rrrrrr4rzIndexMixin._get_arrayXarraycCrr\rr8r<r rrrrDrzIndexMixin._set_intcCrr\rrrrrrKrzIndexMixin._set_arraycCrr\rr8r>r?r rrrrLrzIndexMixin._set_intXintcCrr\rrrrrrQrzIndexMixin._set_arrayXarraycCs4tj||jd}t||\}}||||dS)Nr)rrBrArrrQ)r8r>r?r rRrrrrPsz"IndexMixin._set_arrayXarray_sparseN)__name__ __module__ __qualname____doc__r@rZrrnrrrr!r"r%r(r&r*r,r-r.r/r0r4rDrKrLrQrPrrrrrs4_Xe   rcCst|ds6z!tt|d}t|D]}t|trn tt|d}qWdSWn ty0YdSwt|}|j j dkr>|SdS)zICheck for boolean array or array-like. peek before asarray for array-liker'Nr ) hasattrnextiterrFrboolrr asanyarrayrrl)r< desired_ndimr`rRrrrrhs      rh)rnumpyr_sputilsr_baserrrkintegerrrrrhrrrrs