Ë
    0Ò^jª  ã                   ó¨  — d dl Zd dlmZmZmZ d dl mZmZmZm	Z	m
Z
mZmZmZmZmZmZmZmZ d dlmZ d dlmZ d dlZd„ Zd„ Zd„ Zd	„ Zej8                  j;                  d
eeg«      ej8                  j;                  dddg«      d„ «       «       Zej8                  j;                  d
eeg«      ej8                  j;                  dddg«      d„ «       «       Zy)é    N)Úassert_array_almost_equalÚassert_allcloseÚassert_)ÚarrayÚeyeÚzerosÚ
empty_likeÚemptyÚtril_indices_fromÚtrilÚtriu_indices_fromÚspacingÚfloat32Úfloat64Ú	complex64Ú
complex128)ÚComplexWarning)Úldlc                  ó  — t        d«      } t        j                  t        «      5  t	        | d d …d d…f   «       d d d «       t        j
                  t        «      5  t	        | dz  «       d d d «       y # 1 sw Y   Œ:xY w# 1 sw Y   y xY w)Né   é   ù              ð?)r   ÚpytestÚraisesÚ
ValueErrorr   Úwarnsr   )ÚAs    úm/var/www/ramen.bs-engineer-server.com/venv/lib/python3.12/site-packages/scipy/linalg/tests/test_decomp_ldl.pyÚ	test_argsr      sr   € ÜˆA‹€Aä	�‰”zÓ	"ñ ÜˆAŠa��!�ˆe‰HŒ÷ô 
�‰”nÓ	%ñ ÜˆAˆb‰DŒ	÷ð ÷ð ú÷ð ús   ¥A3ÁA?Á3A<Á?Bc                  óà   — t        dt        ¬«      } t        t        d«      «      \  }}}t        |t	        | «      «       t        |t	        | «      «       t        |t        g t        ¬«      «       y )N)r   r   ©Údtype)r
   Úcomplexr   r   r	   r   Úint)ÚaÚlÚdÚps       r   Útest_empty_arrayr)      sO   € ÜˆfœGÔ$€AÜ”%˜“-Ó �G€A€qˆ!Ü˜a¤¨A£Ô/Ü˜a¤¨A£Ô/Ü˜a¤ r´Ô!5Õ6ó    c                  ó¨  — t        g d¢g d¢g d¢g d¢g«      } t        g d¢g d¢g d¢g d¢g«      }t        g d	¢g d
¢g d¢g d¢g d¢g«      }t        g d¢g d¢g d¢g d¢g d¢g«      }t        g d¢g d¢g d¢g d¢g«      }|||fD ]¢  }t        |«      \  }}}t        |j                  |«      j                  |j                  «      |t        d«      d¬«       t        |d¬«      \  }}}t        |j                  |«      j                  |j                  «      |t        d«      d¬«       Œ¤ t        | d¬«      \  }}}t        |j                  |«      j                  |j                  «      | t        d«      d¬«       t        | dd¬«      \  }}}t        |j                  |«      j                  |j                  «      | t        d«      d¬«       t        |j                  «       j                  d¬«      \  }}}t        t        |j                  |«      j                  |j                  «       j                  «      |z
  «      t        d«      t        d«      d¬«       y )N)yö(\�ÂõØ¿¸…ëQ¸æ¿ù�Âõ(\�@{®Gázä¿ùq=
×£pÀ®Gáz®Àùffffff@q=
×£pí?)r,   y¸…ëQ¸!@ö(\�Âõü?ù)\�Âõ(À�Âõ(\�â?ùHáz®G@q=
×£pù¿)r-   r/   y¤p=
×£À¸…ëQ¸ž¿ù¤p=
×£ø¿áz®GáÀ)r.   r0   r1   yìQ¸…ëá¿¸…ëQ¸¾?)g      @é
   é   é   )g      $@r   é   r3   )g      ð?r5   é   é	   )g      2@r3   r7   r   )g      J@éa   ép   ék   é2   )g     @X@ér   éY   éb   é   )g      \@r=   é@   é!   é   )g     ÀZ@r>   rA   é<   éI   )g      I@r?   rB   rD   éM   )ç       @r   éüÿÿÿr   é   )rF   éþÿÿÿrI   r2   éøÿÿÿ)g      ÀrI   rB   rJ   rG   )g        r2   rJ   rB   éúÿÿÿ)g      @rJ   rG   rK   r2   )yÃõ(\�Âõ¿        ù                rL   rL   )yHáz®Gù?ÍÌÌÌÌÌì¿y=
×£p½!À        rL   rL   )y®Gáz®@áz®GáÊ?yq=
×£pý¿¸…ëQ¸ž?y…ëQ¸…À        rL   )yHáz®G@      ø¿y{®Gázü¿áz®Gáò¿y)\�Âõ(¼?)\�Âõ(¼¿yq=
×£pý¿        g     @�@r   )ÚatolÚrtolF©Úlower©Ú	hermitian©rP   rR   )rH   rH   )	r   r   r   ÚdotÚTr   Úconjr   r   )	r%   ÚbÚcr'   ÚeÚxr&   r(   Úus	            r   Útest_simpler\      sô  € ÜÒAÚ@ÚCÚAðCó 	D€Aô 	ŠÚÚÚðó 	 €Aô 	Ò&Ú%Ú$Ú%Ú#ð	%ó 	&€Aô 	Ò Ú#Ú#Ú"Ú#ð	%ó 	&€Aô
 	Ò.Ú1Ú8ÚAðCó 	D€Að ��AˆYò KˆÜ�a“&‰ˆˆ1ˆaÜ˜Ÿ™˜a›Ÿ™ Q§S¡SÓ)¨1´7¸5³>ÈÕJä�a˜uÔ%‰ˆˆ1ˆaÜ˜Ÿ™˜a›Ÿ™ Q§S¡SÓ)¨1´7¸5³>ÈÖJðKô �!˜uÔ%�G€A€qˆ!Ü�A—E‘E˜!“H—L‘L §¡Ó% q¬w°u«~ÀAÕFä�!˜5¨EÔ2�G€A€qˆ!Ü�A—E‘E˜!“H—L‘L §¡Ó% q¬w°u«~ÀAÕFô �!—&‘&“(—*‘* AÔ&�G€A€qˆ!Ü”D˜Ÿ™˜q›Ÿ™ a§f¡f£h§j¡jÓ1°!Ñ3Ó4´e¸F³mÜ  ›¨aö1r*   c            	      óˆ  — t         j                  j                  d«      } t        d«      D �]  }| j	                  dd«      }| j	                  d«      r| j                  ||f«      dz   n| j                  ||f«      dz  }||j                  «       j                  z   }|t        |«      | j	                  dd	«      z  z  }t        |d
¬«      }t        |d¬«      }t        |d¬«      \  }}}t        t        ||d d …f   |   «       d|› d�«       t        |d¬«      \  }	}}t        t        |	|d d …f   |   «       d|› d�«       �Œ y )NéÒ  r2   r3   éd   r   r   r   é   ç    €„.Aéÿÿÿÿ)ÚkrO   zSpin z failed)ÚnpÚrandomÚdefault_rngÚrangeÚintegersrV   rU   r   r   r   r   r   Úany)
ÚrngÚ_ÚnrZ   Úl_indÚu_indr[   r'   r(   r&   s
             r   Útest_permutationsro   H   s)  € Ü
�)‰)×
Ñ
 Ó
%€CÜ�2‹Yó =ˆØ�L‰L˜˜CÓ ˆà&)§l¡l°1¤oˆC�J‰J˜˜1�vÓ Ò"¸3¿:¹:ÀqÈ!ÀfÓ;MÈbÑ;PˆØ�—‘“—
‘
‰NˆØ	ŒS�‹V�C—L‘L  CÓ(Ñ(Ñ(ˆÜ! ! rÔ*ˆÜ! ! qÔ)ˆô �a˜q”/‰ˆˆ1ˆaä”C˜˜!šQ˜$™ ™Ó'Ð'¨5°°°7Ð);Ô<ä�a˜q”/‰ˆˆ1ˆaä”C˜˜!šQ˜$™ ™Ó'Ð'¨5°°°7Ð);Ö<ñ!=r*   r"   rl   é   é–   c                 ó:  — t         j                  j                  d«      }d| › d|› �}|j                  | | f«      j                  |«      }||j                  z   }|t        | |¬«       ||j                  dd«      «      z  z  }t        |«      \  }}}t        |d¬«      \  }}	}|t        u rd	nd
}
t        |j                  |«      j                  |j                  «      ||
|¬«       t        |j                  |	«      j                  |j                  «      ||
|¬«       y )Nr^   zFailed for size: ú	, dtype: r!   r`   ra   r   rO   g-Cëâ6?ç»½×Ùß|Û=©rN   Úerr_msg)rd   re   rf   ÚastyperU   r   rh   r   r   r   rT   )rl   r"   rj   ÚmsgrZ   r&   Úd1r(   r[   Úd2rN   s              r   Ú$test_ldl_type_size_combinations_realr{   ]   sî   € ô �)‰)×
Ñ
 Ó
%€CØ˜q˜c ¨5¨'Ð2€Cà�
‰
�A�q�6Ó×!Ñ! %Ó(€AØ	ˆA�C‰C‰€AØŒˆQ�eÔ	™U 3§<¡<°°3Ó#7Ó8Ñ	8Ñ8€Aä�1‹v�H€A€rˆ1Ü�1˜AŒ�H€A€rˆ1ØœGÑ#‰4¨€DÜ�A—E‘E˜"“I—M‘M !§#¡#Ó&¨°¸cÕBÜ�A—E‘E˜"“I—M‘M !§#¡#Ó&¨°¸cÖBr*   c                 óÐ  — t         j                  j                  d«      }d| › d|› �}d| › d|› �}|j                  | | f«      d|j                  | | f«      z  z   j                  |«      }||j	                  «       j
                  z   }|t        | |¬«       ||j                  dd«      «      z  z  }t        |«      \  }}}t        |d	¬
«      \  }	}
}|t        u rdnd}t        |j                  |«      j                  |j	                  «       j
                  «      |||¬«       t        |	j                  |
«      j                  |	j	                  «       j
                  «      |||¬«       |j                  | | f«      d|j                  | | f«      z  z   j                  |«      }||j
                  z   }|t        | |¬«       ||j                  dd«      «      z  z  }t        |d	¬«      \  }}}t        |d	d	¬«      \  }	}
}t        |j                  |«      j                  |j
                  «      |||¬«       t        |	j                  |
«      j                  |	j
                  «      |||¬«       y )Nr^   zHer failed for size: rs   zSym failed for size: r   r!   r`   ra   r   rO   g-Cëâ6*?rt   ru   rQ   rS   )rd   re   rf   rw   rV   rU   r   rh   r   r   r   rT   )rl   r"   rj   Úmsg1Úmsg2rZ   r&   ry   r(   r[   rz   rN   s               r   Ú'test_ldl_type_size_combinations_complexr   n   s   € ô �)‰)×
Ñ
 Ó
%€CØ# A 3 i°¨wÐ7€DØ# A 3 i°¨wÐ7€Dð 
�‰�Q˜�FÓ	˜B˜sŸz™z¨1¨a¨&Ó1Ñ1Ñ	1×9Ñ9¸%Ó@€AØ	ˆ!�&‰&‹(�*‰*‰€AØŒˆQ�eÔ	™U 3§<¡<°°3Ó#7Ó8Ñ	8Ñ8€Aä�1‹v�H€A€rˆ1Ü�1˜AŒ�H€A€rˆ1ØœIÑ%‰4¨5€DÜ�A—E‘E˜"“I—M‘M !§&¡&£(§*¡*Ó-¨q°tÀTÕJÜ�A—E‘E˜"“I—M‘M !§&¡&£(§*¡*Ó-¨q°tÀTÕJð 
�‰�Q˜�FÓ	˜B˜sŸz™z¨1¨a¨&Ó1Ñ1Ñ	1×9Ñ9¸%Ó@€AØ	ˆ!�#‰#‰€AØŒˆQ�eÔ	™U 3§<¡<°°3Ó#7Ó8Ñ	8Ñ8€Aä�1 Ô"�H€A€rˆ1Ü�1˜A¨Ô+�H€A€rˆ1Ü�A—E‘E˜"“I—M‘M !§#¡#Ó&¨°¸dÕCÜ�A—E‘E˜"“I—M‘M !§#¡#Ó&¨°¸dÖCr*   ) Únumpyrd   Únumpy.testingr   r   r   r   r   r   r	   r
   r   r   r   r   r   r   r   r   Únumpy.exceptionsr   Úscipy.linalgr   r   r   r)   r\   ro   ÚmarkÚparametrizer{   r   © r*   r   ú<module>r‡      sÔ   ðÛ ß MÑ M÷*÷ *÷ *õ *õ ,Ý Û òò7ò(1òV=ð* ‡�×Ñ˜ 7¨GÐ"4Ó5Ø‡�×Ñ˜˜r 3˜iÓ(ñCó )ó 6ðCð ‡�×Ñ˜ 9¨jÐ"9Ó:Ø‡�×Ñ˜˜r 3˜iÓ(ñDó )ó ;ñDr*   