generalize log_gamma, log_factorial tests

Author: perazz

test/specialfunctions/test_specialfunctions_gamma.fypp

@@ -80,7 +80,7 @@ contains
 
 
     #:for k1, t1 in INT_KINDS_TYPES
-
+    #:set k2, t2 = REAL_KINDS[-2], REAL_TYPES[-2]
     subroutine test_logfact_${t1[0]}$${k1}$(error)
         type(error_type), allocatable, intent(out) :: error
         integer, parameter :: n = 6
@@ -90,32 +90,32 @@ contains
 
         ${t1}$, parameter :: x(n) = [0_${k1}$, 1_${k1}$, 2_${k1}$, 4_${k1}$,   &
                                      5_${k1}$, 100_${k1}$]
-        real(dp), parameter :: ans(n) = [0.0_dp, 0.0_dp, 0.69314718055994_dp, 3.17805383034794_dp,    &
-                                         4.78749174278204_dp, 3.637393755555e2_dp]
+        ${t2}$, parameter :: ans(n) = [0.0_${k2}$, 0.0_${k2}$, 0.69314718055994_${k2}$, 3.17805383034794_${k2}$,    &
+                                         4.78749174278204_${k2}$, 3.637393755555e2_${k2}$]
 
         #:elif k1 == "int16"
 
         ${t1}$, parameter :: x(n) = [0_${k1}$, 1_${k1}$, 2_${k1}$, 4_${k1}$,   &
                                      7_${k1}$, 500_${k1}$]
-        real(dp), parameter :: ans(n) = [0.0_dp, 0.0_dp, 0.69314718055994_dp, 3.17805383034794_dp,    &
-                                         8.52516136106541_dp, 2.611330458460e3_dp]
+        ${t2}$, parameter :: ans(n) = [0.0_${k2}$, 0.0_${k2}$, 0.69314718055994_${k2}$, 3.17805383034794_${k2}$,    &
+                                         8.52516136106541_${k2}$, 2.611330458460e3_${k2}$]
         #:elif k1 == "int32"
 
         ${t1}$, parameter :: x(n) = [0_${k1}$, 1_${k1}$, 2_${k1}$, 4_${k1}$,   &
                                      12_${k1}$, 7000_${k1}$]
-        real(dp), parameter :: ans(n) = [0.0_dp, 0.0_dp, 0.69314718055994_dp, 3.17805383034794_dp,    &
-                                         1.998721449566e1_dp, 5.498100377941e4_dp]
+        ${t2}$, parameter :: ans(n) = [0.0_${k2}$, 0.0_${k2}$, 0.69314718055994_${k2}$, 3.17805383034794_${k2}$,    &
+                                         1.998721449566e1_${k2}$, 5.498100377941e4_${k2}$]
         #:elif k1 == "int64"
 
         ${t1}$, parameter :: x(n) = [0_${k1}$, 1_${k1}$, 2_${k1}$, 4_${k1}$,   &
                                      20_${k1}$, 90000_${k1}$]
-        real(dp), parameter :: ans(n) = [0.0_dp, 0.0_dp, 0.69314718055994_dp, 3.17805383034794_dp,    &
-                                         4.233561646075e1_dp, 9.366874681600e5_dp]
+        ${t2}$, parameter :: ans(n) = [0.0_${k2}$, 0.0_${k2}$, 0.69314718055994_${k2}$, 3.17805383034794_${k2}$,    &
+                                         4.233561646075e1_${k2}$, 9.366874681600e5_${k2}$]
         #:endif
 
         do i = 1, n
             call check(error, log_factorial(x(i)), ans(i), "Integer kind "     &
-                       //"${k1}$ failed", thr = tol_dp, rel = .true.)
+                       //"${k1}$ failed", thr = tol_${k2}$, rel = .true.)
 
         end do
     end subroutine test_logfact_${t1[0]}$${k1}$
@@ -199,6 +199,7 @@ contains
 
 
 
+    #:set k2, t2 = REAL_KINDS[-1], REAL_TYPES[-1]
     subroutine test_loggamma_${t1[0]}$${k1}$(error)
         type(error_type), allocatable, intent(out) :: error
         integer, parameter :: n = 4
@@ -208,24 +209,24 @@ contains
         #:if k1 == "int8"
 
         ${t1}$, parameter :: x(n) = [1_${k1}$, 2_${k1}$, 10_${k1}$, 47_${k1}$]
-        real(sp), parameter :: ans(n) = [0.0, 0.0, 1.28018274e1, 1.32952575e2]
+        ${t2}$, parameter :: ans(n) = [0.0_${k2}$, 0.0_${k2}$, 1.28018274e1_${k2}$, 1.32952575e2_${k2}$]
 
         #:elif k1 == "int16"
 
         ${t1}$, parameter :: x(n) = [1_${k1}$, 2_${k1}$, 111_${k1}$, 541_${k1}$]
-        real(sp), parameter :: ans(n) = [0.0, 0.0, 4.10322777e2, 2.86151221e3]
+        ${t2}$, parameter :: ans(n) = [0.0_${k2}$, 0.0_${k2}$, 4.10322777e2_${k2}$, 2.86151221e3_${k2}$]
 
         #:elif k1 == "int32"
 
         ${t1}$, parameter :: x(n) = [1_${k1}$, 2_${k1}$, 2021_${k1}$,          &
                                      42031_${k1}$]
-        real(sp), parameter :: ans(n) = [0.0, 0.0, 1.33586470e4, 4.05433461e5]
+        ${t2}$, parameter :: ans(n) = [0.0_${k2}$, 0.0_${k2}$, 1.33586470e4_${k2}$, 4.05433461e5_${k2}$]
 
         #:elif k1 == "int64"
 
         ${t1}$, parameter :: x(n) = [1_${k1}$, 2_${k1}$, 2021_${k1}$,          &
                                      42031_${k1}$]
-        real(sp), parameter :: ans(n) = [0.0, 0.0, 1.33586470e4, 4.05433461e5]
+        ${t2}$, parameter :: ans(n) = [0.0_${k2}$, 0.0_${k2}$, 1.33586470e4_${k2}$, 4.05433461e5_${k2}$]
 
         #:elif t1[0] == "c"
 
@@ -247,7 +248,7 @@ contains
         do i = 1, n
 
             call check(error, log_gamma(x(i)), ans(i), "Integer kind ${k1}$ "  &
-              //"failed", thr = tol_sp, rel = .true.)
+              //"failed", thr = tol_${k2}$, rel = .true.)
 
         end do