Commit ff260ddb authored by ulrich_y's avatar ulrich_y

more polylog tests

parent 7592cc50
......@@ -143,43 +143,85 @@ CONTAINS
implicit none
complex(kind=prec) :: ref, z
z = 0.9999_prec * exp((0.,1._prec))
ref = ( 0.324133535392567026495406418480254855764_prec,1.013852051950427381643871737213288702950_prec)
call check( polylog( 2, z), ref)
ref = ( 0.522770725155376044565670599443640187194_prec,0.869361155897677036715122800524769654185_prec)
call check( polylog( 5, z), ref)
ref = ( 0.539824626930596495579386730660957051319_prec,0.842276214103123630233317308533693430239_prec)
z = 0.9_prec * exp((0._prec,1._prec))
ref = ( 0.31714670838088083420005254335964850081059000145548_prec,0.90620928887768788007657902653829197810196658407346_prec)
call check( polylog(2, z), ref)
ref = ( 0.47248076790053852605402172732597392822155906295748_prec,0.78010930007060002834519910647340859560884976749326_prec)
call check( polylog(5, z), ref)
ref = ( 0.48593029219713200876281747476986815689740323314048_prec,0.75804436457088007970152497542119942774917319813883_prec)
call check( polylog(10, z), ref)
z = 0.99999_prec * exp((0.,1._prec))
ref = ( 0.324137319831170505532582434237087014236_prec,1.013948424389723032721958295757400129230_prec)
call check( polylog( 2, z), ref)
ref = ( 0.522815799416934255669781473407124697155_prec,0.869441778136187554517193738781550493266_prec)
call check( polylog( 5, z), ref)
z = 1.00001_prec * exp((0.,1._prec))
ref = ( 0.324138160221287178367712545364630311917_prec,1.013969840316258729274699928974034201521_prec)
call check( polylog( 2, z), ref)
ref = ( 0.522825815862029305845720413414469502401_prec,0.869459694240961134437078365904557794494_prec)
call check( polylog( 5, z), ref)
z = 0.999999_prec * exp((0.,1._prec))
ref = ( 0.324137698033552982036491149170744105494_prec,1.013958061564363933265675662561552413453_prec)
call check( polylog( 2, z), ref)
ref = ( 0.522820306819813365669890959221612196007_prec,0.869449840381005997194835712531309011624_prec)
call check( polylog( 5, z), ref)
z = 0.9999999_prec * exp((0.,1._prec))
ref = ( 0.324137735851376524606585447637048722447_prec,1.013959025281135047042683412228827082409_prec)
call check( polylog( 2, z), ref)
ref = ( 0.522820757559868506601039369765883468501_prec,0.869450646605697511953904345069485267937_prec)
call check( polylog( 5, z), ref)
z = 0.9999999_prec * exp((0.,2._prec))
ref = (-0.496658534687223270187711312384211419797_prec,0.727145993783645319190052450531069102550_prec)
call check( polylog( 2, z), ref)
ref = (-0.432933763639004718041881713834580186629_prec,0.885265176678646954693047624423223092246_prec)
call check( polylog( 5, z), ref)
z = 0.99_prec * exp((0._prec,1._prec))
ref = ( 0.32369036161935845024104448122222955899340391857079_prec,1.00324349284774627379699444852820582088595888924728_prec)
call check( polylog(2, z), ref)
ref = ( 0.51780997596604204699322616691391794080129634794226_prec,0.86049504083221113347441620727610317572983128003060_prec)
call check( polylog(5, z), ref)
ref = ( 0.53448415482661295991549353367903775505364004919109_prec,0.83392812071435332971141463249398007205963468167980_prec)
call check( polylog(10, z), ref)
z = 0.999_prec * exp((0._prec,1._prec))
ref = ( 0.32409544945505588451045078881514007688743827135543_prec,1.01288825835849929172881742864723514784379050415565_prec)
call check( polylog(2, z), ref)
ref = ( 0.52231995926604156766783810033494360274235745715174_prec,0.86855495448337982380397597048602597837498452189533_prec)
call check( polylog(5, z), ref)
ref = ( 0.53933913318980396043645098080232673436372398744714_prec,0.84151728931442886718138858740874188686013946305784_prec)
call check( polylog(10, z), ref)
z = 0.9999_prec * exp((0._prec,1._prec))
ref = ( 0.32413353539256702649540641848025485576423369942320_prec,1.01385205195042738164387173721328870295049412205747_prec)
call check( polylog(2, z), ref)
ref = ( 0.52277072515537604456567059944364018719426776187131_prec,0.86936115589767703671512280052476965418527852729045_prec)
call check( polylog(5, z), ref)
ref = ( 0.53982462693059649557938673066095705131895675203284_prec,0.84227621410312363023331730853369343023972319744457_prec)
call check( polylog(10, z), ref)
z = 0.99999_prec * exp((0._prec,1._prec))
ref = ( 0.32413731983117050553258243423708701423605015249363_prec,1.01394842438972303272195829575740012923040998957878_prec)
call check( polylog(2, z), ref)
ref = ( 0.52281579941693425566978147340712469715561075153075_prec,0.86944177813618755451719373878155049326662610449505_prec)
call check( polylog(5, z), ref)
ref = ( 0.53987317626370393558176605497595292698896795484936_prec,0.84235210666127870713197117713665972968191528367673_prec)
call check( polylog(10, z), ref)
z = 1.00001_prec * exp((0._prec,1._prec))
ref = ( 0.32413816022128717836771254536463031191720755084579_prec,1.01396984031625872927469992897403420152117604112105_prec)
call check( polylog(2, z), ref)
ref = ( 0.52282581586202930584572041341446950240136278312399_prec,0.86945969424096113443707836590455779449482967107881_prec)
call check( polylog(5, z), ref)
ref = ( 0.53988396500338278677493363054405590339711449183908_prec,0.84236897167615972350483622438537793816611895748722_prec)
call check( polylog(10, z), ref)
z = 0.999999_prec * exp((0._prec,1._prec))
ref = ( 0.32413769803355298203649114917074410549435060587269_prec,1.01395806156436393326567566256155241345384396358214_prec)
call check( polylog(2, z), ref)
ref = ( 0.52282030681981336566989095922161219600778245313269_prec,0.86944984038100599719483571253130901162490098836993_prec)
call check( polylog(5, z), ref)
ref = ( 0.53987803119660494489725285096895119223918936255085_prec,0.84235969591788706954609196967489789002490184573847_prec)
call check( polylog(10, z), ref)
z = 0.9_prec * exp((0._prec,2._prec))
ref = (-0.44452028964890283007349990328300389929008760996086_prec,0.66878743447993535281130919066420953251950351295302_prec)
call check( polylog(2, z), ref)
ref = (-0.38839535801823932606323111787411257099116023643543_prec,0.79888940673536917702425494564670039829270801271179_prec)
call check( polylog(5, z), ref)
ref = (-0.37503747568892401565030939135002191537632464837044_prec,0.81776617491393185281903993184823060507449761681675_prec)
call check( polylog(10, z), ref)
z = 0.99_prec * exp((0._prec,2._prec))
ref = (-0.49145214297542995287792244783916287476339763066814_prec,0.72142557030218510441094486726609727591395278662659_prec)
call check( polylog(2, z), ref)
ref = (-0.42847005378062827950272801697834832300896766072327_prec,0.87664867860304046775940447885176343497500132376663_prec)
call check( polylog(5, z), ref)
ref = (-0.41259541382394926678059535801164943627251464211660_prec,0.89947635154294082643326725592105204892079432293970_prec)
call check( polylog(10, z), ref)
z = 0.999999_prec * exp((0._prec,2._prec))
ref = (-0.49665806619812223943709229668650064888627387195443_prec,0.72714548006682757749406286437255854455016840984651_prec)
call check( polylog(2, z), ref)
ref = (-0.43293336180419752186586468054516010134686757646849_prec,0.88526440139529870422603576510716149051864398952584_prec)
call check( polylog(5, z), ref)
ref = (-0.41676869338400351335064877272829801919405402540781_prec,0.90855360420316407253701176433346888385813040933393_prec)
call check( polylog(10, z), ref)
end subroutine
subroutine do_MPL_tests()
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment