U L?hR@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 fgdfkst dS)NcSsntttgtg|gtg|t}|dk s*tt|dddtsDt|dddj t||ffksjtdS)Nr) r&rr2r3r5r0AssertionError isinstanceargumentZ as_coeff_mul)exprr4mer@V/var/www/html/venv/lib/python3.8/site-packages/sympy/integrals/tests/test_meijerint.pyr7s ztest_rewrite_single..tcSs*tttgtgtgtg|tdks&tdSN)r&rr2r3r4r5r0r:r=r@r@rAtn!sztest_rewrite_single..tnr9cSsHddlm}t||}|dd|dDtt}t|||sDtdS)NrAddcSsg|]}|d|dqS)rr9r@).0resr@r@rA -sz2test_rewrite_single..u..)sympy.core.addrGr&replacerr r.r:)r=r0rGrr?r@r@rAu*s  ztest_rewrite_single..ur)rEr@@T)r0r1r rr&rrrrrHalfrrr:)r7rDrNr@r@rAtest_rewrite_singles& $  &rUc CszttdttgtgtgtgtdttddtdtdddttgtgtgtgtdtdfgdfksvtdS)NrQr9rErT) r'r0rr2r3r4r5r1r:r@r@r@rA test_rewrite1Cs86rWcCs~dd}|dt|dt|ddt|dtd|dttd|tdt|dttddttd dS) NcSsttgtgtgtg||}ttdttdttttti}t|t}|dk sZt t | || t |ts|t dS)N ) rr2r3r4r5randcplxrr(r0r:r.subsdiff)Zfacr grZZintegralr@r@rAr7Is  z0test_meijerint_indefinite_numerically..trEr9rV3/2rQrPz7/3)r0rr7r@r@rA%test_meijerint_indefinite_numericallyHs  r_cCsLtttdd\}}|jr |dks$tttttt\}}|jrD|dksHtdS)NrT)r+r0is_zeror:r)vr3r@r@rAtest_meijerint_definiteYsrbc sttdttdtttttttdifdd}|tgtggtgtggtds`t|ttgtggtgtggtdst|tgtggttgtggdtddstdS)NrXcsDddlm}t|||}|t||}t||tdddS)Nr)Mulg?g)r3r5)Zsympy.core.mulrcrr)r.rZr0)r2r3r nrcm1m2rZr@rAr7ds  ztest_inflate..trQr9) r2rYr3rr4r5r1r0r:r^r@rgrA test_inflate`s "$rhcCs$ddlm}|ddd\}}}tt|d tt|d }t|tdtfdd}t|tdtt t td||ddt|d d|||ddd kst ttt|d tt|d t|ttdtfdd}t|tdtt t tdd|d||d dt|d |dd|d||dd d ksxt tttt|||d tdtfddtt ddt |||kst tttt|||d tdtfddtt ddt |||ks t dS) Nrsymbolsza b cTpositiver9rrErP) sympy.core.symbolrjr r0r!rr%expandrrrr:)rjr2r3r4rMr?r@r@rAtest_recursiveps2 $ @@H , ,rqc Cs.ddlm}ddlm}|ddd\}}}ttggdgg||tgg|dg| dg|dd|dtfjsxt|d dd }tt |tgggdgggt t dtft |d ksttt |tgggdgggt t dtfdd t |d kstt tt |tgggdgggt t dtfd d t s6tt tt t tt ksRt|ddd \}}ttt |t d|d||d |d ksttt d dtt t dtdkst|ddd \}}ttt |d|d t dt\}} t|tt|dt|d|ks&t| dks4ttt| t t|t t dt\}} t|d ||ksvtttt t t dtddfkst|ttt t dtdttd kst|ttt t dt dd ttttt  tt ks tttt d t t tttdfks8ttttt  t t tdks\tttdt dd t t tttddfksttttdt d t t tdkstttt ||d dtdt|dt t tdkstttt dt t ttdfks*tttt tt t dttjdfksVtdd} tdD]>} ttt tt t | t dtfdd | | ksftqftttt tt |t dtfdd tdt|tddkst|d\}}}ttgg|dg| dgt dtgg|dg| dgt dt |d t dtddd|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@fks6tttt |tt |t dtfdd t tt |tt |t dtfkstttt d tt t dtfdd tttdtjdkstddlm } |ddd} tttt t | t dd | | d t kst|ddd }tdt |t|t t ddtt|t |d td|dtj|dd ffddtjft!ddff|ddddfkst|d dd \}}ttt |t| t dt t tf|| dtjd|d t |dtjdks*tdS)!Nr)rprizs t muTrealr9rPr6rkrErmFa brQ)Tzsigma mu)r9T)rETcSs(ddtdt|tdd|SNrEr9rO)r0r[rZ)rdr@r@rArIsztest_meijerint..resza b s lowergammard)integeralphar@rOruza s)"sympy.core.functionrprorjr!rr is_Piecewiser:r0rr;r r(r r%r+rrrrr absrrrrTranger rr'sympy.functions.special.gamma_functionsrzr)rprjr6r7mur2r3sigmair4_rIrdrzr|r@r@rAtest_meijerints  $  ( &    *(.(&, ,$  ,. &, (( "  :8B &" $      &   *4rcCsvddlm}m}tt|tt|ttttdtfddddt t tdtdt ttttksrt tt|tt|ttttdtfdddddtkst ttt tt t ddt tj t dtfddddtdt tt ttjt |t tks t tt |dt t ddt |dt ksLt tt |dt t ddt |dt ksxt t|dt t dd|dt  kst t|dt dt t dd|dt d|dt d dkst t|dt dt dt ddd t |dt ddt |dt dd|dt |dt |dt dt ks^t t|dt |dt t dd|dt d dkst tt d|dt |dt t ddt d|dt ddkst t|dt |dt t t ddt |dt dt |dt d|dt |dt ks.invrEr9) rr'sympy.functions.special.delta_functionsrr6rr7r:rr rr,)rrrr@r@rAtest_inversions  ((".(rcCsddlm}ddlm}|ddd}tt|tdt|tt}t|tt }|j r^t |ddd }| ||}t|tt }|j st |j dd| ||kst |j d d st ||tt d}|j d d|ks|j d d|jkst dS) NrSymbol)InverseLaplaceTransformr2Trkrxr3rrrOrE)rorsympy.integrals.transformsrrrr r6r,r7r~r:rZargsZ as_integral)rrr2Frr3f2ZILTr@r@rA!test_inversion_conditional_output.s   $      rcCsddlm}ddlm}|ddd}|ddd }d d t}|d }tt|tttj r^t tt|tttdkszt tt|tttdkst tt|ttt}|j st t |j dd|st dS) Nrr) DiracDeltarMTrrr4F)Z extended_realrEr9r8) rorrrrr,r r6r7r~r:r;r)rrrMr4r2r8rr@r@rA%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| Dst|dd |D} |j| d| j| d} }tt| t|}|dkr>t| |dksZtqNt| ||dksNtqNq.css |]}|jp|t VqdSrB)r~hasr)rHr0r@r@rA gsz$test_lookup_table..cSsg|]\}}||qSr@r@)rHrr0r@r@rArJjsrgg|=)sympy.core.randomrrrKrGsympy.integrals.meijerintr8r*valuessortedrlistZ free_symbolshasattrrr;allr:rdminr)rrrGZz_dummytablelZformulaZtermsZcondhintrZZaiexpandedr2r3rMr@r@rAtest_lookup_tableRs,    rcCsddlm}ddlm}|tttdtddtdddttdtt dddtt d dkspt tttdtdddtttdtt dddtt d ddtt dd|t ddtd dt t tt d dkst dS) Nrry) powdenestrQTrmZpolarr9rVrw) rrzsympy.simplify.powsimprr!rr0r[rrr:rr)rzrr@r@rAtest_branch_bugss  .2>rcCsdddlm}tttdtddtdt ks4tt|dtdtdd|ddt ks`tdS)NrrrETrm)rrr!rr0rr:rr@r@rAtest_linear_subs~s (rc$s.ddlmddlm}m}ddlm}ddlm}|ddd\}|d dd \}|d dd d d ddt t t t t fdddkst t t t t t t fddkst t t dt t t t fddddkst t t dt t t t fdddddks>t t t t ||t t t ft t t fdddks|t t t t t ||t t t ft t t fddkst t t t t ||t t t ft t t fdd|kst t t t t t ||t t t ft t t fdd|ksJt t t t dt t ||t t t ft t t fddd|kst t t t dt t ||t t t ft t t fddd|kst t t dt t ||t t t ft t t fdd}|tr:t t|ddksXt t t dt t ||t t t ft t t fdd|d|dkst t t t dt fdddkst t t t t dt fdddkst t t dt t dt fddddks.t fdd}|ddksVt |t t kspt |t t dddkst ddd}t|t t dd|t t dd|kst t|t t dd|t t dd|kst t|t t d|t t d|ksLt |ddd \} } t | ddt | | t| | t| t| } t | t dt fdddkst t t | t dt fddd}||d|df| | dd| kfkst t t d| t dt fddd} | d| dkks0t || d|dd| | d| | d| ddksrt |ddd \} }t | dt d|dt| |t| t|} tt | t ddfdddkst tt t | t ddfdd| | |ks t tt t d| t ddfdd| | d| || |dksRt tt t t | t ddfddt| |t| t t| t| |t kst |dddd}dd|dt |dtt d dt|d}|t |t dt fdddkst tt t |t dt fddtdt|ddt|dksXt tt t d|t dt fdd|kst d| dt|dt |ddtt d}|t |t dt fdddkst tt t |t dt fdd|k st tt t d|t dt fdd||dk s4t |t t |td|d|t dt fdddtdt|k s~t |ddd \} }}| |t t || |dt | || |d}tt |t dt fdddk st t |}tt |t dt fddd| |tdd| t|dd| | |dt|k sRt tt t |t dt fddd| |dtdd| t|dd| | |dt|k st |ddd \}}t|t ||||t |||t t|dt|dt||d}tt |t dt fdddk sDt tt t |t dt fddd||dk stt tt t d|t dt fddd|d|d||d |dk st |d!dd \}}t|dtt td"dt| t |dt d|d}d#d$}|t |t dt fdk s:t |t t |t dt f|k s\t |t t |d|t dt f|d|k st |t t |d|t dt fd|d%|dk st |d&dd }t t|dtt| dt |t |td't |t fdk s t |ddd \} } | | t | d| | ddt | | | d}tt |t dt fdk st tt t |t dt fdd(t| | tt| k st tt t t |t dt fdd(t| t t | ttt | ks t |ddd }|d)dd }||t ||dtt || }tt |t dt fdksjt tt t ||t dt f||td||kst dd*lm}|d+dd \} }!t |!dtt d| d d|!d|dt | |!d}"t |"t dt fdddks"t |d,dd-}|d.dd }ttt | |d|}#t |#t t t fdddksxt t t |#t t t fdd|kst t t d|#t t t fddd|d|dkst |ddd }|t tt t |dtt t|t dt ftd|ks*t dS)/Nr) expand_mul)rrj) gammasimp)powsimpzmu1 mu2TZnonzeroz sigma1 sigma2rklambdacSs6dtdt|dt||d d|dSNrEr9)rrr )r0rrr@r@rAnormalsz test_probability..normalcSs|t| |SrB)r )r0rater@r@rA exponentialsz%test_probability..exponentialrmrEr9rQrOcst|tttdtftt tfdd}t|tttt tftdtfdd}||ks|t|S)NrTrm)r!r0r1rr:)r=Zres1Zres2rrZmu1rrZsigma1r@rAEs  ztest_probability..Ez alpha betaZseparaterrtk)r{rlza b przd1 d2rPzlamda mucSst|tSrB)r%rewriter rCr@r@rAz"test_probability..rVr4r])rrd)rznu sigmarrrr3)r}rrorrjZsympy.simplify.gammasimprrrr!r0rr:r1rr r%rr rrrrrrrrrr)$rrjrrZmu2Zsigma2rransr|betaZbetadistjr2r3rchiZ chisquaredpZdagumr Zd1Zd2rZlamdardistZmysimpr4ZdistnrdrnurriceZlaplacer@rrAtest_probabilitys    & $$          $  $       $   .66.(.$ :"  . :"" *8"& & 8"* . .  "  " D"2: 6  "   *  D   $  *rc Csddlm}ddlm}ddlm}m}m}m}m }m }t t t tt ttttdtfddd|jdd |ttkstt tt tttdtfddd||dtkstt tt ttd tdtfddd||d t||kstt tt ttd tdtfddd||d t||ks`t|d dd }t tt tt|tfdd||kstt tt tt|tfdd||td kstt ttttdtfdd|tkstt |tttdtfdd|tks(tt tt ttdd|tt|dtksbtt tt td tdd||dttt tttkst|ddd} t t| | | dd| d|| kstt t| | | dd| d|| kstt |dttdd|t|dttt ksZtt |d ttdd|td  |dtd ttt d tt d kstt t t |tttdd|jdd |tdt kstt |ttddt|tttks"tt || | dd| || t| ksRtt |ttddt|tttks~tt || | dd| || || kstt |ttt tdtfddtdkstt |dttttdtfddtd d kstdS)z% Test various exponential integrals. rr)r)ChiCiEiShiSiexpintrETrr)funcr9rQr7rkrmrNrrPN)ror%sympy.functions.elementary.hyperbolicr'sympy.functions.special.error_functionsrrrrrrr%r r!r r8r0r1rrrpr:rrrrZas_independentrr) rrrrrrrrr7rNr@r@rA test_expintKs   "     0  **  $  " "  6  ,0,00rc Csddlm}m}ddlm}m}ddlm}ddlm }m }m }m }ddl m} m} | |tttdd|t td tddfkst| |tttdd|ttt td d kfkst| |tttddttd td d td d tt td d kfkst| |ttttd ddttd kttd k@fksJt| |d ttttd d} | djdd| d ftdtd d tktd d tkBfd tdtd td d dftdkfkstt|t|dttdtfddtd td kstt|t|d ttdtfddttjtd d ksNttd ttd td tddt|d t ttd d kft |d tdfkstdS)Nr)racoth)ratanr)rE1rr)fourier_transformlaplace_transformTr$r9rErxrOF)Znoconds)deeprSrm)!rrr(sympy.functions.elementary.trigonometricrrrrrrrrrrrrr0r6rr:rrr2r factorrprrr!rrTrr) rrrrrrrrrrrrr@r@rA test_messysL   &  $"      2rcCsJttt tdtt tfddt tttttdksFtdS)Nr9TrmrP)r!r rr0rrrr:r@r@r@rAtest_issue_6122s"rcCs>dtttttdd}t|tdd}|tr:tdS)NrErQTrm)r0r2r3rr!rrr:)r=Zantir@r@rAtest_issue_6252srcCsDttttdtdtt tftttdks@tdSrv) r!r rr0rr%rrr:r@r@r@rAtest_issue_6348s. rcCshddlm}m}tttttddt|tks:tttt ttddt|tksdtdS)Nrfresnelcfresnelsr9) rrrrr!rrr0r:rrr@r@rA test_fresnels*rcCstttttdkstdSrB)r(r0r:r@r@r@rAtest_issue_6860srcCsrtttdtdt}|tdtddtdtddksNt|ttjtjt ddksntdS)Nr9rQrV) r(r0rZtogetherr:Z_eval_intervalrZ NegativeOneZOnerrr@r@rAtest_issue_7337s0rcCshttttt ttt dttt dtt tt tdtddksdtdSr)r(rr0r r7r:r@r@r@rAtest_issue_8368s$  rcCstddlm}m}tdtttd|ddtd|ftd|fdtd|d|d|d|ksptdS)NrhwrEr9rQ) sympy.abcrrr!rr1r0r:rr@r@rAtest_issue_10211s2&rcCsrddlm}|ddd\}}tdttd|ddt| |fd||dt|d|dksntdS) Nrrizy LTrkrEr9rQ)rorjr!rr0r:)rjr1Lr@r@rAtest_issue_11806s  ("rcCsddlm}ddlm}m}t|d|d|dd|dd}d|d |d td td dftd df|dtdt t |d}| || d d st dS)Nr)RR)RrMr9g?TrmgUUUUUU?g?rQgrVg-q=)Zsympy.polys.domains.realfieldrrrrMr!rrrrrZalmosteqrdr:)rrrMrr\r@r@rAtest_issue_10681s &*rcCs@ddlm}|ddd}tdtdtt|fd|kssr      $ ,) j0   G 6&