o
    h|4                     @   s  d dl Z d dlZd dlmZ d dlZd dlm	  m
Z dZdd eD Zeg dg dg dgZed dgd	d ggZejg d
g dg dg dgedZdd eD Zdd eD Zdd eD Ze jjdeedZe jjdeedZe jjdeedZedd Zedd Zedd Zedd Zedd Zed d! Z ed"d# Z!ed$d% Z"ed&d' Z#ed(d) Z$ed*d+ Z%ed,d- Z&ed.d/ Z'ed0d1 Z(ed2d3 Z)ed4d5 Z*ed6d7 Z+ed8d9 Z,ed:d; Z-ed<d= Z.ed>d? Z/ed@dA Z0edBdC Z1edDdE Z2dFdG Z3edHdI Z4ee jdJg dKdLdM Z5ee jdJdNdOgdPdQ Z6edRdS Z7edTdU Z8edVdW Z9dXdY Z:dZd[ Z;ed\d] Z<e jd^ej	j=d_d`gfej	j>d_d`gfej	j?dadbgfge jdcej@ejAgddde ZBdfdg ZCdhdi ZDdjdk ZEdldm ZFdndo ZGdpdq ZHdrds ZIdtdu ZJdvdw ZKe jdxej	jLej	jMfdydz ZNd{d| ZOd}d~ ZPdd ZQdd ZRe jddej	jSfdej	jTfdej	jUfdej	jVfdej	jWfdej	jXfdej	jYffdd ZZdS )    N)bsrcoocsccsrdiadoklilc                 C   s   g | ]}t tj| d qS )_array)getattrscipysparse).0T r   \/var/www/html/scripts/venv/lib/python3.10/site-packages/scipy/sparse/tests/test_array_api.py
<listcomp>	   s    r   )r         r   )r   r   r      )r      r   r   r   r   )r   r   r   r   )r   r   r   r   )r   r   r   r   )r   r   r   r   )dtypec                 C      g | ]}|t qS r   Ar   sparrayr   r   r   r           c                 C   r   r   Br   r   r   r   r   !   r   c                 C   r   r   Xr   r   r   r   r   "   r   r   )idsr   r    c                 C   sL   t | jddtjrJ d| jddjdksJ | jddjdks$J d S )Nr   axisExpected array, got matrix)r   r   )r   )
isinstancesumnpmatrixshaper   r   r   r   test_sum/   s
   r*   c                 C   s    t | jddtjrJ dd S )Nr   r"   r$   )r%   meanr'   r(   r   r   r   r   	test_mean7   s   r,   c                 C   s   t | drt| jddtjrJ dt | dr&t| jddtjr&J dt | dr9t| jddtjr9J dt | drLt| jddtjrNJ dd S d S )Nminr   r"   r$   maxargminargmax)hasattrr%   r-   r'   r(   r.   r/   r0   r   r   r   r   test_min_max=   s$   



r2   c                 C   s   t |  tjrJ dd S Nr$   )r%   todenser'   r(   r   r   r   r   test_todenseN   s   r5   c                 C   sn  | j jd d dv rd S tt | dd d f  W d    n1 s$w   Y  tt | d d df  W d    n1 sAw   Y  tt | dddgf  W d    n1 s^w   Y  tt | ddgdf  W d    n1 s{w   Y  | dg jsJ d| dddggf jsJ d| ddggdf jsJ d| d d ddgf jsJ dd S )	Nr   )r   r   r   r   r   r   z(Expected sparse array, got sparse matrixz"Expected ndarray, got sparse arrayz)Expected sparse array, got something else)	__class____name__pytestraisesNotImplementedError	_is_arrayr   r   r   r   test_indexingT   s$   "r<   c                 C   s*   t j| j}t| | t jrJ dd S r3   )r'   randomr)   r%   r(   )r   r    r   r   r   test_dense_additionk   s   r>   c                 C   s   | |  j s	J dd S r3   r;   r   r   r   r   test_sparse_additionq      r@   c                 C   s(   t | |   | d ksJ d S Nr   )r'   allr4   powerr   r   r   r   test_elementwise_mulv   s   (rE   c                 C   s   t t d |   W d    n1 sw   Y  t t tdtjt	d
dd  W d    n1 s;w   Y  td|  |  d ksOJ t|  |  |  d ks`J d S )Nr      r   )r8   r9   	TypeError
ValueErrorr'   eyer   r   	csr_arrayarangereshaperC   r4   r   r   r   r   test_elementwise_rmul{   s   
&&rM   c                 C   s,   t | | j  | | j ksJ d S N)r'   rC   r   r4   dotr   r   r   r   test_matmul   s   ,rP   c                 C   s(   | d j s	J d| d j sJ dd S )Nr   r$   r   r?   r   r   r   r   test_pow   s   rQ   c                 C   s   t | |  tjs
J d S rN   )r%   r'   ndarrayr   r   r   r   test_sparse_divide   s   rS   c                 C   sB   t t | |   jsJ W d    d S 1 sw   Y  d S rN   )r8   warnsRuntimeWarningr4   r;   r   r   r   r   test_sparse_dense_divide   s   "rV   c                 C   s   | d j s	J dd S )Nr   r$   r?   r   r   r   r   test_dense_divide   rA   rW   c                 C   8   t tj | j W d    d S 1 sw   Y  d S rN   )r8   rT   r'   VisibleDeprecationWarningr   r   r   r   r   test_no_A_attr      "rZ   c                 C   rX   rN   )r8   rT   r'   rY   Hr   r   r   r   test_no_H_attr   r[   r]   c                 C   s$   |  djsJ | djsJ d S )Nr   )_getcolr;   _getrowr   r   r   r   test_getrow_getcol   s   r`   c                 C   s*   t | }t|g d | g d  d S )Nr   r   r   r   )splaaslinearoperatornptassert_allclose)r   Lr   r   r   test_as_linearoperator   s   
 rg   c                 C   sL   | j jd d dkrd S t| }|jsJ t| tj	|   d S Nr   r   )
r6   r7   rb   invr;   rd   re   r4   r'   linalgr   Cr   r   r   test_inv   s
   

 rm   c                 C   sV   | j jd d dkrd S tj| }t| }|jsJ t	|
 t|
  d S rh   )r6   r7   r   r   
csc_matrixrb   expmr;   rd   re   r4   )r   Bmatrl   r   r   r   	test_expm   s   

rq   c              	   C   sH   | j jd d dkrd S tt| tddgt| ddg  d S Nr   r   r   r   )	r6   r7   rd   re   rb   expm_multiplyr'   arrayro   r   r   r   r   test_expm_multiply   s   ru   c                 C   s&   t | }t|tj|   d S rN   )rb   normrd   re   r'   rj   r4   )r   rl   r   r   r   	test_norm   s   
rw   c                 C   s(   t | }t|tj|  d d S )Nr   )rb   
onenormestrd   re   r'   rj   rv   r4   rk   r   r   r   test_onenormest   s   
ry   c                 C   sF   | j jd d dvrd S tt| ddgtj| 	 ddg d S )Nr   )r   r   r   r   )
r6   r7   rd   re   rb   spsolver'   rj   solver4   r   r   r   r   test_spsolve   s   r|   c                  C   s8   t jg dg dg dg dg} t| g d d S N)r   r   r   r   )r   r   r   r   )r   r   r   r   )r   r   r   r   ra   )r   r   rJ   rb   spsolve_triangularr   r   r   r   test_spsolve_triangular   s   r   c                 C   sR   | j jd d dkrd S t| }t|tddgtj	| 
 ddg d S rr   )r6   r7   rb   
factorizedrd   re   r'   rt   rj   r{   r4   )r   LUr   r   r   test_factorized  s   
r   solver)
bicgbicgstabcgcgsgmreslgmresminresqmrgcrotmktfqmrc                 C   s`   |dkri }nddi}t t|| tddgfi |\}}|dks$J tj|ddgdd d S )	Nr   atolgh㈵>r   r   r   g?r   )r
   rb   r'   rt   rd   re   )r   r   kwargsxinfor   r   r   test_solvers  s   &r   lsqrlsmrc                 C   s0   t t|| g d^}}t| | g d d S )Nr   r   r   )r
   rb   rd   re   )r   r   r   _r   r   r   test_lstsqr  s   r   c                 C   s.   t j| dd\}}t| | |d |  d S Nr   kr   )rb   eigsrd   re   r    evr   r   r   	test_eigs)  s
   
r   c                 C   s8   | | j  } tj| dd\}}t| | |d |  d S r   )r   rb   eigshrd   re   r   r   r   r   
test_eigsh2  s   

r   c                 C   sZ   t j| dd\}}}tj|  \}}}t|}t|d d }tj||dd d S )Nr   r   gMbP?r   )	rb   svdsr'   rj   svdr4   sortrd   re   )r    usvhu2s2vh2r   r   r   	test_svds<  s
   
r   c                  C   R   t jg dg dg dg dg} t| }t|t	g dg d d S r}   )
r   r   	csc_arrayrb   splurd   re   r{   r'   rt   r    r   r   r   r   	test_spluE     
$r   c                  C   r   r}   )
r   r   r   rb   spilurd   re   r{   r'   rt   r   r   r   r   
test_spiluP  r   r   c                 C   s    t | d  |  d  d S rB   )rd   assert_equalr4   r   r   r   r   test_power_operator[  s    r   zcls,indices_attrsindicesindptrrowcolexpected_dtypec           	      C   s  t jtddd}|j|j||j	|ff}| |}|D ]}t
||j|ks.J q"| |dd}|D ]}t
||j|ksCJ q7t| t jjjr| }|j|j||j|f}| |}|D ]}t
||j|kspJ qd| |dd}|D ]}t
||j|ksJ qyd S d S )N	   r   )r   r   )r)   )r   r   	coo_arrayr'   rK   rL   datar   astyper   r
   r   
issubclass_compressed
_cs_matrixtocsrr   r   )	clsindices_attrsr   input_array	coo_tupleresultattrinput_array_csr	csr_tupler   r   r   test_index_dtype_compresseda  s4   



r   c                  C   s   t jg d} | jrJ d S )N)r   r   r   )r   r   diagsr;   mr   r   r   test_default_is_matrix_diags  s   r   c                  C      t jd} | jrJ d S Nr   )r   r   rI   r;   r   r   r   r   test_default_is_matrix_eye     r   c                  C   s$   t jg dddd} | jrJ d S )Nr   r   r   )r   r   spdiagsr;   r   r   r   r   test_default_is_matrix_spdiags  s   r   c                  C   r   r   )r   r   identityr;   r   r   r   r   test_default_is_matrix_identity  r   r   c                  C   s@   t jtddgddggtddgddgg} | jrJ d S )Nr   r   r   r   r   r   kronr'   rt   r;   r   r   r   r   !test_default_is_matrix_kron_dense     (r   c                  C   s@   t jtddgddggtddgddgg} | jrJ d S )Nr   r   r   r   r   r   r   r   r   r   "test_default_is_matrix_kron_sparse  r   r   c                  C   s@   t jtddgddggtddgddgg} | jrJ d S )Nr   r   )r   r   kronsumr'   rt   r;   r   r   r   r   test_default_is_matrix_kronsum  r   r   c                  C      t jdd} | jrJ d S r   )r   r   r=   r;   r   r   r   r   test_default_is_matrix_random     r   c                  C   r   r   )r   r   randr;   r   r   r   r   test_default_is_matrix_rand  r   r   fnc                 C   sD   t jtd}t jddgddgg}| ||g}|jr J dS )zbSame idea as `test_default_construction_fn_matrices`, but for the
    stacking creation functions.r   r   r   N)r   r   
coo_matrixr'   rI   r;   )r   r   r   r   r   r   r   test_default_is_matrix_stacks  s   r   c                  C   sx   t jtd} t jdgdgg}t jdgg}t j| ||f}|jr)J t j| dgd|gg}|jr:J dS )z]Same idea as `test_default_construction_fn_matrices`, but for the block
    creation functionr   r   r   N)r   r   r   r'   rI   
block_diagr;   bmat)r   r   rl   r   r   r   r   ,test_blocks_default_construction_fn_matrices  s   
r   c               	   C   sz   t D ]8} ttj|  d}|ddgg}|j| ksJ |j| ks"J tt d|_W d    n1 s5w   Y  qd S )Nr	   r   r   qqq)	sparray_typesr
   r   r   format_formatr8   r9   AttributeError)fmtarr_clsMr   r   r   test_format_property  s   r   c                  C   sx   t jd} t j| }| jrJ |jsJ t j|sJ t j| s&J t j| r0J t j|  r:J d S r   )r   r   rI   rJ   r;   issparser4   r   ar   r   r   test_issparse     

r   c                  C   sx   t jd} t j| }| jrJ |jsJ t j|rJ t j| s&J t j| r0J t j|  r:J d S r   )r   r   rI   rJ   r;   
isspmatrixr4   r   r   r   r   test_isspmatrix  r   r   )r   r   r   r   r   r   r   r   r   c                 C   sr   t jjd| d}t j|| }|jrJ |jsJ ||r!J ||s'J || r/J || r7J d S )Nr   )r   )r   r   rI   rJ   asformatr;   r4   )r   r   r   r   r   r   r   test_isspmatrix_format  s   

r   )[r8   numpyr'   numpy.testingtestingrd   scipy.sparser   scipy.sparse.linalgr   rj   rb   r   sparray_classesrt   r   r   floatr    sparrayssquare_sparrayseig_sparraysmarkparametrizeparametrize_sparraysparametrize_square_sparraysparametrize_eig_sparraysr*   r,   r2   r5   r<   r>   r@   rE   rM   rP   rQ   rS   rV   rW   rZ   r]   r`   rg   rm   rq   ru   rw   ry   r|   r   r   r   r   r   r   r   r   r   r   rJ   r   r   int64int32r   r   r   r   r   r   r   r   r   r   hstackvstackr   r   r   r   r   isspmatrix_bsrisspmatrix_cooisspmatrix_cscisspmatrix_csrisspmatrix_diaisspmatrix_dokisspmatrix_lilr   r   r   r   r   <module>   s2   






























	

#
	







