+ nDj!-^RIHt^RIt^RIHtHtHtHtHtH t ^RI t ^RI H t ^RI t^RI HtHtHtHtHtHtHtHt^RIHtHtHtHtHtHtHtHtH t H!t!H"t"^RI#H$t$^RI%H&t&^R I'H(t(H)t)^R I*H+t+^R I%H,t,^R I-H.t.^R I/H0t0R;t1t2]0R,R,R,t1]0R,R,R,t2]Pf]Ph.t5]Pl]Pn.t8]5]8,t9Rt:Rt;Rt<!RR4t=!RR4t>] P~PR]94] P~PR.RO4R44tA!RR4tB!RR4tC!R R!4tD!R"R#4tE] P~PR]94] P~PR$.RO4] P~PR%R&R'.4] P~PR(R)R&.4] P~PR*^^ .4R+44444tFR,tGR-tHR.tIR/tJ!R0R14tK!R2R34tLR4tMR5tNR6tOR7tPR8tQR9tRR:tSR;tTR<tUR=tVR>tWR?tXR@tYRAtZ!RBRC4t[RDt\REt]RFt^RGt_] P~PRHRIRJ7RK4taRLtb] P~PRMRR.4] P~PR]54] P~PRN]c!]d!^444] P~PRO]c!]d!^444] P~PRP]c!]d!^444] P~PRQ^^.4] P~PRR^^.4RRSl4444444te] P~PR]54RT4tf] P~PRURN^ /RO^ /RP^ /RQ^ /RV^ /RR^ /34RW4tg] P~PRX]P!.RO.RO.RO.RO.RO.RO.4]P!.RO4]P!.RO.RO.RO.RO.RO.RO.4]P!.RO.RO.RO.RO.43.4RY4ti] P~PR]94RZ4tj] P~PR[]P!.RO4]P!.RO4]P!.RO4]P!.RO4]P!.RO4]P!.RO4]P!.RO4]P!R\R].RR^.R_R.R`R.R\Ra..4]P!R^.^R.^R.RR.R^..43 ]P!.RO4]P!.RO4]P!.RO4]P!.RO4]P!.RO4]P!.RO4]P!.RO4]P!RR.RR.RR.RR.RR..4]P!RR.RR.RR.RR.RR..43 .4Rb4tk] P~PR]94] P~PR$.RO4] P~PR*^^ .4Rc444tl] P~PR]94] P~PR.RO4Rd44tm] P~PRe]c!]n!]5]5,]9444Rf4to] P~PRe]c!]n!]5]5,]9444Rg4tp] P~PRe]c!]n!]5]5,]9444Rh4tq] P~PRi]P!.RO4]P!.RO4]P!.RO4]P!.RO4]P!^^ .^ ^.^^ .^^A.^^..4]P!Rj^.^R.^R.R^.^R..43]P!.RO4]P!.RO4]P!.RO4]P!.RO4]P!ERER.ERER.ERER.ERER..4]P!RER.RR.RR.RR..43.4Rk4trRlts] P~PRm]c!]n!]9]5]5,444] P~PRn]c!]d!^444Ro44tt] P~PRm]c!]n!]9]5]5,444] P~PRn]c!]d!^444Rp44tu] P~PRm]c!]n!]9]5]5,444] P~PRn]c!]d!^444Rq44tv] P~PRm]c!]n!]9]5]5,444] P~PRn]c!]d!^444Rr44tw] P~PRs^]P!.ER O4]P!.ER O4]P!.ER O4]P!.ER O.ER O.ERO.ERO.43.4Rt4tx] P~PR]94] P~PRu.ERO4Rv44tyRwtz] P~PRx.ERO4] P~PRyRzR{.4R|44t{] P~PRxR}R~.4] P~PRyRzR{.4R44t|] P~PR]94] P~PR.ERO4R44t}] P~PR]94R4t~] P~PR]94] P~PRRRH.4] P~PRRR).4R444t] P~PR]94] P~PR^^.4] P~PRRR).4R444t] P~PR]9^,4] P~PRRRH.4] P~PRRR).4R444t] P~PR]P!.RO4]P!.RO4]P!.RO4]P!R\R].RR^.R_R.R`R.R\Ra..4]P!R^.^R.^R.RR.R^..43]P!.RO4]P!.RO4]P!.RO4]P!RR.RR.RR.RR.RR..4]P!RR.RR.RR.RR.RR..43.4R4t] P~PRm]c!]n!]9]5]5,444] P~PRRR3R)R3.4R44t] P~PRm]c!]n!]9]5]5,444] P~PRRR3R)R3.4R44t] P~PRm]c!]n!]9]5]5,444] P~PRRR3R)R3.4R44t] P~PR]P!.RO4]P!.RO4]P!^^ .^ ^.^^ .^^A.^^..4]P!Rj^.^R.^R.R^.^R..43]P!.RO4]P!.RO4]P!ERER.ERER.ERER.ERER..4]P!RER.RR.RR.RR..43.4R4t] P~PRRRH.4] P~PR]94R44t] P~PR]94R4t] P~PR]P!.ERO.ERO.ERO.ERO.4]P!.ERO.ERO.ERO.ERO.4^^3]P!.ERO.ERO.ERO.ERO.4]P!.ERO.ERO.ER O.ER!O.4^^3.4R4t] P~PR]94R4t] P~PR]94R4t] P~PR]P!.ER"O.ER#O.ER$O.ER%O.4]P!.ER&O.ER'O.ER(O.ER)O.4]P!.ER*O.ER+O.ER,O.ER-O.4]P!.ER.O4RR3]P!.ER/O.ER0O.ER1O.ER2O.4]P!.ER3O.ER4O.ER5O.ER6O.4]P!.ER7O.ER8O.ER9O.ER:O.4]P!.ER.O4RR3.4R4t] P~PR]94R4t] P~PR]P!.ER;O.ERO.4]P!.ER?O.ERnpflapackHAS_LP64division HAS_ILP64print_functionabsolute_importfind_best_lapack_type) r__doc__pytestskipsetsplitlistdir startswithappend)names ignore_listmissingr's r1test_lapack_documentedrJ/s ~~ ,- $$& 'E KfGF $$)@! NN4  b=NNN=r3cV^RIHp^RIHpW8XgQh^RIHp^RIHpW#8XgQhR#)r)r:)r8N)scipy.linalg.blasr:scipy.linalg.lapackr8)blas_has_ilp64lapack_has_ilp64 blas_has_lp64lapack_has_lp64s r1(test_lp64_ilp64_blas_lapack_both_or_nonerRFs&=A  -- -;?  ++ +r3c8a]tRt^PtoRtRtRtRtRtVt R#)TestFlapackSimplec .RO.RO.RO.p.RO.R O.R O.R O.p\Fp\RVR7pVfKV!V4wrVrxp \V '*\V 44\ WQ4\ Wg3^\ V^,4^, 34\ V\P!\ V444V!V^^R7wrVrxp \V '*\V 44K R#) gebalr/N)permutescale)rV)) )rVrrga2U0*3?)r]rrgMb`?)r`rVrr)rrVrr) DTYPESr#rreprrrlenr6r ) selfaa1r/fbalohipivscaleinfos & r1 test_gebalTestFlapackSimple.test_gebalRs  9 -E 6Ay)*1 &BB$ Hd4j ) %b , "As1Q4y{#3 4 %hA @)*2q)B &BB$ Hd4j )r3c.RO.RO.RO.p\P3F;p\RVR7pVfKV!V4wrEp\V'*\ V44K= R#)gehrdrXNikiifii"iii)r6float64r#rrd)rfrgr/rihttaurns& r1 test_gehrdTestFlapackSimple.test_gehrdgsU    jj\E 6AyaDMBT Hd4j ) "r3c \P!^^.^^..4p\P!^^.^^..4p\P!^ ^ .^ ^ ..4pRpREFpVPV4VPV4VPV4rp\R V34wp VP 4'dV^;;,R, uu&RpV !WgV4wrp \ \P !Wj4\P !W4,W,4V !WgWVR7wrp \ \P !VP4PV 4\P !WP4P4,W,^R7V !WgVR R7wrp \ \P !Wj4\P !W4, W,^R7EK R#) rVTfdFDr*C)tranatranbdecimal)isgnN)trsyl) r6arrayr-r#isupperrdot conjugater~) rfrgbctransr/rhb1c1rxrZrns & r1 test_trsylTestFlapackSimple.test_trsylrs~ HHq!fq!f% & HHq!fq!f% & HHq"gBx( )E%!((5/188E?BB%j2%8FE}}1 "22.NAd %bffRmbffQm&C&+j 2#22%HNAd %FF2<<>++Q/"&&LLN"?@C'G<4' ffRVVBZ0!U* ffRVVBFF2JQ%?@!T) ffRVVBFF2JQ%?@ ,1*r3N) __name__ __module__ __qualname____firstlineno__ror{rr__static_attributes____classdictcell__ __classdict__s@r1rTrTPs** *=8--r3rTc8a]tRt^toRtRtRtRtRtVt R#)TestLeastSquaresSolversc \PPR4p\\4Fwr#^ p^p^pVPWE34P V4pVPV4P V4p\ RVR7wr\WWV4p V !WxV R7wr\V ^84V !WxRV,V R7wr\V ^84K \EFp\P!RR.RR.R R ..VR7p\P!.ROVR7p\ RWx34wrpVPwrE\VP4^8XdVP^,pM^p\WWV4p V!WxV R7wppp \VR R\P!RR .VR7^\P!V4P,R 7V!V4wpp \!VV4EK \"EFp\P!RR.RR.RR..VR7p\P!.ROVR7p\ RWx34wrpVPwrE\VP4^8XdVP^,pM^p\WWV4p V!WxV R7wppp \VR R\P!RR.VR7^\P!V4P,R 7V!V4wpp \!VV4EK R #)rXlworkTTCCrr?@@@@ @N-@rtol)gels gels_lwork0@g1@g4@)rrgeqrfr祪,?@@?@@ @ffffff?ry1@@y4@R ?\j,?W?)r6random default_rng enumeratercr-r#rr REAL_DTYPESrr.rerfinfoepsr r+)rfr0indr/mnnrhsrhrglsglslwrr5rnrrrlqrr lqr_truths& r1 test_gels!TestLeastSquaresSolvers.test_gelssii##D)#F+JCAADQF#**51BA%%e,B)*@NJC#5Q5ER51JAq DAI R6#;eDJAq DAI ,![EC:::'.35B,E:B&63bX'? #De88DA288}!xx{#:!:Ee4LCD AcrFBHH.A.@.B38%:"$BHHUO$7$7!7 9"'r Iq!Q sI ./!2$^EHc?$h/$h/18=?B6eDB&63bX'? #De88DA288}!xx{#:!:Ee4LCD AcrFHH&J&J&L+028:288E?;N;N8N P"'r Iq!Q sI ./$r3c \EFcp\P!RR.RR.RR..VR7p\P!.ROVR7p\RW#34wrEVPwrg\ VP4^8XdVP^,pM^pV!WgVR4wrp \ \P!V 44p T p V!W#WRRR4wrpp \VRR\P!RR .VR7^\P!V4P,R 7\V\P!R R .VR7^\P!V4P,R 7EKf \EFqp\P!RR.RR.RR..VR7p\P!.ROVR7p\RW#34wrEVPwrg\ VP4^8XdVP^,pM^pV!WgVR4wp pr\ \P!V 44p \ V4pT p V!W#V VV RRR4wrpp \VRR\P!RR.VR7^\P!V4P,R 7\V\P!R R.VR7^\P!V4P,R 7EKt R#)rrrrrrrXFNrrYN))1)@*@.?U.*@_Y@r)gelsd gelsd_lworkrrrrrrrrrr rr6rr#r.reintrealrrrr+)rfr/rhrrrrrrworkiworkrnr iwork_sizersrankrwork rwork_sizes& r1 test_gelsd"TestLeastSquaresSolvers.test_gelsds [EC:::'.35B,E:B!12J35("< E88DA288}!xx{!,A$ ; D &EJ$RU%' 7 A$ AcrFBHH.A.@.B38%:"$BHHUO$7$7!7 9 Arxx);):)"?"ABQ ]+- .%Eu#%VnnfTrTlVAYU 9K-K,LMNnneDbDk401 6)4q&!A$C LLQ Q1OPQsS  ""%RVVBHHUO,,.CDEBHHRZZ,@,D,D(D "**!5!9!99 ;r3rN)rrrrrFrrrs@r1r'r's;;r3r'c <a]tRtRto]P P R]4R4t]P P R]UUu.uF#pRFpVR8XdV\9dKW3NK K% upp4]P P RRR .4]P P R RR.4R 444t ]P P R .RO.RO.RO.4R 4t Rt ]P P RRR.4R4t RtVtR#uuppi) TestTbtrsir/c XV\9dh\P!.R O.RO.VR7p\P!RR .RR.RR.RR..VR7p\P!^^.RR.^^.^R..VR7pMV\9dk\P!.RO.RO.RO.VR7p\P!RR.RR.RR .R!R"..VR7p\P!RR#.R$R%.R&R'.R(R)..VR7pM\ RV R24h\ RVR7pV!W#RR7wrg\ V^4\Wd^R R 7R #)*zTest real (f07vef) and complex (f07vsf) examples from NAG Examples available from: * https://www.nag.com/numeric/fl/nagdoc_latest/html/f07/f07vef.html * https://www.nag.com/numeric/fl/nagdoc_latest/html/f07/f07vsf.html rXg333333*@y@z Datatype z not understood.tbtrsLabruploh㈵>rr0Np= ף)rRgQ@gHzG@g{Gz?)ggq= ףp@gHzGrgp= ף0g(\+gףp= 0g(\gHzG,gQ#rr)y ףp= Q@y{Gz@GzyQ?HzGy)\(??)yQ Q @yq= ףpGz@yףp= ?{Gzr)yQ?q= ףp @y)\(zGrryQ! ףp= yףp= 8Gzyp= #/)\h7y\(LHzG @yQHz6@yףp= 3@(\=y{Gz-333333yQ+3@GzT5@y?@y?yyt&m=#yi6@Ug$B@y[a^C?b->y-@ji&*!)rr6rr+ ValueErrorr#rr)rfr/rNrx_outrKrrns&& r1test_nag_example_f07vef_f07vsf(TestTbtrs.test_nag_example_f07vef_f07vsfsy K 424 %'B65/!6* %.!5/+ % &A HHq!f!2h !f "g'$) *E n $KF:<!&'B=.9*N;)>:)>:< % &A HHr6l%w/46JK35JKM$) *E y/?@A A '%82-T1qt4r3z dtype,transNrrOUrLr6c \PPR4p^^^rp\RVR7p VR8Hp W,p W, p \ W)^, R4p V Uu.uFq\ V4, NK ppVUu.uFp\ V3W4NK ppVR8Xd\P!WaR7VV &\VV RR7p\P!V^,V3V4p\V 4F9wppVPV4VV\V^4\VV,V413&K; \ Wg3W4pV !VVW2VR7wpp\V^4VR8Xd\VV,VRR 7R #VR 8Xd"\VP V,VRR 7R #VR 8Xd0\VP P#4V,VRR 7R #\%R 4huupiuupi)irKrXrYdia)offsetsformat)rNrrOrr6rXg-C6 ?rr~rzInvalid trans argumentNr)r6r RandomStater#r9rr2r r"r rdiagonalrminrrr~rrT)rfr/rrOr6r0rrkdrKis_upperkukl band_offsetsr band_widthswidthbandsrgrNrowkrrns&&&&& r1test_random_matricesTestTbtrs.test_random_matrices'sii##D) A 6CK ] WRq"- +782F4T1 C< AE14 0 c\ ACC!GQT 2 c\ ACCMMOa/ >56 699+s G?G"zuplo,trans,diagc  \R\PR7p\PP R4pVPR4pVPR4p\ \ WFWqW#4R#)z?Test if invalid values of uplo, trans and diag raise exceptionsrKrXrN)r]r[)r[r])r#r6rxrrrr)rfrOrr6rKr0rNrs&&&& r1&test_invalid_argument_raises_exception0TestTbtrs.test_invalid_argument_raises_exception[sR! ;ii##D) ZZ  JJv iAUAr3c \P!R\R7p\P!^\R7p\R\R7p^VR&V!WRR7wrE\ V^4R#)a(Test if a matrix with a zero diagonal element is singular If the i-th diagonal of A is zero, ?tbtrs should return `i` in `info` indicating the provided matrix is singular. Note that ?tbtrs requires the matrix A to be stored in banded form. In this form the diagonal corresponds to the last row.rXrKrYrMNr\r])rr\)r6r floatr#r)rfrNrrKr5rns& r1test_zero_element_in_diagonal'TestTbtrs.test_zero_element_in_diagonalisQWWV5 ) GGAU # 65 2-T1r3zldab,n,ldb,nrhsc \P!W3\R7p\P!W43\R7p\R\R7p\ \ WuV4R#)z2Test ?tbtrs fails correctly if shapes are invalid.rXrKNr6r rrr#rr)rfldabrldbrrNrrKs&&&&& r1test_invalid_matrix_shapes$TestTbtrs.test_invalid_matrix_shapesys@ WWdYe , GGSKu - 6iA.r3rN)rXr~r)rYrXInvalid)rYr{rX)r{rXrX)r^r^rr^)r^r^r\r^)rrrrr?r$r%rcrVrrkrnrsryrr)r/rrs00@r1rIrIs; [[Wf-+5.+5Z [[]+1N+1%e&+slu 7K-un?N-+1NO [[Vc3Z0 [[Vc3Z0,711 O ,7\ [[.33356B 6B  [[.**1/ /kNs D % D rIrrOrYrLr6rXrc \V V V V V 2P44p\PP V4pVPWD3R7VPWD3R7R,,pVP \P !^VP^^ 4V44pWx, p\P!V\P4'd VPMTpVR8Xd\P!V4M\P!V4pVR8Xd5V\P!V4R\P3,,pVPV4p\!RV34p V !WqW#R7wrVR8Xd\P"P%V\P&R7p \P"P%\P"P)V4\P&R7p ^W,, pMI\P"P%V^R7p\!R V34wppV!V4wpppV!VWR7wr^p\+WVR 7R #) rr*rYr)trcon)rrOr6IordrrNgecongetrf)rCencoder6rrpermutedlogspaceintegers issubdtypefloatingrrrr6r7r-r#linalgrinfrr)r/rrOr6rseedr0Aoffsetr}rEr5norm_A norm_inv_Aranormrrluipvtrnrs&&&&& r1 test_trconrs 5'$vdVA3/668 9D ))   %C  #**1&*"9""<B DA!!,|011a|/=CG IA!QQ7%a|44 AA!QQ!<%a|44 AA!QQ!<%1a|44 AA!QQQQ!L%a|44 AA!RbA>%1a|44 A11!< #==4)g   MM!E( #!"==?AeHbhhqwwqz6JK%%"((1771:*>SI%$5$559r3c\R\PRR7pRp\WVRRR7pVR8Xg VR8XgQhR#R#) gesdd_lwork preferredr/ilp64iA%T)r+r*i`DiDN)r#r6float32r) sgesdd_lworkrrs r1 test_sgesdd_lwork_bug_workaroundrsO $M*57L A T1a(;(;&;#F2771:&3S!!&a 5T1 MM! ) IIaggaj !!$ YYqwwqz!| $"Q$( "d(  FF1u %qsA(A!QqSMAbqEAaDq5)CFRXXa^,CCAq! A //!R(SS\' vvacc266!<( aa(;(;&;#Fr3rNrVr\) rrrrr?r$r%rrrrrrs@r1rr5sh [[Wk2,3,  [[Wk2 [[S&)7G*37Gr3rc a]tRtRto]P P R]4R4t]P P R] !] !] ]444]P P RR 4R44t Rt VtR#) TestHetrdiy complex_dtypecp\P!RVR7p\RV34p\\W24R#)rrXhetrdNrr)rfrrrs&& r1test_hetrd_with_zero_dim_array(TestHetrd.test_hetrd_with_zero_dim_arrayzs, HHV= 1 1$/j%+r3zreal_dtype,complex_dtyperc \P!W3VR7p\R V34wrV\P!^W^,,^,^,VR7R\P!^W^,,^,^,VR7,,V\P!V4&\P !V\P !\P!V444R FpV!WR7wr\V ^4K \Wa4p V!V^V R7wrrp \V ^4\W^\P!V4P,RR7\V \P !\P!V444\V R4\VR4V!WJR7wrrp \V ^4\P!WBR7p\P!VP^,\R7pWVV3&\P!VP^,^, \R7pWV^,V3&WVV^,3&\P !WVR7p\#V^, 4Fp\P!WR7pV RV1V^,3,VRV%RVV&\P !WVR7VV,\P$!V\P&!V44,, p\P(!VV4pK \P*!VR 4p\P&!VP,V,4VV&\P(!\P&!VP,4\P(!VV44p\VV^ \P!V4P,RR7R#) rXr*rrrr/rrN)r hetrd_lworkrrVr)r6r r#rr fill_diagonalrr6rrrrrrr.rr r9rrrrr~)rfr real_dtyperrrrrr5rnrrrrrzr~rjrrrDrrrQHAQs&&&& r1 test_hetrdTestHetrd.test_hetrds HHaV= 1 5t <  IIaaC!AZ 8299QQ3 1 J?? @ "  q !" BGGBGGAJ/0A!!-GA q ! {.!&aq >T1a(<(@(@&@sK2772771:./3S!!&a 5T1 MM! . IIaggaj ,!Q$ YYqwwqz!|3 /"Q$( "bd(  FF1} -qsA0A!QqSMAbqEAaDq=1a&288Arwwqz223Aq! A //!R(WWQSS\*' vvbggacclBFF1aL1  !"RXXj1555C r3rNr)rrrrr?r$r%r+rrCziprrrrrs@r1rrysz [[_n=,>,  [[7!#k>"BCE [[S&)A*EAr3rc \\4EF_wr\RVR7wr#\V^^^R7pV^8d`\P !.RO.R O.R O.R O.R O.R O.VR7p\P !.ROVR7p\P !RR.VR7pMb\P !.RO.RO.RO.RO.RO.RO.4p\P !R.R.R.R.R.R..4p\P !^VR7p\P !.RO.RO.VR7pV!WXWgVR7wrp V^8d\P !.RO4p M\P !.RO4p \W^R7EKb R#)gglserX)rrprrrN)r gglse_lwork)g= ףp=g{Gzg(\ؿ?)zGgHzG?gףp= ӿQ)ffffff@gQ?g?gffffffֿ)rg{Gz?Qg{Gz?)333333?g333333?rg ףp= )g{Gz{Gz?gzG?)grgGz?gHzGgzGg= ףp=?)yQ?QyQQ?yQ{Gz@y= ףp=?)y\(\￮Gz?y333333RQ?yQzG?yQQ?)yףp= ?q= ףpݿy)\(?{Gz?y)\(?(\ſy(\333333?)yGz?RQ?yRQ?HzGy\(\ ףp= ׿y)\(?ɿ)y(\?RQ?y?{Gz?y(\ſq= ףpݿyQ?q= ףp?)yHzG?Qѿy?QyQ뱿Gz?yp= ף?p= ף?yRQ ףp= ?yffffff?GzyzGGzyQ?ffffff @yp= ף)\(@y(\@Q?)rrr)rrrr)^"L?\}?rr)y!f?$_Kdy^gŵ翸F @y!f?}dy61ŵe_ @)rrcr#rr6rr r) rr/func func_lworkrrgrrrr5resultrs r1 test_gglsersj' +,D279zQ!q9 7553434 6 >C DA AOA"b/AMNKKKK MNA ;-%,&-%,&-%, ()A %(A HH'):;5 I!!E:1a 7xx!-.H xx!:;H "&A>>Y(r3c ^\PPR4p\\\ ,4EFnwr^ pV^8d?\ RVR7p\ R VR7wrVVPW334PV4pM[\ RVR7p\ R VR7wrVVPW334VPW334R,,PV4pWwP4P,^, ^\P!W2R7,,p\V^4p\WC4p V!Wy^R7wrp V!WV^R7wr\\^V , \PP!V^R7, 4V ,^84EKq R#) r sytrf_lworkrX hetrf_lworkr*rr)rgipivrrrN)syconsytrf)heconhetrf)r6rrrrcr+r#r-rr~r rrrrrcond)r0rr/rrfunconfunctrfrrrldur r5rconds r1test_sycon_heconrsL ))   %C~ 56   7)-uEJ./AOOFQF#++E2A*-uEJ./AOOFQF#cjj!&8&;;CCEJA^Q 266!#9!9 9Q z-qQ7 1C%qAAeGbiinnQ!n445e;a?@-7r3c\PPR4p\\4EFwr^ p\ RVR7wrErgVPW334P V4pWP,^, pVPW334P V4p WP,^, ^\P!W2R7,,p V!W4wrp \V ^8H4V!V 4wr\V ^8H4V!W4wr\V ^8H4V!V4wrp \V ^8H4\WRR7EK! R#)rrXgiUMu?rN)rsygstsyevdsygvd) r6rrrrr#r-r~r rr)r0rr/rrrrrrBeig_gvdr5rnrrgeigs r1 test_sygstrs  ))   %C ,  %57I(@(@/2Qw!|! D D?(r3c \PPR4p^ ^^r2p\\4EFwrE\ RVR7wrg\ WqV4pV^8d\\VPW34PV44p VPW334PV4p \ RVR7wrM\VPW34VPW34R,,PV44p VPW334VPW334R,,PV4p \ RVR7wr\ WV4p V!WR7wrp\P!\P!WR7VRVR13,34p\P!W%R7p\V4Uu.uFUpVVV,VV.R3,PPVV.R3,P44,, NKW pp\\PV4pV^8dRMRp^ \P !V!R 4P"4,pV !WWR7wpp\%V^8H4\'VVPV 4, \)V 4VR R 7V !WV VV R 7wpp\%V^8H4\'VVP4PPV 4, \)V 4VR R 7V !WV R V R7wpp\%V^8H4\'VV PV4, \)V 4VR R 7V !WV R VV R7wpp\%V^8H4\'VV PVP4P4, \)V 4VR R 7EK R#uupi)a  This test performs a matrix multiplication with an arbitrary m x n matrix C and a unitary matrix Q without explicitly forming the array. The array data is encoded in the rectangular part of A which is obtained from ?TZRZF. Q size is inferred by m, n, side keywords. rrXr*rr)Nr~rrrr/rr)rr)rrrr%)ormrz ormrz_lwork)unmrz unmrz_lwork)r6rr^rrcr#rrr-r4r r9r~rrrr(rrrr)r0qmqncnrr/r&r)lwork_rzrrorun_mrz orun_mrz_lw lwork_mrzr*rzrnr+r,rrrrtolcqs r1test_ormrz_unmrzrFs$ ))   %CRBB' *+C168!(3 7SZZ)0078A B8$++E2A$45M;@%B !Hkcjj"*SZZ-A"-DDLLUSTARH% B8(V4wr\V ^8H4\\WW4V!W>V4wr\V ^8H4\\!WE4TV^,^8Xd^M^R7EK R#)zy Test Cholesky factorization of a positive definite Rectangular Full Packed (RFP) format array and solve a linear system rr*rXrN)pftrsrdr1r2)r6rr^rrcr,r-rr~r r r#rrrrr)r0rr/rrrBf1Bf2rqrdr1r2r4rnrkr5s r1 test_pftrsrtus ))   %C'   7!#((1."33;;EBAFFHJJ3q6)A!''.ACC!CF("A !Qu %AaC85)AaC85)%57@=B &C"e!H  = 1#.  i:;1!,  !%+t/2Qw!|! D D7(r3c \PPR4p\\4EFwr^pV^8dvVP W34VP W34R,,P V4pWDP4P,V\V4,,pMIVP W34P V4pWDP,V\V4,,pV^8dRMRp\RRV R23VR7wrgpV!V4wrVP V^4P V4p V!V^R V ^V 4p V!W<4wr\T \V PV P4P4)^V,,4V^,^8Xd^M^R 7EK R #) zL Test for performing a symmetric rank-k operation for matrix in RFP format. rr*rhr1r2frkrXrNr)r6rr^rrcr,r-rr~r r#rrr)r0rr/rrprefixr1r2shfrkr4r5rAfp_outA_outs r1test_sfrk_hfrkr|sh ))   %C'   7!#((1."33;;EBAFFHJJ3q6)A!''.ACC!CF("AQwC.VHC./Q5:<eq HHQN ! !% (1b!Q,$!%quuQVVXZZ/@.@1Q3.F)G/2Qw!|! D D#(r3c \PPR4p\\4EFwr^ pV^8deVP R ^W334VP R ^W334R,,P V4pWDP4P,pMLVP R ^W334P V4pWDP,V\V4,,p^d\P!V!R4P4,p\RVR7wrgp\W^R7p \V^RR7wrp \W^R7p V!V^V R7wrpV!W^R7wppp\\!VR4\!WR 3,R4VR R 7\V^RR7wprV!V^R7wrpV!W^R7wppp\\#V^4\#VV R 3,^4VR R 7EK R #)zh Test for going back and forth between the returned format of he/sytrf to L and D factors/permutations. rr*rrXrF)r hermitianrr)rr/Ni)syconvr rr)r6rr^rrcrandintr-rr~r r(rr#rrrrr)r0rr/rrrDrtrf trf_lworklwrLDpermrr rnrgrrYs r1 test_syconvrs ))   %C'   7S"qf-S"qf-b0128&- FFHJJA CaV,33E:ACC!CF("A"**U3Z__--!13BIN"PY I 2!u5 d I 2aq34CQ/ 1dQ d11W:r&:2N!u5 1aq/4CQ/ 1dQ D47Q$7cKK5(r3c0a]tRtRtoRtRtRtRtVtR#) TestBlockedQRizX Tests for the blocked QR factorization, namely through geqrt, gemqrt, tpqrt and tpmqr. c \PPR4p\\4EFwr#^pV^8d?VP WD4VP WD4R,,P V4pM VP WD4P V4p^d\P!V!R4P4,p\RVR7wrxV!WE4wrp V ^8XgQh\P!V R4\P!WCR7,p \P!WCR7W,V PP4,, p \P!V 4p\V PP4V ,\P!WCR7VRR7\W,WVRR7V^8dAVP WD4VP WD4R,,P V4pRpM"VP WD4P V4pRpRFpR V3FpV!WVVVR 7wpp V ^8XgQhVV8XdV PP4pMT pVR8Xd VV,pM VV,p\VVVRR7VV3R8XgKsV!WV4wpp V ^8XgQh\!VV4K K \#\$WWR R 7\#\$WWR R 7EK R#)rr*rrXrr/rr~rLrXrrrrr0N)geqrtgemqrtrrLrrLrX)r6rr^rrcr,r-r(rr#rr r~rrrrrr)rfr0rr/rrrDrrrgtrnrrrr transposerrrqqC c_defaults& r1test_geqrt_gemqrtTestBlockedQR.test_geqrt_gemqrtsdii##D)#F+JCAQwXXa^chhqnR&77??FXXa^++E2bjjs11C,-@NMEqJA$19 92!77Aq&);;A A ACCHHJNBFF1,B!# % AE1R 8QwXXa^chhqnR&77??F XXa^++E2 "!9-E$Q14uEGAt19$9 )CCHHJs{UU#Ar"=e} 2*0q/ 4#qy(y$Y2).#0 )V3 ? )VC @ @u,r3c \PPR4p\\4EFGwr#^pV^8d|VP WD4VP WD4R,,P V4pVP WD4VP WD4R,,P V4pM@VP WD4P V4pVP WD4P V4p^d\P!V!R4P4,p\RVR7wr^V^,V3EF0p V!WWV4wrrV^8XgQh\\P!V R4\P!VR44\\P!WV, ^, 4\P!WjV, ^, 44\P!WjV, 4\P!WV, 4pp\P!\P!WCR7V34p\P!^V,VR7VV ,VPP!4,, p\P!\P!V 4\P"!V 434p\%VPP!4V,\P!^V,VR7VRR7\%VV,\P!\P!V4V34VRR7V^8d~VP WD4VP WD4R,,P V4pVP WD4VP WD4R,,P V4pRpMBVP WD4P V4pVP WD4P V4pRpREFpR V3EFpV !WV VVVVR 7wpppV^8XgQhVV8XdVPP!4pMTpVR8Xd?\P!VV3^R 7p\P!VV3^R 7pVV,pM=\P!VV3^R 7p\P!VV3^R 7pVV,p\%VVVRR7VV3R8XgKV !WV VV4wpp pV^8XgQh\VV4\V V4EK EK \'\(WWVVR R 7\'\(WWVVR R7EK3 EKJ R#)rr*rrXrr/rr~rLrXrrrrr0N)tpqrttpmqrtrrr)r6rr^rrcr,r-r(rr#rrrr8r r~rrrrr)!rfr0rr/rrrrDrrlrgrrrnB_pentb_pentrrrrrrrrrrrcdCDqCDr d_defaults!& r1test_tpqrt_tpmqrtTestBlockedQR.test_tpqrt_tpmqrts(ii##D)#F+JCAQwXXa^chhqnR&77??FXXa^chhqnR&77??FXXa^++E2XXa^++E2bjjs11C,-@NMEa^ %aA 1 aqy yRWWQ^RWWQ^<RWWQA 2BGGA1uqy4IJ"$E!2BGGA1u4ENNBFF1$:F#CDFF1q5.Q1CCNNBGGAJ a0@#AB  Qq1uE0J%(r3Ar~~rwwqz66J'K%(r37!#((1.*;;CCEJA!#((1.*;;CCEJA #I!//6A!//6A #I&D"%y!1%+A!Q27&9 1d#qy(y I- ! A !A3;!#AQ!?B!#AQ!?B"#b&C!#AQ!?B!#AQ!?B"$q&C'CcC %=J69?aA9N6Iy$#'19,9(A6(A63"2':iA!QSIiA!QcJJA$,r3rN) rrrrr>rrrrrs@r1rrs <A|QKQKr3rc\PPR4p\\4EFwr^ p^p\ RVR7pV^8d|VP W3V, 4PV4RVP W3V, 4PV4,,pWfP4P,pM9VP W3V, 4PV4pWfP,pV!V4wrxr\V4p RWV, R1W, R13&\V ^4R\P!\P4P,p R\P!\P4P,p VR 9dT MT p\!Wh^, ,RV^, 3,V P4PV ,RVR7V!V^R 7wrxr\#V4pRWV, R1W, R13&\V ^4R\P!\P4P,p R\P!\P4P,p VR 9dT MT p\!Wh^, ,RV^, 3,WP4P,RVR7EK R#) rpstrfrXr*rNr)rQrrr[r6rr^rrcr#r,r-rr~rrrrrrxrr)r0rr/rrrrrpivr_crnrY single_atol double_atolr0rLs r1 test_pstrfro! ))   %C'    6 7aC ''.chhqA#6F6M6Me6T1TTAFFHJJAaC ''.ACCA!!H G "'(CGH T1RXXbjj1555 RXXbjj1555 !Vm{a%CE*AFFHJJN$O!!1- G "'(CGH T1RXXbjj1555 RXXbjj1555 !Vm{a%CE*A N$OOG(r3c\PPR4p\\4EFwr^ p^p\ RVR7pV^8d|VP W3V, 4PV4RVP W3V, 4PV4,,pWfP4P,pM9VP W3V, 4PV4pWfP,pV!V4wrxr\V4p RWV, R1W, R13&\V ^4R\P!\P4P,p R\P!\P4P,p VR 9dT MT p\!Wh^, ,RV^, 3,V P4PV ,RVR7V!V^R 7wrxr\#V4pRWV, R1W, R13&\V ^4R\P!\P4P,p R\P!\P4P,p VR 9dT MT p\!Wh^, ,RV^, 3,WP4P,RVR7EK R#) rpstf2rXr*rNrr)rQrrr)r0rr/rrrrrrrrnrYrrr0rLs r1 test_pstf2rrr3c\P!.RO.RO.R O.R O.4p\P!.R O.R O.R O.4p\\4EF(wr#V^8d5\P!.RO.RO.RO.RO.4pVP V4pM_\P!.RO.RO.RO.VR7pV\P!.RO.RO.RO.4R,, pVP V4p\ RVR7pV!V4wrgrrV^8d6\ VP V4VR,V,V,^RR7K\ VP V4VR,V,V,^RR7EK+ R#)?r*rXgeequNrrQ)rrg1w-!?gd`TRۿ)rgsrr)gs?rg2%䃮g,eX)rgsFg%ug??)y/nҿ&?yDioɴ?Af?y o_[ Acп)ysֿAfҿyPkw?JY8y5;NёCl?)yYڊ?1*?y=yXѿ@a+?yhoſFx)gЈBgtBgffffff@gٓ)@gg#fDgffffff)gHzG?gQg'Vgp= ף)g(\rgS7нr)gq= ףpgAg(\)g333333gBg333333ÿ)gZ9=gQgֽ)gffffff@gtޅBr2)g(\g Zgq= ףp?)gEop=gQ?gZEqҽ)r)N)r6rrrcr-r#r) desired_real desired_cplxrr/rrrrrowcndcolcndamaxrns r1 test_geeqursi88>@@@BCL 88111 23L '  7CEEEGHAA:::>!,AdqBGGAJS ! rd 2bT!V ;< 266"))B*++b01AAdGAdG!>>!((2<<*@JAdqBGGAJ IJr3c 8\PP^*4p^ pVPV4pVPV4VPV4R,,p\ \ 4EF%wrEV^8d=VPW4pVP V4pWb,pVP V4pMXVPW4VPW4R,,pVP V4pWc,pVP V4p\RVR7p\RVR7p V!V^R7wrrV !WW^R7wrV^8d'\VP V4W, ^R7EK\VP V4W, ^R7EK( R#) rr*getc2rXgesc2r,) overwrite_rhsrN) r6rr^r,rrcr-r#r)r0rrrrr/rrrrrr jpivrnrrZs r1test_getc2_gesc2rsR ))   #C A88A;L88A;!R/L'  7AA AA!!22AA AA 6 6$QA6$!< 7 %l&9&9%&@&'gq : : &l&9&9%&@&'gq : :+(r3rjobajobujobvjobrjobpc 6\PP^*4pVwr^d\P!V4P,p \ WV4p \ RVR7p V^8pV^8pV^8H;'dW8Hp\P!V 4pV^8H;'dV'*;'dV'*pV^8H;'dT;'dT'*;'dTpV^8H;'dT;'dT'*;'dTpV'dR pMV'g V'dR pM^pV^8dV^8Xd\\WW#WEWv4 R#V !V VVVVVVR7wpppppp\VV4V'EgV^,V^,, VRV ,p\V\V RR7V R7V^8XdVRRV 13,pV'dPV'dH\V\P!V4,VP4P,WR7V'dB\VP4PV,\P !V 4V R7V'dB\VP4PV,\P !V 4V R7\V^,\P"P%V 44\V^,\P&!V44\V^,^4R#R#) aTest the lapack routine ?gejsv. This function tests that a singular value decomposition can be performed on the random M-by-N matrix A. The test performs the SVD using ?gejsv then performs the following checks: * ?gejsv exist successfully (info == 0) * The returned singular values are correct * `A` can be reconstructed from `u`, `SIGMA`, `v` * Ensure that u.T @ u is the identity matrix * Ensure that v.T @ v is the identity matrix * The reported matrix rank * The reported number of singular values * If denormalized floats are required Notes ----- joba specifies several choices effecting the calculation's accuracy Although all arguments are tested, the tests only check that the correct solution is returned - NOT that the prescribed actions are performed internally. jobt is, as of v3.9.0, still experimental and removed to cut down number of test cases. However keyword itself is tested externally. gejsvrX)rrrrjobtrNF)r+rr)rrS)r6rr^rrr2r#rrrrrrr6rr~identityr matrix_rank count_nonzero)rr/rrrrrrr0rrr0rrlsvecrsvecl2tran is_complexinvalid_real_jobvinvalid_cplx_jobuinvalid_cplx_jobv exit_statussvarrrrrnsigmas&&&&&&&& r1test_gejsv_generalr.spB ))   #C DA %$$ $D#D5A WE 2E 1HE 1HEai % %afF#JHHUHHZMMU-=-=v(>MM:MMU-=-=v(>MM:  /   qtqyi4t4N',Q262626262626 (8$Q4 T;'{!WtAw&#bq'1E E3qU#;$ Gqya!eHBGGEN 2QVVXZZ ?N Q ATJ Q ATJ q299#8#8#; < q2#3#3E#: ; q1 %5r3c\RVR7pV!R4wr#rErg\V^4\VPR4\VPR4\V\P!R.VR74\P !RVR7pV!V4wr#rErg\V^4\VPR4\VPR4\V\P!R.VR74\P !RVR7pV!V4wr#rErg\V^4\VPR4\VPR4\V\P!.VR74\P !\P!^d4P^ ^ 44PV4p\P!WP,4pVPR4p V!V4p \W4R#) z*Test edge arguments return expected statusrrXrrNrVrV)rVr)r#rr.r6rr sinrrMr-asfortranarrayr~rr) r/rrrrrrrnrAcr5s & r1test_gejsv_edge_argumentsrs WE 2E$)9 CAUq&!&!bhht512 E"A#(8 CAUq&!&!bhht512 e$A#(8 CAUq&!&!bhhr/0 ryy~%%b"-.55e Q* Q+1 Q  Q  Q( + Q; z1du, d, dl, du_exp, d_exp, du2_exp, ipiv_exp, b, xg@gffffff@g%@g@rg3@c \RV^,V^,34wrV !W!V4wrrpp\W4\W4\WRR7\W4V !WWW4wpp\VV4R#)rrrNr)r#r)rrrdu_expd_expdu2_expipiv_exprrrrrrrrr rnrs&&&&&&&&& r10test_gttrf_gttrs_NAG_f07cdf_f07cef_f07crf_f07csfrKssd$$6A1GLE$)"$4!CStTC!C BD)D##3T5MGTGQr3c \PPR4pVPV4VPV4R,,pVPV^, 4VPV^, 4R,,pVPV^, 4VPV^, 4R,,p\P!V4\P!VR4,\P!V^4,p\P!V\P 4'd2VP VP VP VP 3wrtrVVPV4VPV4VPV4VPV43wrtrV\PPYqR8Xd\PM^R7p\R V34wrV !WTV4wrTrkrV !WTWkWVR7wr\R V34wppV!V4wppp V!VWR7wpp\P!V4PR,p\VVVR7R#) iTfr*r~rrg?rNr)rgtconr)r6rrr6rrrr-rrrr#rrr)r/rrr0rrrrrrrrr rnrEr5rrrrrrs&&& r1 test_gtconrs ))   )C 1  1 b((A AE SZZA.r1 1B AE SZZA.r1 1B  RWWR_$rwwr1~5A }}UBKK((vvqvvrww7 b88E?AHHUORYYu5EryyQVGWWLA" IINN1CK"&&QN ?E#$6=LE!&rb!1B2D 2"4T :FC#$6=LE51XNBd 2u (FC 88E?   $DC4(r3cV\RVR7pVwr4V!W4R7wrV\V^4R#) geqrfp_lworkrXrNr)r/r.rrrrrns&& r1test_geqrfp_lworkrs-$^EBL DA(KEqr3z ddtype,dtypec\\PP^*4p^d\P!V4P,p^ p\ V3W4^,p\ V^, 3W4p\P !V4\P !VR 4,\P !\P!V4^4,pVP4VP4.p\RVR7p V !WV4wrp \WX^,4\Wh^,4\V ^RV R2R7\P !V R 4\P !\P!V44,p \P !V 4p\W}V,V P4P,VR7\ V3W4pW,p\RVR7pV!WP4V4wpp \V ^RV R2R7\VVVR7R #) rpttrfrXzpttrf: info = z , should be 0)err_msgrpttrszpttrs: info = Nr)r6rr^rrr2r6rrr#r rr rrr~)ddtyper/r0r0rrrrrrr_ernrLrrrr_xs&& r1test_pttrf_pttrsr!s ))   #C rxx"" "D A $QD&6:A#QqSFE7A  RWWQ^#bggbggaj!&<ggaj2771a=(2771b>9 Hrwwy{{ "  (h < '1 > E GGAJA &B 7 GGAJB '"''!Q- 7 !<r3zdtype,realtyper.cR\P!V4P,p\RVR7p^ p\ WWR4wrgrV!WgWR7wrr\ V ^RV R24\ \P!\V4^,4\P!V 4VR7V'd\ V \P!V 4P,\P!V4VR7\ V \P!V 4,\P!V 4P,WR7R#R#) z Tests the ?pteqr lapack routine for all dtypes and compute_z parameters. It generates random SPD matrix diagonals d and e, and then confirms correct eigenvalues with scipy.linalg.eig. With applicable compute_z=2 it tests that z can reform A. rpteqrrXrrr2r.zinfo = z, should be 0.rN) r6rrr#r3rrsortrrr~rr6)r/r-r.r0r5rrrrr2d_pteqre_pteqrz_pteqrrns&&& r1 test_pteqrr;L s ## #D ge 4E A"5AAJA!&+a&O#GgqGD689BGGDGAJ')9E"'''"2"4"44bkk!n! # "'''"22RWWW5E5G5GG & r3c~\RVR7p^ p\WWB4wrVrxV!V^, WhVR7wrrV ^8gQhR#r5rXr2r.Nr#r3 r/r-r.r5rrrrr2r8r9r:rns &&& r1test_pteqr_error_non_spdrAl sI ge 4E A"5AAJA!',AE1Y&O#Gg !8O8r3c \RVR7p^ p\WWB4wrVrx\\W5RRWhVR7\\W5VRRWR7V'd\\W5WhRRVR7R#R#)r5rXNr>r)r#r3rrT) r/r-r.r5rrrrr2s &&& r1"test_pteqr_raise_error_wrong_shaperCz si ge 4E A"5AAJA!*esVQyI*e#2!Ij%A3B9Mr3c\RVR7p^ p\WWB4wrVrx^V^&^V^&V!WVWR7wrrV ^8gQhR#r=r?r@s &&& r1test_pteqr_error_singularrE sS ge 4E A"5AAJA! AaD AaD&+AA&K#Gg !8O8r3zcompute_z,d,e,d_expect,z_expectctRp\RVPR7p\P!V4\P!V^4,\P!VR4,pV!WWpR7wrr\ WVR7\ \P !V 4\P !V4VR7R#)zV Implements real (f08jgf) example from NAG Manual Mark 26. Tests for correct outputs. rr5rXr6rNr)r#r/r6r6rr) r.rrr(z_expectr0r5r2rr_zrns &&&&& r1test_pteqr_NAG_f08jgfrI s| D gagg 6E  RWWQ]"RWWQ^3Aq@BBBt,BFF2Jx 0t "$r3c\P!.4p\RVPR7p\ \ W4R#)rLrXN)r6rr#r/rr)A_emptyrLs r1#test_geqrfp_errors_with_empty_arrayr[ s)hhrlG hgmm \V4F-pVVVV,V,V,VV,V,3&K/ \V4Fp\P!VV,4VVV,VV,3&\P!VV,4VVV,V,VV,V,V,3&\P!VV,4)VVV,VV,V,V,V,3&\P!VV,4VVV,V,VV,3&K VV,V,p\VVR R \P !V4P",R 7R#) rqrrcsdrsrXrlrworkNrrg@rQ)rr rvsr!r#rdictrrr`r6r r9cosrrrr)rorrrr]Xdrvrtrulwvalscs11cs12cs21cs22thetau1u2v1tv2trnrYVHrn11n12n21n22SonerDXcs& r1test_orcsd_uncsdr# s$2s!AK'$TC!Tk }/@/@/CAucK.? @OHC 31 (E!$gu $sG?ra=)qQ-(ra=)qQ-(AvD) GE9 %T 1,UG 46 DJJqMQWWQZ 'N4::a=/J112 +>@ DJJqMQWWQZ 'N4::a=/J112 +>@r3cn\PP^*4p\R VR7wrE^ p\ V^, 3WR7p\ V3WR7p\ V^, 3WR7p \P !VR 4\P !V4,\P !V ^4,p \ V^3WR7p V\ 9dRMRp V'dV P4PMT V ,p VR8Xd V!WxV 4M R.^,wrppppV!WxWW,WVVVR7 pVw ppppppppppVR8Xd1^VR &^VR &V!WxW4pVw ppppppppppV^8gQR 4hR#VR8Xd<^VR &^VR &^VR &V!WxWW.VVVVR 7 pVw ppppppppppV^8gQR 4hR#R#) rrXrr~rrKNrrXz&info should be > 0 for singular matrix)rrrrrr rr) r6rr^r#r2r6rrr~)r/rrr0rrrrrrrrrrrrrrrrrrrrrr rrrrrns&&& r1test_gtsvx_error_singularr s ))   #C#$6eDLE A $acV5 BB#QD?A $acV5 BB B"''!*$rwwr1~5A#QF%AAK'CSE!qA-A!CKbRdVAX)DtUE5bRE7I@I=CS$feT4 s{"2"' DMARdD&%tTaxAAAx BRb "T#%1 DMARdD&%tTaxAAAx r3c\PP^*4p\R VR7wrE^ p\ V^, 3WR7p\ V3WR7p\ V^, 3WR7p \P !VR 4\P !V4,\P !V ^4,p \ V^3WR7p V\ 9dRMRp V'dV P4PMT V ,p VR8Xd V!WxV 4M R.^,wrppppVR8Xdo\\WGRR WV W,WVVVR7 \\WGVRR WW,WVVVR7 \\WGWRR V W,WVVVR7 \\WGWV RR W,WVVVR7 R#\\WGWV W,VRR VVVVR7 \\WGWV W,WRR VVVR7 \\WGWV W,WVRR VVR7 \\WGWV W,WVVRR VR7 R#) rrXrr~rrKNrXrrr) r6rr^r#r2r6rrr~rrTr)r/rrr0rrrrrrrrrrrrrrrrs&&& r1"test_gtsvx_error_incompatible_sizer s ))   #C#$6eDLE A $acV5 BB#QD?A $acV5 BB B"''!*$rwwr1~5A#QF%AAK'CSE!qA-A!CKbRdVAX)DtUE5 s{j%CR!$E 7 j%QsVR$E 7 j%Q3B$E 7 iA1Sb6$E 7 j%QA$s)E 7 j%QA$s8E 7 j%QA$s)U < j%QA$E#2JU r (8?@ww)HHI)Nr3cR#N)NNNrrs&&r1rr -?r3c :\PP^*4p^d\P!V4P,p\ RVR7p^p\ V3W4^,p\ V^, 3W4p \P!V4\P!V R4,\P!\P!V 4^4,p \ V^3WR7p W,p V!W4wrpVP4V P4V P4.pV!WWWR7wrppppp\VV^,4\V V^,4\V V^,4\V^8HRV R24\V V4\P!VR4\P!\P!V44,p\P!V 4p\V VV,\P!V4P,VR7\!VR4'd QR V 24h\VP"R8HR VP" R 24\VP"R8HR VP" R 24R #)a This tests the ?ptsvx lapack routine wrapper to solve a random system Ax = b for all dtypes and input variations. Tests for: unmodified input parameters, fact options, incompatible matrix shapes raise an error, and singular matrices return info of illegal value. ptsvxrXrrrefzinfo should be 0 but is .rrrrz# but should be ({x_soln.shape[1]},)rNr)r[)r6rr^rrr#r2r6rrr rrr rr~rr.)r/r-r df_de_lambdar0r0rrrrrrrrrrnrrrrrrLrs&&&& r1 test_ptsvxr s ))   #C %$$ $D WE 2E A#QD(81$**>8"8: DJJ$ >$**>8"8:r3c<\RVPR7!W4#rrrs&&r1rrS rr3cR#rrrs&&r1rrU rr3c D\PP^*4p\RVR7p^p\ V3W4^,p\ V^, 3W4p\P !V4\P !VR4,\P !\P !V4^4,p \ V^3WR7p W,p V!Wx4wrp\\WWRRWW,V R7\\WWVRRWWR7\\WWWRRW,V R7R#)rrrXrNrr) r6rr^r#r2r6rrrTr)r/r-rrr0rrrrrrrrrrns&&&& r1test_ptsvx_error_raise_errorsrO s ))   #C WE 2E A#QD(81mU1X1Xm>@ B'74($Du Q%(HBq 1 "4 (CB!$/a5)ICq q6$sK;3L+L #L cV\PPR4p\P!V4P^d,p^ p\ W3.WR7p\ RVR7wrVV!RVRR7p\VR ,^4V^,pVR ,p VR ,p V\9d$\V\P!V4^VR7\W,V P4P,V^VR7V!W^^4p\VR ,^4V^,pVR ,p V\9d$\V\P!V4^VR7\W,V P4P,V^VR7\VR ,V ^VR7R#)rrrXcR#Nrrs&r1r!test_gees_trexc.. Dr3FrrQN)geestrexcrrSr_r_rr r6rr^rrr2r#rr+rrrr~) r/r0r0rrgrrrrr2d2s & r1test_gees_trexcr sA ))   %C 88E?  s "D A#QF%AA"#4EBKD .! 7FQq Ar A 4B 2771:AD9AEAFFHJJ&= 1A FQq Ar A 2771:AD9AEAFFHJJ&=AdGRad3r3zt, expect, ifst, ilstcRp\RVPR7pV!WW^R7p\VR,^4V^,p\W0VR7R#)z[ This test implements the example found in the NAG manual, f08qfc, f08qtc, f08qgc, f08quc. rrrX)wantqrNr)r#r/rr)rifstilstexpectr0rrs&&&& r1test_trexc_NAGr sH< D WAGG 4E 11 -FQq AFD)r3c\PPR4p\P!V4P^d,p^ p\ W3.WR7p\ W3.WR7p\ RVR7wrgV!RWERRR7p\VR ,^4V^,p V^,p VR ,p VR ,p V R ,V R ,, p V R ,V R ,, pV\9dG\V \P!V 4^VR7\V \P!V 4^VR7\W,V P4P,V^VR7\W,V P4P,V^VR7V!WW^^4p\VR ,^4V^,p V^,p V^,p V^,p V\9dG\V \P!V 4^VR7\V \P!V 4^VR7\W,V P4P,V^VR7\W,V P4P,V^VR7\V R ,V R ,, V^VR7\V R,V R,, V ^VR7R#)rrrXcR#rrrs&r1r!test_gges_tgexc..< rr3Fr, overwrite_brQN)ggestgexcrrrSrrrr)r/r0r0rrgrrrrrrrr2d1rs& r1test_gges_tgexcr1 s  ))   %C 88E?  s "D A#QF%AA#QF%AA"#4EBKD .!Eu MFQq Aq Ar Ar A 41T7 B 41T7 B 2771:AD92771:AD9AEAFFHJJ&=AEAFFHJJ&= 1q! $FQq Aq Aq Aq A 2771:AD92771:AD9AEAFFHJJ&=AEAFFHJJ&=AdGag%r=AdGag%r=r3c\PPR4p\P!V4P^d,p^ p\ W3.WR7p\ R VR7wrVpV!RVRR7p\VR ,^4V^,p VR ,p V R ,p V\9d$\V \P!V 4^VR7\W,V P4P,V^VR7\P!V4p ^V ^&\W|V 4p V\9d V!WWR7pMV!WWV ^,R7p\VR ,^4V^,p V^,p V\9d$\V \P!V 4^VR7\W,V P4P,V^VR7\V R,V ^VR7R #)rrrXcR#rrrs&r1r!test_gees_trsen..l rr3FrrQrrliworkN)rtrsen trsen_lworkrrSrrr6rr^rrr2r#rr+rrrr~r r)r/r0r0rrgrrrrrr2rselectrs& r1test_gees_trsenra s ))   %C 88E?  s "D A#QF%AA/( 7D.! 7FQq Ar A 4B 2771:AD9AEAFFHJJ&= XXa[FF1I ; 2E v!1v!qBQq Aq A 2771:AD9AEAFFHJJ&=AdGRad3r3z*t, q, expect, select, expect_s, expect_sepg?g(\ @gRQ?gK?c4RpRp\RVPR7wr\WV4p VP\9d V!W WR7p MV!W WV ^,R7p \ V R ,^4V ^,pV ^,pVP\9dV ^,p V ^,p MV ^,p V ^,p \ W1V,VP 4P,VR7\ V^V , VR7\ V^V , VR7R#) zK This test implements the example found in the NAG manual, f08qgc, f08quc. r{Gz?rXrrrN)rrr)r#r/rr+rrrr~)rrrrexpect_s expect_sepr0atol2rrrrrseps&&&&&& r1test_trsen_NAGr  s` D E) 1E ; 2Eww. v!1v!qBQq Aq Aww. 1IQi 1IQiFEAFFHJJ.T:Ha!e%0JCe4r3c N\PPR4p\P!V4P^d,p^ p\ W3.WR7p\ W3.WR7p\ R VR7wrgpV!RWERRR7p \V R ,^4V ^,p V ^,p V R ,p V R ,p V R ,V R ,, pV R,V R,, pV\9dG\V \P!V 4^VR7\V \P!V 4^VR7\W,V P4P,V^VR7\W,V P4P,V^VR7\P!V4p^V^&\VVW4pV^,^,V^,3pV!VWWVR7p \V R ,^4V ^,p V ^,p V R,p V R,p V\9dG\V \P!V 4^VR7\V \P!V 4^VR7\W,V P4P,V^VR7\W,V P4P,V^VR7\V R ,V R ,, V^VR7\V R,V R,, V^VR7R#)rrrXcR#rrrs&r1r!test_gges_tgsen.. rr3FrrQrN)rtgsen tgsen_lworkrrrSrrirr)r/r0r0rrgrrrrrrrrr2rrrrs& r1test_gges_tgsenr sc ))   %C 88E?  s "D A#QF%AA#QF%AA/( 7D.!Eu MFQq Aq Ar Ar A 41T7 B 41T7 B 2771:AD92771:AD9AEAFFHJJ&=AEAFFHJJ&= XXa[FF1I ; 5E1XaZq "E 61U 3FQq Aq Ar Ar A 2771:AD92771:AD9AEAFFHJJ&=AEAFFHJJ&=AdGag%r=AdGag%r=r3za, b, c, d, e, f, rans, lansc dRp \RVR7p V !WW#WE4wrrp\V^4\V R^\P!V4P ^d,RR7\VR^\P!V4P ^d,RR7\WV RR 7\WV R R 7R #) rtgsylrXrzSCALE must be 1.0rr0rrzDIF must be nearly 0zSolution for R is incorrect)r0rzSolution for L is incorrectN)r#rrr6rr)rgrrrrriranslansr/r0rroutloutrZdifrns&&&&&&&&& r1test_tgsyl_NAGr sH D ge 4E#(qQ#: DDqE3QRXXe_-@-@-D/1C1288E?+>+>s+B24DT9;DT9;r3rijobc xV\P8XdRMRp\PPR4p^ ^re\ VP R^ WU.4P V4VP R^ WU.4P V4RR7vrxp \ VP R^ Wf.4P V4VP R^ Wf.4P V4RR7vrp VP R^WV.4P V4p VP R^WV.4P V4p \RVR7pV!WzWWWR7wpppppV^8XgQR4hVR 8gQR 4hV^8Xd7\VR ^\P!V4P^d,R R 7MVR 8gQR 4hV^8:EdVR8Xd@W,VV ,, pVV ,pW,VV ,, pVV ,pMVR8Xd\P!V4V,\P!V4V,,pVV ,pV\P!V 4,V\P!V 4,,pRV,V ,p\XXVR RR7\XXVR RR7R#R#)gMbP?g|=lOElt/r)outputrrX)rrzINFO is non-zerorzSCALE must be non-negativezDIF must be 0 for ijob =0rzDIF must be non-negativerXr~zlhs1 and rhs1 do not match)r0rrzlhs2 and rhs2 do not matchNirr) r6rrrruniformr-r#rrrr)r/rrr0r0rrrgrr5rrrrirrrrZrrnlhs1rhs1lhs2rhs2s&&& r1 test_tgsylr#G sl BJJ&4ED ))   0 1C rq#++c2v.55e<++c2v.55e<!HA1#++c2v.55e<++c2v.55e<!HA1 BA6"))%0A BA6"))%0A ge 4E#(qQ/4$A D$sD 19(((9 C<555< qySqrxx/B/B3/F ; =cz555z qy C<8dQh&D7D8dQh&D7D c\<<?T)BLLOd,BBD7D",,q/)D2<<?,BBD:adB < >r3mtypecXVR8Xd"V\9d\P!R4\PP R4p^^rTV\ 9dCVPWD3R7VPWD3R7R,,PV4pM"VPWD3R7PV4pVR8XdWfP,MWfP4P,pVPWE3R7PV4pV R2V R2V R23p\WR 7wrp V !WBR 7p V !WlVR 7wrpV^8XgQhV !WWrR 7wppV^8XgQh\P!V4Pp\VV,V^dV,V,R 7R#)r_zhetrs not for real dtypes.*M/t|0rr*r^rrtrsrXrr )rgr rrrN)rr?r@r6rrr+rr-r~rr#rrr)r$r/rr0rrrrrGrrr'rrr rnrrs&&& r1 test_sy_hetrsr(z sk }+- 01 ))   0 1C!t  [[qf[ % ! (@(C C K KE R KKaVK $ + +E 2D=CCa&&(**nA ! #**51AgSMeWI.5' =E*5>CC a %E!6OCt 199Ca5GAt 199 ((5/  CAaCQs+r3cVR8Xd"V\9d\P!R4\P!4R8Xd"V\ 9d\P!R4\ PPR4p^pVPWD34VPWD34R,,p\ P!V\ P4'd VPpVPV4pVR8XdWUP,MWUP4P,pV R2V R23p\W`R 7wrxV!WQR 7wrp V ^8XgQhV!WVR 7wrV ^8XgQh\ P !V4P"p \ P$P'V4pV'dG\)\ P*!V 4\ P*!V4^dV,V ,R 7R #\)\ P,!V 4\ P,!V4^dV,V ,R 7R #)r_zhetri not for real dtypes.z win-arm64z1Test segfaulting on win-arm64 in CI, see gh-23133r&r*r^rr1rXr)rgr rrN)rr?r@ sysconfig get_platformr+r6rrrrrr-r~rr#rrrrrrr)r/rr$r0rrrGrr1rr rnA_invrrs&&& r1test_sy_he_trir- s }+- 01;.5N3J GH ))   0 1C A A6SZZ/22A }}UBKK(( FF AD=CCa&&(**nAgSMeWC= (E3HC!)OCt 199e4KE 199 ((5/  C ))-- C  3q59E 3q59Er3rz uplo, m, nc\PPR4pVPW#3R7PV4p\ RV34wrV!WWR7p VR8Xd\P !V4M\P !V4pVR8Xd&\P!\W#44p ^W{V 3&V !W4p \WRR7R#)l8#q9 r)rOr6rY>rN)lantrr) r6rrr-r#rrrr`r) rrOrrr6r/r0rr0rrErDrs &&&&&& r1 test_lantrr1 s ))   1 2C  &&u-A#$6=LE d .Cck rwwqzA s{ IIc!i Q$ .CC4(r3c  \PPR4p^ ^^rTpVPW334VPW334R,,pVP\P!^VP ^^ 4V44pWg, p\P !V\P4'd VPpVPV4p^V\P!W4^,4&^V\P!W5)^, 4&\WTV4p^V,V,^,p \P!W3VR7p WVR1R3&\PPYaR8Xd\P M^R7p \#R V 34wrV !WV4wrpV !WWNVV R7^,p\#R V34wppV!V4^,pV!VWR7^,p\%VV^R 7R#) 2 r*rXNr)r~r)rrdrcrNr rrr)gbcongbtrfr)r6rrrrrrrrr-rrrr rrrr#r)r/rr0rrcrdrrtmpLDABrNrr4r5lu_bandr r5rErrrrs&& r1 test_gbconr9 s ))   ,CAq2A A6SZZ/22A \\"++aa)LERR(G1 TR$  C$$6=LE5 q!B E %a (CC1%r3c2\PPR4p^ ^^rTpVPW334VPW334R,,p\P!V\P4'd VP pVP V4p^V\P!W4^,4&^V\P!W5)^, 4&\WTV4p\RV34wrV !W4p V!WWG4p \WRR7R#)r3r*r/rN)langbr) r6rrrrrr-rrrr#r) r/rr0rrcrdrrNr;rrrEs && r1 test_langbr< s ))   ,CAq2A A6SZZ/22A }}UBKK(( FF A$%Abooaa !%&Abooaq!" BA B#$6=LE .C " !CC4(r3 compute_vc\PPR4p^ pVPW0R7pVPV^, VR7p\P!VR4\P!V4,\P!V^4,p\PP V4p\ R4pV!WEVR7wrp V ^8XgQh\\P!V 4\P!V44V'dc\P!V4Pp \V \P!V 4,V P,WlR,R7R#R#)lJ*Os{UrXstevd)r=rrNr) r6rrr6reigvalshr#rr7rrr~) r/r=r0rrrrrr?rYr+rnrs && r1 test_stevdrAs ))   2 3C A 1 "A 1q5 &A 2#bggam3A ))  Q C W %Eqy1JA$ 199BGGAJ -hhuo!!BGGAJ,a3h?r3)rq)r^r[rC)r~1O)r_r^rr3)g)\(@gp= ףgffffff?g ףp= )gQ?gQgGz?g(\)gQ޿gQgGz?gzGʿ)gQ?gQ?gHzG?g)\(?)ggq= ףp@g333333r)ףp= ?g(\rg(\)g cZB#@gI.!v @g?ܵ?r)gC?g=yX5gc=yXga4?)gB`"?g:pΈҞgʡE?gn4@?)g[B>٬?g٬\m?gJ{/L?g Oe?)gc]Fgꕲ q׿g\m?fc]F)g؁sFڿgZB>?g0L F%?gq= ףp)g?gR!u?guVſg&Sٿ)gǘ?gV-g ^)p?g( )gFx $g6[ ٿgUN@giq?)g1Zd?gOnӿgΈ ?g_vO?)g}?5^Iؿg58EGr?gio?g7[ A)g@rffffff?r)r1rffffff@)333333 @ @rr)rrrJrK)rLrMr`rgC>r)rrHra)r[r\r]r^r^g gffffff&rrSr@@???)rPrQrRrS)?rTffffff @333333ӿ333333@ffffff ????)rXrY@rY)rTrUrVrW)rXrYrZrYy~:pffffff?)rQrRrSy333333@y@@y333333 @3333332@y333333yffffff-ffffff#@y333333yfffff?@y333333"@y𿚙?yffffff(@y@y?@y@@yy@)rDr~rE))r\r`)r`r\rB)r]r#r^)rrra)r]rbr\rV)rOgK=Urr)r^).)y0@0@y2@"?)r^rbrVr])rYrPrbyP@0@y0@y@W@O@yN@PyS@TyQ@Ry,@;yA@.@y)gp= ף@rgq= ףp?rF)g\(\ @g ףp= g?)gŏ1w- @gR'?g/n?g&䃞ͪ?)g cZB>?gCl?g:pΈڿg??)gaTR'?gSۿg}гY?g%uο)g\mg٬\m?gAf?gL F%u)gǘgŏ1w-!?g333333?gz6?)rq)r`r_r)evevdevrevx)rVr#dr)rg)\(rgQ?)r皙rffffff?)rgrhg?)rrrrh)rhlV}gV_?g|?5^?)g?rhgV/?g;On?)rrrggj+)y ףp= ? ףp= ׿yRQȿQ?y)\(?п)r@yQ ףp= yq= ףpͿp= ף?)rr @yGz?(\?)rrr@)rly1%Ŀ q?ys??ܵ|ȿyHzG??ܵ?)rrmyV/?ݓ?yjt?vտ)rrrnyB>٬?=U?)rrrrk)g/$?g QIg~jtx?gJ4?)r58EGrgGr?gyX5;?)rg?߾rogt?)rrrg yǹ)g؁sF?g_L?gGz?gUN@?)goT?g0*g' gz6>W)g(g&䃞ͪӿgbX9ҿg-!lV?)gb=y?gۊe?rjg8EGr?)rig?gQg(\ſ)g ףp= ?gQ?rGri)g)\(ܿgQտgQg(\?)rg{GzԿgp= ףg)\(?)rVrrrV)yqhyfc]F?ڊe׿yMbȿ&S?y&1??п)ry?5^I @yo0*yZd;OͿ~:p?)rryx $(@4@y[ A?&?)rrry?ܵ@St$)y?ܵ꿽R!uy2U0*6[?yV-?=yXy8m4?1%̿)ySt$?\mҿyʡE?S㥛?y~:p cڿyK7A`?[ A?)y:pΈ~jtԿyH}?9#J{yH}? cZy+eXw? -ٿ)y"u? c?y?տN@ayRQȿ{GzĿyh"lxz?EGrǿ)y47 )yS!uq F%u @y yտGx $(?y3ı.n?rh|)yv? F%uyd`TR?I&†ۿyN@?ݓy4@ @ ^)?)ys{ @ o_yH.@|Pk @y 0*?*:Hy]m{?Gz)y) 0[rsH //*!!!+.2$0,### +,V4V< *+F3I> zz2::& ,, . ~ %*O.,]-]-@p9p9f&)"BCD*--*;;(;;DK/K/^&)1#s,#s,q"g& )'--2* )F).'4T'6T48AGAGHLL^.?bA63>1>FD&DR7!t7:t&8R36D@"DJD4 LFUKUKp%PP%PP*/Z <DKMKMK$:>&&!12+.eAh0eAh0eAh0!Q(!Q(b&))111/3b&J+. / F(!1+!1+!1+!1+!1+!1+ '11988%?%?%@%=%?%? %AB 88$DE88%G%F%E%H%G%G %IJ 88%G%H%F%G%IJKLM$ ,%M$ , &)XK*XKvL88$9:XX&@AXX&;<XX&9:XX&BCXXl3XXo6XXSz(,d|(+T{(+U|(+T{ ')*  XXAw()2w()2w(*Bx(*Aw '012& XX&GHXX'45 XX&FGXX'ABXX':;XX&?@XXo6XX*'=(3\'B(4i'@(3['A(1:'> '@A  XX'7(.'7(.'8(/'8(.'7 '9:%%)* V W* V (&)1q"g&)'2*)2&)"BCD*c+ ";VDEG,&G,&^c+ ";VDEG0G0S{!:FCDFMFM&A((#56((?3((#45((#89((QGaVaVb"X%&G$-.((S!Hq"g2wQ%&G$-. ((#34((#>?((=1((#56((VW$57G%+V$4vv6F$HI((T5MD$<$%)4L$23 E&1'&1.@)c&+ *CDEGd58n5&6G&:)c&+ *CDEGd58n56G)c&+ *CDEGd58n5N6GN)c&+ *CDEGd58n56G:88$;<88$5688$DE88%F%G%G%F%HI J  =  =&)(@A=$B*=$@.#>?t - L.@ LD%=1t - L.2 L6*01 '2+ '6*+C++C\&)t}5#s,7@-6*7@t&)1v.#s,-B-/*-B`&(+t}5#s,-<-6,-<`88$9:88$?@88$9:88c3Z#tsDk&(%[3+%?@88b!Wq"g2wR&(!W%./ 088$EF88%>?88$DE88h %;&1;%?&2H%=&1;%?&/%< %>? 88ff%57G&,g%6&8I&,f%5%78 9:;());())4FK1<==1>,?@,NO ?@BC 1: C@1:h)4FK1<==1>,?@,NO ?@BC K C@K()4FK1<==1>,?@,NO ?@BC " C@"J88$6788O488aWq!fq!fr2h&'W%./88c1X2wB&(!Wq"g%78 9 88$4588$?@88h %:&.%9&.%9&.%9%;<88fg%6&,f%5&,f%5&,f%5%78 9 :;" *#;" *5$-0&)+=*1+=\&)!4*!4Hhh))*)+, hh211134  hhJKHHJK hh5556 78   0 *10 * &),>*,>^&)*4**4Z0hh112134 hh033134 hh213245 hh|  hh6755 67 hh5656 78 hh6555 67 hh|3')T5U)T5B&)5>*5>p"hh((((*+ hh((((*+ hh((((*+ hh((((*+ hh((((*+ hh((((*+ hh((((*+ hh((((*+9,-!.D+.;/E!.F;"+.*-1->2./->`4,/&)&),**0,.4,/&)&)F**0F6m!45(EF#s,&) )*-F6 ) &)1!&2*!&Hm!45&))*6)&+.f-@./@r3