Commit 7170d56e authored by ulrich_y's avatar ulrich_y

Added more digits to MPL tests

The numbers used in ac315e3e are not reproducible,
used something close instead and re-computed
reference
parent 13945806
......@@ -18,7 +18,7 @@ PROGRAM TEST
call parse_cmd_args()
#endif
tol = 8e-10
tol = zero * 1.e5_prec
call do_MPL_tests()
call do_GPL_tests()
! call do_shuffle_tests() ! put this somewhere else
......@@ -68,19 +68,20 @@ CONTAINS
complex(kind=prec) :: ref
print*, 'doing MPL tests...'
ref = (0.022696600480693277651633_prec,0.)
call test_one_MPL((/ 1,1 /),(/ (0.3156498673740053_prec,0.), (0.3431255827785649_prec,0.) /),ref, '1.1')
ref = (0.0226974036143118403367048874190211851348_prec,0.)
call test_one_MPL((/ 1,1 /),cmplx([ 119._prec/377._prec, 35._prec / 102._prec ],kind=prec),ref, '1.1')
ref = (2.31346156303083313941538625891818239608127E-4_prec,0.)
call test_one_MPL((/ 1,1 /),cmplx([0.03_prec, 100._prec-30._prec*sqrt(11._prec) ], kind=prec),ref, '1.2')
ref = (0.00023134615630308335448329926098409_prec,0.)
call test_one_MPL((/ 1,1 /),(/ (0.03_prec,0.), (0.5012562893380046_prec,0.) /),ref, '1.2')
ref = (2.34617486908178608355021488772249510275409004880846912378356e-5_prec,0.)
call test_one_MPL((/ 2,1,2 /),cmplx([0.03_prec, 100._prec-30._prec*sqrt(11._prec), 250._prec/9._prec + sqrt(6875._prec)/3._prec ], kind=prec),ref, '1.3')
ref = (0.000023446106415452030937059124671151_prec,0.)
call test_one_MPL((/ 2,1,2 /),(/ (0.03_prec,0.), (0.5012562893380046_prec,0.), (55.3832_prec,0.) /),ref, '1.3')
ref = (-0.06565799418838388900030511939327151809905663800733251335678_prec,0.)
call test_one_MPL((/1, 1/), (/(-0.25_prec,0.),(-2._prec,0.) /), ref, '1.4')
ref = (-0.0319989639656484125354488353677231989818697987217232462808_prec,0.)
call test_one_MPL((/2, 1/), (/(-0.25_prec,0.),(-2._prec,0.) /), ref, '1.4')
call test_one_MPL((/2, 1/), (/(-0.25_prec,0.),(-2._prec,0.) /), ref, '1.5')
end subroutine do_MPL_tests
subroutine test_one_condensed(m,z,y,k,ref,test_id)
......
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