+
    iA                     t    ^ RI Ht ^ RI Ht ]P                  P                  tR	R ltR	R ltR	R ltR tR t	R# )
    )mp)libmpc           	      $   V^8  d   \        RV 4       V P                  p\        P                  ! V 4      w  r4V^8  d   \        RV4       \        RV4       W4,          VP	                  4       ,          p\        P
                  ! R\        P                  ! \        P                  4      ,          4      p^ p\        V4       F>  p\        V4       F,  p	V\        W	V3,          WYV3,          ,
          4      ,          pK.  	  K@  	  WrV,          ,          p^ p
\        V4       F7  p\        V^,           V4       F  p	V
\        WIV3,          4      ,          p
K  	  K9  	  V^ 8  d   \        RWz4       V^8  d   \        RV4       Wv8  g   Q hV
^ 8X  g   Q hR# )   zoriginal matrix (hessenberg):
Q:
zH:
皙?zdifference (H):B:
N)
printrowsr   
hessenbergtranspose_conjexplogepsxrangeabs)AverbosenQHBr   err0xyerr1s   &&         U/var/www/html/photoedit/myenv/lib/python3.14/site-packages/mpmath/tests/test_eigen.pyrun_hessenbergr   	   s>   {/3	A==DA{fQfQ	  ""A
&&rvvbff~%
&CDAYACA#Q3((D   	EMDDAYAq!ACA#KD "  {,{fa::199    c           	      8   V^8  d   \        RV 4       V P                  p\        P                  ! V 4      w  r4V^8  d   \        RV4       \        RV4       W4,          VP	                  4       ,          pW3P	                  4       ,          p\        P
                  ! R\        P                  ! \        P                  4      ,          4      p^ p\        V4       F>  p	\        V4       F,  p
V\        W
V	3,          WZV	3,          ,
          4      ,          pK.  	  K@  	  WV,          ,          p^ p\        V4       FK  p	\        V4       F9  p
W8X  d   WjV	3;;,          ^,          uu&   V\        WjV	3,          4      ,          pK;  	  KM  	  WV,          ,          p^ p\        V4       F7  p	\        V	^,           V4       F  p
V\        WJV	3,          4      ,          pK  	  K9  	  V^ 8  d   \        RWV4       V^8  d   \        RV4       W8  g   Q hW8  g   Q hV^ 8X  g   Q hR# )r   zoriginal matrix (schur):
r   zR:
r   zdifference (S):r	   N)
r
   r   r   schurr   r   r   r   r   r   )r   r   r   r   Rr   Cr   r   r   r   r   err2s   &&           r   	run_schurr%   .   s   {*A.	A88A;DA{fafa	  ""A	A
&&rvvbff~%
&CDAYACA#Q3((D   	EMDDAYAvA#!CA#KD  
 	EMDDAYAq!ACA#KD "  {T2{fa::::199r   c                    V^8  d   \        RV 4       V P                  p\        P                  ! V RRR7      w  r4pV^8  d%   \        RV4       \        RV4       \        RV4       \        P                  ! R\        P
                  ! \        P                  4      ,          4      p^ p\        V4       F  pWRV3,          ,          W8,          VRV3,          ,          ,
          p	\        V\        P                  ! V	4      4      pWHR3,          V ,          WHR3,          W8,          ,          ,
          p	\        V\        P                  ! V	4      4      pK  	  WrV,          ,          pV^ 8  d   \        R	V4       Wv8  g   Q hR
# )r   zoriginal matrix (eig):
T)leftrightzE:
zEL:
zER:
r   :NNNzdifference (E):N)
r
   r   r   eigr   r   r   r   maxmnorm)
r   r   r   EELERr   r   ir   s
   &&        r   run_eigr0   \   s   {(!,	Aqt4IA2{fagrgr
&&rvvbff~%
&CDAY1Q3K!$AaC.(4!%sGaK"qS'AD.(4!%  	EMD{&::r   c                     ^ p \        ^4       EF  p^\        \        P                  ! 4       ^,          4      ,           p\        P                  ! 4       R8  d   ^\        P                  ! W"4      ,          ^,
          p\        P                  ! 4       R8  dE   V^
,          p\        V4       F,  p\        V4       F  p\        W4V3,          4      W4V3&   K  	  K.  	  M^\        P                  ! W"4      ,          ^,
          R^\        P                  ! W"4      ,          ^,
          ,          ,           p\        P                  ! 4       R8  d   V^
,          p\        V4       Fs  p\        V4       Fa  p\        \        P
                  ! W4V3,          4      4      R\        \        P                  ! W4V3,          4      4      ,          ,           W4V3&   Kc  	  Ku  	  \        W0R7       \        W0R7       \        W0R7       EK  	  R# )r   g      ?y              ?r   N)
r   intr   rand
randmatrixreimr   r%   r0   )vr/   r   r   r   r   s         r   test_eig_dynr9   |   s]   	AAYBGGIM""779s?BMM!''!+Awwy3RA#AY!$QsVA# ' # R]]1((1,a"--:M6MPQ6Q0RRAwwy3RA#AY!$RUU1qS6]!3b3ruuQsV};M6M!MA# ' # 	q&!!) r   c                     ^ p . p\         P                  ! . RO. RO. RO.4      pVP                  V4       VP                  VP                  4       4       \         P                  ! . RO. RO. RO.4      pVP                  V4       VP                  VP                  4       4       \         P                  ! . RO. RO. RO.4      pVP                  V4       VP                  VP                  4       4       \         P                  ! . R	O. R
O. RO.4      pVP                  V4       VP                  VP                  4       4       V F'  p\	        W R7       \        W R7       \        W R7       K)  	  R# )r   r2   N)   r   r   )r   r;   r   )r   r   r;   )r;   r   r   )r;   r   r   )r   r;   r   )r   r   r   )r   r   r   )r   r   r   )r   matrixappend	transposer   r%   r0   )r8   ASr   s      r   test_eigr@      s   	A	B
		9 	A IIaLIIakkm
		9 	A IIaLIIakkm
		9 	A IIaLIIakkm
		9 	A IIaLIIakkmq&!! r   N)r   )
mpmathr   r   backendr   r   r%   r0   r9   r@    r   r   <module>rD      s8     			"J,\@ 0 r   