-iip(dgZddlZddlmZddlmZddlmZej ej ej ej dZ ejejejejdZdddZGd d Z ddZdS)_svdpN)aslinearoperator) LinAlgError)_propack)fdFD)LMSMcJeZdZdZdZdZedZedZdS)_AProdz Wrapper class for linear operator The call signature of the __call__ method matches the callback of the PROPACK routines. c t||_dS#t$r*ttj||_YdSwxYwN)rA TypeErrornpasarray)selfrs _/var/www/html/bet.cuttalo.com/ml/venv/lib/python3.11/site-packages/scipy/sparse/linalg/_svdp.py__init__z_AProd.__init__(sR 5%a((DFFF 5 5 5%bjmm44DFFFF 5s0A  A c|dkr!|j||dd<dS|j||dd<dS)Nr)rmatvecrmatvec)rtransamnxys r__call__z_AProd.__call__.sI Q;;6==##AaaaDDD6>>!$$AaaaDDDc|jjSr)rshapers rr$z _AProd.shape4s v|r"c |jjS#t$rD|jt j|jjdjcYSwxYw)Nr)rdtypeAttributeErrorrrzerosr$r%s rr'z _AProd.dtype8s` B6<  B B B6==$&,q/!:!:;;A A A A Bs A AAN) __name__ __module__ __qualname____doc__rr!propertyr$r'r"rrr!sw 555 %%% XBBXBBBr"rr TFMb`?c |td|}|dvrtd|s|dkrtdt|}|jj} t |}t |}nW#t$rJtj tj d|rd }nd }t |}t |}YnwxYw|j \}}|d ks|t||krtd |d |z}|d}t|d z|d z|}||krtd|d|d|rd nd}|rd nd}tj ||d zfd|}tj ||fd|}tj |d|}tj |d|}|q|||dddf<tj tj d|r+|dddfxxd||zz cc<n+ ||dddf<n #t$rtd|wxYw| +tjtj|j} | tj|jdz} |rptj| | | |f|}|||z }|t||z ||krtd|dkrtdn+tj| | | f|}tjt)t+| t)t+|fd} d}!|s|rP||zd|dzzzd|zzdz}"|"t-d |dzzd|zzdz|!t-||zz }"d!|z}#n9||zd|zzd|dzzzdzt-||zd|zdzz}"d|zdz}#tj |"|}$tj |#tj}%tj d |}&tj d tj}'|r"tj ||z|z|}(|$|(|%f})n|$|%f})|dtjtjjdtj"}*|r,|t:||||||||| |||||g|)|| |&|'|*R}+n|||||||| |||||g |)|| |&|'|*R}+|+dkrt=d#|+d$|+dkrt=d%|d&|d'|ddd|f||ddd|fj |fS)(ax Compute the singular value decomposition of a linear operator using PROPACK Parameters ---------- A : array_like, sparse matrix, or LinearOperator Operator for which SVD will be computed. If `A` is a LinearOperator object, it must define both ``matvec`` and ``rmatvec`` methods. k : int Number of singular values/vectors to compute which : {"LM", "SM"} Which singular triplets to compute: - 'LM': compute triplets corresponding to the `k` largest singular values - 'SM': compute triplets corresponding to the `k` smallest singular values `which='SM'` requires `irl_mode=True`. Computes largest singular values by default. irl_mode : bool, optional If `True`, then compute SVD using IRL (implicitly restarted Lanczos) mode. Default is `True`. kmax : int, optional Maximal number of iterations / maximal dimension of the Krylov subspace. Default is ``10 * k``. compute_u : bool, optional If `True` (default) then compute left singular vectors, `u`. compute_v : bool, optional If `True` (default) then compute right singular vectors, `v`. tol : float, optional The desired relative accuracy for computed singular values. If not specified, it will be set based on machine precision. v0 : array_like, optional Starting vector for iterations: must be of length ``A.shape[0]``. If not specified, PROPACK will generate a starting vector. full_output : bool, optional If `True`, then return sigma_bound. Default is `False`. delta : float, optional Level of orthogonality to maintain between Lanczos vectors. Default is set based on machine precision. eta : float, optional Orthogonality cutoff. During reorthogonalization, vectors with component larger than `eta` along the Lanczos vector will be purged. Default is set based on machine precision. anorm : float, optional Estimate of ``||A||``. Default is ``0``. cgs : bool, optional If `True`, reorthogonalization is done using classical Gram-Schmidt. If `False` (default), it is done using modified Gram-Schmidt. elr : bool, optional If `True` (default), then extended local orthogonality is enforced when obtaining singular vectors. min_relgap : float, optional The smallest relative gap allowed between any shift in IRL mode. Default is ``0.001``. Accessed only if ``irl_mode=True``. shifts : int, optional Number of shifts per restart in IRL mode. Default is determined to satisfy ``k <= min(kmax-shifts, m, n)``. Must be >= 0, but choosing 0 might lead to performance degradation. Accessed only if ``irl_mode=True``. maxiter : int, optional Maximum number of restarts in IRL mode. Default is ``1000``. Accessed only if ``irl_mode=True``. rng : `numpy.random.Generator`, optional Pseudorandom number generator state. When `rng` is None, a new `numpy.random.Generator` is created using entropy from the operating system. Types other than `numpy.random.Generator` are passed to `numpy.random.default_rng` to instantiate a ``Generator``. Returns ------- u : ndarray The `k` largest (``which="LM"``) or smallest (``which="SM"``) left singular vectors, ``shape == (A.shape[0], 3)``, returned only if ``compute_u=True``. sigma : ndarray The top `k` singular values, ``shape == (k,)`` vt : ndarray The `k` largest (``which="LM"``) or smallest (``which="SM"``) right singular vectors, ``shape == (3, A.shape[1])``, returned only if ``compute_v=True``. sigma_bound : ndarray the error bounds on the singular values sigma, returned only if ``full_output=True``. Nz:`rng` must be a normalized numpy.random.Generator instance>r r z#`which` must be either 'LM' or 'SM'r z#`which`='SM' requires irl_mode=Truer)r'r r rz.k must be positive and not greater than m or n iz3kmax must be greater than or equal to k, but kmax (z) < k ()r )orderr')sizey?zv0 must be of length g?z0shifts must satisfy k <= min(kmax-shifts, m, n)!zshifts must be >= 0!i  )lowhighr5r'z#An invariant subspace of dimension z was found.zk=z0 singular triplets did not converge within kmax=z iterations)! ValueErrorupperrr'char_lansvd_irl_dict _lansvd_dictKeyErrorr iscomplexobjemptyr$minr)loweruniformsqrtfinfoepsarrayintboolmaxint32isupperintegersiinfoint64uint64_which_converterrconjT),rkwhichirl_modekmax compute_u compute_vv0 full_outputtoldeltaetaanormcgselr min_relgapshiftsmaxiterrngaprodtyp lansvd_irllansvdrrjobujobvuvsigmabnddoptionioptionNBlworkliworkworkiworkdparmiparmzworkworks rng_stateinfos, rrr@sAr {UVVV KKMME L  >??? @ >??? 1IIE + C #%c* c" ### ?28AS111 2 2 CCC%c* c"# ;DAq A1s1ayy==IJJJ |!t q1ua!eT " "D axx + + +&' + + +,, ,  11qD 11qD !TAXc555A !T#S111A HQc 5 5 5E (1Csyy{{ 3 3 3C  z++1+%%!!!Q$ ?28AS111 2 2 0 aaadGGGrCKKQK/// /GGG :AaaadGG : : :8Q8899 9 : } )** {hsmm4' C(E3z:#))++NNN >AXF s4&=!Q'' ' '<== = aZZ344 4(E3.ciikkBBBhDIIDII7sCCCG BIA$' !AdF*Q. QtQwY4'!+RAq \:::4A$47*Q.QUAdFQJ1G1GG4! 8E - - -D HV28 , , ,E HQciikk * * *E HQbh ' ' 'E {{}}QS111eU"e  "(););)?"#2966IIz*5141fa#E3<49<;B<"<$)<+0<2;<<< vdD!Q4eQsAII%I'.I05I7<I>GIII axx C$ C C CEE E  & & & & & &'' ' QQQU8UAaaa!eHMMOO-s 22s,BACC I I=)r TNTTNFrNNrFTr0NNN)__all__numpyrscipy.sparse.linalgr scipy.linalgrrslansvddlansvdclansvdzlansvdrD slansvd_irl dlansvd_irl clansvd_irl zlansvd_irlrCrXrrr/r"rrs )000000$$$$$$                     BBBBBBBB>15JK8<;?c3c3c3c3c3c3r"