From 5d80ae02f10cb7350925d449e9b5bfe1dc9a712e Mon Sep 17 00:00:00 2001 From: Yannick Ulrich Date: Mon, 30 Sep 2019 11:41:12 +0200 Subject: [PATCH] More tests --- checks/test.f90 | 33 +++++++++++++++++++++++++++++++++ 1 file changed, 33 insertions(+) diff --git a/checks/test.f90 b/checks/test.f90 index 8932562..6206696 100644 --- a/checks/test.f90 +++ b/checks/test.f90 @@ -388,6 +388,39 @@ CONTAINS 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') + + ! 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 -- GitLab