Cpj*.S/rSSKrSSKJr SSKJr SSKJr SSK J r \ R\ R\ R\ RS.r\ R \ R"\ R$\ R&S.rSSS .r"S S 5rS S jrg)_svdpN) HAS_ILP64)aslinearoperator) LinAlgError)_propack)fdFD)LMSMcD\rSrSrSrSrSr\S5r\S5r Sr g) _AProd"z| Wrapper class for linear operator The call signature of the __call__ method matches the callback of the PROPACK routines. c[U5Ulg![a' [[R"U55Ulgf=fN)rA TypeErrornpasarray)selfrs U/var/www/html/pdf-tiff/venv/lib/python3.13/site-packages/scipy/sparse/linalg/_svdp.py__init___AProd.__init__)s6 5%a(DF 5%bjjm4DF 5s.AAcUS:XaURRU5USS&gURRU5USS&g)Nr)rmatvecrmatvec)rtransamnxys r__call___AProd.__call__/s5 Q;66==#AaD66>>!$AaDc.URR$r)rshapers rr( _AProd.shape5svv||r&cURR$![aS URR[R "URR S55Rs$f=f)Nr)rdtypeAttributeErrorrrzerosr(r)s rr, _AProd.dtype9sU B66<<  B66==$&&,,q/!:;AA A BsAA54A5)rN) __name__ __module__ __qualname____firstlineno____doc__rr$propertyr(r,__static_attributes__r&rrr"s; 5 % BBr&rc^ Uc [S5eUR5nUS;a [S5eU(dUS:Xa [S5e[U5nURRn[ Un[ UnURunnUS :dU[UU5:a [S 5eUcS U-nUcSn[US -US -U5nXA:a[SUSUS35eU(aS OSnU(aS OSn[R"UUS -4SUS9n[R"UU4SUS9n[R"USUR5S9n[R"USUR5S9nUcfURUS9USS2S4'[R"[R"SUS95(a!USS2S4==SURUS9-- ss'O UUSS2S4'U c4[R "[R""U5R$5n U c#[R""U5R$S-n U(a][R&"XX4UR5S9nUcXA- nU[UU- UU5:a [S5eUS:a [S5eO%[R&"XU 4UR5S9n[((a[R*O[R,n [R&"[/[1U 55[/[1U554U S9n!Sn"U(dU(aHUU-SUS---SU--S-n#U#[3SUS--SU--S-U"[3UU5-5- n#S U-n$O5UU-SU--SUS---S-[3UU-SU-S-5-n#SU-S-n$[R"U#UR5S9n%[R"U$U S9n&[R"S UR5S9n'[R"S U S9n(UR55(a![R"UU-U-US9n)U%U)U&4n*OU%U&4n*[7US!5(aMUR9S[R:"[R*5R2S[R<S"9n+OLUR?S[R:"[R*5R2S[R<S"9n+U(a(U"[@UUUUUUUUU UUUUU/U*QUPU!PU'PU(PU+P76n,OU"UUUUXU UUUUU/ U*QUPU!PU'PU(PU+P76n,U,S:a[CS#U,S$35eU,S:a[CS%US&US'35eUSS2SU24UUSS2SU24RE5RFU4$![aI [R"[R"SUS95(aS nOS n[ Un[ UnGNf=f![a [SU35ef=f)(a4 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 rz#`which` must be either 'LM' or 'SM'rz#`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!  integers)lowhighr<r,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.loweruniformsqrtfinfoepsarrayrint64int32intboolmaxisupperhasattrrDiinfouint64randint_which_converterrconjT)-rkwhichirl_modekmax compute_u compute_vv0 full_outputtoldeltaetaanormcgselr min_relgapshiftsmaxiterrngaprodtyp lansvd_irllansvdr r!jobujobvuvsigmabnddoption int_dtypeioptionNBlworkliworkworkiworkdparmiparmzworkworks rng_stateinfos- rrrAsNr {UVV KKME L >??  >?? 1IE ++  C #%c* c" ;;DAq A1s1ay=IJJ |!t q1ua!eT "D x gaS +, , 1qD1qD !TAXc5A !T#S1A HHQc 5E ((1Csyy{ 3C  z++1+%!Q$ ??288AS1 2 2 adGrCKKQK// /G :AadG  } ))* {hhsm4'((E:#))+N >XF s4&=!Q' '<= = aZ34 4((E.ciikB%I288IhhDIDI7yIG BIA$' !AdF*Q. QtQwY4'!+RAq \::4A$47*Q.QUAdFQJ1GG4! 88E -D HHV9 -E HHQciik *E HHQi (E {{}}QS1eU"e sJLL*..Qbii! KK*..Qbii  *5141fa#E3<49<;B<"<$)<+0<2;< dD!QeQsAII%I'.I05I7<I>GI ax1${ CE E 6 &' ' QU8UAa!eHMMO--s 22O # ??288AS1 2 2CC%c* c"#Z :4QC89 9 :s*V= X=AXXX,)r TNTTNFrNNrFTgMb`?NNN)__all__numpyrscipy.linalg.lapackrscipy.sparse.linalgr scipy.linalgrrslansvddlansvdclansvdzlansvdrK slansvd_irl dlansvd_irl clansvd_irl zlansvd_irlrJr`rrr7r&rrs ))0$                             BB>15JK8<;?o3r&