a khR@sddlmZddlmZmZmZmZddlmZddl m Z ddl m Z m Z mZmZddlmZmZmZddlmZmZmZddlmZdd lmZmZdd lmZmZm Z m!Z!dd l"m#Z#m$Z$dd l%m&Z&m'Z'dd l(m)Z)m*Z*ddl+m,Z,m-Z-ddl.m/Z/ddl0m1Z1ddl2m3Z3m4Z4m5Z5m6Z6m7Z7m8Z8m9Z9ddl:m;Z;ddlZ?ddl@mAZAmBZBmCZCmDZDmEZEmFZFmGZGmHZHmIZIddZJddZKddZLddZMddZNdd ZOe;d!d"ZPd#d$ZQd%d&ZRd'd(ZSd)d*ZTe;d+d,ZUd-d.ZVd/d0ZWe;d1d2ZXe;d3d4ZYd5d6ZZd7d8Z[d9d:Z\d;d<Z]d=d>Z^d?d@Z_dAdBZ`dCdDZadEdFZbdGdHZcdIdJZddKdLZedMdNZfdOdPZgdQdRZhdSdTZidUdVZjdWS)X) expand_func)IRationaloopi)S)default_sort_key)Absargre unpolarify)exp exp_polarlog)coshacoshsinh)sqrt) Piecewisepiecewise_fold)cossinsincasin)erferfc)gamma polygamma)hypermeijerg)Integral integrate hyperexpandsimplify)_rewrite_single _rewrite1meijerint_indefinite _inflate_g_create_lookup_tablemeijerint_definitemeijerint_inversion)slow)verify_numericallyrandom_complex_number) xyabcdstzc Cs"dd}dd}|tdt|tddtd|tdttdtdtd|tdt|ttdd}|tt ttttttttttd dttd ttd dd tdd t j td d fd fdtd dd ffdt dt ttd fgdfksJdS)NcSsntttgtg|gtg|t}|dus*Jt|dddtsDJ|dddjt||ffksjJdS)Nr) r&rr2r3r5r0 isinstanceargumentZ as_coeff_mul)exprr4mer?R/var/www/auris/lib/python3.9/site-packages/sympy/integrals/tests/test_meijerint.pyr7s ztest_rewrite_single..tcSs*tttgtgtgtg|tdus&JdSN)r&rr2r3r4r5r0r<r?r?r@tn!sztest_rewrite_single..tnr9cSsHddlm}t||}|dd|dDtt}t|||sDJdS)NrAddcSsg|]}|d|dqS)rr9r?).0resr?r?r@ -z2test_rewrite_single..u..)sympy.core.addrFr&replacerr r.)r<r0rFrr>r?r?r@u*s   ztest_rewrite_single..ur)rDr?@T) r0r1r rr&rrrrrHalfrr)r7rCrNr?r?r@test_rewrite_singles$ $  &*rUc CszttdttgtgtgtgtdttddtdtdddttgtgtgtgtdtdfgdfksvJdS)NrQr9rDrT)r'r0rr2r3r4r5r1r?r?r?r@ test_rewrite1Cs86rWcCs~dd}|dt|dt|ddt|dtd|dttd|tdt|dttddttd dS) NcSsttgtgtgtg||}ttdttdttttti}t|t}|dusZJt | || t |ts|JdS)N ) rr2r3r4r5randcplxrr(r0r.subsdiff)Zfacr grZZintegralr?r?r@r7Is   z0test_meijerint_indefinite_numerically..trDr9rV3/2rQrPz7/3)r0rr7r?r?r@%test_meijerint_indefinite_numericallyHs  r_cCsLtttdd\}}|jr |dus$Jttttt\}}|jrD|dusHJdS)NrT)r+r0is_zeror)vr3r?r?r@test_meijerint_definiteYsrbc sttdttdtttttttdifdd}|tgtggtgtggtds`J|ttgtggtgtggtdsJ|tgtggttgtggdtddsJdS)NrXcsDddlm}t|||}|t||}t||tdddS)Nr)Mulg?g)r3r5)Zsympy.core.mulrcrr)r.rZr0)r2r3r nrcm1m2rZr?r@r7ds  ztest_inflate..trQr9)r2rYr3rr4r5r1r0r^r?rgr@ test_inflate`s "$rhcCs$ddlm}|ddd\}}}tt|d tt|d }t|tdtfdd}t|tdtt t td||ddt|d d|||ddd ksJttt|d tt|d t|ttdtfdd}t|tdtt t tdd|d||d dt|d |dd|d||dd d ksxJtttt|||d tdtfddtt ddt |||ksJtttt|||d tdtfddtt ddt |||ks JdS) Nrsymbolsza b cTpositiver9rrDrP) sympy.core.symbolrjr r0r!rr%expandrrr)rjr2r3r4rMr>r?r?r@test_recursiveps2 $ @@H , ,rqc Cs.ddlm}ddlm}|ddd\}}}ttggdgg||tgg|dg| dg|dd|dtfjsxJ|d dd }tt|tgggdgggttdtft |d ksJtt|tgggdgggttdtfdd t |d ksJt tt|tgggdgggttdtfd d t s6Jt t ttt tksRJ|ddd \}}ttt|td|d||d |d ksJttd dt t tdtdksJ|ddd \}}tt t|d|d tdt\}} t|tt|dt|d|ks&J| dks4Jtt | tt |ttdt\}} t|d ||ksvJtt ttt dt ddfksJ|tt ttdtdt td ksJ|tt t tdtdd t t ttt ttks Jtt td tt tttdfks8Jtt tt tt tdks\Jtt dtdd tt tttddfksJtt tdtd tt tdksJtt t||d dtdt|dtt tdksJtttdtt ttdfks*Jtt t tttdttjdfksVJdd} tdD]>} tt t ttt| tdtfdd | | ksfJqfttt t tt|tdtfdd tdt|tddksJ|d\}}}ttgg|dg| dgtdtgg|dg| dgtdt|d tdtddd|dt d|d t |d|d|t | d|d|d t |d|d|d t |d|d|d t|d kt|td dk@t|dt|dt|dk@fks6Jttt|tt|tdtfdd t tt|tt|tdtfksJtt td tttdtfdd tttdtjdksJddlm} |ddd} ttt t t| tdd | | d tksJ|ddd }tdt|t|ttddtt|t |d td|dtj|dd ffddtjft ddff|ddddfksJ|d dd \}}ttt|t | tdtt tf|| dtjd|d t |dtjdks*JdS)!Nr)rprizs t muTrealr9rPr6rkrDrmFa brQ)Tzsigma mu)r9T)rDTcSs(ddtdt|tdd|SNrDr9rO)r0r[rZ)rdr?r?r@rHsztest_meijerint..resza b s lowergammard)integeralphar?rOruza s)!sympy.core.functionrprorjr!rr is_Piecewiser0rr:r r(r r%r+rrrrr absrrrrTranger rr'sympy.functions.special.gamma_functionsrzr)rprjr6r7mur2r3sigmair4_rHrdrzr|r?r?r@test_meijerints  $  ( &     *(.(&, ,$  ,. &, (( "  :8B &" $      &   *4rcCsvddlm}m}tt|tt|ttttdtfddddt t tdtdt ttttksrJtt|tt|ttttdtfdddddtksJttt tt t ddt t j t dtfddddtdt tt tt jt |t tks Jtt |dt t ddt |dt ksLJtt |dt t ddt |dt ksxJt|dt t dd|dt  ksJt|dt dt t dd|dt d|dt d dksJt|dt dt dt ddd t |dt ddt |dt dd|dt |dt |dt dt ks^Jt|dt |dt t dd|dt d dksJtt d|dt |dt t ddt d|dt ddksJt|dt |dt t t ddt |dt dt |dt d|dt |dt ks.invrDr9) rr'sympy.functions.special.delta_functionsrr6rr7rr rr,)rrrr?r?r@test_inversions  ((".(rcCsddlm}ddlm}|ddd}tt|tdt|tt}t|tt }|j r^J|ddd }| ||}t|tt }|j sJ|j dd| ||ksJ|j d d sJ||tt d}|j d d|ks|j d d|j ksJdS) NrSymbol)InverseLaplaceTransformr2Trkrxr3rrrOrD)rorsympy.integrals.transformsrrrr r6r,r7r~rZargsZ as_integral)rrr2Frr3f2ZILTr?r?r@!test_inversion_conditional_output.s   $      rcCsddlm}ddlm}|ddd}|ddd }d d t}|d }tt|tttj r^Jtt|tttduszJtt|tttdusJtt|ttt}|j sJt |j dd|sJdS) Nrr) DiracDeltarMTrrr4F)Z extended_realrDr9r8) rorrrrr,r r6r7r~r:r)rrrMr4r2r8rr?r?r@%test_inversion_exp_real_nonreal_shiftCs      rcCsbddlm}m}ddlm}ddlm}i}t||D]}t |t dD] \}}}} i} t |j |gD]2} t | dr| jr|dd| | <qp|d d | | <qpt|t s|| }d d |D} td d| DsJ|dd |D} |j| d| j| d} }tt| t|}|dkr>t| |dksZJqNt| ||dksNJqNq.css |]}|jp|t VqdSrA)r~hasr)rGr0r?r?r@ grJz$test_lookup_table..cSsg|]\}}||qSr?r?)rGrr0r?r?r@rIjrJrgg|=)sympy.core.randomrrrKrFsympy.integrals.meijerintr8r*valuessortedrlistZ free_symbolshasattrrr:allrdminr)rrrFZz_dummytablelZformulaZtermsZcondhintrZZaiexpandedr2r3rMr?r?r@test_lookup_tableRs,    rcCsddlm}ddlm}|tttdtddtdddttdtt dddtt d dkspJtttdtdddtttdtt dddtt d ddtt dd|t ddtd dt t tt d dksJdS) Nrry) powdenestrQTrmZpolarr9rVrw) rrzsympy.simplify.powsimprr!rr0r[rrrr)rzrr?r?r@test_branch_bugss  .2>rcCsdddlm}tttdtddtdt ks4Jt|dtdtdd|ddt ks`JdS)NrrrDTrm)rrr!rr0rrr?r?r@test_linear_subs~s (rc$s.ddlmddlm}m}ddlm}ddlm}|ddd\}|d dd \}|d dd d d ddt t t t t fdddksJt t t t t t fddksJt t dt t t t fddddksJt t dt t t t fdddddks>Jt t t ||t t t ft t t fdddks|Jt t t t ||t t t ft t t fddksJt t t t ||t t t ft t t fdd|ksJt t t t t ||t t t ft t t fdd|ksJJt t t dt t ||t t t ft t t fddd|ksJt t t dt t ||t t t ft t t fddd|ksJt t dt t ||t t t ft t t fdd}| tr:Jt|ddksXJt t dt t ||t t t ft t t fdd|d|dksJt t t dt fdddksJt t t t dt fdddksJt t dt t dt fddddks.Jfdd}|ddksVJ|t t kspJ|t t dddksJddd}t|t t dd|t t dd|ksJt|t t dd|t t dd|ksJt|t t d|t t d|ksLJ|ddd \} } t | ddt | | t| | t| t| } t | t dt fdddksJt t | t dt fddd}||d|df| | dd| kfksJt t d| t dt fddd} | d| dkks0J|| d|dd| | d| | d| ddksrJ|ddd \} }t | dt d|dt| |t| t|} tt | t ddfdddksJtt t | t ddfdd| | |ks Jtt t d| t ddfdd| | d| || |dksRJtt t t | t ddfddt| |t| t t| t| |t ksJ|dddd}dd|dt |dtt d dt|d}|t |t dt fdddksJtt t |t dt fddtdt|ddt|dksXJtt t d|t dt fdd|ksJd| dt|dt |ddtt d}|t |t dt fdddksJtt t |t dt fdd|k sJtt t d|t dt fdd||dk s4J|t t |td|d|t dt fdddtdt|k s~J|ddd \} }}| |t t || |dt | || |d}tt |t dt fdddk sJt |}tt |t dt fddd| |tdd| t|dd| | |dt|k sRJtt t |t dt fddd| |dtdd| t|dd| | |dt|k sJ|ddd \}}t|t ||||t |||t t|dt|dt||d}tt |t dt fdddk sDJtt t |t dt fddd||dk stJtt t d|t dt fddd|d|d||d |dk sJ|d!dd \}}t|dtt td"dt| t |dt d|d}d#d$}|t |t dt fdk s:J|t t |t dt f|k s\J|t t |d|t dt f|d|k sJ|t t |d|t dt fd|d%|dk sJ|d&dd }t t|dtt| dt |t |td't |t fdk s J|ddd \} } | | t | d| | ddt | | | d}tt |t dt fdk sJtt t |t dt fdd(t| | tt| k sJtt t t |t dt fdd(t| t t | ttt | ks J|ddd }|d)dd }||t ||dtt || }tt |t dt fdksjJtt t ||t dt f||td||ksJdd*lm}|d+dd \} }!t |!dtt d| d d|!d|dt | |!d}"t |"t dt fdddks"J|d,dd-}|d.dd }ttt | |d|}#t |#t t t fdddksxJt t |#t t t fdd|ksJt t d|#t t t fddd|d|dksJ|ddd }|t tt t |dtt t|t dt ftd|ks*JdS)/Nr) expand_mul)rrj) gammasimp)powsimpzmu1 mu2TZnonzeroz sigma1 sigma2rklambdacSs6dtdt|dt||d d|dSNrDr9)rrr )r0rrr?r?r@normalsz test_probability..normalcSs|t| |SrA)r )r0rater?r?r@ exponentialsz%test_probability..exponentialrmrDr9rQrOcst|tttdtftt tfdd}t|tttt tftdtfdd}||ks|J|S)NrTrm)r!r0r1r)r<Zres1Zres2rrZmu1rrZsigma1r?r@Esztest_probability..Ez alpha betaZseparaterrtk)r{rlza b przd1 d2rPzlamda mucSst|tSrA)r%rewriter rBr?r?r@rJz"test_probability..rVr4r])rrd)rznu sigmarrrr3)r}rrorrjZsympy.simplify.gammasimprrrr!r0rr1rr r%rr rrrrrrrrrr)$rrjrrZmu2Zsigma2rransr|betaZbetadistjr2r3rchiZ chisquaredpZdagumr d1Zd2rZlamdardistZmysimpr4ZdistnrdrnurZriceZlaplacer?rr@test_probabilitys    & $$      $ $    $   .66.(.$ :"  . :"" *8"& & 8"* . .  "  " D"2: 6  "   *  D   $  *rc Csddlm}ddlm}ddlm}m}m}m}m }m }t t t tt ttttdtfddd|jdd |ttksJt tt tttdtfddd||dtksJt tt ttd tdtfddd||d t||ksJt tt ttd tdtfddd||d t||ks`J|d dd }t tt tt|tfdd||ksJt tt tt|tfdd||td ksJt ttttdtfdd|tksJt |tttdtfdd|tks(Jt tt ttdd|tt|dtksbJt tt td tdd||dttt tttksJ|ddd} t t| | | dd| d|| ksJt t| | | dd| d|| ksJt |dttdd|t|dttt ksZJt |d ttdd|td  |dtd ttt d tt d ksJt t t |tttdd|jdd |tdt ksJt |ttddt|tttks"Jt || | dd| || t| ksRJt |ttddt|tttks~Jt || | dd| || || ksJt |ttt tdtfddtdksJt |dttttdtfddtd d ksJdS)z% Test various exponential integrals. rr)r)ChiCiEiShiSiexpintrDTrr)funcr9rQr7rkrmrNrrPN)ror%sympy.functions.elementary.hyperbolicr'sympy.functions.special.error_functionsrrrrrrr%r r!r r8r0r1rrrprrrrZas_independentrr) rrrrrrrrr7rNr?r?r@ test_expintKs   "     0  **  $  " "  6  ,0,00rc Csddlm}m}ddlm}m}ddlm}ddlm }m }m }m }ddl m} m} | |tttdd|t td tddfksJ| |tttdd|ttt td d kfksJ| |tttddttd td d td d tt td d kfksJ| |ttttd ddttd kttd k@fksJJ| |d ttttd d} | djdd| d ftdtd d tktd d tkBfd tdtd td d dftdkfksJt|t|dttdtfddtd td ksJt|t|d ttdtfddttjtd d ksNJtd ttd td tddt|d t ttd d kft|d tdfksJdS)Nr)racoth)ratanr)rE1rr)fourier_transformlaplace_transformTr$r9rDrxrOF)Znoconds)deeprSrm) rrr(sympy.functions.elementary.trigonometricrrrrrrrrrrrrr0r6rrrr2r Zfactorrprrr!rrTrr) rrrrrrrrrrrrr?r?r@ test_messysD  4  $"      2rcCsJttt tdtt tfddt tttttdksFJdS)Nr9TrmrP)r!r rr0rrrr?r?r?r@test_issue_6122s"rcCs>dtttttdd}t|tdd}|tr:JdS)NrDrQTrm)r0r2r3rr!rr)r<Zantir?r?r@test_issue_6252srcCsDttttdtdtt tftttdks@JdSrv)r!r rr0rr%rrr?r?r?r@test_issue_6348s. rcCshddlm}m}tttttddt|tks:Jtttttddt|tksdJdS)Nrfresnelcfresnelsr9) rrrrr!rrr0rrr?r?r@ test_fresnels*rcCstttttdusJdSrA)r(r0r?r?r?r@test_issue_6860srcCsrtttdtdt}|tdtddtdtddksNJ|ttjtjtddksnJdS)Nr9rQrV) r(r0rZtogetherZ_eval_intervalrZ NegativeOneZOnerrr?r?r@test_issue_7337s0rcCshttttt ttt dttt dtt tt tdtddksdJdSr)r(rr0r r7r?r?r?r@test_issue_8368s0 rcCstddlm}m}tdtttd|ddtd|ftd|fdtd|d|d|d|kspJdS)NrhwrDr9rQ) sympy.abcrrr!rr1r0rr?r?r@test_issue_10211s2&rcCsrddlm}|ddd\}}tdttd|ddt| |fd||dt|d|dksnJdS) Nrrizy LTrkrDr9rQ)rorjr!rr0)rjr1Lr?r?r@test_issue_11806s  ("rcCsddlm}ddlm}m}t|d|d|dd|dd}d|d |d td td dftd df|dtdt t |d}| || d d sJdS)Nr)RR)RrMr9g?TrmgUUUUUU?g?rQgrVg-q=) Zsympy.polys.domains.realfieldrrrrMr!rrrrrZalmosteqrd)rrrMrr\r?r?r@test_issue_10681s &*rcCs@ddlm}|ddd}tdtdtt|fd|kssr      $ ,) j0   G 6&