U Zʹj<?@sddlTddlmZddlZdejkrBeejejddZndZdejkZda dd Z d d Z gZ e edd e ed ddD]Ze dedgqedd D](Ze dedge dedgqe ddgddgddgddgddgddgddgddgddgddgddgddgd dgd!dgd"dgd#dgd$dgd%dgd&dgd'dgd$dgd(dgd)dgd*dgd&dgd+dgd,dgd-d.gd/d.gd0d1gd2d3gd4d5gd6d7gd8d9gd:d;gdgd?dgd@dgdAdgdBdgdCdgdDdgdEdgdFdGgdHdIgdJdKgdLdMgdNdgdOdgdPdgdQdgdRdgdSdgdRdgdTdgdUdgdVdgdWdgdVdgdXdYgg=7Z dZD]Ze d[edge d\edge d]edge d^edge d_edge d`edge daedge dbedge dcedge ddedge deedge dfedgqLdgD]jZdhD]^Ze dieefdge djeefdge dkeefdge dleefdgq8q0edmkr ddnlmZdoZe D]:ZeZe eeZedpeeeeee7Zqedqeedrt dS)s)*)ifacNz-dpsz-raisec Cs^z ||}WnRtd7atr"tt|tdddl}|tdtYdSX|j}|j}|j}|j} t d} d} t ||t || krd} ttd|t |dt |dt t ||ftd7atrt t || t | | krDd} ttd|t |dt | dt t || ftd7atrDt | sZtj d |dS) NrZ EXCEPTIONrz&Error! %s (re = %s, wanted %s, err=%s) z&Error! %s (im = %s, wanted %s, err=%s)z%s ok; )errcountraise_print tracebackprint_tbsysexc_inforealimagepsabsZnstr SystemExitstdoutwrite) namefunczyxr ZxreZximZyreZyimZtolerrr[/var/www/html/bwave/python/venv/lib/python3.8/site-packages/mpmath/tests/extratest_gamma.pychecksF ,,rcCs|\}}td|dt_t|}tdt_|dkrTt|}t|}t|}t|}n t|}t |}||}d|}dD]f}|tkrq|t_td|t dt||t dt||t d t||t d t||td t_qxdS) Nz Testing z =i2r) r(<ZxiXriizdps = %sgammargammaloggamma factorialr") r mpdpsevalmaxdpsr)r+r,r*expr)caserresultZ gamma_valZ loggamma_valZ factorial_valZ rgamma_valr.rrrtestcase4s2    r4iz%s8z%s+0.5z%s+0.37z(0.1+1j)z (-0.1+1j)z(0.1-1j)z (-0.1-1j)Z10jz-10jZ100jZ10000jz -10000000jz (10**100)*jz125+(10**100)*jz-125+(10**100)*jz(10**10)*(1+j)z(10**10)*(-1+j)z(10**100)*(1+j)z(10**100)*(-1+j)z(1.5-1j)z(6+4j)z(4+1j)z(3.5+2j)z(-6-4j)z(-2-3j)z (-2.5-2j)z(3+3j)z(2-2j)1023zlog(2)4zlog(6)5zlog(24)z0.5z log(pi)/2z1.5zlog(sqrt(pi)/2)z2.5zlog(3*sqrt(pi)/4)z mpf('0.37')z0.25z)log(sqrt(2*sqrt(2*pi**3)/agm(1,sqrt(2))))z-0.4z mpf('-1.9')z mpf('12.8')z mpf('33.7')z mpf('95.2')z mpf('160.3')z mpf('2057.8')Z25z log(ifac(24))80z log(ifac(79))500zlog(ifac(500-1))Z8000zlog(ifac(8000-1))z8000.5z mpf('8000.1')zmpf('1.37e10')zmpf('1.37e10')*(1+j)zmpf('1.37e10')*(-1+j)zmpf('1.37e10')*(-1-j)zmpf('1.37e100')zmpf('1.37e100')*(1+j)zmpf('1.37e100')*(-1+j)zmpf('1.37e100')*(-1-j)z3+4jampc('-1.7566267846037841105306041816232757851567066070613445016197619371316057169472361826396083480461846305298860734828967253578064447068977111523651210600259708734715632405373076389685095561916961679704883904239638670319343338908388009531786948197210025029725361069435208930363494971027388382086721660805397916323064321605458016797620170995150951921863546031736733861250062671478363174983174780484475256740163443225458588326103258610863362045919438223029718235161814175530618223688296232894588415495615809337292518431903058265147109853171056894218498782764388681620045286085387381541336752982963143014622747051765799672222008686321794822143126731612769761171322046332838061619713895191371243359764435612951384238091232760634271570950240717650166551484551654327989936028503008171693413044615024511055703811707517257682549003543406938864812466781522545540015867361207626414225907787661003765157377139385212757490499491284143906816424244705094759339932733567910991920631339597278805393743140853391550313363278558195609260225928','4.74266443803465792819488940755002274088830335171164611359052405215840070271590681300937317113976705186354250813687568855081767037900279030487082277549828099966758775645041925653923672591196104389515931289826469459903721798606134294436498090428077839141927485901735557543641049637962003652638924845391650954629013775555010722490760652938524839066763429718336190205584222879898420095911804502117983417158744776290996876098194664579906420307070808945181689246805549314043258530272479246115112769957368212585759640878745385160943755234939803677490810820437032389675754312185365002552976365531236035424489891346371159557028288389233931136182050741628120897320644145308130874835332032440560546577484241423134079056537777170351934430586103623577814746004431994179990531852293907799261385520580149820193022197572124649872089512234542069845198000512157973103058858459643347618317513706729969847568154109777507997488135638784405288158432214886648743541773208808731479748217023665577802702269468013673719173759245720489020315779001'))"@z (2+j)*%s/3z (-2+j)*%s/3z (1+2*j)*%s/3z (2-j)*%s/3z (20+j)*%s/3z (-20+j)*%s/3z (1+20*j)*%s/3z (20-j)*%s/3z (200+j)*%s/3z (-200+j)*%s/3z(1+200*j)*%s/3z (200-j)*%s/3)rrrrAr#iiiir8i7ii߱)z1e-5z1e-20z1e-100z1e-10000zfadd(%s,'%s',exact=True)zfsub(%s,'%s',exact=True)zfadd(%s,'%sj',exact=True)zfsub(%s,'%sj',exact=True)__main__) default_timergz Test time:z Total time:zErrors:)ZmpmathZ mpmath.libmprrargvintindexr0r r rr4Z testcaseslistrangenappendrt__name__ZtimeitrKclockZtot_timer2t1t2r rrrrs   &"]