ۓ,j3 LSr/SQrSSKrSSKJr SSKJrJ r SSK J r SSK J r SSK Jr SS K Jr SS KJr SS KJr SS KJrJr \ r SSKr\"S R;\R<"\R>R@5S55 \RCS5 S%Sjr"S&Sjr#S'Sjr$"SS\5r%"SS\5r&"SS\5r'"SS\5r("SS\5r)"SS\5r*"SS\5r+"S S!\5r,\)r-\$r.\(r/\%r0S"r1S#r2\3S$:Xa\2"5 gg! N=f)(zR Author(s): Matthew Loper See LICENCE.txt for licensing and contact information. )invsvddetslogdetpinvlstsqnormN)Ch depends_on) NanDivide)asarray)sqrt)sum) concatenate)randn)rowcol tensorinvcRUcUb [S5e[[US-55$)Nz('ord' and 'axis' should be None for now.)NotImplementedErrorch_sqrtch_sum)xordaxiss ڛC:\Users\johnn\AppData\Local\Packages\PythonSoftwareFoundation.Python.3.13_qbz5n2kfra8p0\LocalCache\local-packages\Python313\site-packages\chumpy/linalg.pyrr s+ $*!"LMM 6!Q$<  cUS:wa [S5e[S5nXlXl[ URR U5UR- S-SS9n[ n[ nX4XV4$)Nz%non-default rcond not yet implementedc6[U5RU5$N)rdot)abs rlstsq..3saQr rr )r) Exceptionr r&r'rr%r)r&r'rcondr residualsrankss rrr/se {?@@ '(A C C!##''!*qss*a/q9I DA   r c(US:wa [S5eUS:wa [S5eURSURS:nU(a URn[US9n[ XS9n[ XUS9nU(aXTUR4$XdUR4$)Nr zfull_matrices must be 0r zcompute_uv must be 1r)rsvd_d)rr1svd_v)r*shapeTSvdDSvdVSvdU)r full_matrices compute_uvneed_transposer1r2svd_us rSvdr<<s122Q.//WWQZ!''!*,N CC 1IE 1 "E 1 /EUWW$$UWW$$r c*\rSrSrSrSrSrSrSrg)PinvRmtxc^URnURSURS:a>URR[ URUR555nO=[ URRU55RUR5nX0lg)Nr r )r@r3r4r%Inv_result)selfwhichr@results r on_changedPinv.on_changedUsphh 99Q<#))A, &UUYYs377355>23F3(,,SUU3F r c.URR$r$)rCrrDs r compute_rPinv.compute_r]s||~~r cjXRLa%URRUR5$gr$)r@rCdr_wrtrDwrts rcompute_dr_wrtPinv.compute_dr_wrt`s) ((?<<&&txx0 0 r )rCN __name__ __module__ __qualname____firstlineno__dtermsrGrLrR__static_attributes__r rr>r>Rs F1r r>c*\rSrSrSrSrSrSrSrg) LogAbsDetwrc[RRURR5uUlUlgr$)nplinalgrrrJsign)rDrEs rrGLogAbsDet.on_changedzs'"$))"3"3DFFHH"= 42r r]c$\rSrSrSrSrSrSrg) SignLogAbsDet) logabsdetcZURRnURR$r$)rkrJrb)rD_s rrLSignLogAbsDet.compute_rs! NN  ~~"""r cgr$r[rPs rrRSignLogAbsDet.compute_dr_wrtsr r[NrUrVrWrXrYrLrRrZr[r rriris F#r ric$\rSrSrSrSrSrSrg)Detrch[RRURR5$r$)r`rarrrJrKs rrL Det.compute_ryy}}TVVXX&&r cXRLaS[UR[RR URR5R -5$gr$)rrrJr`rarr4rPs rrRDet.compute_dr_wrts? &&=tvv dffhh 7 9 99: : r r[Nrqr[r rrsrss F';r rsc$\rSrSrSrSrSrSrg)rBr&ch[RRURR5$r$)r`rarr&rJrKs rrL Inv.compute_rrwr c XRLagURnURS::a![R"X"R 5*$[R "USURSURS45n[R"USS5n[R "USURSURS45n[R"[URS5Vs/sH)n[R"X$X45R *PM+ sn5R n[R"U5nU$s snf)Nrr"r ) r&rJndimr`kronr4reshaper3rollaxisdstackrangesp block_diag)rDrQAinvAinvTirFs rrRInv.compute_dr_wrts ff vv 99>GGD&&)) )::db$**R.$**R.%IJDKKb"-EJJur5;;r?EKKO&LMEYYuTZZXY]G[\G[!%(!;!=!= =G[\]__F]]6*F  ]s70Er[Nrqr[r rrBrBs F'r rBc@\rSrSrSr\"S5S5rSrSrSr g)r5rc[RRURRSS9nUSUSUSR /nSUS[R "US5[R"S5:'U$)NF)r8r r rg)r`rarrrJr4absspacing)rDrFs rUDVSvdD.UDVsitvvxxu=)VAYq 479q "&&#bjjm34 r c URS$)Nr )rrKs rrLSvdD.compute_rsxx{r cBXRLagURup#nURRRnUSUS2SUS24nUSUS2SUR24n[ R "SX$5nURURSS45nU$)Nr r z ik,jk->kijr")rrrJr3sizer`einsumr)rDrQudvshprFs rrRSvdD.compute_dr_wrts ff  ((affhhnn gs1vgwAw  gs1vgww <.a" 56 r r[N) rUrVrWrXrYr rrLrRrZr[r rr5r5s' F_  r r5c(\rSrSrSrSrSrSrSrg)r6r1rc4URRS$)Nrr1rrKs rrLSvdV.compute_rzz~~a  r c XRLagURRup#nURRRRn[ USUS5n[ R"U5nX7SUR&[ R"USUSUSUS45nUSn USn X:de[US5Hn [US5Hn [U 5Hn [U S-U 5Hn[ R"X~X}/X}X~//5n[ R"X+U 4XLU4-X+U4*XLU 4-/5n[ RRUU5nUSXXU4'US*XXU 4'M M M M URSURS:Xde[ R"SU*U5RURRURR45$)Nr r rz ak,ijkl->alij)rr1rrJr3maxr`zerosrrarrayrasolverr)rDrQU_DVrmxszDomegaMNrjklr@rhsrFs rrRSvdV.compute_dr_wrts ff  Qjjllnn""3q63q6" HHTN(277 #a&#a&#a&#a&9: F Fvvs1vA3q6]qA"1Q3] hhT!$KT!$K()*!hhA#q1v A#wq1v~'FG!#c!:*0)Ag*0)Ag+"# wwqzU[[^+,+yy1"e4<>#  $DLL)M "))  CEE-/ 0 3bii6G6G6Nq6Q QR '#,q/  % ) )$ /0 "))  CEE "1 %& '#,q/r __main__)NN)r")r r )4__doc____all__numpyr` scipy.sparsesparserchr r ch_opsr r ch_asarrayrrrr reorderingrr ch_randomrrutilsrrinspectexecjoingetsourcelinesrarappendrrr<r>r]rirsrBr5r6r7rrrrrrrUr[r rrs% D)#!5/ G'' (;(;