
    tKg                         d Z ddlZ	 ddlZddlmZ  ej                  d      Zd ZdZd Z	e
dk(  r e	        yy# e$ r Y w xY w)	zPrecompute the polynomials for the asymptotic expansion of the
generalized exponential integral.

Sources
-------
[1] NIST, Digital Library of Mathematical Functions,
    https://dlmf.nist.gov/8.20#ii

    N)Polyxc           	         t        dt              g}t        |       D ]b  }|j                  t        dd|z  t        z  z
  t              ||   z  t        t        t        dz   z        ||   j	                         z  z          d |S )N      )r   r   rangeappenddiff)KAks      f/home/alanp/www/video.onchill/myenv/lib/python3.12/site-packages/scipy/special/_precompute/expn_asy.py
generate_Ar      sj    	aA1X	a!A#a%i#AaD(41q5	?1Q499;+FFG H    ze/* This file was automatically generated by _precompute/expn_asy.py.
 * Do not edit it manually!
 */
c                  b   t        t               t        j                  j	                  ddd      } d}t        |      }t        | dz   d      5 }|j                  t               |j                  dt        |       d       t        |      D ]_  \  }}d	j	                  |j                         D cg c]  }t        |j                  d
             c}       |j                  d| d       a d	j	                  t        |dz         D cg c]  }d| 	 c}       |j                  d       d	j	                  |D cg c]  }t        |j                                c}       |j                  d       d d d        t        j                   | dz   |        y c c}w c c}w c c}w # 1 sw Y   2xY w)Nz..cepheszexpn.h   z.newwz#define nA 
z,    zstatic const double Az[] = {tmp};
r   r   z$static const double *A[] = {{tmp}};
z$static const int Adegs[] = {{tmp}};
)print__doc__ospathjoinr   openwriteWARNINGlen	enumeratecoeffsstrevalfr   degreerename)fnr   r   fr   Akr   s          r   mainr)   "   sP   	'N	dHh	/B
A1A	b6k3	1		+c!fXR()q\EArII=As1772;'=>GG+A3o>? " 			E!a%L1LqQqc7L12	78		a0a3ryy{#a01	78 
  IIb6k2 >10 
 	s7   A"F%/!F9F%	F
'F%< F 
F%F%%F.__main__)r   r   sympyr   symbolsr   ImportErrorr   r   r)   __name__ r   r   <module>r0      sd    
	cA
& zF K  		s   : AA