ó <Zc@sßddlmZmZed„ƒZed„ƒZed„ƒZed„ƒZeded„ƒZed„ƒZ ed„ƒZ ed dd „ƒZ ed dd „ƒZ ed ded „ƒZ ed ded „ƒZeed„ƒZeed„ƒZeed„ƒZed„ƒZeed„ƒZed„ƒZed„ƒZed„ƒZed„ƒZed„ƒZed„ƒZed„ƒZed„ƒZed„ƒZed„ƒZed„ƒZdS(i(tdefunt defun_wrappedcCsk|j|dƒ}d|j|jƒ||jdd|ƒ}|j|ƒsg|j|ƒ|j}n|S(Niÿÿÿÿiii(ii(ii(tsquare_exp_argtsqrttpithyp1f1t_ret_imtj(tctxtztz2tv((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyt _erf_complexs -cCs±|j|ƒdkrn|j|ƒ}|j|dtƒ}|j|ƒ|j|jƒ|jdd|ƒ}nd|j|ƒ}|j |ƒs­d|j |ƒ|j }n|S(Nitexacti(ii(ii( treRtfnegtTruetexpRRthyperuR RRR(R R R tnz2R ((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyt _erfc_complex s5cCs£|j|ƒ}|j|ƒrFy|j|ƒSWqFtk rBqFXn|j|ƒr–|j r–y t|ƒ|j|jƒƒSWq–tk r’q–Xn|j|ƒS(N( tconvertt _is_real_typet_erftNotImplementedErrort_is_complex_typetimagttypetrealR (R R ((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyterfs   cCs£|j|ƒ}|j|ƒrFy|j|ƒSWqFtk rBqFXn|j|ƒr–|j r–y t|ƒ|j|jƒƒSWq–tk r’q–Xn|j|ƒS(N( RRt_erfcRRRRRR(R R ((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyterfc's   cCs|jdd}|rM|j||d|ƒ}|j|j|d|ƒ}n|j||d|ƒ}|dkrŒ|j||dtƒ}n|S(NiitpreciR(R!tfmultfdivtoneR(R R tmultt reciprocalR!R ((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyR6s cCsr|s |S|j|ƒ}d|j|jƒ||jdd|ƒ}|j|ƒsn|j|ƒ|j}n|S(Niii(ii(ii(RRRRRRR(R R R R ((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyterfiBs-cs*ˆjˆƒ}|ˆks3|dks3|dkr@ˆjdƒS|‰ˆsPˆSˆdkrcˆjSˆdkrvˆjStˆƒdkr¡dˆddˆ}n^ˆjdˆjtˆƒddƒ}ˆjˆƒˆj|ˆj|ƒƒˆjdƒ}ˆj d 7_ ˆj ‡‡fd †|ƒS( Niÿÿÿÿis*erfinv(x) is defined only for -1 <= x <= 1gÍÌÌÌÌÌì?g¨o™Óe1á?ig†åÏ·ê?ii csˆj|ƒˆS(N(R(tt(R tx(s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyt]s( Rt bad_domaintinftninftabstlnRtsignRR!tfindroot(R R)txretatu((R R)s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyterfinvLs $   (6icCsI|j|ƒ}|j||d d|dƒ||jd|jƒS(Ni(RRRR(R R)tmutsigma((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pytnpdf_scCsR||||jdƒ}|dkr9|j| ƒdSd|j|ƒdSdS(Niii(RR R(R R)R6R7R3((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pytncdfds c Csp||krd}n8|sw|dkrH|dkrH|j||ƒ}qM|||j|d||d|ƒ|}nÖ|j|ƒ\}}|dkrï||j krÎ|j } |jd9_|| 7}qï|dkrï|j|8_qïn|||j|d||d|ƒ} |||j|d||d|ƒ} | | |}|rl||j||ƒ:}n|S(Niiiiüÿÿÿ(tbetathyp2f1t nint_distanceR!teps( R R3tbtx1tx2t regularizedR tmtdthts1ts2((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pytbetaincls(  /    ((cCs•t|ƒ}|j|ƒ}|dkr9|j}t}n|j|ƒ}||jk}|dkru|j}t}n|j|ƒ}||jk}|pœ|sô|rç|j|ƒdkrÁ|jS|j|ƒdkrÝ|jS|jSn|j |ƒS||kr|jS|j|ƒ|j|ƒkr<|j ||||ƒ S|r_|r_|j ||||ƒ S|rx|j |||ƒS|r‘|j |||ƒSdS(Ni(tboolRtNonetzerotFalseR,RR$tnantgammatgammainct_gamma3t _upper_gammat _lower_gamma(R R R3R>RAtlower_modifiedtupper_modified((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyRN…s<            csoˆj|ƒr"t|ƒˆjƒS|g|‰ˆjˆdtƒ‰‡‡‡‡fd†}ˆj||gƒS(NRcsGˆjˆƒˆ|gd|dggˆdgd|gˆf}|fS(Niiÿÿÿÿ(R(R tT1(tGR>R tnegb(s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyRDµs@(tisnpintRR,RRt hypercomb(R R R>RARD((RUR>R RVs<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyRQ®s  cs™ˆj|ƒr«y…|rƒˆj|ƒr:t|ƒˆjƒSˆj}z0ˆjd7_ˆj|ˆƒˆj|ƒSWd|ˆ_Xnˆj|ˆƒSWq«tk r§q«Xn|dkrψdkrÏ|ˆdS|dkrÿˆd ksóˆd krÿ|ˆdSˆjˆdt ƒ‰|g|‰y5‡‡‡‡fd†}ˆj ||gd t ƒSWn<ˆj k r”‡‡‡‡fd †}ˆj ||gƒSXdS( Ni iiÿÿÿÿiiyð?RcsF|d}ˆjˆƒˆgd|ggˆd| ggdˆfgS(Ni(R(R tr(RUR3R tnega(s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyRDÖs t force_seriescssgd|dg|gˆggdf}ˆjˆƒ ˆ|gd|dggˆdgd|gˆf}||fS(Niiiÿÿÿÿ(R(R RTtT2(RUR3R RZ(s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyRDÛs(Ayð¿ð¿yð¿ð?( tisintRWRRJR!t_gamma_upper_intRMRRRRXt NoConvergence(R R R3RAtorigRD((RUR3R RZs<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyRPºs2 !   $  c Cs|j|ƒ}|r"|r"|jSzÏ|jd7_|j||d|ƒ}|j||d|ƒ}||}|j|ƒt|j|ƒ|j|ƒƒdkr¦|S|sð|j|d|d|ƒ}|j|d|d|ƒ}||}|SWd|jd8_Xt‚dS(NiRAiöÿÿÿi(RWRJR!RNtmagtmaxR( R R R3R>RAtpoleRTR\tR((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyROás"  4 cCs4|j|ƒrI|j|ƒrIy|j||ƒSWqItk rEqIXn|j|ƒsg|j|ƒro||S||jkr†d|S|dkrÌ|j|ƒdkrºt|ƒ|jƒS|j|dSn|dkrê|j | ƒ|S|dkr|j | ƒ|d|dS||d|j d||ƒS(Niiiÿÿÿÿi( R]Rt _expint_intRtisnanR,RRR$RRN(R tnR ((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pytexpintûs$    cCsr|r?|dkr|jS|j|j|ƒƒ|j|jƒS|sI|S|dkr\|jS|j|j|ƒƒS(Nii(RJteiR/tln2R-(R R toffset((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pytlis & cCs3y|j|ƒSWntk r.|j|ƒSXdS(N(t_eiRt _ei_generic(R R ((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyRis c Cs}||jkr|S||jkr)|jS|j|ƒdkrýy¥|j|}|j|ƒ|jddgg|d|jdtƒ|}|j |ƒ}|dkr»||j |j 7}n|dkrÞ||j |j 8}n|SWqý|j k rùqýXn||j dddd|ƒ|j}|j |ƒr`|d|j|ƒ|j|j|ƒ7}n||jt|ƒƒ7}|S(NitmaxtermsR[iigà?(R,R-RJRaR$RthyperR!RRRRR_thyp2f2teulertlogR.(R R RYR tim((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyRn&s* ;  &.cCs6y|j|ƒSWntk r1|jd|ƒSXdS(Ni(t_e1RRh(R R ((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyte1Bs cCs3y|j|ƒSWntk r.|j|ƒSXdS(N(t_ciRt _ci_generic(R R ((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pytciIs cCs|j|ƒrB||jkr%|jS||jkrB|jdSn|j|j|dtƒ}|j|dtƒ}d|j |ƒ|j |ƒ}|j |ƒ}|j |ƒ}|dkrÿ|dkrÜ||jd7}n|dkrÿ||jd8}qÿn|dkrN|dkr+||jd7}n|dkrN||jd8}qNn|j |ƒr{|dkr{|j |ƒ}n|S(Nyð?Rgà?iyà?( tisinfR,RJR-RR"RRRRiRRR(R R tjztnjzR tzrealtzimag((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyRxPs.       cCs3y|j|ƒSWntk r.|j|ƒSXdS(N(t_siRt _si_generic(R R ((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pytsids cCsE|j|ƒrF||jkr)d|jS||jkrFd|jSn|j|ƒdkr|j|j|dtƒ}|j|dtƒ}d|j |ƒ|j |ƒ}|j |ƒ}|dkrÚ|d|j8}n|dkrú|d|j7}n|j |ƒr|j |ƒ}n|S||j d d d d ||ƒSdS(Ngà?gà¿iÿÿÿÿRyà¿iiiigп(ii(ii(ii( RzR,RR-RaR"RRRRiRRthyp1f2(R R R{R|R R}((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pyR€ks$    cCs·|j|dtƒ}d|j|ƒ|j|ƒ}|j|ƒ}|j|ƒ}|dkrs||jd7}n@|dkr“||jd8}n |dkr³||jd7}n|S(NRgà?iyà?yð?(RRRiRRR(R R tnzR R}R~((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pytchi€s    cCsÃ|j|ƒdkr|j|dtƒ}d|j|ƒ|j|ƒ}|j|ƒ}|dkry|d|j8}n|dkr™|d|j7}n|S||jd d d d ||ƒSdS( NiÿÿÿÿRgà?iyà?iiigÐ?(ii(ii(ii(RaRRRiRRR‚(R R RƒR R~((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pytshiŽs   cCsu||jkr|jdƒS||jkr8|jdƒS|j|dd|jd d d |jd |ddƒS( Ngà?gà¿iiiiii(ii(ii(ii(R,tmpfR-RR‚(R R ((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pytfresnels›s   cCsf||jkr|jdƒS||jkr8|jdƒS||jdd d |jd |ddƒS( Ngà?gà¿iiiii(ii(ii(ii(R,R†R-R‚R(R R ((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pytfresnelc£s   N( t functionsRRR RRR RKRR'R5R8R9RGRIRNRQRPRORhRlRiRnRvRyRxRR€R„R…R‡Rˆ(((s<lib/python2.7/site-packages/mpmath/functions/expintegrals.pytsH    ( &