+
    iJ"                         ^ RI Ht ^ RI Ht ]P                  P                  tRR ltRR ltRR ltRR ltR t	RR lt
R	 tR
 tR tR tR tR tR tR tR tRR ltR# )    )mp)libmpc                    V'       d   \        R \        V 4      4       \        P                  ! V 4      w  r#V\        P                  ! V4      ,          VP                  4       ,          pW,
          pW3P                  4       ,          \        P                  ! V P                  4      ,
          pV'       d   \        RV4       \        RV4       \        P                  ! V4      p\        P                  ! V4      pV'       d   \        RVRVR4       \        RVRVR4       \        P                  ! R\        P                  ! \        P                  4      ,          4      p	Wy8  g   Q hW8  g   Q hV# original matrix:
zeigenvalues:
zeigenvectors:
zdifference:
皙?)printstrr   eigsydiag	transposeeyerowsmnormexplogeps
AverboseDQBCENCNEr   s
   &&        _/var/www/html/photoedit/myenv/lib/python3.14/site-packages/mpmath/tests/test_eigen_symmetric.py	run_eigsyr    	   s    "CF+88A;DA	BGGAJ&A	A	KKMBFF166N*A"#	!B	!BmRq$/mRq$/
&&#rvv&
'C8O88O8I    c                    V'       d   \        R \        V 4      4       \        P                  ! V 4      w  r#V\        P                  ! V4      ,          VP                  4       ,          pW,
          pW3P                  4       ,          \        P                  ! V P                  4      ,
          pV'       d   \        RV4       \        RV4       \        P                  ! V4      p\        P                  ! V4      pV'       d   \        RVRVR4       \        RVRVR4       \        P                  ! R\        P                  ! \        P                  4      ,          4      p	Wy8  g   Q hW8  g   Q hV# r   )r
   r   r   eigher   transpose_conjr   r   r   r   r   r   r   s
   &&        r   	run_eigher%   $   s    "CF+88A;DA	BGGAJ))++A	A	/A"#	!B	!BmRq$/mRq$/
&&#rvv&
'C8O88O8Ir!   c                    V P                   V P                  rC\        P                  ! R \        P                  ! \        P
                  4      ,          4      pV'       d"   \        R\        V 4      4       \        RV4       \        P                  ! WR7      w  rgp\        P                  ! VP                  VP                   4      p	\        \        W44      4       F  p
Wz,          WV
3&   K  	  V'       d@   \        R\        V4      4       \        R\        V4      4       \        R\        V4      4       Wi,          V,          V ,
          p\        P                  ! V4      pV'       d   \        R\        V4      RV4       W8  g   Q hWP                  4       ,          \        P                  ! VP                   4      ,
          p\        P                  ! V4      pV'       d   \        R	\        V4      RV4       W8  g   Q hVP                  4       V,          \        P                  ! VP                  4      ,
          p\        P                  ! V4      pV'       d   \        R
\        V4      RV4       W8  g   Q hR# r	   r   full)full_matriceszU:
zS:
zV:
zC
r   zD:
zE:
N)r   colsr   r   r   r   r
   r   svd_rzerosxrangeminr   r   r   r   r)   r   mnr   US0VSjr   errr   r   s   &&&            r   	run_svd_rr8   ?   s   66166q
&&rvvbff~%
&C"CF+fm$xx9HA1
 ACIA#  fc!ffc"gfc!f		AA
((1+CeSVT3'99	KKMBFF166N*A
((1+Cfc!fdC(99	BFF166N*A
((1+Cfc!fdC(99r!   c                    V P                   V P                  rC\        P                  ! R \        P                  ! \        P
                  4      ,          4      pV'       d"   \        R\        V 4      4       \        RV4       \        P                  ! WR7      w  rgp\        P                  ! VP                  VP                   4      p	\        \        W44      4       F  p
Wz,          WV
3&   K  	  V'       d@   \        R\        V4      4       \        R\        V4      4       \        R\        V4      4       Wi,          V,          V ,
          p\        P                  ! V4      pV'       d   \        R\        V4      RV4       W8  g   Q hWP                  4       ,          \        P                  ! VP                   4      ,
          p\        P                  ! V4      pV'       d   \        R	\        V4      RV4       W8  g   Q hVP                  4       V,          \        P                  ! VP                  4      ,
          p\        P                  ! V4      pV'       d   \        R
\        V4      RV4       W8  g   Q hR# r'   )r   r*   r   r   r   r   r
   r   svd_cr,   r-   r.   r   r$   r   r/   s   &&&            r   	run_svd_cr;   f   s   66166q
&&rvvbff~%
&C"CF+fm$xx9HA1
 ACIA#  fc!ffc"gfc!f		AA
((1+CeSVT3'::	/A
((1+Cfc!fdC(99	Q/A
((1+Cfc!fdC(99r!   c                 6   R p\         P                  ! \        V4      V 4      w  rEV\         P                  ! V4      ,          pV\         P                  ! V4      ,          p\         P                  ! V4      V8  g   Q h\         P                  ! V4      V8  g   Q hR# )h㈵>N)r   gauss_quadraturelenmatrixr   )qtypeabr   des   &&&   r   	run_gaussrF      so    
Cs1vu-DA1A1A88A;88A;r!   c                    \         P                  ! W 4      p\        V 4       FL  p\        V 4       F:  p\        ^\         P                  ! 4       ,          ^,
          V,          4      W#V3&   K<  	  KN  	  V# )z$
random matrix with integer entries
)r   r@   r-   intrand)r1   ranger   ir6   s   &&   r   irandmatrixrL      sX     			!AAYARWWY*e34AcF   Hr!   c                     \         P                  ! ^^.^^..4      p \        V 4       \        V 4       \         P                  ! ^R.R^..4      p \        V 4       \        V 4       \         P                  ! . RO. RO. RO.4      p \        V 4       \        V 4       \         P                  ! . RO. RO. RO.4      p \        V 4       \        V 4       \         P                  ! ^R	.R
^..4      p \        V 4       \         P                  ! . RO. RO. RO.4      p \        V 4       \         P                  ! . RO. RO. RO.4      p \        V 4       R# )   N)rN         )rP         )rQ   rS      )rN   r   rQ   )r   rR      )rQ   rU   rT   y      @      @y      @      )rN   y             &r   )y              &@rR   y              =@)r   y             =rT   )rN   y      &@      1@y      @      3@)y      &@      1rR   y      *      7@)y      @      3y      *      7rT   )r   r@   r    r%   )r   s    r   test_eighe_fixed_matrixrV      s    
		Aq6Aq6"#AaLaL
		As8c2Y'(AaLaL
		:{J78AaLaL
		9i34AaLaL 			At9tQi()AaL
		<>?AaL
		)+CE\]^AaLr!   c                     ^p \        ^
4       Fm  p^\        P                  ! W 4      ,          ^,
          p\        ^ V 4       F+  p\        V^,           V 4       F  pW#V3,          W$V3&   K  	  K-  	  \        V4       Ko  	  R# )rT   N)r-   r   
randmatrixr    NrB   r   rK   r6   s        r   test_eigsy_randmatrixr[      sl    	ABZa##a'1AAE1%Q3A# &  	! r!   c                     ^p \        ^
4       F  p^\        P                  ! W 4      ,          ^,
          R^\        P                  ! W 4      ,          ^,
          ,          ,           p\        ^ V 4       Fa  p\        P                  ! W#V3,          4      W#V3&   \        V^,           V 4       F%  p\        P                  ! W#V3,          4      W$V3&   K'  	  Kc  	  \        V4       K  	  R# )rT                 ?N)r-   r   rX   reconjr%   rY   s        r   test_eighe_randmatrixr`      s    	ABZq$$q(B!bmmA6I2IA2M,NN1AUU1qS6]AcFAE1%Q3A# & 
 	! r!   c                      ^p ^p\        ^
4       FT  p\        W4      p\        ^ V 4       F+  p\        V^,           V 4       F  pW4V3,          W5V3&   K  	  K-  	  \        V4       KV  	  R# )   N)r-   rL   r    rZ   RrB   r   rK   r6   s         r   test_eigsy_irandmatrixre      sc    	A	ABZ
a
1AAE1%Q3A# &  	! r!   c                  n   ^p ^p\        ^
4       F  p\        W4      R\        W4      ,          ,           p\        ^ V 4       Fa  p\        P                  ! W4V3,          4      W4V3&   \        V^,           V 4       F%  p\        P                  ! W4V3,          4      W5V3&   K'  	  Kc  	  \        V4       K  	  R# )rb   r]   N)r-   rL   r   r^   r_   r%   rc   s         r   test_eighe_irandmatrixrg      s    	A	ABZ
a
b;q#44
41AUU1qS6]AcFAE1%Q3A# & 
 	! r!   c                  ,   \        ^4       EF  p \        P                  ! 4       R8  p^\        \        P                  ! 4       ^
,          4      ,           p^\        \        P                  ! 4       ^
,          4      ,           p^\        P                  ! W#4      ,          ^,
          p\        P                  ! 4       R8  dE   V^
,          p\        V4       F,  p\        V4       F  p\        WEV3,          4      WEV3&   K  	  K.  	  \        WARR7       EK  	  R# )rT         ?Fr)   r   N)r-   r   rI   rH   rX   r8   rK   r(   r0   r1   r   xys          r   test_svd_r_randrn      s    AYwwy3BGGIN##BGGIN##a##a'779s?GAAYAq1v;AcF #  	!U;; r!   c                     \        ^4       EFz  p \        P                  ! 4       R8  p^\        \        P                  ! 4       ^
,          4      ,           p^\        \        P                  ! 4       ^
,          4      ,           p^\        P                  ! W#4      ,          ^,
          R^\        P                  ! W#4      ,          ^,
          ,          ,           p\        P                  ! 4       R8  d   V^
,          p\        V4       Fs  p\        V4       Fa  p\        \        P
                  ! WEV3,          4      4      R\        \        P                  ! WEV3,          4      4      ,          ,           WEV3&   Kc  	  Ku  	  \        WARR7       EK}  	  R# )rT   ri   r]   Frj   N)r-   r   rI   rH   rX   r^   imr;   rk   s          r   test_svd_c_randrq     s   AYwwy3BGGIN##BGGIN##q$$q(B!bmmA6I2IA2M,NN779s?GAAYAruuQsV}-Sq1v5G0GGAcF #  	!77 r!   c            	      ^   \         P                  ! R \         P                  ! \         P                  4      ,          4      p . RO. RO. RO. R	O. R
O. RO. RO. RO.p\         P                  ! V4      p\         P                  ! \         P
                  ! R4      ^\         P
                  ! R4      ^ ^ .4      p\         P                  ! VRR7      pW2,          p\         P                  ! V4      V 8  g   Q h\         P                  ! VRR7      pW2,          p\         P                  ! V4      V 8  g   Q hR# )r	   i  i  F)
compute_uvN)   
   rN   rR   rQ   )   rQ   ru   r      )rS   rx   rO   rR   )rS   rz   rb   )	   rw   rU   rz   rb   )r{   rU   irT   rx   )rN   i   rT   rU   )rb   rT   r   rz   rN   )	r   r   r   r   r@   sqrtr+   r   r:   )r   rB   rC   r5   s       r   test_svd_test_caser~     s     &&rvvbff~%
&C									 A 			!A
		2774="bggclAq9:A
'AFA88A;
'AFA88A;r!   c                      RR .p ^^.p\        RW4       . R	Op . R
Op\        RW4       . ROp . ROp\        RW4       RR.p RR.p\        RW4       . ROp . ROp\        RW4       . ROp . ROp\        RW4       R# )gEy?legendre
legendre01g'e?gE%[?hermitelaguerreNgEy)g{[lg0#;r   g0#;?g{[l?)ΞS?$%ڡ?gr3yV4?r   r   )gHƱ?g?g
p?gü6?)>c@C?:Ά_?r   r   g'e)g"U) g䬀)r   g䬀)?g"U) @)8Dn?R1?g??r   r   )g\m?gpPZ@g\(@)g F?g1m<?goF?)rF   )rB   rC   s     r   test_gauss_quadrature_staticr   *  s    	z"A
aAj! QAPAj! <A<AlA"	z"A
z"AiOAOAi/A/Aj! r!   c           
      \  a aaa ^o\         P                  ! ^S,          ^4      oV3R loRVVV 3R llpV! RR R^.4       V! RR ^ ^.4       V! RR \         P                  ) \         P                  .4       V! R	R
 ^ \         P                  .4       V! RR ^ \         P                  .^\         P                  ! ^4      ,          R7       V! RR R^.4       V! RR R^.4       V! RR R^.^\         P                  ! ^4      ,          ^\         P                  ! ^4      ,          R7       R# )rT   c                    < ^ p\        \        S4      ^,
          RR4       F  pW,          SV,          ,           pK  	  V# )r   rx   )r-   r?   )rl   rrK   r   s   &  r   F(test_gauss_quadrature_dynamic.<locals>.FH  s7    A
B+A!A ,r!   c                 h  <a \         P                  ! SWVR 7      w  rV^ p\        \        V4      4       F%  pWvV,          S! WX,          4      ,          ,          pK'  	  \         P                  ! VV3R lV4      p	\         P
                  ! Wy,
          4      p
S'       d   \        W
Wy4       V
R8  g   Q hR# )alphabetac                 .   < S! V 4      S! V 4      ,          # N )rl   r   FWs   &r   <lambda><test_gauss_quadrature_dynamic.<locals>.run.<locals>.<lambda>U  s    beadlr!   r=   N)r   r>   r-   r?   quadfabsr
   )rA   r   rd   r   r   XWrB   rK   rC   cr   r1   r   s   &f&&&      r   run*test_gauss_quadrature_dynamic.<locals>.runN  s    ""1e4HAA1!$A   GG*A.GGAEN%A!4xxr!   r   c                     ^# rU   r   rl   s   &r   r   /test_gauss_quadrature_dynamic.<locals>.<lambda>^  s    ar!   r   c                     ^# r   r   r   s   &r   r   r   _  s    r!   r   c                 >    \         P                  ! V ) V ,          4      # r   r   r   r   s   &r   r   r   `  s    RVVQBqD\r!   r   c                 0    \         P                  ! V ) 4      # r   r   r   s   &r   r   r   a  s    bffaRjr!   	glaguerrec                 f    \         P                  ! V 4      \         P                  ! V ) 4      ,          # r   )r   r}   r   r   s   &r   r   r   b  s    rwwqz"&&!*4r!   )r   
chebyshev1c                 V    ^\         P                  ! ^W ,          ,
          4      ,          # r   r   r}   r   s   &r   r   r   c  s    "''!AC%. 0r!   
chebyshev2c                 H    \         P                  ! ^W ,          ,
          4      # r   r   r   s   &r   r   r   d  s    !#r!   jacobic                     ^V ,
          ^\         P                  ! ^4      ,          ,          ^V ,           ^\         P                  ! ^4      ,          ,          ,          # r   )r   mpfr   s   &r   r   r   e  s3    QqSAbffQiK0AaC1RVVAY;3GGr!   r   N)r   r   rx   )r   rX   infr   )r   r   r   r   r1   s   f @@@r   test_gauss_quadrature_dynamicr   C  s    	A
a!eQA   
K"a)kAq6*	)RVVGRVV+<=
(1bff+64q"&&k1rvvVWy=Y02q':.Q8G"aZ[^`^d^def^gZgpqtvtztz{|t}p}r!   N)F)FT)ru   )mpmathr   r   backendr-   r    r%   r8   r;   rF   rL   rV   r[   r`   re   rg   rn   rq   r~   r   r   r   r!   r   <module>r      so     			66%N%N8
<86!2"@r!   