+
    ‰°³jJ"  ã                   óÂ   € ^ 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W‰8  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
   &&        Úb/home/luca/screen-checking/.venv/lib/python3.14/site-packages/mpmath/tests/test_eigen_symmetric.pyÚ	run_eigsyr    	   sñ   € ßÜÐ"¤C¨£FÔ+ä�8Š8�A‹;�D€AØ	ŒB�GŠG�A‹J�˜Ÿ™›Õ&€AØ	�€AØ	�K‰K‹MÕœBŸFšF 1§6¡6›NÕ*€AçÜÐ Ô"ÜÐ Ô#ä	�Š�!‹€BÜ	�Š�!‹€BçÜˆm˜R  q¨$Ô/Üˆm˜R  q¨$Ô/ä
�&Š&�#œŸšœrŸv™v›Õ&Ó
'€CàŒ8€Oˆ8ØŒ8€Oˆ8à€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W‰8  g   Q hV# r   )r
   r   r   Úeigher   Útranspose_conjr   r   r   r   r   r   r   s
   &&        r   Ú	run_eigher%   $   s÷   € ßÜÐ"¤C¨£FÔ+ä�8Š8�A‹;�D€AØ	ŒB�GŠG�A‹J�˜×)Ñ)Ó+Õ+€AØ	�€AØ	×ÑÓÕ¤§¢¨¯©£Õ/€AçÜÐ Ô"ÜÐ Ô#ä	�Š�!‹€BÜ	�Š�!‹€BçÜˆm˜R  q¨$Ô/Üˆm˜R  q¨$Ô/ä
�&Š&�#œŸšœrŸv™v›Õ&Ó
'€CàŒ8€Oˆ8ØŒ8€Oˆ8à€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,          WšV
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       WÅ8  g   Q hWˆP                  4       ,          \        P                  ! VP                   4      ,
          p\        P                  ! V4      pV'       d   \        R	\        V4      RV4       WÅ8  g   Q hVP                  4       V,          \        P                  ! VP                  4      ,
          p\        P                  ! V4      pV'       d   \        R
\        V4      RV4       WÅ8  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˜  € à�6‰6�1—6‘6€qä
�&Š&�”r—v’vœbŸf™f“~Õ%Ó
&€CçÜÐ"¤C¨£FÔ+Üˆf�mÔ$ä�xŠx˜Ô9�H€Aˆ1ä
�Š�—‘˜Ÿ™Ó €AÜ”C˜“IÖˆØ•ˆˆAˆ#‹ñ ÷ Üˆf”c˜!“fÔÜˆf”c˜"“gÔÜˆf”c˜!“fÔà	���	�A�€AÜ
�(Š(�1‹+€CßÜˆe”S˜“V˜T 3Ô'ØŒ9Ðˆ9à	�K‰K‹MÕœBŸFšF 1§6¡6›NÕ*€AÜ
�(Š(�1‹+€CßÜˆf”c˜!“f˜d CÔ(ØŒ9Ðˆ9à	�‰‹˜ÕœBŸFšF 1§6¡6›NÕ*€AÜ
�(Š(�1‹+€CßÜˆf”c˜!“f˜d CÔ(ØŒ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,          WšV
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       WÅ8  g   Q hWˆP                  4       ,          \        P                  ! VP                   4      ,
          p\        P                  ! V4      pV'       d   \        R	\        V4      RV4       WÅ8  g   Q hVP                  4       V,          \        P                  ! VP                  4      ,
          p\        P                  ! V4      pV'       d   \        R
\        V4      RV4       WÅ8  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ž  € à�6‰6�1—6‘6€qä
�&Š&�”r—v’vœbŸf™f“~Õ%Ó
&€CçÜÐ"¤C¨£FÔ+Üˆf�mÔ$ä�xŠx˜Ô9�H€Aˆ1ä
�Š�—‘˜Ÿ™Ó €AÜ”C˜“IÖˆØ•ˆˆAˆ#‹ñ ÷ Üˆf”c˜!“fÔÜˆf”c˜"“gÔÜˆf”c˜!“fÔà	���	�A�€AÜ
�(Š(�1‹+€CßÜˆe”S˜“V˜T 3Ô'ØŒ:Ðˆ:à	×ÑÓÕ¤§¢¨¯©£Õ/€AÜ
�(Š(�1‹+€CßÜˆf”c˜!“f˜d CÔ(ØŒ9Ðˆ9à	×ÑÓ˜QÕ¤§¢¨¯©£Õ/€AÜ
�(Š(�1‹+€CßÜˆf”c˜!“f˜d CÔ(ØŒ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ä×Òœs 1›v uÓ-�D€AØŒ�Š�1‹Õ€AØŒ�Š�1‹Õ€Aä�8Š8�A‹;˜ÔÐÐÜ�8Š8�A‹;˜ÔÐÒ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Ü�AŽYˆÜ˜–ˆAÜ˜œRŸWšW›Y�¨Õ*¨eÕ3Ó4ˆA�ˆc‹Fó ñ ð €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      &@      1ÀrR   y      *À      7@)y      @      3Ày      *À      7ÀrT   )r   r@   r    r%   )r   s    r   Útest_eighe_fixed_matrixrV   £   sù   € Ü
�	Š	�A�q�6˜A˜q˜6Ð"Ó#€AÜˆa„LÜˆa„Lä
�	Š	�A�s�8˜c 2˜YÐ'Ó(€AÜˆa„LÜˆa„Lä
�	Š	’:š{ªJÐ7Ó8€AÜˆa„LÜˆa„Lä
�	Š	’9šiªÐ3Ó4€AÜˆa„LÜˆa„Lô 	�	Š	�A�t�9˜t Q˜iÐ(Ó)€AÜˆa„Lä
�	Š	’<¢²Ð>Ó?€AÜˆa„Lä
�	Š	Ò)Ò+CÒE\Ð]Ó^€AÜˆa†Lr!   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ä�BŽZˆØ”—’˜aÓ#Õ# aÕ'ˆä˜˜1–ˆAÜ˜A �E 1Ö%�Ø˜Q˜3���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ä�BŽZˆØ”—’˜qÓ$Õ$ qÕ(¨B°!´b·m²mÀAÓ6IÕ2IÈAÕ2MÕ,NÕNˆä˜˜1–ˆAÜ—U’U˜1˜q˜S�6“]ˆA�ˆc‰FÜ˜A �E 1Ö%�ÜŸš  Q 3¥›��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ä�BŽZˆÜ
�aÓ
ˆä˜˜1–ˆAÜ˜A �E 1Ö%�Ø˜Q˜3���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ä�BŽZˆÜ
�aÓ
˜b¤;¨qÓ#4Õ4Õ
4ˆä˜˜1–ˆAÜ—U’U˜1˜q˜S�6“]ˆA�ˆc‰FÜ˜A �E 1Ö%�ÜŸš  Q 3¥›��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¿   € Ü�A�YˆÜ�wŠw‹y˜3‰ˆØ””B—G’G“I •NÓ#Õ#ˆØ””B—G’G“I •NÓ#Õ#ˆØ”—’˜aÓ#Õ# aÕ'ˆÜ�7Š7‹9�sŒ?Ø��GˆAÜ˜A–Y�Ü ž�AÜ˜q 1 �v›;�A˜�c“Fó #ñ ô 	�!°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  € Ü�A�YˆÜ�wŠw‹y˜3‰ˆØ””B—G’G“I •NÓ#Õ#ˆØ””B—G’G“I •NÓ#Õ#ˆØ”—’˜qÓ$Õ$ qÕ(¨B°!´b·m²mÀAÓ6IÕ2IÈAÕ2MÕ,NÕNˆÜ�7Š7‹9�sŒ?Ø��GˆAÜ˜A–Y�Ü ž�AÜœrŸušu Q¨ s¥V›}Ó-°´S¼¿º¸qÀ1À½v»Ó5GÕ0GÕG�A˜�c“Fó #ñ ô 	�!°×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   iùÿÿÿrT   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Ô   € ô �&Š&�”r—v’vœbŸf™f“~Õ%Ó
&€Câ	Ú	Ú	Ú	Ú	Ú	Ú	Ú	ð	 €Aô 	�	Š	�!‹€AÜ
�	Š	”2—7’7˜4“= "¤b§g¢g¨c£l°A°qÐ9Ó:€Aä
�Š� Ô'€AØ…F€AÜ�8Š8�A‹;˜ÔÐÐä
�Š� Ô'€AØ…F€AÜ�8Š8�A‹;˜ÔÐÒ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# )gŸ“ÿE§yâ?ÚlegendreÚ
legendre01gß'Üež æ?gEç%”ø[ì?ÚhermiteÚlaguerreNgŸ“ÿE§yâ¿)g{¨[àlÿì¿g 0‹ý#;á¿r   g 0‹ý#;á?g{¨[àlÿì?)ççÎžSÎ?ç�ý$%Ú¡Þ?gžr3yV4â?r…   r„   )gþáúHÆ±?gíþàÕ?g
ö€õ�på?gÀÃ¼à6Çí?)ç>Šcò@CÆ?çá:Î†_ÞÔ?r‡   r†   gß'Üež æ¿)g"·ô¢U) Àgä¬€) ¬î¿r   gä¬€) ¬î?g"·ô¢U) @)ç…8ÿDŸn”?çRž1Ù?g�Á ø?î?r‰   rˆ   )gÐ\mæœÚ?gÜpP©¯Z@gÅñð\ç(@)g‰§Ë FÁæ?g²1m<ÓÑ?gšçoßöF…?)rF   )rB   rC   s     r   Útest_gauss_quadrature_staticrŠ   *  s�   € Ø	�zÐ"€AØ
ˆaˆ€AÜˆj˜!Ô âQ€AÚP€AÜˆj˜!Ô â<€AÚ<€AÜˆl˜AÔ"à	�zÐ"€AØ
�zÐ"€AÜˆi˜ÔâO€AÚO€AÜˆ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 ›e¡a¨£džlr!   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€   ù€ Ü×"Ò" 1 eÀ4ÔH‰ˆàˆÜœ˜A›–ˆAØ�1•™˜!�$›•ÕŠAñ  ô �GŠGÕ*¨AÓ.ˆä�GŠG�A•E‹NˆçÜ�%˜AÔ!à�4ŒxˆŠ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   € œRŸVšV Q B q¥Dœ\r!   rƒ   c                 ó0   € \         P                  ! V ) 4      # r•   r§   r£   s   &r   r˜   r¤   a  s   € œbŸfšf a Rœjr!   Ú	glaguerrec                 óf   € \         P                  ! V 4      \         P                  ! V ) 4      ,          # r•   )r   r}   r   r£   s   &r   r˜   r¤   b  s   € œrŸwšw q›z¬"¯&ª&°!°«*Ö4r!   )r’   Ú
chebyshev1c                 óV   € ^\         P                  ! ^W ,          ,
          4      ,          # r¢   ©r   r}   r£   s   &r   r˜   r¤   c  s   €  ¤"§'¢'¨!¨A­C­%£.Ö 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   € ˜Q˜q�S A¤b§f¢f¨Q£i¥KÕ0°A°aµC¸1¼R¿VºVÀA»Y½;Õ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˜!•e˜QÓ€Aõ÷ñ ñ  ˆ
‘K " a Ô)Ùˆ‘k A q 6Ô*Ùˆ	Ñ)¬R¯V©V¨G´R·V±VÐ+<Ô=Ùˆ
Ñ(¨1¬b¯f©f¨+Ô6ÙˆÑ4°q¼"¿&¹&°kÈ1ÌrÏvÊvÐVWËyÍ=ÕYÙˆÑ0°2°q°'Ô:ÙˆÑ.°°Q°Ô8ÙˆÑGÈ"ÈaÈÐZ[Ô^`×^dÒ^dÐefÓ^gÕZgÐpqÔtv×tzÒtzÐ{|Ó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!   