a am/@sdZddlZddlmZzeeejfZWne yDeejfZYn0ddZ GdddZ d d Z d d Z d dZdddZddZddZdS)z*Indexing mixin for sparse matrix classes. N) 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 c/Users/vegardjervell/Documents/master/model/venv/lib/python3.9/site-packages/scipy/sparse/_index.py_broadcast_arrays s   r c@seZdZdZddZddZddZdd Zd d Zd d Z ddZ ddZ ddZ ddZ ddZddZddZddZddZd d!Zd"d#Zd$d%Zd&d'Zd(S)) IndexMixinzS This class provides common dispatching and validation logic for indexing. cCs||\}}t|trft|tr.|||St|trD|||S|jdkrZ|||Stdn&t|trt|tr| ||St|tr|tdkr||kr| S| ||S|jdkr| ||Stdn|jdkrt|tr| ||St|tr|||Snvt|tr.| ||St|trDtdnH|jddkr|jdksp|jddkr||dddf|St||\}}|j|jkrtd|jdkr|jt|j|jdS|||S)Nrzindex results in >2 dimensionsr'number of row and column indices differdtype)_validate_indices isinstance INT_TYPES _get_intXintslice_get_intXslicendim_get_intXarray IndexError_get_sliceXintcopy_get_sliceXslice_get_sliceXarray_get_arrayXint_get_arrayXsliceshape_get_columnXarrayZravelr size __class__rZ atleast_2dr_get_arrayXarray)selfkeyrowcolr r r __getitem__ sH                        , zIndexMixin.__getitem__c Csb||\}}t|tr\t|tr\tj||jd}|jdkrDtd||||j ddSt|t rtj | |j ddddf}n t|}t|t rtj | |j ddddf}|jdkr|dddf}n t|}t||\}}|j |j kr tdddlm}||r|jdkr<|d}|d}|j ddkoX|j ddk}|j ddkov|j ddk} |s|j d|j dkr| s|j d|j dkstd|j ddks|j ddkrdS|jdd }|||||n^tj||jd}|j |j kr4t||j }|jdkrDdS||j }||||dS) Nrrz&Trying to assign a sequence to an itemrr) isspmatrixzshape mismatch in assignmentT)r)rrrrasarrayrr" ValueError _set_intXintZflatrZarangeindicesr Z atleast_1drr rbaser*ZtocooZsum_duplicates_set_arrayXarray_sparseZsqueezeZ broadcast_toZreshape_set_arrayXarray) r%r&rr'r(ijr*Z broadcast_rowZ broadcast_colr r r __setitem__KsV  $  "         zIndexMixin.__setitem__cCs|j\}}t|\}}t|rVt|}|| ks8||krDtd||dkrl||7}nt|tsl|||}t|rt|}|| ks||krtd||dkr||7}nt|ts|||}||fS)Nzrow index (%d) out of rangerzcolumn index (%d) out of range)r _unpack_indexrintrrr _asindices)r%r&MNr'r(r r r rs&          zIndexMixin._validate_indicesc Cszt|}Wn4tttfyB}ztd|WYd}~n d}~00|jdvrVtd|jdkrd|S|}||krtd|| }|dkr|| krtd|||us|j j s| }||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)rzIndex dimension must be <= 2rindex (%d) out of range) rr+r, TypeError MemoryErrorrrr"maxminrZowndatar)r%idxlengthreZmax_indxZmin_indxr r r r7s&      zIndexMixin._asindicescCsP|j\}}t|}|| ks$||kr0td||dkr@||7}||tdS)zGReturn a copy of row i of the matrix, as a (1 x n) row vector. r;rN)r r6rrrr%r2r8r9r r r getrows  zIndexMixin.getrowcCsP|j\}}t|}|| ks$||kr0td||dkr@||7}|td|S)zMReturn a copy of column i of the matrix, as a (m x 1) column vector. r;rN)r r6rrrrCr r r getcols  zIndexMixin.getcolcCs tdSNNotImplementedErrorr%r'r(r r r rszIndexMixin._get_intXintcCs tdSrFrGrIr r r rszIndexMixin._get_intXarraycCs tdSrFrGrIr r r rszIndexMixin._get_intXslicecCs tdSrFrGrIr r r rszIndexMixin._get_sliceXintcCs tdSrFrGrIr r r rszIndexMixin._get_sliceXslicecCs tdSrFrGrIr r r rszIndexMixin._get_sliceXarraycCs tdSrFrGrIr r r rszIndexMixin._get_arrayXintcCs tdSrFrGrIr r r rszIndexMixin._get_arrayXslicecCs tdSrFrGrIr r r r!szIndexMixin._get_columnXarraycCs tdSrFrGrIr r r r$szIndexMixin._get_arrayXarraycCs tdSrFrGr%r'r(rr r r r-szIndexMixin._set_intXintcCs tdSrFrGrJr r r r1szIndexMixin._set_arrayXarraycCs4tj||jd}t||\}}||||dS)Nr)rr+Ztoarrayrr r1)r%r'r(r_r r r r0sz"IndexMixin._set_arrayXarray_sparseN)__name__ __module__ __qualname____doc__r)r4rr7rDrErrrrrrrrr!r$r-r1r0r r r r r s(+4  r cCs0ddlm}m}t||tjfr>|jdkr>|jjdkr>| St |}t|t rt |dkrf|\}}qt |dkr|dt d}}qtdnLt|}|dur|t d}}n,|jdkrt|t dfS|jdkr| S||s||rtdt|}t|}|durt|}|dur(t|}||fS) z Parse index. Always return a tuple of the form (row, col). Valid type for row/col is integer, slice, or array of integers. r)spmatrixr*r:rrNzinvalid number of indiceszoIndexing with sparse matrices is not supported except boolean indexing where matrix and index are equal shapes.)r/rPr*rrZndarrayrrkindZnonzero_check_ellipsistuplelenrr_compatible_boolean_index_boolean_index_to_array)indexrPr*r'r(r@Zbool_rowZbool_colr r r r5s@          r5cCs |turtdtdfSt|ts&|St|D]\}}|tur.|}qLq.|St|dkrhtdtdfSt|dkr|dkr|dturtdtdfStd|dfS|dtdfSg}||ddD]}|tur||q|t|}tdd|}|d|tdf|t|S)z6Process indices with Ellipsis. Returns modified index.Nrr:r)EllipsisrrrS enumeraterTappendr>)rWr3vZfirst_ellipsistailndZnslicer r r rR$s0      rRcCst|}|jjdkr|SdS)z8Returns a compatible array if elements are boolean. rN)rZ asanyarrayrrQr@r r r _maybe_bool_ndarrayKs  r_r:cCsP|dkr dSztt|d}Wnty2YdS0t|trBdSt||dS)zQReturns True if first element of the incompatible array type is boolean. rNT)nextiterr<rbool_first_element_bool)r@Zmax_dimfirstr r r rcTs  rccCst|dst|rt|SdS)ztReturns a boolean index array that can be converted to integer array. Returns None if no such array exists. rN)hasattrrcr_r^r r r rUcsrUcCs |jdkrtdt|dS)Nrzinvalid index shaper)rrrwherer^r r r rVms rV)r:)rOZnumpyrZsputilsrr6longintegerr NameErrorr r r5rRr_rcrUrVr r r r s  [.'