Commit 4b18fe2b authored by ulrich_y's avatar ulrich_y

Restructed and added to F tests

parent c830b6d1
......@@ -330,50 +330,59 @@ CONTAINS
! Thanks to Roman Zwicky and Ben Pullin for these tests
ref = (0.392707112217551702879328061598355445_prec, - 1.274969448494380061814943180491890080_prec)
print*, ' ', 'testing GPL ', 'F.1', ' ...'
res = G( (/ (0._prec,0.), (1._prec,0.), (0._prec, 1.5_prec) /) )
call check(res,ref)
call test_one_flat((/ (0._prec,0.), (1._prec,0.), (0._prec, 1.5_prec) /),ref,'F.1')
ref = (2.09624324167194065961839363660174566_prec, 0.62644605179087454819421905065313983_prec)
print*, ' ', 'testing GPL ', 'F.2', ' ...'
res = G( (/ (1._prec,0.), (-1._prec,0.), (2.5_prec, -2.4_prec) /) )
call check(res,ref)
call test_one_flat( (/ (1._prec,0.), (-1._prec,0.), (2.5_prec, -2.4_prec) /) , ref, 'F.2')
ref = conjg(ref)
print*, ' ', 'testing GPL ', 'F.3', ' ...'
res = G( (/ (1._prec,0.), (-1._prec,0.), (2.5_prec, +2.4_prec) /) )
call check(res,ref)
call test_one_flat( (/ (1._prec,0.), (-1._prec,0.), (2.5_prec, +2.4_prec) /) , ref, 'F.3')
ref = (0.6874619495289224183167486785286777066_prec, -1.8261934916106546308783576818077830287_prec)
print*, ' ', 'testing GPL ', 'F.4', ' ...'
res = G( (/ (-1._prec,0.), (1._prec,0.), (2.5_prec, +2.4_prec) /) )
call check(res,ref)
call test_one_flat( (/ (-1._prec,0.), (1._prec,0.), (2.5_prec, +2.4_prec) /) , ref, 'F.4')
ref = (0.320016770069038023050391296154549_prec , 0.064263450879286017225353602844105_prec)
print*, ' ', 'testing GPL ', 'F.5', ' ...'
res = G( (/ (0._prec,0.), (1._prec,0.), (-1._prec,0.), (0.3_prec, -1.2_prec) /) )
call check(res,ref)
call test_one_flat( (/ (0._prec,0.), (1._prec,0.), (-1._prec,0.), (0.3_prec, -1.2_prec) /) , ref, 'F.5')
ref = conjg(ref)
print*, ' ', 'testing GPL ', 'F.6', ' ...'
res = G( (/ (0._prec,0.), (1._prec,0.), (-1._prec,0.), (0.3_prec, +1.2_prec) /) )
call check(res,ref)
call test_one_flat( (/ (0._prec,0.), (1._prec,0.), (-1._prec,0.), (0.3_prec, +1.2_prec) /) , ref, 'F.6')
ref = (0.8382358254435272068734922352_prec, - 0.3702062785327198487992149341_prec)
print*, ' ', 'testing GPL ', 'F.7', ' ...'
res = G( (/ (0._prec,0.), (1._prec,0.), (-1._prec,0.), (1.1_prec,2._prec) /) )
call check(res,ref)
call test_one_flat( (/ (0._prec,0.), (1._prec,0.), (-1._prec,0.), (1.1_prec,2._prec) /) , ref, 'F.7')
ref = (0.185156872427485220072774923047908301422_prec,-0.249989197161429773744237773045427322847_prec)
print*, ' ', 'testing GPL ', 'F.8', ' ...'
res = G( (/ (1._prec,0.), (-1._prec,0.), (2.3_prec,-1._prec), (1.1_prec,-0.1_prec) /) )
call check(res,ref)
call test_one_flat( (/ (1._prec,0.), (-1._prec,0.), (2.3_prec,-1._prec), (1.1_prec,-0.1_prec) /) , ref, 'F.8')
ref = (-1.11161035333074623447215317094270858897_prec,-0.875225967273157437459269290416981432294_prec)
print*, ' ', 'testing GPL ', 'F.9', ' ...'
res = G( (/ (1._prec,0.), (1.1_prec,-0.1_prec), (1._prec,0.), (2.3_prec,-1.2_prec), (3.4_prec,9.7_prec) /) )
call check(res,ref)
call test_one_flat( (/ (1._prec,0.), (1.1_prec,-0.1_prec), (1._prec,0.), (2.3_prec,-1.2_prec), (3.4_prec,9.7_prec) /) , ref, 'F.9')
ref = (-0.1109203962012759970680116512959180346924_prec,0.114779444698808046532782358409267479097_prec)
call test_one_flat( (/ (0._prec, 0.) , (-0.1_prec, - 0.3_prec), (1._prec,0.), (0.5_prec, 0._prec) /) , ref, 'F.10')
ref = (-0.108930339927322068475622176344923224046_prec,0.1122893584726586134006078500824852031925_prec)
call test_one_flat( (/ (-0.005_prec, 0.) , (-0.1_prec, - 0.3_prec), (1._prec,0.), (0.5_prec, 0._prec)/) , ref, 'F.11')
ref = (-1.00642162671475190543135817048062920163_prec,0.75972801214517468120067717498723544972_prec)
call test_one_flat( (/ (-0.005_prec, 0.) , (-0.1_prec, - 0.3_prec), (1._prec,0.), (1.5_prec, 0._prec)/) , ref, 'F.12')
ref = (-0.764543307991553066235740993950606803628_prec, 0.54861507469010653128533842086292369104_prec)
call test_one_flat( (/ (-0.1_prec, - 0.3_prec), (1._prec,0.), (-0.005_prec, 0.) , (0.5_prec, 0._prec)/) , ref, 'F.13')
ref = (4.52004218277686921073832347986485146573_prec,-1.812384329980915889835338149845575221705_prec)
call test_one_flat( (/ (1.005_prec,0.), (0._prec,0.), (1.1_prec,0.3_prec), (1._prec,0.) /) , ref, 'F.14')
ref = (1.6706939608077788063504210821631008732_prec,-0.631820518538505729899318406678510198266_prec)
call test_one_flat( (/ (1.1_prec,0.), (0._prec,0.), (1.1_prec,0.3_prec), (1._prec,0.) /) , ref, 'F.10')
ref = (0.198693810038777868937967610639933071867_prec,+0.72807783717224553483870820314753060868_prec)
call test_one_flat( (/ (0._prec,0.), (1.0_prec,0._prec), (1.1_prec,0.), (1.1_prec,+0.3_prec) /) , ref, 'F.11')
ref = (1.6706939608077788063504210821631008732_prec,-0.631820518538505729899318406678510198266_prec)
call test_one_flat( (/ (1.1_prec,0.), (0._prec,0._prec), (1.1_prec,+0.3_prec), (1._prec,0.) /) , ref, 'F.12')
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')
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