U Zʹj|-@sjddlmZmZeddZeddZeddZedd Zed;d d Zed dZeddZ edddZ ed?ddZ ed@ddZedAddZedBdd Zed!d"ZedCd#d$Zed%d&Zed'd(Zed)d*Zed+d,Zed-d.Zed/d0Zed1d2Zed3d4Zed5d6Zed7d8Zed9d:ZdS)D)defun defun_wrappedcCsL||d}d||j||dd|}||sH|||j}|S)Nrrrsquare_exp_argsqrtpiZhyp1f1_re_imjctxzz2vr\/var/www/html/bwave/python/venv/lib/python3.8/site-packages/mpmath/functions/expintegrals.py _erf_complexs  " rcCs|||dkrL||}|j|dd}||||j|dd|}nd||}||sxd| ||j }|S)NrTexactrr) rer fnegexpr r Zhyperurr rr)rrrZnz2rrrr _erfc_complex s & rcCs|||}||r6z ||WStk r4YnX||rr|jsrzt|||jWStk rpYnX||SN) convert _is_real_typeZ_erfNotImplementedError_is_complex_typeimagtyperealrrrrrrerfs   r'cCs|||}||r6z ||WStk r4YnX||rr|jsrzt|||jWStk rpYnX||Sr) rr Z_erfcr!r"r#r$r%rr&rrrerfc's   r(FcCsb|jdd}|r6|j|||d}|j|j||d}n|j|||d}|dkr^|j||dd}|S)N)precrTr)r+fmulZfdivone)rrZmultZ reciprocalr+rrrrr 6sr cCsR|s|S||}d||j||dd|}||sN|||j}|S)Nrrrr rrrrerfiBs " r.cs}|ks"|dks"|dkr,dS|s8SdkrFjSdkrTjStdkrvddd}nHdjtdd}||d}j d 7_ fd d |S) Nrrz*erfinv(x) is defined only for -1 <= x <= 1g?goe1?rgϷ?r cs|Sr)r')trxrr]zerfinv..) r Z bad_domaininfninfabslnr signr r+Zfindroot)rr2Zxreaurr1rerfinvLs     (r<cCs>||}|||d d|d||d|jS)Nr)rrr r )rr2musigmarrrnpdf_s r@cCsD||||d}|dkr.|| dSd||dSdS)Nrr=r)r r(r')rr2r>r?r:rrrncdfdsrAc Cs||krd}n|sV|dkr0|dkr0|||}q||||d||d||}n||\}}|dkr||j kr|j } |jd9_|| 7}n|dkr|j|8_||||d||d|} ||||d||d|} | | |}|r||||}|S)Nr=rr)betaZhyp2f1Z nint_distancer+eps) rr:bx1Zx2 regularizedrmdhs1s2rrrbetaincls(&     rMNcCs t|}||}|dkr&|j}d}n||}||jk}|dkrN|j}d}n||}||jk}|s|s|r||dkr|jS||dkr|jS|jS||S||kr|jS||||kr||||| S|r|r| |||| S|r| |||S|r| |||SdS)NFr=) boolrzeror5rr-nangammagammainc_gamma3 _upper_gamma _lower_gamma)rrr:rErGZlower_modifiedZupper_modifiedrrrrRs<      rRcsP|rt|jS|g|jddfdd}||gS)NTrcs2|gd|dggdgd|gf}|fS)Nrrr)rT1GrErZnegbrrrJs,z_lower_gamma..h)isnpintr$r5r hypercomb)rrrErGrJrrXrrUs   rUcsH|rzl|rf|r*t|jWSj}z*jd7_||WWS|_Xn|WSWntk rYnX|dkrdkr|dS|dkr̈dksdkr|dSjdd |g|z&fd d }j ||gdd WSj k rBfd d } ||gYSXdS)Nr/rrr=ryy?Trcs4|d}gd|ggd| ggdfgSNrrV)rrrYr:rZnegarrrJsz_upper_gamma..h) force_seriescsRgd|dg|gggdf} |gd|dggdgd|gf}||fS)Nrr=rrV)rrWT2r^rrrJs.) isintrZr$rOr+Z_gamma_upper_intrQr!rr[ NoConvergence)rrr:rGorigrJrr^rrTs2      rTc Cs||}|r|r|jSz|jd7_|j|||d}|j|||d}||}||t||||dkr||W>S|s|j|d||d}|j|d||d}||}|WSW5|jd8_XtdS)N)rGir=)rZrOr+rRmagmaxr!) rrr:rErGZpolerWr`RrrrrSs" $ rScCs||r8||r8z|||WStk r6YnX||sL||rT||S||jkrfd|S|dkr||dkrt||jS|j|dS|dkr| | |S|dkr| | |d|dS||d| d||S)Nrr=rr) rar Z _expint_intr!isnanr5rr$r-rrR)rnrrrrexpints$ rjcCsT|r.|dkr|jS|||||jS|s6|S|dkrD|jS|||S)Nrr)rOeir8Zln2r6)rroffsetrrrlisrmcCs0z ||WStk r*||YSXdSr)Z_eir! _ei_genericr&rrrrks rkcCs||jkr|S||jkr|jS||dkrzr|j|}|||jddgg||jdd|}||}|dkr||j |j 7}|dkr||j |j 8}|WS|j k rYnX|| dddd||j }||r|d||||j|7}n||t|7}|S)NrT)Zmaxtermsr_r=r?)r5r6rOrer-rZhyperr+rr rrbZhyp2f2Zeulerlogr7)rrr]rZimrrrrn&s4     $rncCs2z ||WStk r,|d|YSXdSr\)Z_e1r!rjr&rrre1Bs rqcCs0z ||WStk r*||YSXdSr)Z_cir! _ci_genericr&rrrciIs rscCs||r.||jkr|jS||jkr.|jdS|j|j|dd}|j|dd}d||||}| |}| |}|dkr|dkr||jd7}|dkr||jd8}|dkr|dkr||jd7}|dkr||jd8}| |r|dkr| |}|S)N?Trror=?) isinfr5rOr6r r,rrrkr rr )rrjznjzrzrealzimagrrrrrPs.      rrcCs0z ||WStk r*||YSXdSr)Z_sir! _si_genericr&rrrsids r|cCs||r2||jkrd|jS||jkr2d|jS||dkr|j|j|dd}|j|dd}d||||}| |}|dkr|d|j8}|dkr|d|j7}| |r| |}|S|| dd d d ||SdS) NrorTryr=rrgп) rvr5r r6rer,rrrkr r hyp1f2)rrrwrxrryrrrr{ks$      r{cCs|j|dd}d||||}||}||}|dkrR||jd7}n.|dkrj||jd8}n|dkr||jd7}|S)NTrror=rurt)rrkr rr )rrnzrryrzrrrchis  rcCs||dkrn|j|dd}d||||}||}|dkrT|d|j8}|dkrj|d|j7}|S||dddd ||SdS) NrTrror=rurrg?)rerrkrr r~)rrrrrzrrrshis rc Cs^||jkr|dS||jkr(|dS|j|dd|ddd|jd |d d S) Nror}r)rr)r)r)rr))r5mpfr6r r~r&rrrfresnelss     rc CsP||jkr|dS||jkr(|dS||ddd|jd |ddS) Nror})rr)r)r)rr)r)r5rr6r~r r&rrrfresnelcs     r)rF)r=r)r=r)r=rF)r=NF)F)F)F)F)Z functionsrrrrr'r(r r.r<r@rArMrRrUrTrSrjrmrkrnrqrsrrr|r{rrrrrrrrsl           (  &