§
    fŠtjñs  ã                   ó   — d dl Z d dlZd dlmZmZmZ d dlmZ d dl	m
Z d dlmc mZ d dl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mZmZ d dlm Z  d dl!m"Z"m#Z#m$Z$ d d	l%m&Z&m'Z' e j(        j)        Z) e"e¦  «         G d
„ d¦  «        ¦   «         Z* e"e¦  «         G d„ d¦  «        ¦   «         Z+ e"e¦  «         G d„ d¦  «        ¦   «         Z, e"e¦  «         G d„ d¦  «        ¦   «         Z- e"e¦  «         G d„ d¦  «        ¦   «         Z. e"e¦  «         G d„ d¦  «        ¦   «         Z/ e"e¦  «         G d„ d¦  «        ¦   «         Z0dddddœd„Z1 e"e¦  «         G d„ d¦  «        ¦   «         Z2 e"ej3        ¦  «         G d„ d¦  «        ¦   «         Z4dS )é    N)Úassert_equalÚassert_almost_equalÚassert_allclose)Úgiven)	ÚrombÚnewton_cotesÚcumulative_trapezoidÚ	trapezoidÚquadÚsimpsonÚ
fixed_quadÚqmc_quadÚcumulative_simpson)Ú%_cumulative_simpson_unequal_intervals)ÚstatsÚspecialÚ	integrate)Úskip_xp_invalid_arg)Úmake_xp_test_caseÚxp_default_dtypeÚis_numpy)Úxp_assert_closeÚxp_assert_equalc                   ó   — e Zd Zd„ Zd„ ZdS )ÚTestFixedQuadc                 ór   ‡— dŠdd‰z  z  }t          ˆfd„dd‰¬¦  «        \  }}t          ||d¬¦  «         d S )	Né   é   é   c                 ó   •— | d‰z  dz
  z  S )Nr   r   © )ÚxÚns    €úc/var/www/html/CA-Chatbot/venv/lib/python3.11/site-packages/scipy/integrate/tests/test_quadrature.pyú<lambda>z+TestFixedQuad.test_scalar.<locals>.<lambda>   s   ø€  a¨!¨A©#°©'¡l€ ó    r   ©r#   çê-�™—q=©Úrtol)r   r   )ÚselfÚexpectedÚgotÚ_r#   s       @r$   Útest_scalarzTestFixedQuad.test_scalar   sT   ø€ ØˆØ�a˜‘c‘7ˆÝÐ2Ð2Ð2Ð2°A°q¸AÐ>Ñ>Ô>‰ˆˆQå˜˜X¨EÐ2Ñ2Ô2Ð2Ð2Ð2r&   c                 ó¢   ‡— d}t          j        dd|z  ¦  «        Šd‰dz   z  }t          ˆfd„dd|¬¦  «        \  }}t          ||d¬¦  «         d S )	Nr   r   r   c                 ó"   •— | ‰d d …d f         z  S ©Nr!   )r"   Úps    €r$   r%   z+TestFixedQuad.test_vector.<locals>.<lambda>#   s   ø€  a¨¨1¨1¨1¨d¨7¬¡m€ r&   r   r'   r(   r)   )ÚnpÚaranger   r   )r+   r#   r,   r-   r.   r3   s        @r$   Útest_vectorzTestFixedQuad.test_vector   si   ø€ ØˆÝŒI�a˜˜1™ÑÔˆØ�a˜!‘e‘9ˆÝÐ3Ð3Ð3Ð3°Q¸¸QÐ?Ñ?Ô?‰ˆˆQÝ˜˜X¨EÐ2Ñ2Ô2Ð2Ð2Ð2r&   N)Ú__name__Ú
__module__Ú__qualname__r/   r6   r!   r&   r$   r   r      s2   € € € € € ð3ð 3ð 3ð3ð 3ð 3ð 3ð 3r&   r   c                   ó   — e Zd Zd„ Zd„ ZdS )ÚTestRombc                 óš   — t          t          |                     d¦  «        ¦  «        |                     d|j        ¬¦  «        ¦  «         d S )Ng      1@g      `@©Údtype)r   r   r5   ÚasarrayÚfloat64©r+   Úxps     r$   Ú	test_rombzTestRomb.test_romb)   s<   € Ý�˜RŸYšY t™_œ_Ñ-Ô-¨r¯zªz¸%ÀrÄz¨zÑ/RÔ/RÑSÔSÐSÐSÐSr&   c                 ól  — t          j        d¦  «        }t          j        d|z  ¦  «        }t          |                     |¦  «        ¦  «        }t          d„ t          j        |¦  «        t          j        |¦  «        ¦  «        \  }}t          ||                     ||j	        ¬¦  «        dd¬¦  «         d S )Né   çš™™™™™É?c                 óT   — t          j        t          j        d| z  ¦  «        ¦  «        S )NrF   )r4   ÚcosÚarray©r"   s    r$   r%   z,TestRomb.test_romb_gh_3731.<locals>.<lambda>1   s   € ¥R¤V­B¬H°S¸±U©O¬OÑ%<Ô%<€ r&   r=   g:Œ0âŽyE>r   ©r*   Úatol)
r4   r5   rH   r   r?   r   ÚminÚmaxr   r@   )r+   rB   r"   ÚyÚvalr,   r.   s          r$   Útest_romb_gh_3731zTestRomb.test_romb_gh_3731,   s“   € åŒI�fÑÔˆÝŒF�3�q‘5‰MŒMˆÝ�2—:’:˜a‘=”=Ñ!Ô!ˆÝÐ<Ð<½b¼fÀQ¹i¼iÍÌÐPQÉÌÑSÔS‰ˆ�!Ý˜˜RŸZšZ¨¸¼
˜ZÑCÔCÈ$ÐUVÐWÑWÔWÐWÐWÐWr&   N)r7   r8   r9   rC   rQ   r!   r&   r$   r;   r;   '   s:   € € € € € ðTð Tð TðXð Xð Xð Xð Xr&   r;   c                   ó   — e Zd Zd„ Zd„ ZdS )ÚTestNewtonCotesc                 ó°  — d}t          |d¦  «        \  }}t          ||t          j        ddg¦  «        z  ¦  «         t	          ||dz   dz  ¦  «         d}t          |d¦  «        \  }}t	          ||t          j        g d¢¦  «        z  dz  ¦  «         t	          ||dz   d	z  ¦  «         d}t          |d¦  «        \  }}t	          ||t          j        g d
¢¦  «        z  dz  ¦  «         t	          ||dz   dz  ¦  «         d}t          |d¦  «        \  }}t	          ||t          j        g d¢¦  «        z  dz  ¦  «         t	          ||dz   dz  ¦  «         dS )z5Test the first few degrees, for evenly spaced points.r   ç      à?é   ç      (@r   )ç      ð?ç      @rX   ç      @é   g     €¦@)rX   ç      @r\   rX   ç       @g     P¹@r   )ç      @ç      @@rW   r_   r^   g     €V@é   g     ˆ=AN)r   r   r4   rI   r   )r+   r#   ÚwtsÚerrcoffs       r$   Útest_newton_cotesz!TestNewtonCotes.test_newton_cotes7   sq  € àˆÝ# A qÑ)Ô)‰ˆˆWÝ�S˜!�BœH c¨3 ZÑ0Ô0Ñ0Ñ1Ô1Ð1Ý˜G a¨¡d U¨4¡ZÑ0Ô0Ð0àˆÝ# A qÑ)Ô)‰ˆˆWÝ˜C ¥2¤8¨O¨O¨OÑ#<Ô#<Ñ!<¸SÑ!@ÑAÔAÐAÝ˜G a¨¡d U¨6¡\Ñ2Ô2Ð2àˆÝ# A qÑ)Ô)‰ˆˆWÝ˜C ¥2¤8Ð,@Ð,@Ð,@Ñ#AÔ#AÑ!AÀ#Ñ!EÑFÔFÐFÝ˜G a¨¡d U¨6¡\Ñ2Ô2Ð2àˆÝ# A qÑ)Ô)‰ˆˆWÝ˜C ¥2¤8Ð,HÐ,HÐ,HÑ#IÔ#IÑ!IÈ$Ñ!NÑOÔOÐOÝ˜G a¨¡d U¨9¡_Ñ5Ô5Ð5Ð5Ð5r&   c                 óV  — t          j        g d¢¦  «        }|dz  }t          |¦  «        \  }}d}t          j        ||¦  «        }t	          ||¦  «         t          j        g d¢¦  «        }|dz  }t          |¦  «        \  }}d}t          j        ||¦  «        }t	          ||¦  «         dS )z9Test newton_cotes with points that are not evenly spaced.)ç        g      ø?ç       @r   gUUUUUU@)re   gffffffö?gÍÌÌÌÌÌ @r\   ç      "@N)r4   rI   r   Údotr   )r+   r"   rO   ra   rb   Úexact_integralÚnumeric_integrals          r$   Útest_newton_cotes2z"TestNewtonCotes.test_newton_cotes2M   s¯   € õ ŒH�_�_�_Ñ%Ô%ˆØˆq‰DˆÝ# A‘”‰ˆˆWØˆÝœ6 # q™>œ>ÐÝÐ,¨nÑ=Ô=Ð=åŒHÐ)Ð)Ð)Ñ*Ô*ˆØˆq‰DˆÝ# A‘”‰ˆˆWØˆÝœ6 # q™>œ>ÐÝÐ,¨nÑ=Ô=Ð=Ð=Ð=r&   N)r7   r8   r9   rc   rk   r!   r&   r$   rS   rS   5   s2   € € € € € ð6ð 6ð 6ð,>ð >ð >ð >ð >r&   rS   c                   ó`   — e Zd Zd„ Zej                             dddg¦  «        d„ ¦   «         ZdS )ÚTestSimpsonc           
      óz  — t          j        d¦  «        }t          t          |¦  «        d¦  «         t          t          |d¬¦  «        d¦  «         t          t          |t          j        ddd¦  «        ¬¦  «        d	¦  «         t          j        d
dd¦  «        }d„ }t          t           ||¦  «        |¬¦  «        d¦  «         t          j        d
dd¦  «        }t          t           ||¦  «        d¬¦  «        d¦  «         t          j        d¦  «                             dd¦  «        }t          j        d¦  «                             ddd¦  «        } ||¦  «        }t          d¦  «        D ]Å}t          |||¬¦  «        }t          j        |dg¬¦  «        }|D ]—}t          |j
        ¦  «        }	|	                     |t          d ¦  «        ¦  «         |t          |	¦  «                 d         dz  dz  |t          |	¦  «                 d         dz  dz  z
  }
t          ||j
                 |
¦  «         Œ˜ŒÆt          j        d¦  «                             dd¦  «        } ||¦  «        }t          ||d¬¦  «        }d|d d …d
f         |d d …df         z   z  |d d …d
f         |d d …df         z
  z  }
t          ||
¦  «         t          j        d¦  «                             dd¦  «        }t          j        d¦  «                             ddd¦  «        } ||¦  «        }t          d¦  «        D ]Å}t          |||¬¦  «        }t          j        |dg¬¦  «        }|D ]—}t          |j
        ¦  «        }	|	                     |t          d ¦  «        ¦  «         |t          |	¦  «                 d         dz  dz  |t          |	¦  «                 d         dz  dz  z
  }
t          ||j
                 |
¦  «         Œ˜ŒÆt          j        dg¦  «        }t          j        |d¦  «        }t          t          ||d¬¦  «        d¦  «         t          t          ||d¬¦  «        d¦  «         t          j        g d¢¦  «        }t          j        |d¦  «        }t          t          ||d¬¦  «        d¦  «         t          t          ||d¬¦  «        d¦  «         t          j        g d¢g d¢g d¢g¦  «        }t          j        |d¦  «        }g d¢}d gdz  }t          t          ||d¬¦  «        |¦  «         t          t          ||d¬¦  «        |¦  «         t          j        g d¢g d¢g d!¢g¦  «        }t          j        |d¦  «        }g d"¢}g d#¢}t          t          ||d¬¦  «        |¦  «         t          t          ||d¬¦  «        |¦  «         d S )$NrE   é€   rU   ©Údxé@   r   r   rJ   é    r   c                 ó   — | dz  S )Nr   r!   rJ   s    r$   Úfz#TestSimpson.test_simpson.<locals>.fi   s   € Ø�a‘4ˆKr&   ç      5@r`   rf   ér   é   ç      P@rV   ©r"   ÚaxisÚmulti_index)Úflagséÿÿÿÿé   r   é   r[   é}   re   )rV   rV   rV   rV   )r   r   r   r   )re   re   re   re   ç«ªªªªJe@)r   r   rx   rs   )re   g      a@g      ‘@g      Á@)r‚   r‚   gªªªª*UÅ@)r4   r5   r   r   Úlinspacer   ÚreshapeÚrangeÚnditerÚlistr|   ÚinsertÚsliceÚtuplerI   Úpower)r+   rO   r"   ru   ÚaÚiÚrÚitr.   ÚidxÚintegralÚ	zero_axisÚdefault_axiss                r$   Útest_simpsonzTestSimpson.test_simpsona   s“  € ÝŒI�b‰MŒMˆÝ•W˜Q‘Z”Z Ñ%Ô%Ð%Ý•W˜Q 3Ð'Ñ'Ô'¨Ñ,Ô,Ð,Ý•W˜Q¥"¤+¨a°°BÑ"7Ô"7Ð8Ñ8Ô8¸"Ñ=Ô=Ð=õ ŒK˜˜1˜aÑ Ô ˆð	ð 	ð 	õ 	�   !¡¤¨Ð*Ñ*Ô*¨DÑ1Ô1Ð1õ ŒK˜˜1˜aÑ Ô ˆÝ�   !¡¤¨Ð-Ñ-Ô-¨sÑ3Ô3Ð3õ ŒI�b‰MŒM×!Ò! ! QÑ'Ô'ˆÝŒI�c‰NŒN×"Ò" 1 a¨Ñ+Ô+ˆØˆAˆa‰DŒDˆÝ�q‘”ð 	=ð 	=ˆAÝ˜˜Q QÐ'Ñ'Ô'ˆAÝ”˜1 ] OÐ4Ñ4Ô4ˆBØð =ð =�Ý˜2œ>Ñ*Ô*�Ø—
’
˜1�e D™kœkÑ*Ô*Ð*Ø�U 3™ZœZœ=¨Ô,¨aÑ/°!Ñ3°a½¸c¹
¼
´mÀAÔ6FÈÑ6IÈAÑ6MÑM�Ý  "¤.Ô 1°8Ñ<Ô<Ð<Ð<ð	=õ ŒI�b‰MŒM×!Ò! ! QÑ'Ô'ˆØˆAˆa‰DŒDˆÝ�A˜ Ð$Ñ$Ô$ˆà˜!˜A˜A˜A˜q˜Dœ' A a a a¨ d¤GÑ+Ñ,°°!°!°!°Q°$´¸!¸A¸A¸A¸q¸D¼'Ñ0AÑBˆÝ˜˜8Ñ$Ô$Ð$õ ŒI�b‰MŒM×!Ò! ! QÑ'Ô'ˆÝŒI�c‰NŒN×"Ò" 1 a¨Ñ+Ô+ˆØˆAˆa‰DŒDˆÝ�q‘”ð 	=ð 	=ˆAÝ˜˜Q QÐ'Ñ'Ô'ˆAÝ”˜1 ] OÐ4Ñ4Ô4ˆBØð =ð =�Ý˜2œ>Ñ*Ô*�Ø—
’
˜1�e D™kœkÑ*Ô*Ð*Ø�U 3™ZœZœ=¨Ô,¨aÑ/°!Ñ3°a½¸c¹
¼
´mÀAÔ6FÈÑ6IÈAÑ6MÑM�Ý  "¤.Ô 1°8Ñ<Ô<Ð<Ð<ð	=õ ŒH�a�S‰MŒMˆÝŒH�Q˜‰NŒNˆÝ�  Q¨QÐ/Ñ/Ô/°Ñ5Ô5Ð5Ý�  Q¨RÐ0Ñ0Ô0°#Ñ6Ô6Ð6åŒH�\�\�\Ñ"Ô"ˆÝŒH�Q˜‰NŒNˆÝ�  Q¨QÐ/Ñ/Ô/°Ñ5Ô5Ð5Ý�  Q¨RÐ0Ñ0Ô0°#Ñ6Ô6Ð6åŒH�l�l�l L L L°,°,°,Ð?Ñ@Ô@ˆÝŒH�Q˜‰NŒNˆØ(Ð(Ð(ˆ	Ø!�{ Q‘ˆÝ�  Q¨QÐ/Ñ/Ô/°Ñ;Ô;Ð;å�  Q¨RÐ0Ñ0Ô0°,Ñ?Ô?Ð?åŒH�l�l�l L L L°.°.°.ÐAÑBÔBˆÝŒH�Q˜‰NŒNˆØ0Ð0Ð0ˆ	Ø>Ð>Ð>ˆÝ�  Q¨QÐ/Ñ/Ô/°Ñ;Ô;Ð;Ý�  Q¨RÐ0Ñ0Ô0°,Ñ?Ô?Ð?Ð?Ð?r&   ÚdroplastFTc                 óú   — t          j        g d¢g d¢g¦  «        }|r|d d …d d…f         }t          |d¬¦  «        }t          t          j        |t           j        ¬¦  «        d¬¦  «        }t	          ||¦  «         d S )N)r   r   r   r   r   r   éüÿÿÿr[   )r   r   r   r—   é
   é   éþÿÿÿr˜   r~   ©r{   r=   )r4   rI   r   r@   r   )r+   r•   rO   Úresultr,   s        r$   Útest_simpson_2d_integer_no_xz(TestSimpson.test_simpson_2d_integer_no_x¯   s”   € õ ŒHÐ/Ð/Ð/Ø3Ð3Ð3ð5ñ 6ô 6ˆàð 	Ø�!�!�!�S�b�S�&”	ˆAÝ˜ Ð$Ñ$Ô$ˆÝ�2œ8 A­R¬ZÐ8Ñ8Ô8¸rÐBÑBÔBˆÝ�V˜XÑ&Ô&Ð&Ð&Ð&r&   N)r7   r8   r9   r”   ÚpytestÚmarkÚparametrizer�   r!   r&   r$   rm   rm   _   s\   € € € € € ðK@ð K@ð K@ð\ „[×Ò˜Z¨%°¨Ñ7Ô7ð	'ð 	'ñ 8Ô7ð	'ð 	'ð 	'r&   rm   c                   óx   — e Zd Zd„ Zd„ Zd„ Zd„ Zej         	                    dddg¦  «        d„ ¦   «         Z
d	„ Zd
S )ÚTestCumulative_trapezoidc                 ó   — |                      ddd¬¦  «        }|}t          ||d¬¦  «        }|                     g d¢¦  «        }t          ||¦  «         t          ||d ¬¦  «        }t          ||dd …         ¦  «         d S )	Nrš   r   r[   ©Únumr   ©Úinitial)re   ç      ø¿ç       Àr¨   re   r   ©rƒ   r	   r?   r   )r+   rB   r"   rO   Úy_intÚ
y_expecteds         r$   Útest_1dz TestCumulative_trapezoid.test_1d¾   sŒ   € Ø�KŠK˜˜A 1ˆKÑ%Ô%ˆØˆÝ$ Q¨°1Ð5Ñ5Ô5ˆØ—Z’ZÐ 9Ð 9Ð 9Ñ:Ô:ˆ
Ý˜˜zÑ*Ô*Ð*å$ Q¨°4Ð8Ñ8Ô8ˆÝ˜˜z¨!¨"¨"œ~Ñ.Ô.Ð.Ð.Ð.r&   c                 óÊ  — |                      |                     dt          |¦  «        ¬¦  «        d¦  «        }|}t          ||d¬¦  «        }|                     g d¢g d¢gg d¢g d	¢gg d
¢g d¢gg¦  «        }t          ||¦  «         g d¢}t          g d¢|¦  «        D ]E\  }}t          ||d|¬¦  «        }|j        dk    sJ ‚t          ||d |¬¦  «        }|j        |k    sJ ‚ŒFd S )Né   r=   ©rV   r   r   r   r¦   )re   rU   rf   ç      @)re   r±   ç      $@ç     €0@)re   ç      !@ç      2@g     €<@)re   ç      )@g      :@ç     @D@)re   r³   g      A@g     @J@)re   ç     €4@g      E@g      P@))r   r   r   )rV   r   r   )rV   r   rV   ©r   r   r   )r§   r{   )r„   r5   r   r	   r?   r   ÚzipÚshape)	r+   rB   r"   rO   r«   r¬   Úshapesr{   r»   s	            r$   Útest_y_nd_x_ndz'TestCumulative_trapezoid.test_y_nd_x_ndÈ   sC  € Ø�JŠJ�r—y’y Õ2BÀ2Ñ2FÔ2F�yÑGÔGÈÑSÔSˆØˆÝ$ Q¨°1Ð5Ñ5Ô5ˆØ—Z’ZÐ"4Ð"4Ð"4Ø"6Ð"6Ð"6ð"8à"6Ð"6Ð"6Ø"7Ð"7Ð"7ð"9à"7Ð"7Ð"7Ø"7Ð"7Ð"7ð"9ð	!:ñ ;ô ;ˆ
õ 	˜˜zÑ*Ô*Ð*ð 3Ð2Ð2ˆÝ˜y˜y˜y¨&Ñ1Ô1ð 	(ð 	(‰KˆD�%Ý(¨¨A°q¸tÐDÑDÔDˆEØ”; )Ò+Ð+Ð+Ð+Ý(¨¨A°tÀ$ÐGÑGÔGˆEØ”; %Ò'Ð'Ð'Ð'Ð'ð		(ð 	(r&   c           	      ó:  — |                      |                     dt          |¦  «        ¬¦  «        d¦  «        }|                     dt          |¦  «        ¬¦  «        dz  }|                     g d¢g d¢gg d¢g d	¢gg¦  «        |                     g d
¢gg d¢gg d¢gg¦  «        |                     g d¢g d¢gg d¢g d¢gg d¢g d¢gg¦  «        f}t	          g d¢|¦  «        D ];\  }}t          ||d |j        |         …         |d ¬¦  «        }t          ||¦  «         Œ<d S )Nr¯   r=   r°   r   r   )rY   ç      @rZ   r^   )r]   rg   r²   ç      &@)g      D@g      F@g      H@g      J@)g      L@g      N@ry   g      Q@)rf   r\   rY   r¿   )r²   rÀ   rW   ç      *@)rµ   g      3@ç      4@rv   )rU   r¿   g     €1@)r±   rv   g     ÀJ@)r´   g     €B@g     `V@)r¶   g     €J@g     `_@)r³   g     @Q@g     0d@)r¸   g     @U@g     °h@r¹   )r"   r{   r§   )r„   r5   r   r?   rº   r	   r»   r   )r+   rB   rO   r"   Úys_expectedr{   r¬   r«   s           r$   Útest_y_nd_x_1dz'TestCumulative_trapezoid.test_y_nd_x_1dÝ   s¡  € Ø�JŠJ�r—y’y Õ2BÀ2Ñ2FÔ2F�yÑGÔGÈÑSÔSˆØ�IŠI�aÕ/°Ñ3Ô3ˆIÑ4Ô4°aÑ7ˆð �JŠJÐ)Ð)Ð)Ø+Ð+Ð+ð-à-Ð-Ð-Ø-Ð-Ð-ð/ð0ñ 1ô 1ð �JŠJÐ)Ð)Ð)Ð*Ø-Ð-Ð-Ð.Ø-Ð-Ð-Ð.ð0ñ 1ô 1ð �JŠJ˜˜˜Ø)Ð)Ð)ð+à)Ð)Ð)Ø+Ð+Ð+ð-à+Ð+Ð+Ø+Ð+Ð+ð-ð	.ñ /ô /ð0ˆõ !$ I I I¨{Ñ ;Ô ;ð 	/ð 	/ÑˆD�*Ý(¨¨a°°´¸´°Ô.?ÀdØ15ð7ñ 7ô 7ˆEå˜E :Ñ.Ô.Ð.Ð.ð	/ð 	/r&   c                 óô  — |                      ddd¬¦  «        }t          |¦  «        }|                     g d¢¦  «        }t          ||¦  «         t          |d¬¦  «        }|                     g d¢¦  «        }t          ||¦  «         t          |d	¬
¦  «        }|                     g d¢¦  «        }t          ||¦  «         t          |d	d¬¦  «        }|                     g d¢¦  «        }t          ||¦  «         d S )Nrš   r   r[   r¤   )r¨   r©   r¨   re   r   r¦   )r   r¨   r©   r¨   re   rV   rp   )ç      Àç      ÀrÆ   re   )rq   r§   )r   rÆ   rÇ   rÆ   re   rª   )r+   rB   rO   r«   r¬   s        r$   Útest_x_nonez$TestCumulative_trapezoid.test_x_noneõ   s  € Ø�KŠK˜˜A 1ˆKÑ%Ô%ˆå$ QÑ'Ô'ˆØ—Z’ZÐ 5Ð 5Ð 5Ñ6Ô6ˆ
Ý˜˜zÑ*Ô*Ð*å$ Q°Ð2Ñ2Ô2ˆØ—Z’ZÐ 8Ð 8Ð 8Ñ9Ô9ˆ
Ý˜˜zÑ*Ô*Ð*å$ Q¨1Ð-Ñ-Ô-ˆØ—Z’ZÐ 5Ð 5Ð 5Ñ6Ô6ˆ
Ý˜˜zÑ*Ô*Ð*å$ Q¨1°aÐ8Ñ8Ô8ˆØ—Z’ZÐ 8Ð 8Ð 8Ñ9Ô9ˆ
Ý˜˜zÑ*Ô*Ð*Ð*Ð*r&   r§   r   rU   c                 óÀ   — |                      ddd¬¦  «        }t          j        t          d¬¦  «        5  t	          ||¬¦  «         ddd¦  «         dS # 1 swxY w Y   dS )z4If initial is not None or 0, a ValueError is raised.r   r˜   r¤   z	`initial`©Úmatchr¦   N)rƒ   rž   ÚraisesÚ
ValueErrorr	   )r+   r§   rB   rO   s       r$   Útest_initial_errorz+TestCumulative_trapezoid.test_initial_error  s§   € ð
 �KŠK˜˜2 2ˆKÑ&Ô&ˆÝŒ]�:¨[Ð9Ñ9Ô9ð 	5ð 	5Ý  ¨GÐ4Ñ4Ô4Ð4ð	5ð 	5ð 	5ñ 	5ô 	5ð 	5ð 	5ð 	5ð 	5ð 	5ð 	5ð 	5øøøð 	5ð 	5ð 	5ð 	5ð 	5ð 	5s   ´AÁAÁAc                 ó´   — t          j        t          d¬¦  «        5  t          |                     g ¦  «        ¬¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )NúAt least one point is requiredrÊ   ©rO   )rž   rÌ   rÍ   r	   r?   rA   s     r$   Útest_zero_len_yz(TestCumulative_trapezoid.test_zero_len_y  s—   € ÝŒ]�:Ð-MÐNÑNÔNð 	3ð 	3Ý  2§:¢:¨b¡>¤>Ð2Ñ2Ô2Ð2ð	3ð 	3ð 	3ñ 	3ô 	3ð 	3ð 	3ð 	3ð 	3ð 	3ð 	3ð 	3øøøð 	3ð 	3ð 	3ð 	3ð 	3ð 	3s   œ$AÁAÁAN)r7   r8   r9   r­   r½   rÄ   rÈ   rž   rŸ   r    rÎ   rÒ   r!   r&   r$   r¢   r¢   ¼   s˜   € € € € € ð/ð /ð /ð(ð (ð (ð*/ð /ð /ð0+ð +ð +ð& „[×ÒØ�A�s�8ñô ð5ð 5ñô ð5ð3ð 3ð 3ð 3ð 3r&   r¢   c                   ó<   — e Zd Zd„ Zd„ Zd„ Zed„ ¦   «         Zd„ ZdS )ÚTestTrapezoidc           	      ó0  — |                      ddd¦  «        }t          |                     d|dz  z  ¦  «        |                     d|                     |j        ¦  «        z  ¦  «        z  d¬¦  «        }t          ||                     d¦  «        ¦  «         d S )Néöÿÿÿr˜   gš™™™™™¹?ç      à¿r   rp   rX   )r5   r
   ÚexpÚsqrtr?   Úpir   )r+   rB   r"   rŽ   s       r$   Útest_simplezTestTrapezoid.test_simple  s‚   € Ø�IŠI�c˜2˜rÑ"Ô"ˆÝ�b—f’f˜S 1¨¡6™\Ñ*Ô*¨R¯WªW°Q¸¿ºÀBÄEÑ9JÔ9JÑ5JÑ-KÔ-KÑKÐPSÐTÑTÔTˆå˜˜2Ÿ:š: c™?œ?Ñ+Ô+Ð+Ð+Ð+r&   c                 óD  — |                      ddd¦  «        }|                      ddd¦  «        }|                      ddd¦  «        }|                     |¦  «        |d         |d         z
  z  }|dxx         dz  cc<   |dxx         dz  cc<   |                     |¦  «        |d         |d         z
  z  }|dxx         dz  cc<   |dxx         dz  cc<   |                     |¦  «        |d         |d         z
  z  }|dxx         dz  cc<   |dxx         dz  cc<   |d d …d d f         |d d d …d f         z   |d d d d …f         z   }|                     ||d d …d d f         z  d¬¦  «        }	|                     ||d d d …d f         z  d¬¦  «        }
|                     ||d d d d …f         z  d¬¦  «        }t          ||d d …d d f         d¬	¦  «        }t	          ||	¦  «         t          ||d d d …d f         d¬	¦  «        }t	          ||
¦  «         t          ||d d d d …f         d¬	¦  «        }t	          ||¦  «         t          ||d¬	¦  «        }t	          ||	¦  «         t          ||d¬	¦  «        }t	          ||
¦  «         t          ||d¬	¦  «        }t	          ||¦  «         d S )
Nr   r   rV   r   r   é   r~   r›   rz   )rƒ   Ú	ones_likeÚsumr
   r   )r+   rB   r"   rO   ÚzÚwxÚwyÚwzÚqÚqxÚqyÚqzrŽ   s                r$   Ú	test_ndimzTestTrapezoid.test_ndim  sá  € Ø�KŠK˜˜1˜aÑ Ô ˆØ�KŠK˜˜1˜aÑ Ô ˆØ�KŠK˜˜1˜bÑ!Ô!ˆà�\Š\˜!‰_Œ_  !¤ q¨¤t¡Ñ,ˆØ
ˆ1ˆˆŒ�‰
ˆˆ‰Ø
ˆ2ˆˆŒ�!‰ˆˆ‰Ø�\Š\˜!‰_Œ_  !¤ q¨¤t¡Ñ,ˆØ
ˆ1ˆˆŒ�‰
ˆˆ‰Ø
ˆ2ˆˆŒ�!‰ˆˆ‰Ø�\Š\˜!‰_Œ_  !¤ q¨¤t¡Ñ,ˆØ
ˆ1ˆˆŒ�‰
ˆˆ‰Ø
ˆ2ˆˆŒ�!‰ˆˆ‰àˆaˆaˆa��tˆmÔ˜q  a a a¨ œÑ.°°4¸¸a¸a¸a°<´Ñ@ˆà�VŠV�A˜˜1˜1˜1˜d D˜=Ô)Ñ)°ˆVÑ2Ô2ˆØ�VŠV�A˜˜4    D˜=Ô)Ñ)°ˆVÑ2Ô2ˆØ�VŠV�A˜˜4  q q q˜=Ô)Ñ)°ˆVÑ2Ô2ˆõ �a˜1˜Q˜Q˜Q  d˜]Ô+°!Ð4Ñ4Ô4ˆÝ˜˜2ÑÔÐÝ�a˜1˜T ! ! ! T˜\œ?°Ð3Ñ3Ô3ˆÝ˜˜2ÑÔÐÝ�a˜1˜T 4¨¨¨˜\œ?°Ð3Ñ3Ô3ˆÝ˜˜2ÑÔÐõ �a˜1 1Ð%Ñ%Ô%ˆÝ˜˜2ÑÔÐÝ�a˜1 1Ð%Ñ%Ô%ˆÝ˜˜2ÑÔÐÝ�a˜1 1Ð%Ñ%Ô%ˆÝ˜˜2ÑÔÐÐÐr&   c           
      ó¤  — |                      |                     ddd¦  «        d¦  «        }|                      |                     ddd¦  «        d¦  «        }|                     ddd¦  «        }t          t          ||d¬¦  «        |¦  «         t          t          ||                     g d	¢¦  «        d¬¦  «        |¦  «         t          t          ||                     g d	¢¦  «        d d …d f         d¬¦  «        |¦  «         t          j        t          ¦  «        5  t          ||                     g d	¢¦  «        d d d …f         d¬¦  «         d d d ¦  «         n# 1 swxY w Y   |                     g d
¢¦  «        }t          t          ||d¬¦  «        |¦  «         t          t          ||                     ddd¦  «        d¬¦  «        |¦  «         d S )Nr   é   é   )rV   r˜   éÈ   i|  r˜   rz   )r   r²   rÂ   )r·   g     P`@g     �k@r   é	   )r„   rƒ   r   r
   r?   rž   rÌ   Ú	Exception)r+   rB   r"   rO   Úout0Úout1s         r$   Útest_gh21908zTestTrapezoid.test_gh21908C  s  € à�JŠJ�r—{’{ 1 b¨"Ñ-Ô-¨wÑ7Ô7ˆØ�JŠJ�r—{’{ 1 b¨"Ñ-Ô-¨wÑ7Ô7ˆà�{Š{˜3  RÑ(Ô(ˆÝ�	 ! q¨qÐ1Ñ1Ô1°4Ñ8Ô8Ð8Ý�	 ! r§z¢z°-°-°-Ñ'@Ô'@ÀqÐIÑIÔIÈ4ÑPÔPÐPåÝ�a˜2Ÿ:š: m m mÑ4Ô4°Q°Q°Q¸°WÔ=ÀAÐFÑFÔFØñ	
ô 	
ð 	
õ Œ]�9Ñ%Ô%ð 	Gð 	Gå�a˜2Ÿ:š: m m mÑ4Ô4°T¸1¸1¸1°WÔ=ÀAÐFÑFÔFÐFð	Gð 	Gð 	Gñ 	Gô 	Gð 	Gð 	Gð 	Gð 	Gð 	Gð 	Gøøøð 	Gð 	Gð 	Gð 	Gð �zŠzÐ/Ð/Ð/Ñ0Ô0ˆÝ�	 ! q¨qÐ1Ñ1Ô1°4Ñ8Ô8Ð8ÝÝ�a˜2Ÿ;š; q¨!¨RÑ0Ô0°qÐ9Ñ9Ô9Øñ	
ô 	
ð 	
ð 	
ð 	
s   Ä2EÅEÅ Ec                 óÂ  — t          j        d¦  «        }||z  }|dk    }t           j                             ||¬¦  «        }d}t	          t          ||¦  «        |¦  «         t           j                             ||¬¦  «        }t	          t          ||¦  «        |¦  «         t           j                             ||¬¦  «        }t	          t          ||¦  «        |¦  «         d S )Nr[   r   )ÚmaskrÁ   )r4   r5   ÚmarI   r   r
   )r+   rB   r"   rO   ró   ÚymrŽ   Úxms           r$   Útest_maskedzTestTrapezoid.test_masked[  sÁ   € õ ŒI�a‰LŒLˆØ�‰EˆØ�AŠvˆÝŒU�[Š[˜ ˆ[Ñ&Ô&ˆØˆÝ�	 " aÑ(Ô(¨!Ñ,Ô,Ð,åŒU�[Š[˜ ˆ[Ñ&Ô&ˆÝ�	 " bÑ)Ô)¨1Ñ-Ô-Ð-åŒU�[Š[˜ ˆ[Ñ&Ô&ˆÝ�	 ! RÑ(Ô(¨!Ñ,Ô,Ð,Ð,Ð,r&   c                 ó6  — t          t          d¦  «        ¦  «        }d„ |D ¦   «         }t          j        |t          j        ¬¦  «        }t          j        |t          j        ¬¦  «        }t          ||¦  «        }t          ||¦  «        }t          ||¦  «         d S )Nr[   c                 ó   — g | ]}||z  ‘ŒS r!   r!   )Ú.0Úts     r$   ú
<listcomp>z1TestTrapezoid.test_array_like.<locals>.<listcomp>n  s   € ÐÐÐ�qˆQ�‰UÐÐÐr&   r=   )r‡   r…   r4   r?   r@   r
   r   )r+   r"   rO   ÚxarrÚyarrÚresÚresarrs          r$   Útest_array_likezTestTrapezoid.test_array_likel  s„   € Ý•�q‘”‰NŒNˆØÐ˜AÐÑÔˆÝŒz˜!¥2¤:Ð.Ñ.Ô.ˆÝŒz˜!¥2¤:Ð.Ñ.Ô.ˆÝ˜˜1‰oŒoˆÝ˜4 Ñ&Ô&ˆÝ˜˜VÑ$Ô$Ð$Ð$Ð$r&   N)	r7   r8   r9   rÛ   rè   rñ   r   r÷   r  r!   r&   r$   rÔ   rÔ     sm   € € € € € ð,ð ,ð ,ð#ð #ð #ðJ
ð 
ð 
ð0 ð-ð -ñ Ôð-ð %ð %ð %ð %ð %r&   rÔ   c                   ó`  — e Zd Zd„ Zdd„Zej                             dddg¦  «        ej                             ddd	g¦  «        d
„ ¦   «         ¦   «         Zej                             dddgddgddgddgg¦  «        d„ ¦   «         Z	ej                             dddg¦  «        d„ ¦   «         Z
d„ ZdS )ÚTestQMCQuadc           	      óÒ  ‡— ‰                      ddg¦  «        }‰                      ddg¦  «        }d}t          j        t          |¬¦  «        5  t	          d||¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          j        t
          |¬¦  «        5  t	          d„ ||¦  «         d d d ¦  «         n# 1 swxY w Y   ˆfd„}d	}t          ‰¦  «        rEt          j        t          |¬¦  «        5  t	          |||¦  «         d d d ¦  «         n# 1 swxY w Y   nDt          j        t
          |¬¦  «        5  t	          |||¦  «         d d d ¦  «         n# 1 swxY w Y   d
}t          j        t          |¬¦  «        5  t	          d„ ||d¬¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          j        t          |¬¦  «        5  t	          d„ ||d¬¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          j        t          |¬¦  «        5  t	          d„ ||d¬¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          j        t
          |¬¦  «        5  t	          d„ ||t          j
                             d¦  «        ¬¦  «         d d d ¦  «         n# 1 swxY w Y   d}t          j        t          |¬¦  «        5  t	          d„ ||d¬¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )Nre   rX   z`func` must be callable.rÊ   za duckz/`func` must evaluate the integrand at points...c                  ó   — dS ©Nr   r!   r!   r&   r$   r%   z3TestQMCQuad.test_input_validation.<locals>.<lambda>‚  s   € ˜Q€ r&   c                 óH   •— | j         dk    sJ ‚‰                     | ¦  «        S r  )Úndimrß   ©r"   rB   s    €r$   Úfuncz/TestQMCQuad.test_input_validation.<locals>.func„  s"   ø€ Ø”6˜Q’;�;�;�;Ø—6’6˜!‘9”9Ðr&   z8Exception encountered when attempting vectorized call...z`n_points` must be an integer.c                 ó   — dS r  r!   rJ   s    r$   r%   z3TestQMCQuad.test_input_validation.<locals>.<lambda>‘  ó   € ˜q€ r&   g     �@)Ún_pointsz!`n_estimates` must be an integer.c                 ó   — dS r  r!   rJ   s    r$   r%   z3TestQMCQuad.test_input_validation.<locals>.<lambda>•  r  r&   r´   )Ún_estimatesz8`qrng` must be an instance of scipy.stats.qmc.QMCEngine.c                 ó   — dS r  r!   rJ   s    r$   r%   z3TestQMCQuad.test_input_validation.<locals>.<lambda>™  r  r&   )Úqrngz8`qrng` must be initialized with dimensionality equal to c                 ó   — dS r  r!   rJ   s    r$   r%   z3TestQMCQuad.test_input_validation.<locals>.<lambda>�  r  r&   r   z,`log` must be boolean \(`True` or `False`\).c                 ó   — dS r  r!   rJ   s    r$   r%   z3TestQMCQuad.test_input_validation.<locals>.<lambda>¡  r  r&   r˜   ©Úlog)r?   rž   rÌ   Ú	TypeErrorr   rÍ   r   ÚwarnsÚUserWarningr   ÚqmcÚSobol)r+   rB   rŒ   ÚbÚmessager
  s    `    r$   Útest_input_validationz!TestQMCQuad.test_input_validationx  sA  ø€ Ø�JŠJ˜˜B�xÑ Ô ˆØ�JŠJ˜˜B�xÑ Ô ˆà,ˆÝŒ]�9¨GÐ4Ñ4Ô4ð 	%ð 	%Ý�X˜q !Ñ$Ô$Ð$ð	%ð 	%ð 	%ñ 	%ô 	%ð 	%ð 	%ð 	%ð 	%ð 	%ð 	%øøøð 	%ð 	%ð 	%ð 	%ð DˆÝŒ]�:¨WÐ5Ñ5Ô5ð 	&ð 	&Ý�Y�Y  1Ñ%Ô%Ð%ð	&ð 	&ð 	&ñ 	&ô 	&ð 	&ð 	&ð 	&ð 	&ð 	&ð 	&øøøð 	&ð 	&ð 	&ð 	&ð	ð 	ð 	ð 	ð 	ð MˆÝ�B‰<Œ<ð 	%Ý”�k°Ð9Ñ9Ô9ð %ð %Ý˜˜q !Ñ$Ô$Ð$ð%ð %ð %ñ %ô %ð %ð %ð %ð %ð %ð %øøøð %ð %ð %ð %øõ ”�z°Ð9Ñ9Ô9ð %ð %Ý˜˜q !Ñ$Ô$Ð$ð%ð %ð %ñ %ô %ð %ð %ð %ð %ð %ð %øøøð %ð %ð %ð %ð 3ˆÝŒ]�9¨GÐ4Ñ4Ô4ð 	9ð 	9Ý�[�[ ! Q°Ð8Ñ8Ô8Ð8ð	9ð 	9ð 	9ñ 	9ô 	9ð 	9ð 	9ð 	9ð 	9ð 	9ð 	9øøøð 	9ð 	9ð 	9ð 	9ð 6ˆÝŒ]�9¨GÐ4Ñ4Ô4ð 	9ð 	9Ý�[�[ ! Q°CÐ8Ñ8Ô8Ð8ð	9ð 	9ð 	9ñ 	9ô 	9ð 	9ð 	9ð 	9ð 	9ð 	9ð 	9øøøð 	9ð 	9ð 	9ð 	9ð MˆÝŒ]�9¨GÐ4Ñ4Ô4ð 	7ð 	7Ý�[�[ ! Q¨XÐ6Ñ6Ô6Ð6ð	7ð 	7ð 	7ñ 	7ô 	7ð 	7ð 	7ð 	7ð 	7ð 	7ð 	7øøøð 	7ð 	7ð 	7ð 	7ð MˆÝŒ]�:¨WÐ5Ñ5Ô5ð 	Að 	AÝ�[�[ ! Q­U¬Y¯_ª_¸QÑ-?Ô-?Ð@Ñ@Ô@Ð@ð	Að 	Að 	Añ 	Aô 	Að 	Að 	Að 	Að 	Að 	Að 	Aøøøð 	Að 	Að 	Að 	Að BˆÝŒ]�9¨GÐ4Ñ4Ô4ð 	0ð 	0Ý�[�[ ! Q¨BÐ/Ñ/Ô/Ð/ð	0ð 	0ð 	0ñ 	0ô 	0ð 	0ð 	0ð 	0ð 	0ð 	0ð 	0ð 	0øøøð 	0ð 	0ð 	0ð 	0ð 	0ð 	0s¢   ÁA+Á+A/Á2A/ÂB2Â2B6Â9B6Ã.DÄDÄDÄ3EÅEÅEÅ9FÆFÆ!FÇG#Ç#G'Ç*G'ÈH,È,H0È3H0É2JÊJÊJÊ:KËK Ë#K é   r   Nc           	      óÐ  ‡‡— t          ‰¦  «        }|€t          j        d¦  «        }d}t          j        |¦  «        }t          j        |¦  «        }ˆfd„Št          j                             d¦  «        }	t          j         	                    ||	¬¦  «        }
t          j        |¦  «        }t          j        |¦  «        |z  }t          ‰‰                     ||¬¦  «        ‰                     ||¬¦  «        |||
¬¦  «        }t          j                             ||||¬¦  «        }t          j        |dz
  d	¦  «        |j        z  }t#          |j        ‰                     ||¬¦  «        |¬
¦  «         t          j        |¦  «        |j        z  dk    sJ ‚t          j                             d¦  «        }	t          j         	                    ||	¬¦  «        }
t          ˆˆfd„‰                     ||¬¦  «        ‰                     ||¬¦  «        ||d|
¬¦  «        }|j        j        ‰j        k    rdnd}t#          ‰                     ‰                     |j        ¦  «        ¦  «        |j        |¬¦  «         ‰                     |j        dz   ¦  «        t          j        |¦  «        dk     r‰j        ndk    sJ ‚t#          ‰                     |j        ¦  «        |j        ||dz  ¬¦  «         d S )Nr   c                 óz   •— ‰                      d‰                     | | z  d¬¦  «        z  ¦  «        d‰j        z  z  S )Nr×   r   r›   r   )rØ   rß   rÚ   r	  s    €r$   r
  z$TestQMCQuad.basic_test.<locals>.func«  s9   ø€ à—6’6˜$ §¢¨¨!©°! Ñ!4Ô!4Ñ4Ñ5Ô5¸¸R¼U¹ÑCÐCr&   l   é^^@E9®? )Úseedr=   )r  r  r  )Úlower_limitr   g×£p=
×ï?©rL   r   c                  ó4   •— ‰                       ‰| Ž ¦  «        S r2   r  )Úargsr
  rB   s    €€r$   r%   z(TestQMCQuad.basic_test.<locals>.<lambda>¼  s   ø€ ¨¯ª¨t¨t°T¨{Ñ(;Ô(;€ r&   T)r  r  r  r  ç›+¡†›„=g�íµ ÷ÆÀ>r)   y                éd   rK   )r   r4   ÚonesÚzerosÚeyeÚrandomÚdefault_rngr   r  r  r   r?   Úmultivariate_normalÚcdfr   ÚstdtritÚstandard_errorr   r‘   Úprodr>   r@   ÚrealrØ   ÚimagrÚ   )r+   r  r  ÚsignsrB   r>   r  ÚmeanÚcovÚrngr  rŒ   r  rÿ   ÚrefrL   Úlogresr*   r
  s       `             @r$   Ú
basic_testzTestQMCQuad.basic_test£  sÄ  øø€ Ý  Ñ$Ô$ˆØˆ=Ý”G˜A‘J”JˆEØˆÝŒx˜‰~Œ~ˆÝŒf�T‰lŒlˆð	Dð 	Dð 	Dð 	Dð 	Dõ Œi×#Ò#Ð$7Ñ8Ô8ˆÝŒy�Š˜t¨#ˆÑ.Ô.ˆÝŒH�T‰NŒNˆÝŒG�D‰MŒM˜EÑ!ˆÝ�t˜RŸZšZ¨°˜ZÑ7Ô7¸¿ºÀAÈU¸Ñ9SÔ9SØ (°kÈðNñ Nô NˆåÔ'×+Ò+¨A¨t°SÀaÐ+ÑHÔHˆÝŒ˜{¨1™}¨eÑ4Ô4°sÔ7IÑIˆÝ˜œ b§j¢j°¸E jÑ&BÔ&BÈÐNÑNÔNÐNÝŒw�u‰~Œ~˜cœlÑ*¨QÒ.Ð.Ð.Ð.åŒi×#Ò#Ð$7Ñ8Ô8ˆÝŒy�Š˜t¨#ˆÑ.Ô.ˆÝÐ;Ð;Ð;Ð;Ð;ØŸ*š* Q¨e˜*Ñ4Ô4°b·j²jÀÈ%°jÑ6PÔ6PØ#+¸Ø"¨ð/ñ /ô /ˆð œÔ*¨b¬jÒ8Ð8ˆuˆu¸dˆÝ˜Ÿš §¢ v¤Ñ 7Ô 7Ñ8Ô8¸#¼,ÈTÐRÑRÔRÐRØ�wŠw�v”¨Ñ+Ñ,Ô,½"¼'À%¹.¼.È1Ò:LÐ:L°´°ÐRSÒTÐTÐTÐTÝ˜Ÿš˜vÔ4Ñ5Ô5ØÔ*°¸DÀ¹Hð	Fñ 	Fô 	Fð 	Fð 	Fð 	Fr&   r  i   r  rx   c                 ó6   — |                       |||¬¦  «         d S )N)rB   ©r:  )r+   r  r  rB   s       r$   Ú
test_basiczTestQMCQuad.test_basicÆ  s"   € ð 	�Š˜ +°"ˆÑ5Ô5Ð5Ð5Ð5r&   r4  rX   g      ð¿c                 ó4   — |                       ||¬¦  «         d S )N)r4  rB   r<  )r+   r4  rB   s      r$   Ú	test_signzTestQMCQuad.test_signË  s   € à�Š˜e¨ˆÑ+Ô+Ð+Ð+Ð+r&   r  FTc           	      ó4  — d}t          j        t          |¬¦  «        5  t          d„ |                     ddg¦  «        |                     ddg¦  «        |¬¦  «        }d d d ¦  «         n# 1 swxY w Y   |j        |r|j         ndk    sJ ‚|j        dk    sJ ‚d S )Nz-A lower limit was equal to an upper limit, sorÊ   c                 ó   — dS r  r!   rJ   s    r$   r%   z'TestQMCQuad.test_zero.<locals>.<lambda>Ó  s   €  Q€ r&   r   r   r  )rž   r  r  r   r?   r‘   Úinfr0  )r+   r  rB   r  rÿ   s        r$   Ú	test_zerozTestQMCQuad.test_zeroÏ  sõ   € àAˆÝŒ\�+¨WÐ5Ñ5Ô5ð 	Yð 	YÝ˜;˜;¨¯
ª
°A°q°6Ñ(:Ô(:¸B¿JºJÈÈ1ÀvÑ<NÔ<NÐTWÐXÑXÔXˆCð	Yð 	Yð 	Yñ 	Yô 	Yð 	Yð 	Yð 	Yð 	Yð 	Yð 	Yøøøð 	Yð 	Yð 	Yð 	YàŒ|¨3Ð 5 ¤  °AÒ6Ð6Ð6Ð6ØÔ! QÒ&Ð&Ð&Ð&Ð&Ð&s   ž?A)Á)A-Á0A-c                 óà   — d„ }t          |dd¦  «        }t          j                             dd¬¦  «        t          j                             dd¬¦  «        z
  }t	          |j        |d¦  «         d S )Nc                 óD   — t           j                             | d¬¦  «        S )Nr   ©Úscale)r   ÚnormÚpdfrJ   s    r$   r
  z-TestQMCQuad.test_flexible_input.<locals>.funcÚ  s   € Ý”:—>’> !¨1�>Ñ-Ô-Ð-r&   r   r   r   rF  g{®Gáz„?)r   r   rH  r.  r   r‘   )r+   r
  rÿ   r8  s       r$   Útest_flexible_inputzTestQMCQuad.test_flexible_input×  sl   € ð	.ð 	.ð 	.õ �t˜Q Ñ"Ô"ˆÝŒj�nŠn˜Q aˆnÑ(Ô(­5¬:¯>ª>¸!À1¨>Ñ+EÔ+EÑEˆÝ˜œ c¨4Ñ0Ô0Ð0Ð0Ð0r&   )r  r   NN)r7   r8   r9   r  r:  rž   rŸ   r    r=  r?  rC  rJ  r!   r&   r$   r  r  v  s  € € € € € ð)0ð )0ð )0ðV!Fð !Fð !Fð !FðF „[×Ò˜Z¨$°¨Ñ7Ô7Ø„[×Ò˜]¨Q°¨GÑ4Ô4ð6ð 6ñ 5Ô4ñ 8Ô7ð6ð „[×Ò˜W¨¨B x°#°s°¸cÀ2¸YÈÈSÈ	Ð&RÑSÔSð,ð ,ñ TÔSð,ð „[×Ò˜U U¨D MÑ2Ô2ð'ð 'ñ 3Ô2ð'ð1ð 1ð 1ð 1ð 1r&   r  r~   ©r"   rq   r§   r{   c                óf  ‡	— | j         |         dk     r-|€t          | |||d ¬¦  «        S |t          | |||d¬¦  «        z   S t          j        | |d¦  «        } t          j        |¦  «        dk    rt          j        ||d¦  «        n|}t          j        |¦  «        dk    rt          j        ||d¦  «        n|}t          j        |¦  «        dk    rt          j        ||d¦  «        n|}| j         d         Š	|�|t          j        ‰	¦  «        z  n|}|d u }|rdn|}t          j        || j         ¦  «        }t          j        || j         d d…         dz   ¦  «        }t          j        | ||fd¬¦  «        }ˆ	fd„}t          j        |d|¦  «        }|r|d	dd …f         n|}t          j        |d|¦  «        }|S )
NrV   )r"   rq   r{   r§   r   r~   r   )r   r›   c                 ód   •— t          | d ‰…         | ‰d‰z  …         | d‰z  d …         ¬¦  «        S )Nr   )r"   r§   )r   )rà   r#   s    €r$   ru   z*cumulative_simpson_nd_reference.<locals>.fý  s8   ø€ Ý! ! B Q B¤%¨1¨Q¨q°©s¨U¬8¸Q¸qÀ¹s¸t¸t¼WÐEÑEÔEÐEr&   .)	r»   r	   r4   Úmoveaxisr  r5   Úbroadcast_toÚconcatenateÚapply_along_axis)
rO   r"   rq   r§   r{   Úinitial_was_nonerà   ru   rÿ   r#   s
            @r$   Úcumulative_simpson_nd_referencerS  â  sØ  ø€ à„wˆt„}�qÒÐØˆ?Ý'¨¨Q°2¸DÈ$ÐOÑOÔOÐOàÕ1°!°q¸RÀdÐTUÐVÑVÔVÑVÐVõ 	Œ�A�t˜RÑ Ô €AÝ$&¤G¨A¡J¤J°¢N N�Œ�A�t˜RÑ Ô Ð ¸€AÝ&(¤g¨b¡k¤k°A¢o o�Œ�R˜˜rÑ	"Ô	"Ð	"¸2€BÝ02´¸Ñ0@Ô0@À1Ò0DÐ0D�bŒk˜' 4¨Ñ,Ô,Ð,È'€Gð 	
Œ�Œ€AØ˜^ˆ�RŒY�q‰\Œ\ÑÐ°€Aà $�ÐØ#Ð0ˆaˆa¨€Gõ 	Œ˜˜1œ7Ñ#Ô#€AÝŒo˜g q¤w¨s°¨s¤|°dÑ':Ñ;Ô;€GÝ
Œ˜˜1˜g�¨RÐ0Ñ0Ô0€AðFð Fð Fð Fð Få
Ô
˜a  QÑ
'Ô
'€Cð +Ð
3ˆ#ˆc�1�2�2ˆgŒ,ˆ,°€CÝ
Œ+�c˜2˜tÑ
$Ô
$€CØ€Jr&   c                   ó>  — e Zd Z ej        d¦  «        Zedz  Zej         	                    dd¦  «        ej         	                    dd¦  «        d„ ¦   «         ¦   «         Z
 ed¬¦  «        ej         	                    d	 ej        d
d¦  «        ¦  «        ej         	                    dd¦  «        ej         	                    dd¦  «        ej         	                    dd¦  «        ej         	                    dd¦  «        d„ ¦   «         ¦   «         ¦   «         ¦   «         ¦   «         ¦   «         Zej         	                    dd eg d¢¬¦  «        fd eeg d¢geeg¬¦  «        fd eeeegeeegd¬¦  «        fd eg g ¬¦  «        fd ed¬¦  «        fd  e ej        d!¦  «        ¬¦  «        fd" e ej        d!¦  «        ¬#¦  «        fd$ ed% ej        d!¦  «        ¬&¦  «        fg¦  «        d'„ ¦   «         Zd(„ Zd)„ Zej                             d*¦  «        ej        j         e ej        ej         ej        ddd*¬+¦  «         ej        d,d*d-¬.¦  «                             d/„ ¦  «        ¬0¦  «        ¬1¦  «        d2„ ¦   «         ¦   «         ¦   «         Zej                             d*¦  «        ej        j         e ej        ej         ej        ddd*¬+¦  «         ej        d,d*d-¬.¦  «                             d3„ ¦  «        ¬0¦  «        ¬1¦  «        d4„ ¦   «         ¦   «         ¦   «         Zd%S )5ÚTestCumulativeSimpsonr   r   Úuse_dx)FTÚuse_initialc                 ó(  — t           j                             d¦  «        }d}|rdnd}|                     |                     ¦   «         ¦  «        }|dk    r;|                     t          j        |                     |¦  «        ¦  «        ¦  «        }nG|                     ||j        ¬¦  «        |z  |                     |                     ¦   «         ¦  «        z   }|                     |dz   |j        ¬¦  «        d d …|j        f         }	|                     |                     |dz   ¦  «        ¦  «        d d …|j        f         }
|                     |
||	z  z  d¬¦  «        }|                     |
||	dz   z  z  |	dz   z  d¬¦  «        }|r|n||d         z
  dd …         }|r|d         nd }|rd	|ind
|i}t          |fi |¤d|i¤Ž}|st          ||d¬¦  «         d S |rdnd}t          ||d¬¦  «         t          ||d d…         ||d d…         d¬¦  «         d S )Nì   F»*ŒBG r˜   rV   r   r=   r   r   r›   rq   r"   r§   gVçž¯â<r)   g{®Gázd?)r4   r+  r,  r?   Úsortr5   r@   Únewaxisrß   r   r   )r+   rV  rW  rB   r7  r#   Úorderrq   r"   r�   ÚcrO   ÚYr8  r§   Úkwargrÿ   Úi0s                     r$   r­   zTestCumulativeSimpson.test_1d  s.  € õ
 Œi×#Ò#Ð$8Ñ9Ô9ˆØˆð Ð"�� ˆØ�ZŠZ˜Ÿ
š
™œÑ%Ô%ˆØ�AŠ:ˆ:Ø—
’
�2œ7 3§:¢:¨a¡=¤=Ñ1Ô1Ñ2Ô2ˆAˆAà—	’	˜! 2¤:�	Ñ.Ô.¨rÑ1°B·J²J¸s¿zºz¹|¼|Ñ4LÔ4LÑLˆAØ�IŠI�e˜a‘i r¤zˆIÑ2Ô2°1°1°1°b´j°=ÔAˆØ�JŠJ�s—z’z %¨!¡)Ñ,Ô,Ñ-Ô-¨a¨a¨a°´¨mÔ<ˆØ�FŠF�1�Q˜‘T‘6 ˆFÑ"Ô"ˆØ�FŠF�1�Q˜˜Q™‘Z‘<  Q¡Ñ'¨aˆFÑ0Ô0ˆØÐ0ˆaˆa Q q¨¤t¡V¨Q¨R¨R¤Lˆð &Ð/�!�A”$�$¨4ˆØ$Ð2��r�
�
¨3°¨(ˆÝ  Ð=Ð= eÐ=Ð=°WÐ=Ð=Ð=ˆð ð 	@Ý˜C ¨5Ð1Ñ1Ô1Ð1Ð1Ð1à!Ð(�� qˆBå˜C ¨6Ð2Ñ2Ô2Ð2å˜C   A œJ¨¨B¨E°¨E¬
¸Ð?Ñ?Ô?Ð?Ð?Ð?r&   T)Úcpu_onlyr{   éýÿÿÿrV   Úx_ndim)r   rV   Úx_len)r   r   r`   Úi_ndim)Nr   rV   rq   )NTc                 ó(  — t           j                             d¦  «        }dd|g}|d         ||         c||<   |d<   |                     ¦   «         }	d|	|<   |dk    r|	nd}
|                     |                     |¬¦  «        ¦  «        }d	\  }}|rF|dk    r|                     |	¬¦  «        n|                     ¦   «         }|                     |¦  «        }ns|dk    r*t          j        |                     |¬¦  «        |¬
¦  «        n-t          j        |                     ||         ¬¦  «        ¦  «        }|                     |¦  «        }|€d n(|                     |                     |
¬¦  «        ¦  «        }t          |||||¬¦  «        }t          t          j        |¦  «        t          j        |¦  «        |€d nt          j        |¦  «        |€d nt          j        |¦  «        |¬¦  «        }t          ||                     |¦  «        d¬¦  «         d S )NrY  r[   é   r~   r   rV   r!   ©Úsize)NNr›   rK  çVçž¯Ò<r)   )	r4   r+  r,  Úcopyr?   rZ  r   rS  r   )r+   r{   rc  rd  re  rq   rB   r7  r»   Úshape_len_1Úi_shaperO   r"   r§   rÿ   r8  s                   r$   Útest_ndzTestCumulativeSimpson.test_nd1  sñ  € õ Œi×#Ò#Ð$8Ñ9Ô9ˆð �A�u�ˆØ!& r¤¨E°$¬KÐˆˆd‰�U˜2‘YØ—j’j‘l”lˆØˆ�DÑØ!'¨1¢ �+�+°"ˆð �JŠJ�s—z’z u�zÑ-Ô-Ñ.Ô.ˆØ‰ˆˆ2Øð 	Ø17¸!²°�—’ �Ñ-Ô-Ð-ÀÇÂÁÄˆBØ—’˜B‘”ˆBˆBà?EÈºz¸z•”˜Ÿš¨˜Ñ/Ô/°dÐ;Ñ;Ô;Ð;Ý”g˜cŸjšj¨e°D¬k˜jÑ:Ô:Ñ;Ô;ð à—
’
˜1‘”ˆAØ ˜.�$�$¨b¯jªj¸¿ºÈ¸Ñ9QÔ9QÑ.RÔ.Rˆõ !  a¨B¸ÀdÐKÑKÔKˆõ .ÝŒJ�q‰MŒM�RœZ¨™]œ]°r°z¨t¨tÅrÄzÐRTÁ~Ä~Ø#˜O�D�Dµ´¸GÑ1DÔ1DÈ4ð
ñ 
ô 
ˆõ 	˜˜RŸZšZ¨™_œ_°5Ð9Ñ9Ô9Ð9Ð9Ð9r&   )r  Úkwarg_updatezx must be strictly increasing)r   r   rV   r   rJ   )r   r   r   r   )r"   rO   r   )r"   rO   r{   rÐ   z-`axis=4` is not valid for `y` with `y.ndim=1`r›   z+shape of `x` must be the same as `y` or 1-Dr[   z'`initial` must either be a scalar or...r¦   z"`dx` must either be a scalar or...N)r"   rq   c           
      óf  ‡— t          ‰                     | j        ¦  «        ‰                     | j        ¦  «        d d d¬¦  «        }ˆfd„|                     ¦   «         D ¦   «         }t          j        t          |¬¦  «        5  t          di t          |fi |¤Ž¤Ž d d d ¦  «         d S # 1 swxY w Y   d S )Nr~   )rO   r"   rq   r§   r{   c                 ó”   •— i | ]D\  }}|t          |t          ¦  «        r'‰                     t          j        |¦  «        ¦  «        n|“ŒES r!   )Ú
isinstancer‡   r?   r4   )rú   ÚkÚvrB   s      €r$   ú
<dictcomp>zATestCumulativeSimpson.test_simpson_exceptions.<locals>.<dictcomp>e  sX   ø€ ð :ð :ð :Ù ˜A˜qð ½
À1ÅdÑ8KÔ8KÐR˜2Ÿ:š:¥b¤j°¡m¤mÑ4Ô4Ð4ÐQRð :ð :ð :r&   rÊ   r!   )	Údictr?   Úy0Úx0Úitemsrž   rÌ   rÍ   r   )r+   r  ro  rB   Úkwargs0s      ` r$   Útest_simpson_exceptionsz-TestCumulativeSimpson.test_simpson_exceptionsX  s%  ø€ õ ˜Ÿš D¤GÑ,Ô,°·
²
¸4¼7Ñ0CÔ0CÈØ#¨"ð.ñ .ô .ˆð:ð :ð :ð :Ø$0×$6Ò$6Ñ$8Ô$8ð:ñ :ô :ˆåŒ]�:¨WÐ5Ñ5Ô5ð 	@ð 	@ÝÐ?Ð?¥ gÐ!>Ð!>°Ð!>Ð!>Ð?Ð?Ð?ð	@ð 	@ð 	@ñ 	@ô 	@ð 	@ð 	@ð 	@ð 	@ð 	@ð 	@ð 	@øøøð 	@ð 	@ð 	@ð 	@ð 	@ð 	@s   ÂB&Â&B*Â-B*c                 ó  — t           j                             d¦  «        }|                     |                     d¬¦  «        ¦  «        }t	          |d¬¦  «        }t          ||                     d|j        ¬¦  «        ¦  «         d S )NrY  r˜   rh  re   rp   rí   r=   )r4   r+  r,  r?   r   r   r)  r@   )r+   rB   r7  rO   rÿ   s        r$   Útest_special_casesz(TestCumulativeSimpson.test_special_casesj  sq   € åŒi×#Ò#Ð$8Ñ9Ô9ˆØ�JŠJ�s—z’z r�zÑ*Ô*Ñ+Ô+ˆÝ   rÐ*Ñ*Ô*ˆÝ˜˜RŸXšX a¨r¬z˜XÑ:Ô:Ñ;Ô;Ð;Ð;Ð;r&   c                 ó˜  — t          j        |d¬¦  «        }t          ||¦  «        }t          |dddd…f         |dddd…f         ¦  «        dddd…f         }g |j        dd…         ¢d‘R }t          j        t          j        |¦  «        |ddd…f         |ddd…f         z
  t          j        |¦  «        gd¬¦  «        }d|dddd…f<   |S )aj  `cumulative_simpson` and `simpson` can be tested against other to verify
        they give consistent results. `simpson` will iteratively be called with
        successively higher upper limits of integration. This function calculates
        the theoretical correction required to `simpson` at even intervals to match
        with `cumulative_simpson`.
        r~   r›   .Nr   re   r   )r4   Údiffr   r»   rP  r)  )r+   rO   r"   ÚdÚsub_integrals_h1Úsub_integrals_h2Úzeros_shapeÚtheoretical_differences           r$   Ú1_get_theoretical_diff_between_simps_and_cum_simpszGTestCumulativeSimpson._get_theoretical_diff_between_simps_and_cum_simpsu  s  € õ ŒG�A˜BÐÑÔˆÝ@ÀÀAÑFÔFÐÝ@Øˆc�4�4�R�4ˆiŒL˜!˜C   2 ˜Iœ,ñ
ô 
à
ˆtˆt�ˆtˆ)ôÐð
 )˜œ   œÐ( aÐ(Ð(ˆÝ!#¤å”˜Ñ%Ô%Ø! # q r r 'Ô*Ð-=¸cÀ3ÀBÀ3¸hÔ-GÑGÝ”˜Ñ%Ô%ðð
 ð"
ñ "
ô "
Ðð -0Ð˜s A D q D˜yÑ)ð &Ð%r&   r˜   )Úmax_dimsÚmin_sideÚmax_siderÖ   F)Ú	allow_nanc                 ó(   — t          | ¦  «        dk    S ©NgH¯¼šò×z>©ÚabsrJ   s    r$   r%   zTestCumulativeSimpson.<lambda>™  ó   € Í#ÈaÉ&Ì&ÐSWÊ-€ r&   )ÚelementsrÑ   c           
      ó^  — d„ }t          |                     |¦  «        d¬¦  «        } ||¦  «        }|                      |t          j        |j        d         ¦  «        ¬¦  «        }t          |ddd…f         |                     |ddd…f         |ddd…f         z   ¦  «        d	¬
¦  «         dS )áe  Theoretically, the output of `cumulative_simpson` will be identical
        to `simpson` at all even indices and in the last index. The first index
        will not match as `simpson` uses the trapezoidal rule when there are only two
        data points. Odd indices after the first index are shown to match with
        a mathematically-derived correction.c                 ó€   ‡ — t          j        ˆ fd„t          d‰ j        d         dz   ¦  «        D ¦   «         d¬¦  «        S )Nc                 óF   •— g | ]}t          ‰d d|…f         d¬¦  «        ‘ŒS ).NrX   rp   ©r   )rú   r�   rO   s     €r$   rü   z|TestCumulativeSimpson.test_cumulative_simpson_against_simpson_with_default_dx.<locals>.simpson_reference.<locals>.<listcomp>¦  s2   ø€ ÐNÐNÐN°•˜˜3   ˜7œ¨Ð,Ñ,Ô,ÐNÐNÐNr&   r   r~   r   r›   ©r4   Ústackr…   r»   rÑ   s   `r$   Úsimpson_referencezhTestCumulativeSimpson.test_cumulative_simpson_against_simpson_with_default_dx.<locals>.simpson_reference¤  sJ   ø€ Ý”8ØNÐNÐNÐNµe¸A¸q¼wÀr¼{È1¹}Ñ6MÔ6MÐNÑNÔNÐUWðñ ô ð r&   rX   rp   r~   rJ   .r   Ng¼‰Ø—²Òœ<r#  )r   r?   r…  r4   r5   r»   r   )r+   rO   rB   r—  rÿ   r8  r„  s          r$   Ú7test_cumulative_simpson_against_simpson_with_default_dxzMTestCumulativeSimpson.test_cumulative_simpson_against_simpson_with_default_dx“  sÕ   € ð"	ð 	ð 	õ
 ! §¢¨A¡¤°3Ð7Ñ7Ô7ˆØÐ Ñ"Ô"ˆØ!%×!WÒ!WØ•”˜1œ7 2œ;Ñ'Ô'ð "Xñ "
ô "
Ðõ 	Ø��Q�R�R�ŒL˜"Ÿ*š* S¨¨a¨b¨b¨¤\Ð4JÈ3ÐPQÐPRÐPRÈ7Ô4SÑ%SÑTÔTØð	
ñ 	
ô 	
ð 	
ð 	
ð 	
r&   c                 ó(   — t          | ¦  «        dk    S r‹  rŒ  rJ   s    r$   r%   zTestCumulativeSimpson.<lambda>¹  rŽ  r&   c           
      ó>  — d|j         d         dz
  z  }t          j        dd|j         d         ¬¦  «        }|dd…         d|z  t          j                             ddt          |¦  «        dz
  ¦  «        z  z   |dd…<   d„ }t          |                     |¦  «        |                     |¦  «        ¬	¦  «        } |||¦  «        }|                      ||¦  «        }t          |d
dd…f         |                     |d
dd…f         |d
dd…f         z   ¦  «        ¦  «         dS )r‘  r˜   r~   r   r   r¤   NrF   c                 ó„   ‡ ‡— t          j        ˆˆ fd„t          d‰ j        d         dz   ¦  «        D ¦   «         d¬¦  «        S )Nc           	      óZ   •— g | ]'}t          ‰d d|…f         ‰d d|…f         ¬¦  «        ‘Œ(S ).NrJ   r”  )rú   r�   r"   rO   s     €€r$   rü   zlTestCumulativeSimpson.test_cumulative_simpson_against_simpson.<locals>.simpson_reference.<locals>.<listcomp>Ê  s>   ø€ ÐTÐTÐT°q•˜˜3   ˜7œ q¨¨b¨q¨b¨¤zÐ2Ñ2Ô2ÐTÐTÐTr&   r   r~   r   r›   r•  )rO   r"   s   ``r$   r—  zXTestCumulativeSimpson.test_cumulative_simpson_against_simpson.<locals>.simpson_referenceÈ  sO   øø€ Ý”8ØTÐTÐTÐTÐT½EÀ!ÀQÄWÈRÄ[ÐQRÁ]Ñ<SÔ<SÐTÑTÔTØðñ ô ð r&   rJ   .)
r»   r4   rƒ   r+  ÚuniformÚlenr   r?   r…  r   )	r+   rO   rB   Úintervalr"   r—  rÿ   r8  r„  s	            r$   Ú'test_cumulative_simpson_against_simpsonz=TestCumulativeSimpson.test_cumulative_simpson_against_simpson³  s6  € ð" �q”w˜r”{ Q‘Ñ'ˆÝŒK˜˜2 1¤7¨2¤;Ð/Ñ/Ô/ˆØ�!�"�"”˜˜H™¥R¤Y×%6Ò%6°r¸1½cÀ!¹f¼fÀq¹jÑ%IÔ%IÑIÑIˆˆ!ˆ"ˆ"‰ð	ð 	ð 	õ ! §¢¨A¡¤°"·*²*¸Q±-´-Ð@Ñ@Ô@ˆØÐ  1Ñ%Ô%ˆØ!%×!WÒ!WØˆqñ"
ô "
Ðõ 	Ø��Q�R�R�ŒL˜"Ÿ*š* S¨¨a¨b¨b¨¤\Ð4JÈ3ÐPQÐPRÐPRÈ7Ô4SÑ%SÑTÔTñ	
ô 	
ð 	
ð 	
ð 	
r&   )r7   r8   r9   r4   r5   rx  rw  rž   rŸ   r    r­   Úskip_xp_backendsrn  rv  r{  r}  r…  Ú	fail_slowÚslowr   Úhyp_numÚarraysr@   Úarray_shapesÚstÚfloatsÚfilterr˜  r   r!   r&   r$   rU  rU    sï  € € € € € à	ˆŒ�1‰Œ€BØ	ˆQ‰€Bà„[×Ò˜X }Ñ5Ô5Ø„[×Ò˜]¨MÑ:Ô:ð!@ð !@ñ ;Ô:ñ 6Ô5ð!@ðF Ð˜tÐ$Ñ$Ô$Ø„[×Ò˜V Y R¤Y¨r°1Ñ%5Ô%5Ñ6Ô6Ø„[×Ò˜X vÑ.Ô.Ø„[×Ò˜W iÑ0Ô0Ø„[×Ò˜X }Ñ5Ô5Ø„[×Ò˜T <Ñ0Ô0ð:ð :ñ 1Ô0ñ 6Ô5ñ 1Ô0ñ /Ô.ñ 7Ô6ñ %Ô$ð:ðB „[×ÒÐ8Ø	(¨$¨$°°°Ð*>Ñ*>Ô*>Ð?Ø	(¨$¨$°"°l°l°lÐ1CÈÈBÀxÐ*PÑ*PÔ*PÐQØ	(¨$¨$°"°b¸"°À"ÀbÈ"ÀÐTUÐ*VÑ*VÔ*VÐWØ	)¨4¨4°"¸Ð+;Ñ+;Ô+;Ð<Ø	8¸$¸$ÀA¸,¹,¼,ÐGØ	6¸¸¸y¸r¼yÈ¹|¼|Ð8LÑ8LÔ8LÐMØ	2°D°DÀÀÄÈ1ÁÄÐ4NÑ4NÔ4NÐOØ	-¨t¨t°d¸y¸r¼yÈ¹|¼|Ð/LÑ/LÔ/LÐMð	;ñ 	ô 	ð@ð @ñ	ô 	ð@ð<ð <ð <ð&ð &ð &ð< „[×Ò˜2ÑÔØ„[ÔØ
€UØ
ˆ'Œ.ØŒJØ ˆGÔ ¨!°aÀ"ÐEÑEÔEØ�R”Y˜s B°%Ð8Ñ8Ô8×?Ò?Ð@WÐ@WÑXÔXð
ñ 
ô 
ðñ ô ð
ð 
ñô ñ Ôñ Ôð
ð. „[×Ò˜2ÑÔØ„[ÔØ
€UØ
ˆ'Œ.ØŒJØ ˆGÔ ¨!°aÀ"ÐEÑEÔEØ�R”Y˜s B°%Ð8Ñ8Ô8×?Ò?Ð@WÐ@WÑXÔXð
ñ 
ô 
ðñ ô ð
ð 
ñô ñ Ôñ Ôð
ð 
ð 
r&   rU  c                   ó¸   — e Zd Zd„ Zd„ Zej                             d e e	ddd¦  «        ¦  «         e e	ddd	¦  «        ¦  «        z   ¦  «        d
„ ¦   «         Z
dS )ÚTestLebedevc                 óš   — d}t          j        t          |¬¦  «        5  t          j        d¦  «         d d d ¦  «         d S # 1 swxY w Y   d S )NzOrder n=-1 not available...rÊ   r~   )rž   rÌ   ÚNotImplementedErrorr   Úlebedev_rule)r+   r  s     r$   r  z!TestLebedev.test_input_validationÚ  s’   € à/ˆÝŒ]Õ.°gÐ>Ñ>Ô>ð 	'ð 	'ÝÔ" 2Ñ&Ô&Ð&ð	'ð 	'ð 	'ñ 	'ô 	'ð 	'ð 	'ð 	'ð 	'ð 	'ð 	'ð 	'øøøð 	'ð 	'ð 	'ð 	'ð 	'ð 	's   žA Á AÁAc                 óÜ   — d„ }t          j        d¦  «        \  }}| ||¦  «        z  }d}t          ||d¬¦  «         t          t          j        |¦  «        dt          j        z  ¦  «         d S )Nc                 ó6   — t          j        | d         ¦  «        S )Nr   )r4   rØ   rJ   s    r$   ru   z&TestLebedev.test_quadrature.<locals>.fã  s   € Ý”6˜!˜Aœ$‘<”<Ðr&   é   gá9‰-@r&  r)   r   )r   r®  r   r4   rß   rÚ   )r+   ru   r"   Úwrÿ   r8  s         r$   Útest_quadraturezTestLebedev.test_quadratureà  sv   € ð	 ð 	 ð 	 õ Ô% bÑ)Ô)‰ˆˆ1Ø�!�!�A‘$”$‰hˆØˆÝ˜˜S uÐ-Ñ-Ô-Ð-Ý�œ˜q™	œ	 1¥r¤u¡9Ñ-Ô-Ð-Ð-Ð-r&   r\  rV   rs   r   é#   é„   rg  c                 óR  — t          j        |¦  «        \  }}t          j        dd¬¦  «        5  t	          j        |j        d¬¦  «        }t          |j        dd¬¦  «         d d d ¦  «         n# 1 swxY w Y   t          t          j	        |¦  «        dt          j
        z  ¦  «         d S )NÚignore)ÚdivideÚinvalidr   r›   rj  r#  r   )r   r®  r4   Úerrstater   Údirectional_statsÚTr   Úmean_resultant_lengthrß   rÚ   )r+   r\  r"   r²  rÿ   s        r$   Útest_propertieszTestLebedev.test_propertiesì  sæ   € åÔ% eÑ,Ô,‰ˆˆ1åŒ[ °(Ð;Ñ;Ô;ð 	Fð 	FÝÔ)¨!¬#°AÐ6Ñ6Ô6ˆCÝ˜CÔ5°q¸uÐEÑEÔEÐEð	Fð 	Fð 	Fñ 	Fô 	Fð 	Fð 	Fð 	Fð 	Fð 	Fð 	Føøøð 	Fð 	Fð 	Fð 	Fõ 	�œ˜q™	œ	 1¥R¤U¡7Ñ+Ô+Ð+Ð+Ð+s   ®3A-Á-A1Á4A1N)r7   r8   r9   r  r³  rž   rŸ   r    r‡   r…   r¾  r!   r&   r$   r«  r«  Ø  s•   € € € € € ð'ð 'ð 'ð
.ð 
.ð 
.ð „[×Ò˜W d d¨5¨5°°B¸©?¬?Ñ&;Ô&;¸d¸dÀ5À5ÈÈSÐRSÑCTÔCTÑ>UÔ>UÑ&UÑVÔVð,ð ,ñ WÔVð,ð ,ð ,r&   r«  )5rž   Únumpyr4   Únumpy.testingr   r   r   Ú
hypothesisr   Úhypothesis.strategiesÚ
strategiesr§  Úhypothesis.extra.numpyÚextrar¤  Úscipy.integrater   r   r	   r
   r   r   r   r   r   Úscipy.integrate._quadraturer   Úscipyr   r   r   Úscipy.conftestr   Úscipy._lib._array_apir   r   r   Úscipy._lib._array_api_no_0dr   r   rŸ   r¡  r   r;   rS   rm   r¢   rÔ   r  rS  rU  r®  r«  r!   r&   r$   ú<module>rÌ     s  ðà €€€Ø Ð Ð Ð Ø LÐ LÐ LÐ LÐ LÐ LÐ LÐ LÐ LÐ LØ Ð Ð Ð Ð Ð Ø "Ð "Ð "Ð "Ð "Ð "Ø (Ð (Ð (Ð (Ð (Ð (Ð (Ð (Ð (ð;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð ;ð NÐ MÐ MÐ MÐ MÐ Mà +Ð +Ð +Ð +Ð +Ð +Ð +Ð +Ð +Ð +Ø .Ð .Ð .Ð .Ð .Ð .Ø OÐ OÐ OÐ OÐ OÐ OÐ OÐ OÐ OÐ OØ HÐ HÐ HÐ HÐ HÐ HÐ HÐ Hà”;Ô/Ð àÐ�:ÑÔð3ð 3ð 3ð 3ð 3ñ 3ô 3ñ Ôð3ð  Ð�4ÑÔð
Xð 
Xð 
Xð 
Xð 
Xñ 
Xô 
Xñ Ôð
Xð Ð�<Ñ Ô ð&>ð &>ð &>ð &>ð &>ñ &>ô &>ñ !Ô ð&>ðR Ð�7ÑÔðY'ð Y'ð Y'ð Y'ð Y'ñ Y'ô Y'ñ ÔðY'ðx ÐÐ'Ñ(Ô(ðV3ð V3ð V3ð V3ð V3ñ V3ô V3ñ )Ô(ðV3ðr Ð�9ÑÔð\%ð \%ð \%ð \%ð \%ñ \%ô \%ñ Ôð\%ð~ Ð�8ÑÔðh1ð h1ð h1ð h1ð h1ñ h1ô h1ñ Ôðh1ðV -1°TÀ4Èbð "ð "ð "ð "ð "ðJ ÐÐ%Ñ&Ô&ðM
ð M
ð M
ð M
ð M
ñ M
ô M
ñ 'Ô&ðM
ð` Ð�9Ô)Ñ*Ô*ð,ð ,ð ,ð ,ð ,ñ ,ô ,ñ +Ô*ð,ð ,ð ,r&   