Commit 5d80ae02 authored by ulrich_y's avatar ulrich_y

More tests

parent ab33be23
...@@ -388,6 +388,39 @@ CONTAINS ...@@ -388,6 +388,39 @@ CONTAINS
ref = (0.1810339553848393655117844582129810543006_prec,+0.82543851493400141056822586864697825094_prec) ref = (0.1810339553848393655117844582129810543006_prec,+0.82543851493400141056822586864697825094_prec)
call test_one_flat( (/ (0._prec,0.), (1.0_prec,0._prec), (1._prec,0.), (1.1_prec,+0.3_prec) /) , ref, 'F.13') call test_one_flat( (/ (0._prec,0.), (1.0_prec,0._prec), (1._prec,0.), (1.1_prec,+0.3_prec) /) , ref, 'F.13')
! This one converges really slowly!
! TODO
call set_options(liinf=200000)
tol = 2e-9_prec
ref = (-3.028608056170828558746818459366566689807_prec,-0.52999686156911157999452896083600882192_prec)
call test_one_flat( (/ (1._prec,0._prec), (0.0_prec,0._prec), (0._prec,1._prec), (1.1_prec,0.0_prec) /) , ref, 'F.14')
! Here the branch cut matters, in Mathematica this is entered as G(1_-, 0, -I, -1.1)
ref = (0.0538179677556747874824671812943783830465_prec,0.340519960719077282936661573786600733497_prec)
print*, ' ', 'testing GPL ', 'F.15', ' ...'
res = G((/ inum(1._prec,-1), izero, inum((0._prec,-1._prec), -1) /), inum(-1.1_prec, +1))
call check(res,ref)
call set_options(liinf=1000)
tol = zero * 1.e5_prec
ref = (0.0075181720252360930529649927100295277234_prec,0.009279206321129807482628652573319115651_prec)
call test_one_flat( (/ (0._prec,0.), (0._prec,0.), (1.0_prec,0._prec), (0.1_prec,-0.1_prec), (0.11_prec,0._prec) /) , ref, 'F.16')
ref = (0.0135493735561310633082925053080300174678_prec,+0.01851696361979639851211857931403696163_prec)
call test_one_flat( (/ (0._prec,0.), (0._prec,0.), (1.0_prec,0._prec), (0.1_prec,-0.1_prec), (0.15_prec,0._prec) /) , ref, 'F.17')
ref = (-2.41176903997917647290181947143140323970805e-4_prec,0.00129690184788114442363777169655122)
call test_one_flat( (/ (0._prec,0.), (1._prec,0.), (0.5_prec,-0.1_prec), (-1.3_prec,10.2_prec), (0.4_prec,0._prec) /) , ref, 'F.18')
ref = (-0.003113755848404291707649322614093090379815_prec,-4.6914973893242872720303973859400498154E-4_prec)
call test_one_flat( (/ (0._prec,0.), (1._prec,0.), (0.1_prec,-0.1_prec), (-1.3_prec,10.2_prec), (0.4_prec,0._prec) /) , ref, 'F.19')
end subroutine do_GPL_tests end subroutine do_GPL_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