+ LV-jp(*R.t^RIt^RIHt^RIHt^RIHtR]PR]PR]PR]P/t R]PR]PR]P R]P"/tR ^R ^/t!R R 4tRR ltR#)_svdpN)aslinearoperator) LinAlgError)_propackfdFDLMSMcPa]tRt^!toRtRtRt]R4t]R4t Rt Vt R#)_AProdz| Wrapper class for linear operator The call signature of the __call__ method matches the callback of the PROPACK routines. c\V4VnR# \d)\\P!T44TnR#i;iN)rA TypeErrornpasarray)selfrs&&j/Volumes/fast/ai/experiments/ui-tars-smoke/.venv/lib/python3.14/site-packages/scipy/sparse/linalg/_svdp.py__init___AProd.__init__(s6 5%a(DF 5%bjjm4DF 5s/AAcV^8Xd!VPPV4VR&R#VPPV4VR&R#)NNNN)rmatvecrmatvec)rtransamnxys&&&&&&r__call___AProd.__call__.s5 Q;66==#AaD66>>!$AaDc.VPP#r)rshapers&rr& _AProd.shape4svv||r$cVPP# \dXTPP\P !TPP ^,44Pu#i;i))rdtypeAttributeErrorrrzerosr&r's&rr+ _AProd.dtype8sU B66<<  B66==$&&,,q/!:;AA A BsAA:9A:)rN) __name__ __module__ __qualname____firstlineno____doc__rr"propertyr&r+__static_attributes____classdictcell__) __classdict__s@rr r !s> 5 % BBr$r c~ Vf \R4hVP4pVR9d \R4hV'gVR8Xd \R4h\V4pVPPp\ V,p\ V,pVPwppV^8gV\VV48d \R 4hVf ^ V,pVfR p\V^,V^,V4pWA8d\R V R V R 24hV'd^M^pV'd^M^p\P!VV^,3RVR7p\P!VV3RVR7p\P!VRVP4R7p\P!VRVP4R7pVfuVPVR7VR&\P!\P!^VR74'd-VR;;,RVPVR7,, uu&MVVR&V f5\P !\P"!V4P$4p V f(\P"!V4P$R,p V'dl\P&!WW3VP4R7pVf WA, pV\VV, VV48d \R4hV^8d \R4hM(\P&!WV 3VP4R7p\P&!\)\+V 44\)\+V443RR7p ^ p!V'g V'dVV,^V^,,,^ V,,^,p"V"\-^V^,,^V,,^,V!\-VV4,4, p"^V,p#MiVV,^ V,,^V^,,,^,\-VV,^V,^,4,p"^V,^,p#\P!V"VP4R7p$\P!V#\P.R7p%\P!^VP4R7p&\P!^\P.R7p'VP14'd-\P!VV,V,VR7p(V$V(V%3p)MV$V%3p)VP3^\P4!\P64P,^\P8R7p*V'd.V!\:V,VVVVVVVV VVVVV.V)OVNV NV&NV'NV*N5!p+MV!VVVVWV VVVVV. V)OVNV NV&NV'NV*N5!p+V+^8d\=RV+ R24hV+^8d\=RV RV R24hVRRV13,VVRRV13,P?4P@V3# \dV\P!\P!^TR74'dRpMRp\ T,p\ T,pELi;i \d\RT 24hi;i) 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 instancer z#`which` must be either 'LM' or 'SM'z#`which`='SM' requires irl_mode=True)r+r rz.k must be positive and not greater than m or niz3kmax must be greater than or equal to k, but kmax (z) < k ()r)orderr+)sizery?zv0 must be of length g?z0shifts must satisfy k <= min(kmax-shifts, m, n)!zshifts must be >= 0!i)lowhighr;r+z#An invariant subspace of dimension z was found.zk=z0 singular triplets did not converge within kmax=z iterations>r r )rr)! ValueErrorupperr r+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@sr {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+%$ ??288AS1 2 2 dGrCKKQK// /G :AdG  } ))* {hhsm4'((E:#))+N >XF s4&=!Q' '<= = aZ34 4((E.ciikBhhDIDI7sCG BIA$' !AdF*Q. QtQwY4'!+RAq \::4A$47*Q.QUAdFQJ1GG4! 88E -D HHV288 ,E HHQciik *E HHQbhh 'E {{}}QS1eU"e  "(();)?)?"#2996I*5141fa#E3<49<;B<"<$)<+0<2;< dD!QeQsAII%I'.I05I7<I>GI ax1${ CE E 6 &' ' QU8UAa!eHMMO--s 22w # ??288AS1 2 2CC%c* c"#Z :4QC89 9 :s$.X?Z"?rs )0$                    !! BB>c3r$