o EOf3@sdZddlmZddlmZddlZddlmZer ddl mZ e ej fZ ddZGd d d Zdd dZddZddZdddZddZddZdS)z0Indexing mixin for sparse array/matrix classes. ) annotations) TYPE_CHECKINGN) isintlikecCs0t||\}}|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. )npZbroadcast_arraysflagsZ writeable)abxyr 3lib/python3.10/site-packages/scipy/sparse/_index.py_broadcast_arrayss  rc@seZdZdZddZddZddZdd Zd d Zd d Z ddZ ddZ ddZ ddZ ddZddZddZddZddZd d!Zd"d#Zd$d%Zd&d'Zd(d)Zd*S)+ IndexMixinzS This class provides common dispatching and validation logic for indexing. cCs"ddlm}t||rtddS)zWe do not currently support 1D sparse arrays. This function is called each time that a 1D array would result, raising an error instead. Once 1D sparse arrays are implemented, it should be removed. r)sparrayzcWe have not yet implemented 1D sparse slices; please index using explicit indices, e.g. `x[:, [0]]`N)Z scipy.sparser isinstanceNotImplementedError)selfrr r r _raise_on_1d_array_slice#s  z#IndexMixin._raise_on_1d_array_slicecCs||\}}t|trDt|tr|||St|tr&||||S|jdkr5||||S|jdkr@|||St dt|trt|trX|| ||St|trq|tdkrk||krk| S| ||S|jdkr|| ||St d|jdkrt|tr||||St|tr|||Sn5t|tr|||St|trt d|jddkr|jdks|jddkr||dddf|St||\}}|j|jkrt d|jdkr|jt|j|jdS|||S)Nrzindex results in >2 dimensionsr'number of row and column indices differdtype)_validate_indicesr INT_TYPES _get_intXintslicer_get_intXslicendim_get_intXarray IndexError_get_sliceXintcopy_get_sliceXslice_get_sliceXarray_get_arrayXint_get_arrayXsliceshape_get_columnXarrayZravelrsize __class__rZ atleast_2dr_get_arrayXarray)rkeyrowcolr r r __getitem__3sV                        &   zIndexMixin.__getitem__c CsN||\}}t|tr.t|tr.tj||jd}|jdkr"td||||j ddSt|t rEtj | |j ddddf}nt|}t|t rntj | |j ddddf}|jdkrm|dddf}nt|}t||\}}|j |j krtdddlm}||r|jdkr|d}|d}|j ddko|j ddk}|j ddko|j ddk} |s|j d|j dkr| s|j d|j dkstd|j ddks|j ddkrdS|jdd }|||||dStj||jd}|j |j krt||j }|jdkrdS||j }||||dS) Nrrz&Trying to assign a sequence to an itemrrissparsezshape mismatch in assignmentT)r")rrrrasarrayrr) ValueError _set_intXintZflatrZarangeindicesr'Z atleast_1drrr _baser1ZtocooZsum_duplicates_set_arrayXarray_sparseZsqueezeZ broadcast_toZreshape_set_arrayXarray) rr,r r-r.ijr1Z broadcast_rowZ broadcast_colr r r __setitem__esV  $  "        zIndexMixin.__setitem__c Cspddlm}t||tjfr*|jdkr*|jjdkr*|j|jkr#t d| \}}nt |\}}|j\}}ddd}t |rZt |}|| ksK||krQt d||dkrY||7}nt|}duri|||d}n t|tst|||}t |rt |}|| ks||krt d||dkr||7}||fSt|} dur|| |d}||fSt|ts|||}||fS)Nr)_spbaserr z.boolean index shape does not match array shapeidxnpt.NDArray[np.bool_] axis_sizeint axis_namestrreturnnpt.NDArray[np.int_]cSs2t||krtd|dt|d|t|S)Nzboolean z index has incorrect length: z instead of )lenr _boolean_index_to_array)r=r?rAr r r _validate_bool_idxs z8IndexMixin._validate_indices.._validate_bool_idxzrow index (%d) out of rangerr-zcolumn index (%d) out of rangecolumn)r=r>r?r@rArBrCrD)r6r<rrZndarrayrrkindr'r nonzero _unpack_indexrr@_compatible_boolean_indexr _asindices) rr,r<r-r.MNrGZbool_rowZbool_colr r r rsB            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) rr2r3 TypeError MemoryErrorr rr)maxminrZowndatar")rr=lengthr eZmax_indxZmin_indxr r r rMs*      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. rPrN)r'r@r rrrr9rNrOr r r _getrow  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. rPrN)r'r@r r!rrWr r r _getcolrYzIndexMixin._getcolcCtNrrr-r.r r r rzIndexMixin._get_intXintcCr[r\r]r^r r r rr_zIndexMixin._get_intXarraycCr[r\r]r^r r r rr_zIndexMixin._get_intXslicecCr[r\r]r^r r r r!r_zIndexMixin._get_sliceXintcCr[r\r]r^r r r r# r_zIndexMixin._get_sliceXslicecCr[r\r]r^r r r r$ r_zIndexMixin._get_sliceXarraycCr[r\r]r^r r r r%r_zIndexMixin._get_arrayXintcCr[r\r]r^r r r r&r_zIndexMixin._get_arrayXslicecCr[r\r]r^r r r r(r_zIndexMixin._get_columnXarraycCr[r\r]r^r r r r+r_zIndexMixin._get_arrayXarraycCr[r\r]rr-r.r r r r r4r_zIndexMixin._set_intXintcCr[r\r]r`r r r r8r_zIndexMixin._set_arrayXarraycCs4tj||jd}t||\}}||||dS)Nr)rr2Ztoarrayrrr8)rr-r.r _r r r r7!sz"IndexMixin._set_arrayXarray_sparseN)__name__ __module__ __qualname____doc__rr/r;rrMrXrZrrrr!r#r$r%r&r(r+r4r8r7r r r r rs,240   rrCctuple[int | slice | npt.NDArray[np.bool_ | np.int_], int | slice | npt.NDArray[np.bool_ | np.int_]]cCst|}t|tr(t|dkr|\}}n8t|dkr$|dtd}}n(tdt|}|dur8|td}}n|jdkrC|tdfS|jdkrL|Sddl m }||sZ||r^td||fS)z Parse index. Always return a tuple of the form (row, col). Valid type for row/col is integer, slice, array of bool, or array of integers. rrrNzinvalid number of indicesr0zoIndexing with sparse matrices is not supported except boolean indexing where matrix and index are equal shapes.) _check_ellipsisrtuplerErr rLrrJr6r1)indexr-r.r=r1r r r rK(s(        rKcCs|tur tdtdfSt|ts|Sddt|D}|s |St|dkr*td|\}tddt|}|d|tdf|||ddS)z6Process indices with Ellipsis. Returns modified index.NcSsg|] \}}|tur|qSr )Ellipsis).0r9vr r r Wsz#_check_ellipsis..rz0an index can only have a single ellipsis ('...')r)rjrrrh enumeraterEr rS)riZellipsis_indicesr9Z num_slicesr r r rgNs  *rgcCst|}|jjdkr |SdS)z8Returns a compatible array if elements are boolean. r N)rZ asanyarrayrrIr=r r r _maybe_bool_ndarraycs  rqrcCsP|dkrdSz tt|d}Wn tyYdSwt|tr!dSt||dS)zQReturns True if first element of the incompatible array type is boolean. rNT)nextiterrQrbool_first_element_bool)r=Zmax_dimfirstr r r ruls  rucCst|ds t|r t|SdS)ztReturns a boolean index array that can be converted to integer array. Returns None if no such array exists. rN)hasattrrurqrpr r r rL{srLcCs |jdkr tdt|dS)Nrzinvalid index shaper)rr rwhererpr r r rFs rF)rCrf)r)reZ __future__rtypingrZnumpyrZ_sputilsrZ numpy.typingZnptr@ZintegerrrrrKrgrqrurLrFr r r r s$      &