Commit 1c6d8bd1 authored by ulrich_y's avatar ulrich_y

Added tests for new plogs

parent e85716b5
......@@ -243,8 +243,85 @@ CONTAINS
res = G((/ 0._prec, 1.3_prec, 0._prec, 0._prec, 4._prec, 1.1_prec /))
call check(res,ref)
print*, ' ', 'testing GPL ', '5.13', ' ...'
ref = (-1.229846637025798984989235266565630713100_prec,-1.05705870670612913651142424798975000765_prec)
res = G([inum(2._prec, +1), inum(7._prec, +1)], inum(5._prec, +1))
call check(res,ref)
print*, ' ', 'testing GPL ', '5.14', ' ...'
ref = (-1.229846637025798984989235266565630713100_prec,+1.05705870670612913651142424798975000765_prec)
res = G([inum(2._prec, -1), inum(7._prec, +1)], inum(5._prec, +1))
call check(res,ref)
print*, ' ', 'testing GPL ', '5.15', ' ...'
ref = (0.2982254208688088675254638762780704094718_prec,0.)
res = G([inum(2._prec, -1), inum(7._prec, +1)], inum(-5._prec, -1))
call check(res,ref)
print*, ' ', 'testing GPL ', '5.17', ' ...'
ref = (0.190800137777535619036913153766083992418_prec, 0.)
res = G((/ 6, 1, 1 /))
call check(res,ref)
print*, ' ', 'testing GPL ', '5.18', ' ...'
ref = (0.058192342415778512650117048874978455691_prec, 0.)
res = G((/ 6, 1, -1 /))
call check(res,ref)
ref = cmplx(log(0.5_prec), pi,kind=prec)
print*, ' ', 'testing GPL ', '5.19a', ' ...'
res = G((/ inum(2._prec, +1) /), inum(3._prec, +1))
call check(res,ref)
print*, ' ', 'testing GPL ', '5.19b', ' ...'
res = G((/ inum(2._prec, -1) /), inum(3._prec, +1))
call check(res,conjg(ref))
print*, ' ', 'testing GPL ', '5.19c', ' ...'
res = G((/ inum(2._prec, +1) /), inum(3._prec, -1))
call check(res,ref)
print*, ' ', 'testing GPL ', '5.19d', ' ...'
res = G((/ inum(2._prec, -1) /), inum(3._prec, -1))
call check(res,conjg(ref))
print*, ' ', 'testing GPL ', '5.19e', ' ...'
res = G((/ inum(-2._prec, +1) /), inum(-3._prec, +1))
call check(res,conjg(ref))
print*, ' ', 'testing GPL ', '5.19f', ' ...'
res = G((/ inum(-2._prec, -1) /), inum(-3._prec, +1))
call check(res,ref)
print*, ' ', 'testing GPL ', '5.19g', ' ...'
res = G((/ inum(-2._prec, +1) /), inum(-3._prec, -1))
call check(res,conjg(ref))
print*, ' ', 'testing GPL ', '5.19h', ' ...'
res = G((/ inum(-2._prec, -1) /), inum(-3._prec, -1))
call check(res,ref)
ref = -(2.374395270272480200677499763071638424_prec, - 1.273806204919600530933131685580471698_prec)
print*, ' ', 'testing GPL ', '5.20a', ' ...'
res = G((/ izero, inum(2._prec, +1) /), inum(3._prec, +1))
call check(res,ref)
print*, ' ', 'testing GPL ', '5.20b', ' ...'
res = G((/ izero, inum(2._prec, -1) /), inum(3._prec, +1))
call check(res,conjg(ref))
print*, ' ', 'testing GPL ', '5.20c', ' ...'
res = G((/ izero, inum(2._prec, +1) /), inum(3._prec, -1))
call check(res,ref)
print*, ' ', 'testing GPL ', '5.20d', ' ...'
res = G((/ izero, inum(2._prec, -1) /), inum(3._prec, -1))
call check(res,conjg(ref))
print*, ' ', 'testing GPL ', '5.20e', ' ...'
res = G((/ izero, inum(-2._prec, +1) /), inum(-3._prec, +1))
call check(res,conjg(ref))
print*, ' ', 'testing GPL ', '5.20f', ' ...'
res = G((/ izero, inum(-2._prec, -1) /), inum(-3._prec, +1))
call check(res,ref)
print*, ' ', 'testing GPL ', '5.20g', ' ...'
res = G((/ izero, inum(-2._prec, +1) /), inum(-3._prec, -1))
call check(res,conjg(ref))
print*, ' ', 'testing GPL ', '5.20h', ' ...'
res = G((/ izero, inum(-2._prec, -1) /), inum(-3._prec, -1))
call check(res,ref)
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