+
    i(                        ^ RI Ht ^ RIt^ RI5 ]P
                  P                  t]P                  t]P                  t	]P                  t
]P                  t]P                  t]! . RO. RO. RO.4      t. ROt]! . RO. RO. RO. RO.4      t. ROt]! . RO. RO. R O. R!O. R"O.4      t. R#Ot]! . R$O. R%O. R&O. R'O. R(O.4      t. R)Ot]! . R*O. R+O. R,O.4      t]! . R-O. R.O. R/O.4      t. R0Ot]! ^R1.^^.R^..4      t. R2Ot]! . R3O. R4O. R5O. R6O.4      t. R7Ot]! . R8O. R9O. R:O.4      t. R;Ot]! . R<O. R=O. R>O.4      t. R?Ot R t!R t"R t#R t$R t%R	 t&R
 t'R t(R t)R t*]+! ^24      R 4       t,R t-R t.R t/R t0R t1R# )@    )divisionN)*c                  2   \         P                  4       p \        p\        V 4      w  r\	        WV4      p\        W4      pV. RO8X  g   Q hV Uu. uF  p\        V^4      NK  	  up. RO8X  g   Q h\        P                  4       p \        p\        V 4      w  r\	        WV4      p\        W4      pV. RO8X  g   Q hV Uu. uF  p\        V^4      NK  	  up. RO8X  g   Q h\        ^4      p V P                  4       p\        V ^R7       W8w  g   Q hR# u upi u upi )   )	overwriteN)r      r      )g-P@grhA@g(mUg^ \^v@g+6CV@)r   r	      r	   )gw]@g_DP@gIϘ}X?gRxgE9)
A3copyb3	LU_decompL_solveU_solveroundA4b4
randmatrix)Abpyxibaks          x/Users/tonyclaw/.openclaw/workspace/skills/math-calculator/venv/lib/python3.14/site-packages/mpmath/tests/test_linalg.pytest_LU_decompr   H   s   
	A
AQ<DAaAA"#$!QE!RL!$ )K K K K
	A
AQ<DAaAA"#$!QE!RL!$ )K K K K1A
&&(Ca18O8 % %s   D<Dc                      \         \        \        3 FA  p \        V 4      p\	        W,          \        V P                  4      ,
          ^4      R8  d   KA  Q h	  R# )r   g+=N)A1A2A5inversemnormeyerows)r   invs     r   test_inverser'   ^   s?    "b\ajQUS[(!,v555     c                  b   ^\         n        \        \        r\	        \        W4      4      w  r#rEV\        \        R4      \        R4      \        R4      ^ .R\        R4      \        R4      \        R4      .R\        R4      \        R	4      \        R
4      .R\        R4      \        R4      \        R4      ..4      8X  g   Q hVRR\        R4      .8X  g   Q h\        \        V^4      ^
4      R8X  g   Q h. ROpR p. p. p	\        ^^4       Fx  p
V! V
4      p \	        \        W4      4      w  r#rE\        V4      p\        V4      pVP                  \        V^4      4       V	P                  \        \        WV4      ^4      4       Kz  	  V Uu. uF  p\        V^
4      NK  	  up. RO8X  g   Q h\        \        V4      \        V	4      ,
          \        4      R8  g   Q hR p. p. p\        ^^
4       F  p
V! V
4      p \	        V P                  4       4      w  r#rEVP                  \        V^4      4       VP                  \        \        V RRV
^,
          13,          W@RV
^,
          3,          4      ^4      4       K  	  \        \        V4      \        V4      ,
          \        4      R8  g   Q hR# u upi )   z3.0z-2.0z-1.0      ?z3.333333333333333z-2.9999999999999991z2.0z-0.66666666666666674z2.8142135623730948z-2.8284271247461898z-1.3333333333333333z-0.20000000000000018z4.2426406871192857z-1.4142135623730949gv@c           	          . p\        ^^4       FI  pTP                  \        ^V ^,           4       Uu. uF  pRW#,           ^,
          ,          NK  	  up4       KK  	  \        V4      # u upi )r   r+   )rangeappendmatrix)nr   r   js   &   r   coefftest_householder.<locals>.coefft   sZ    q"AHHaQ@1bAEAI&&@A ay As   A'
vIh%<=c                     \        ^V ,          V 4      p\        ^V ,          ^R^R7      pVP                  R 4      p\        VR,          4      V,          \        VRV 1^3,          4      ,          pV# )r   minmaxc                 D    \        R \        4       ,          V ,          4      # )y              ?)exppir   s   &r   <lambda>9test_householder.<locals>.hilbert_cmplx.<locals>.<lambda>   s    c"RT'!)nr(   N)NNNr   )hilbertr   applydiag)r0   r   vs   &  r   hilbert_cmplx'test_householder.<locals>.hilbert_cmplx   s_    AaCNqsA21-GG,-3LN4"1"Q$=(r(   r@   Ng      )gJ+Y@gy&1,M@g%C;B@g(\O8@g1Z1@gn(@gMb"@g(\@gq=
ףp@gʡE@gy&1,@g/$@)gRmX.@g8hi\?g5Żڌ^?g?g A8C?)mpdpsA8b8householderextendr/   mpfr   normr-   r.   residualinfr   )r   r   Hr   r   rr   r2   	residualsrefresr0   resrE   residuals_cmplxrefres_cmplxs                  r   test_householderrY   c   s{   BFrqVA\*JA!	%j#f+s6{A.
3"#C(=$>s5zJ
C&',@(A	
 "	3$%c*@&A	
!	"# # # # R234444aR $6666	A IF1a[!H .
a1I1Ia$d8A!,a01  '00isE#rNi0 5# # # # y!F6N2C86AAA OL1b\! *
atAqz*D!AdqsdF)Q!AaC%!A1EF	 
 '&*>>DvMMM) 1s   	J,c                      \        ^4      p \        V 4      w  rp\        W,          W#,          ,
          ^4      R8  g   Q hR# )   gV瞯<N)r   lur#   )r   PLUs       r   test_factorizationr`      s5    1AeGA!qsA'''r(   c            	         \        \        \        \        \        \        4      \        4      \
        4      R 8  g   Q h\        \        \        \        \        \        4      \        4      \
        4      R8  g   Q h\        \        \        \        \        \        4      \        4      \
        4      R8:  g   Q h\        \        \        \        \        \        4      ^ ,          \        4      \
        4      R 8  g   Q h\        \        \        \        \        \        4      ^ ,          \        4      \
        4      R8  g   Q h\        \        \        \        \        \        4      ^ ,          \        4      ^4      R8:  g   Q h\        \        \        \        \        \        4      \        4      ^4      R 8  g   Q h\        \        \        \        \        \        4      ^ ,          \        4      ^4      R 8  g   Q hR# )g|=      ?g333333@Ngo   @)rO   rP   A6lu_solveb6rQ   A7b7rJ   rK   qr_solveA10b10 r(   r   
test_solverl      s>   Xb"-r2C86AAAXb"-r2C83>>>Xb"-r2C8JFFFXb"-a0"5s;fDDDXb"-a0"5s;cAAAXb"-a0"5q9S@@@hsC0#6:VCCChsC03S91=FFFr(   c            	          \        ^R.^R.^^..4      p \        ^\        ,           ^\        ) .4      p\        \        V \	        W4      V4      4      R8  g   Q hR# )r   y               @y              @g|a2U?N)r/   r1   rO   rP   rd   )r   r   s     r   test_solve_overdet_complexrn      sU    B!R1a&)*AAq1"~AHQNA./&888r(   c                  
  a ^\         n        RR.R
R..p \        \        ^4      4      p^^.pRV : RV: R2RV: RV: R2RV : RV: R2RV: RV: R23 F(  o\        P
                  ! \        \        3V3R l4       K*  	  R	# )r*   ffffff@g333333?g?z	lu_solve(z, )z	qr_solve(c                     < \        S 4      # N)eval)r   s   r   r=   test_singular.<locals>.<lambda>   s	    tAwr(   Ng?)rH   rI   reprzerospytestraisesZeroDivisionError
ValueError)r   Br   r   s      @r   test_singularr}      sr    BF
seR[!AU1XA	
AA$%q)A+F$%q)A+FH(*5GHr(   c                  p   \         P                  \         P                  \        4      4      \         P                  . RO. RO. RO.4      8X  g   Q h\         P	                  \        \
        4      p \         P                  \         P                  \        V \
        4      \         P                  4      ^ 8X  g   Q hR# )r   N)r   r   r   )r   r   r   )r?   g      rb   )	fpcholeskyr/   A9cholesky_solveb9rO   rP   rQ   r<   s    r   test_choleskyr      sm    ;;ryy}%Iy/3Z)[[[[
"b!A772;;r1b)2662a777r(   c                     \        \        4      ^8X  g   Q h\        \        \        4      ^4      ^8X  g   Q h\        \        \        4      4      R8X  g   Q h\        \        \
        4      4      R8X  g   Q h\        \        4      ^8X  g   Q h\        \        \        4      4      R8X  g   Q h\        \        ^4      4      ^ 8X  g   Q hR# )r   i*  iC iN)	detr   r   r    r   r   r!   rc   rw   rk   r(   r   test_detr      s    r7a<<R""""R>T!!!R>W$$$r7a<<R>X%%%uQx=Ar(   c                      ^\         n        \        RR.RR..4      p \        V R 4      \	        R4      8X  g   Q h\        V R 4      \	        R4      8X  g   Q h\        V R 4      \	        R	4      8X  g   Q hR
# )r*   g6?gꕲq?ga2U0*?gBiq?c                     \        V ^4      # )r   r#   r<   s   &r   r=   test_cond.<locals>.<lambda>   s
    U1QZr(   z327065209.73817754c                 "    \        V \        4      # rs   )r#   rQ   r<   s   &r   r=   r      s
    U1S\r(   c                     \        V R 4      # )Fr   r<   s   &r   r=   r      s
    U1S\r(   z249729266.80008656N)rH   rI   r/   condrN   r   s    r   	test_condr      sv    BF 66"234A'(C0D,EEEE)*c2F.GGGG)*c2F.GGGGr(   c                  v    \        ^
^
4      p \        \        \        V 4      4      V ,
          ^4      R8  g   Q hR# )
   gj7֖6N)r   r#   r"   r   s    r   test_precisionr      s1    2rA$q(!,v555r(   c                     ^\         n        ^\        n        \        P                  . RO. RO. RO.4      p \        P                  . RO4      p\        P	                  W4      pV^ ,          P
                  R8  g   Q hV^,          P
                  R8  g   Q hV^,          P
                  R8  g   Q hRV^ ,          9   g   Q hR	V^,          9   g   Q hRV^,          9   g   Q hR# )
r*   r4   giE'n@g'ɮ@N)z0.1z0.3z1.0)z7.1z5.5z4.8)z3.2z4.4z5.6)4z0.6z0.5gFxlbO*)rH   rI   ivr/   rd   delta)ar   cs      r   test_interval_matrixr      s    BFBF
		&':;NOPA
		#$A
AAQ4::Q4::Q4::$!,,,$!,,,$!,,,r(   c                      \        ^4      p \        V 4      pV P                  \        V 4      8X  g   Q hRV R&   V P                  e   Q hR# )r	   N)r   r   )r   r   _LU)r   LUs     r   test_LU_cacher      sA    1A	1B55IaL   AcF55==r(   c                      \        ^RRR7      p \        ^^RRR7      p\        W4      \        ^^RRR7      ,           p\        WV4      p\        \	        WV4      ^4      \        \	        WV4      ^4      8  g   Q hR# )	r[   g#B;g@xDr6   i  gh㈵>Nr   gh㈵)r   rd   improve_solutionrO   rP   )r   r   x1x2s       r   test_improve_solutionr      so    1%T*A1aU-A	!*Qu%@	@B	!	#B"A&hqa.@!)DDDDr(   c            	         \        ^4       EF!  p ^p^pW,           \        n        ^ p^p\        ^V^,           4      p\        V4      p\        V Uu. uF  p\	        V4      NK  	  up4      p\        V4      R8  d   \        V4      p	\        V	4      pK(  X	V,          p	V	R,          V,          V	,          p
V	R,          V,          V	,          pV\        n        \        V
RR7      pW,           \        n        W,
          pV\        n        \        V\        4      P                  ^ 4      '       d   EK"  Q h	  ^\        n        R# u upi )r	   g{Gz?pade)methodNr?   )r-   rH   rI   rC   r:   absr   r   expmrO   rQ   ae)r   rI   extradmNdgr   r   expama1e2e1ds                 r   test_exp_pader      s    1X1QqS\HR(RSVR() "gn1AQBbDUQY]UT\A"V$GAs|q!!!!1 2 BF# )s   D?
c            	         ^\         n        R
p ^	p^pRp\        R4      p\        ^ ^
4       EF  p^ pRp\	        V^,          4      p\        V\        4       ,          4      p\        V\        4       ,          4      p	\        \        W4      4      p
\        \        W4      4      p\         P                  W4      pV'       d   RpRp\        ^ V4       Ft  p\        ^ V
4       Fa  p\        WV ,
          \        4       ,          ,           4      p\        WV ,
          \        4       ,          ,           4      p\        VV4      WV3&   Kc  	  Kv  	  M`RpRp\        ^ V4       FL  p\        ^ V
4       F9  p\        WV ,
          \        4       ,          ,           4      p\        V4      WV3&   K;  	  KN  	  \        WVR7      w  pp\        R4      p\        VVV,          ,
          4      pVV8:  g   Q hVR8X  dk   \        \        V
4      VP                  V,          ,
          4      pVV8:  g   Q h\        \        V
4      VVP                  ,          ,
          4      pVV8:  g   Q hVR8X  g   EK4  \        \        V
4      VP                  VP!                  4       ,          ,
          4      pVV8:  g   Q h\        \        V
4      VP!                  4       VP                  ,          ,
          4      pVV8:  d   EK  Q h	  R	# )r*   Fz0.0fullcomplexTreal)edpsz1.0E-11Ni)rH   rI   rN   xrangeboolnintrandintr8   r7   r/   mpcqrrO   r$   T	conjugate)lowlimituplimitmaxmflgzerokexdpsmodenum1num2r   r0   r   dtyper1   r   valval2QRmaxnormn1s                         r   test_qrr     sT   BFHGD
Cu:DAb\1q5k DK DK D D  IIaN CEAa[!Ax8+;TV*CCDCH,<df+D DED d^AcF % ! CEAa[!Ax8+;TV*CCDC XAcF % ! !%(1 i.!a!e)_W}}F?c!fqssQw&'B= =c!fq133w&'B= =Ic!fqssQ[[]223B= =c!fq{{}qss223B= =y r(   )r	   r      )r   r   r	   )r   r   r   )r      r
   )r   r?   r?   r   rG   )r   rG   r	   r?   )r   r	   rG   )r   r   r
   )r	   r   rG   r?   )r   r   r?   r?   r   )r   r   r   r   r?   )r
   r   r   r   )r   r   rG   	   i)r   r[   r   r      )r   r   r   r   2   )gQx$@g=
ףp=g        gQ롿gGz@)gNbX9gRQ?g(\@g(\B!gHzG?)gRQgQ?g rпgp=
ף"@gE)g)\#@g(\@g     @1g;On?g\(8@)gQ"g33333sA@gp=
ףS@gQ?g33333sP@)gfffff!@g
ףp=4@gGz@rp   gףp=
7M@)r   r   r   )rG   r   r[   )r	   r[   i)g	?gZB>@gr3܀o@)gsFZ*o@g#GGge,i9@)g}.AgR#@g(~Ƅ.)g     7@gQ/gQ@g      )r	   r   r   )r   r   r	   )r?   r   r   )r?   rG   r?   )r   r   r?   )r   r   r	   r
   )r
   r   rG   )r   r[   r   )rG   r   g      @)r      g      /)      ?      ?       @r   )g      @g      @g      @)g      @g       @g      "@)r+   r   r+   )2
__future__r   rx   mpmathlibmpbackendr   rH   r   r   r   rL   r   r/   r   b1r    b2r   r   r   r   r!   rc   re   rf   rg   rJ   rK   r   r   ri   rj   r   r'   rY   r`   rl   rn   r}   r   r   r   extradpsr   r   r   r   r   r   rk   r(   r   <module>r      s      			 LL	
**
**nn&& Y     !!!!!	# $
 33333	5 6
 &\  //9; < 'aYV!W  Y  \  $  ,6
5Nn(
G9
H8
H 
"6 6-E8D!r(   