§
    ¢éUg˜
  ã                   óÐ  — d dl Zd dlmZ d dlmZmZmZ d dlm	Z	m
Z
  ej        dd gd dgg¦  «        Z ej        g d¢g¦  «        Zd ggdd	ggd gggd ggd
d	ggd gggdggdd	ggdgggdggd
d	ggdggged ej        d¦  «        eged
 ej        d¦  «        egede eed          ¦  «        z  ej                             e¦  «         eed          ¦  «        z  ggZd	dgddggg d¢gd	gdgdggg d¢g d¢gd	dgddgddggd	dgddggg d¢gd	gdgdggg d¢g d¢gd	dgddgddggg d¢g d¢g d¢g ej        d¦  «         ej        d¦  «         ej        d¦  «        gZd„ Zd„ Zd„ Zd „ ZdS )!é    N)Únorm)Úassert_Úassert_allcloseÚassert_equal)ÚpolarÚeighé   é   )é   r	   r	   Úrightr   Úlefté	   é   )r   r	   r
   )r
   r   r   é   y      @      @)r   r	   ù              @r   )r   r	   y      @       @)r
   y      @      ð¿y       €      Ày      @       Ày      @      à?)i'  é
   r   )éÿÿÿÿr	   r   )r   r   r	   )r   r   )r   r	   )r	   r   c                 óv   — t          | |¬¦  «        \  }}t          ||d¬¦  «         t          ||d¬¦  «         d S )N©ÚsideçVçž¯Ò<©Úatol)r   r   )Úar   Ú
expected_uÚ
expected_pÚuÚps         úc/var/www/surfInsights/venv3-11/lib/python3.11/site-packages/scipy/linalg/tests/test_decomp_polar.pyÚcheck_precomputed_polarr    &   sJ   € õ �˜ÐÑÔ�D€A€qÝ�A�z¨Ð.Ñ.Ô.Ð.Ý�A�z¨Ð.Ñ.Ô.Ð.Ð.Ð.ó    c                 ó¢  — t          j        t          j        t          ¦  «        j        ¦  «        }t          j        | ¦  «        }|j        \  }}t          | d¬¦  «        \  }}t          |j        ||f¦  «         t          |j        ||f¦  «         t          | 
                    |¦  «        | |¬¦  «         ||k    rOt          |                     ¦   «         j         
                    |¦  «        t          j        |¦  «        d¬¦  «         nNt          | 
                    |                     ¦   «         j        ¦  «        t          j        |¦  «        d¬¦  «         t          |                     ¦   «         j        |¦  «         t          |d¬¦  «        }|t          |¦  «        dk             }t!          |dk                         ¦   «         ¦  «         t          | d	¬¦  «        \  }}t          |j        ||f¦  «         t          |j        ||f¦  «         t          | 
                    |¦  «        | |¬¦  «         ||k    rOt          |                     ¦   «         j         
                    |¦  «        t          j        |¦  «        d¬¦  «         nNt          | 
                    |                     ¦   «         j        ¦  «        t          j        |¦  «        d¬¦  «         t          |                     ¦   «         j        |¦  «         t          |d¬¦  «        }|t          |¦  «        dk             }t!          |dk                         ¦   «         ¦  «         d S )
Nr   r   r   r   T)Úeigvals_onlyg›+¡†›„=r   r   )ÚnpÚsqrtÚfinfoÚfloatÚepsÚasarrayÚshaper   r   r   ÚdotÚconjÚTÚeyer   Úabsr   Úall)	r   Úproduct_atolÚaaÚmÚnr   r   ÚevalsÚnonzero_evalss	            r   Úverify_polarr7   .   s¢  € õ ”7�2œ8¥E™?œ?Ô.Ñ/Ô/€Lå	Œ�A‰Œ€BØŒ8�D€A€qå�˜Ð!Ñ!Ô!�D€A€qÝ�”˜1˜a˜&Ñ!Ô!Ð!Ý�”˜1˜a˜&Ñ!Ô!Ð!å�A—E’E˜!‘H”H˜a lÐ3Ñ3Ô3Ð3ØˆA‚v€vÝ˜Ÿš™œœ
Ÿš qÑ)Ô)­2¬6°!©9¬9¸5ÐAÑAÔAÐAÐAå˜Ÿš˜aŸfšf™hœhœjÑ)Ô)­2¬6°!©9¬9¸5ÐAÑAÔAÐAå�A—F’F‘H”H”J Ñ"Ô"Ð"Ý� Ð&Ñ&Ô&€EØ�#˜e™*œ* uÒ,Ô-€MÝˆ]˜aÒ×$Ò$Ñ&Ô&Ñ'Ô'Ð'å�˜Ð Ñ Ô �D€A€qÝ�”˜1˜a˜&Ñ!Ô!Ð!Ý�”˜1˜a˜&Ñ!Ô!Ð!å�A—E’E˜!‘H”H˜a lÐ3Ñ3Ô3Ð3ØˆA‚v€vÝ˜Ÿš™œœ
Ÿš qÑ)Ô)­2¬6°!©9¬9¸5ÐAÑAÔAÐAÐAå˜Ÿš˜aŸfšf™hœhœjÑ)Ô)­2¬6°!©9¬9¸5ÐAÑAÔAÐAå�A—F’F‘H”H”J Ñ"Ô"Ð"Ý� Ð&Ñ&Ô&€EØ�#˜e™*œ* uÒ,Ô-€MÝˆ]˜aÒ×$Ò$Ñ&Ô&Ñ'Ô'Ð'Ð'Ð'r!   c                  óH   — t           D ]\  } }}}t          | |||¦  «         Œd S ©N)Úprecomputed_casesr    )r   r   r   r   s       r   Útest_precomputed_casesr;   U   s?   € Ý+<ð Að AÑ'ˆˆ4�˜ZÝ  4¨°ZÑ@Ô@Ð@Ð@ðAð Ar!   c                  ó8   — t           D ]} t          | ¦  «         Œd S r9   )Úverify_casesr7   )r   s    r   Útest_verify_casesr>   Z   s(   € Ýð ð ˆÝ�Q‰Œˆˆðð r!   )Únumpyr$   Únumpy.linalgr   Únumpy.testingr   r   r   Úscipy.linalgr   r   ÚarrayÚdiag2Úa13r.   r-   r+   r:   Úemptyr=   r    r7   r;   r>   © r!   r   ú<module>rH      s“  ðØ Ð Ð Ð Ø Ð Ð Ð Ð Ð Ø BÐ BÐ BÐ BÐ BÐ BÐ BÐ BÐ BÐ BØ $Ð $Ð $Ð $Ð $Ð $Ð $Ð $ð 	ˆŒ�1�a�&˜1˜a˜&Ð!Ñ"Ô"€Ø€b„h�	�	�	ˆ{ÑÔ€ð 	€c€UˆG�q�c�U˜a˜S˜EÐ"Ø€c€UˆF�a�S�E˜Q˜C˜5Ð!Ø€c€UˆG�q�c�U˜a˜S˜EÐ"Ø€c€UˆF�a�S�E˜Q˜C˜5Ð!Ø
ˆG�V�R”V˜A‘Y”Y Ð&Ø
ˆF�F�B”F˜1‘I”I˜uÐ%Øˆ'�3�t�t˜C œF‘|”|Ñ# S¤U§Y¢Y¨s¡^¤^°D°D¸¸Q¼±L´LÑ%@ÐAðÐ ð ˆ€Vˆa�ˆVÐØ€Y€Y€KØ€Sˆ1ˆ#�ˆs€OØ€Y€Y�	�	�	ÐØˆ€Vˆa�ˆV�a˜�VÐØˆ€Vˆa�ˆYÐØ€Z€Z€LØ€Sˆ1ˆ#�ˆtÐØ€\€\�>�>�>Ð"Øˆ€Vˆd�Fˆ^˜a ˜VÐ$Ø€^€^�[�[�[ ) ) )Ð,Ø€B„HˆVÑÔØ€B„HˆVÑÔØ€B„HˆVÑÔð€ð$/ð /ð /ð$(ð $(ð $(ðNAð Að Að
ð ð ð ð r!   