Commit 01ef395b authored by Luca Naterop's avatar Luca Naterop

some more tests

parent a3fbb738
......@@ -92,19 +92,29 @@ CONTAINS
complex(kind=prec), parameter :: epsilon = 1E-14
print*, 'doing GPL tests...'
ref = dcmplx(0.0819393734128676)
ref = cmplx(0.0819393734128676)
call test_one_condensed((/ 1,1 /),cmplx((/ 1.3d0, 1.1d0 /)),cmplx(0.4),2,ref,'2.1')
ref = dcmplx(0.01592795952537145)
ref = cmplx(0.01592795952537145)
call test_one_condensed((/ 3,2 /),cmplx((/ 1.3d0, 1.1d0 /)),cmplx(0.4),2,ref,'2.2')
ref = dcmplx(0.0020332632172573974)
ref = cmplx(0.0020332632172573974)
call test_one_condensed((/ 4 /),cmplx((/ 0 /)),cmplx(1.6),1,ref,'2.3')
! requires making convergent
ref = cmplx((0.09593041677639341, -0.8829351795197851))
call test_one_flat(cmplx([0,1,3,2]),ref,'2.4')
end subroutine do_GPL_tests
ref = cmplx((0.009947947789928621,0.0))
call test_one_flat(cmplx([0.0, 0.0, -3.3333333333333335, -3.3333333333333335, 1.0]),ref,'2.5')
! requires hoelder convolution
ref = cmplx((-0.012709942828250949,0.0))
call test_one_flat(cmplx([0.0, 3.3333333333333335, 1.0, 3.3333333333333335, 1.0]),ref,'2.5')
end subroutine do_GPL_tests
subroutine do_shuffle_tests()
complex(kind=prec) :: v(2) = cmplx((/1,2/))
complex(kind=prec) :: w(2) = cmplx((/3,4/))
......
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