
    ^Mhk                        d dl Z d dlmZmZ d dlmZmZmZmZ d dl	m
Z
mZmZ d dlmZmZ 	 d dlZn# e$ r  ed          ZY nw xY w	 d dlZn# e$ r  ed          ZY nw xY w eed          d	             Ze j        j         eed           eed
          e j                            d          d                                                 Ze j        j         eed           eed
          d                                     Z eed          d             Ze j        j         eed          d                         ZdS )    N)MissingModulecheck_version)ArgIntArgmp_assert_allcloseassert_mpmath_equal)	compute_gcompute_alpha	compute_d)gammainc	gammainccsympympmathz0.19c            
         t          j        d          5  t          j        d          t          j        d          dz  t          j        d          dz  t          j        d           dz  t          j        d           dz  t          j        d	          d
z  t          j        d          dz  g} t          t	          d          |            d d d            d S # 1 swxY w Y   d S )N         i      i  i;  i % i' i `uicP l     
F    )mpworkdpsmpfr   r	   )gs    l/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/scipy/special/tests/test_precompute_gammainc.pytest_gr      s     
B , ,VAYYq		"bfQiimfSkk\%"&++g!5VF^^I%rvg{'BD 	9Q<<+++	, , , , , , , , , , , , , , , , , ,s   B9CC"Cz0.7zrtol only 2e-11, see gh-6938c                  "   t          j        d          5  t          j        d          t          j        d          t          j        d          dz  t          j        d          dz  t          j        d           dz  t          j        d          dz  t          j        d          dz  t          j        d	           d
z  t          j        d          dz  g	} t          t	          d          |            d d d            d S # 1 swxY w Y   d S )Nr   r   r      $   i  i  irB  r   iS iX 	   )r   r   r   r   r
   )alphas    r   
test_alphar!      s    
B 4 4BF1IIrvayy{BF1IIbL&))C45&++g%rvayy'79 	=++U333	4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4s   C"DDDc                  ^   t          j        d          5  ddt          j        d           fddt          j        d          fddt          j        d           fddt          j        d          fd	dt          j        d
          fd	dt          j        d           fddt          j        d          fddt          j        d           fddt          j        d           fddt          j        d          fddt          j        d           fddt          j        d          fddt          j        d          fddt          j        d           fddt          j        d          fddt          j        d           fddt          j        d           fddt          j        d          fddt          j        d           fddt          j        d           fg} t          d!d"          fd#| D             }d$ | D             }t	          ||           d d d            d S # 1 swxY w Y   d S )%N2   r   z 0.333333333333333333333333333333r   z#0.102618097842403080425739573227e-7r   z#0.185185185185185185185185185185e-2z#0.119516285997781473243076536700e-7   z#0.413359788359788359788359788360e-2z#0.140925299108675210532930244154e-7r   z#0.649434156378600823045267489712e-3z#0.191111684859736540606728140873e-7   z#0.861888290916711698604702719929e-3z#0.288658297427087836297341274604e-7   z#0.336798553366358150308767592718e-3z#0.482409670378941807563762631739e-7   z#0.531307936463992223165748542978e-3z#0.882860074633048352505085243179e-7r   z#0.344367606892377671254279625109e-3z#0.175629733590604619378669693914e-6   z#0.652623918595309418922034919727e-3z#0.377358774161109793380344937299e-6r   z#0.596761290192746250124390067179e-3z#0.870823417786464116761231237189e-6
      c                 4    g | ]\  }}}|         |         S  r,   ).0knstdds       r   
<listcomp>ztest_d.<locals>.<listcomp>N   s'    00091aqtAw000    c                     g | ]
}|d          S )r$   r,   )r-   xs     r   r2   ztest_d.<locals>.<listcomp>O   s    %%%qt%%%r3   )r   r   r   r   r   )datasetresr0   r1   s      @r   test_dr8   ,   s    
B % %q26"DEEEFr26"GHHIq26"GHHHIr26"GHHIq"&!FGGHrBF#HIIIJq"&!FGGHrBF#HIIIJq26"GHHHIr26"GHHIq26"GHHHIr26"GHHIq"&!FGGHrBF#HIIIJq"&!FGGHrBF#HIIIJq26"GHHHIr26"GHHIq26"GHHHIr26"GHHI'K( b"0000000%%W%%%3$$$1% % % % % % % % % % % % % % % % % %s   G?H""H&)H&c            	      ~    t          t          d t          ddd          t          dd          gdddd           d S )	Nc                 0    t          j        | |d          S )NT)bregularizedr   r   ar5   s     r   <lambda>ztest_gammainc.<locals>.<lambda>Y       R[aT%J%J%J r3   r   d   F)inclusive_aFFg<r#   nan_okrtolr/   dps)r   r   r   r,   r3   r   test_gammaincrI   S   s[    
 JJQ777QE$5BB@ @ @ @ @ @r3   c            	          t          d d t          dd          t          dd          gdddd	           t          t          d
 t          dd          t          dd          gdddd	           d S )Nc                 &    t          | |d          S )N  )rH   )r   r>   s     r   r@   z test_gammaincc.<locals>.<lambda>c   s    Yq!%>%>%> r3   c                 0    t          j        | |d          S NT)r?   r<   r=   r>   s     r   r@   z test_gammaincc.<locals>.<lambda>d   rA   r3      rB   FrD   r#   rL   rE   c                 0    t          j        | |d          S rN   r=   r>   s     r   r@   z test_gammaincc.<locals>.<lambda>j   rA   r3   r   r   )r   r   r   r   r,   r3   r   test_gammainccrQ   ^   s    
 >>JJRs2s||4$5BDB B B B 	JJ3Q5$5BB@ @ @ @ @ @r3   )pytestscipy.special._testutilsr   r   scipy.special._mptestutilsr   r   r   r   &scipy.special._precompute.gammainc_asyr	   r
   r   'scipy.special._precompute.gammainc_datar   r   r   ImportErrorr   r   r   markslowxfail_on_32bitr!   xslowr8   rI   rQ   r,   r3   r   <module>r\      s    A A A A A A A A: : : : : : : : : : : :) ) ) ) ) ) ) ) ) ) G G G G G G G G#LLLL # # #M'""EEE#! ! ! !	x	 	 BBB! r6, , , r6ue:;;4 4 <;   4 r6ue!% !%   !%H r6@ @ @ r6@ @  @ @ @s!   1 AAA A A 