§
    fŠtjÆ	  ã                   ó�   — d Z ddlZddlmZ 	 ddlZn# e$ r Y nw xY wd„ Zd„ Zd„ Z	d„ Z
d„ Zd	Zd
Zd„ Zedk    r e¦   «          dS dS )zì
Precompute coefficients of Temme's asymptotic expansion for gammainc.

This takes about 8 hours to run on a 2.3 GHz Macbook Pro with 4GB ram.

Sources:
[1] NIST, "Digital Library of Mathematical Functions",
    https://dlmf.nist.gov/

é    N)Úlagrange_inversionc                 ód  — t          j        d¦  «        dz  g}t          d| ¦  «        D ]†}|d         |z  }t          dt          |¦  «        ¦  «        D ]}|||         ||          z  |dz   z  z  }Œ||d         dt          j        d¦  «        |dz   z  z   z  z  }|                     |¦  «         Œ‡|S )za_k from DLMF 5.11.6é   é   éÿÿÿÿr   )ÚmpÚsqrtÚrangeÚlenÚmpfÚappend)ÚnÚaÚkÚakÚjs        úd/var/www/html/CA-Chatbot/venv/lib/python3.11/site-packages/scipy/special/_precompute/gammainc_asy.pyÚ	compute_ar      sº   € å	Œ�‰Œ�A‰ˆ€AÝ�1�a‰[Œ[ð ð ˆØˆrŒU�1‰WˆÝ�q�#˜a™&œ&Ñ!Ô!ð 	%ð 	%ˆAØ�!�A”$�q˜!˜”u‘*˜a !™eÑ$Ñ$ˆBˆBØ
ˆa�Œd�A�œ˜q™	œ	 1 q¡5Ñ)Ñ)Ñ*Ñ*ˆØ	�Š�‰ŒˆˆØ€Hó    c                 ób   ‡— t          d| z  ¦  «        Šˆfd„t          | ¦  «        D ¦   «         }|S )zg_k from DLMF 5.11.3/5.11.5r   c                 ó|   •— g | ]8}t          j        d ¦  «        t          j        d|¦  «        z  ‰d |z           z  ‘Œ9S )r   g      à?)r   r	   Úrf)Ú.0r   r   s     €r   ú
<listcomp>zcompute_g.<locals>.<listcomp>#   s=   ø€ Ð;Ð;Ð;¨Q�Œ�‰Œ•B”E˜#˜q‘M”MÑ	! ! A a¡C¤&Ñ	(Ð;Ð;Ð;r   )r   r
   )r   Úgr   s     @r   Ú	compute_gr       s6   ø€ å�!�A‘#‰Œ€AØ;Ð;Ð;Ð;µ%¸±(´(Ð;Ñ;Ô;€AØ€Hr   c                 óÜ   — | dk    r/t          j        d| t          j        | dz   ¦  «        z
  z  ¦  «        S | dk     r0t          j        d| t          j        | dz   ¦  «        z
  z  ¦  «         S dS )z6Function from DLMF 8.12.1 shifted to be centered at 0.r   r   r   )r   r	   Úlog)Úlams    r   Úetar    '   sh   € à
ˆQ‚w€wÝŒw�q˜#¥¤ s¨Q¡w¡¤Ñ/Ñ0Ñ1Ô1Ð1Ø	ˆqŠˆÝ”˜˜3¥¤¨¨a©¡¤Ñ0Ñ1Ñ2Ô2Ð2Ð2àˆqr   c                 ó\   — t          j        t          d| dz
  ¦  «        }t          |¦  «        S )zalpha_n from DLMF 8.12.13r   r   )r   Útaylorr    r   )r   Úcoeffss     r   Úcompute_alphar$   1   s'   € åŒY•s˜A˜q 1™uÑ%Ô%€FÝ˜fÑ%Ô%Ð%r   c           	      ób  — |d| z  z   }t          j        d¦  «         dz  g}t          |dz   ¦  «        }t          d|¦  «        D ]&}|                     |dz   ||dz            z  ¦  «         Œ'|g}t          | ¦  «        }t          d| ¦  «        D ]|}g }	t          |d|z  z
  ¦  «        D ]M}|	                     d|z  ||         z  |d         |         z  |dz   ||dz
           |dz            z  z   ¦  «         ŒN|                     |	¦  «         Œ}t          | ¦  «        D ]}||         d|…         ||<   Œ|S )zd_{k, n} from DLMF 8.12.12r   r   é   r   r   N)r   r   r$   r
   r   r   )
ÚKÚNÚMÚd0Úalphar   Údr   r   Údks
             r   Ú	compute_dr.   7   sM  € à	ˆAˆa‰C‰€AÝŒ6�!‰9Œ9ˆ*�Q‰,ˆ€BÝ˜!˜a™%Ñ Ô €EÝ�1�a‰[Œ[ð &ð &ˆØ
�	Š	�1�q‘5˜%  !¡œ*Ñ$Ñ%Ô%Ð%Ð%Ø	ˆ€AÝ�!‰Œ€AÝ�1�a‰[Œ[ð ð ˆØˆÝ�q˜1˜Q™3‘w‘”ð 	Bð 	BˆAØ�IŠI�r˜A‘g˜a œd‘l 1 Q¤4¨¤7Ñ*¨a°!©e°Q°q¸±s´V¸A¸a¹C´[Ñ-@Ñ@ÑAÔAÐAÐAØ	�Š�‰ŒˆˆÝ�1‰XŒXð ð ˆØ�Œt�B�Q�BŒxˆˆ!‰ˆØ€Hr   z½/* This file was automatically generated by _precomp/gammainc.py.
 * Do not edit it manually!
 */

#ifndef IGAM_H
#define IGAM_H

#define K {}
#define N {}

static const double d[K][N] =
{{z
#endif
c                  óp  — t          t          ¦  «         d} d}t          j        d¦  «        5  t	          | |¦  «        }d d d ¦  «         n# 1 swxY w Y   t
          j                             t
          j                             t          ¦  «        ddd¦  «        }t          |dz   d¦  «        5 }|                     t                               | |¦  «        ¦  «         t          |¦  «        D ]‚\  }}d„ |D ¦   «         }|                     d	¦  «         |                     d
                     |¦  «        ¦  «         || dz
  k     r|                     d¦  «         Œm|                     d¦  «         Œƒ|                     t          ¦  «         d d d ¦  «         n# 1 swxY w Y   t          j        |dz   |¦  «         d S )Né   é2   z..Úcepheszigam.hz.newÚwc                 ó>   — g | ]}t          j        |d dd¬¦  «        ‘ŒS )é   r   )Ú	min_fixedÚ	max_fixed)r   Únstr)r   Úxs     r   r   zmain.<locals>.<listcomp>h   s+   € ÐIÐIÐIÀ•2”7˜1˜b¨A¸Ð;Ñ;Ô;ÐIÐIÐIr   Ú{z, r   z},
z}};
)ÚprintÚ__doc__r   Úworkdpsr.   ÚosÚpathÚjoinÚdirnameÚ__file__ÚopenÚwriteÚheaderÚformatÚ	enumerateÚfooterÚrename)r'   r(   r,   ÚfnÚfr   Úrows          r   ÚmainrM   ^   së  € Ý	�'�N„N€NØ
€AØ
€AÝ	Œ�B‰Œð ð Ý�a˜‰OŒOˆðð ð ñ ô ð ð ð ð ð ð øøøð ð ð ð å	Œ�Š•b”g—o’o¥hÑ/Ô/°°xÀÑ	JÔ	J€BÝ	ˆb�6‰k˜3Ñ	Ô	ð 
 1Ø	�Š•—’˜a Ñ#Ô#Ñ$Ô$Ð$Ý ‘l”lð 	!ð 	!‰FˆAˆsØIÐIÀSÐIÑIÔIˆCØ�GŠG�C‰LŒLˆLØ�GŠG�D—I’I˜c‘N”NÑ#Ô#Ð#Ø�1�q‘5ŠyˆyØ—’˜‘”��à—’˜Ñ Ô Ð Ð Ø	�Š•‰Œˆð
ð 
ð 
ñ 
ô 
ð 
ð 
ð 
ð 
ð 
ð 
øøøð 
ð 
ð 
ð 
õ „Iˆb�6‰k˜2ÑÔÐÐÐs$   ­A
Á
AÁAÂ,CFÆFÆFÚ__main__)r<   r>   Úscipy.special._precompute.utilsr   Úmpmathr   ÚImportErrorr   r   r    r$   r.   rE   rH   rM   Ú__name__© r   r   ú<module>rT      sò   ðð	ð 	ð 
€	€	€	Ø >Ð >Ð >Ð >Ð >Ð >ð	ØÐÐÐÐøØð 	ð 	ð 	Ø€Dð	øøøð	ð 	ð 	ðð ð ðð ð ð&ð &ð &ðð ð ð(ð ðð ðð ð ð* ˆzÒÐØ€D�F„F€F€F€Fð Ðs   Ž “š