@@ -89,36 +89,36 @@ contains
89
89
90
90
${t1}$, parameter :: x(n) = [0_${k1}$, 1_${k1}$, 2_${k1}$, 4_${k1}$, &
91
91
5_${k1}$, 100_${k1}$]
92
- real(sp ), parameter :: ans(n) = [0.0, 0.0, 0.693147180, 3.17805383, &
92
+ real(dp ), parameter :: ans(n) = [0.0, 0.0, 0.693147180, 3.17805383, &
93
93
4.78749174, 3.63739376e2]
94
94
95
95
#:elif k1 == "int16"
96
96
97
97
${t1}$, parameter :: x(n) = [0_${k1}$, 1_${k1}$, 2_${k1}$, 4_${k1}$, &
98
98
7_${k1}$, 500_${k1}$]
99
- real(sp ), parameter :: ans(n) = [0.0, 0.0, 0.693147180, 3.17805383, &
99
+ real(dp ), parameter :: ans(n) = [0.0, 0.0, 0.693147180, 3.17805383, &
100
100
8.52516136, 2.61133046e3]
101
101
102
102
#:elif k1 == "int32"
103
103
104
104
${t1}$, parameter :: x(n) = [0_${k1}$, 1_${k1}$, 2_${k1}$, 4_${k1}$, &
105
105
12_${k1}$, 7000_${k1}$]
106
- real(sp ), parameter :: ans(n) = [0.0, 0.0, 0.693147180, 3.17805383, &
106
+ real(dp ), parameter :: ans(n) = [0.0, 0.0, 0.693147180, 3.17805383, &
107
107
1.99872145e1, 5.49810038e4]
108
108
109
109
#:elif k1 == "int64"
110
110
111
111
${t1}$, parameter :: x(n) = [0_${k1}$, 1_${k1}$, 2_${k1}$, 4_${k1}$, &
112
112
20_${k1}$, 90000_${k1}$]
113
- real(sp ), parameter :: ans(n) = [0.0, 0.0, 0.693147180, 3.17805383, &
113
+ real(dp ), parameter :: ans(n) = [0.0, 0.0, 0.693147180, 3.17805383, &
114
114
4.23356165e1, 9.36687468e5]
115
115
116
116
#:endif
117
117
118
118
do i = 1, n
119
119
120
120
call check(error, log_factorial(x(i)), ans(i), "Integer kind " &
121
- //"${k1}$ failed", thr = tol_sp , rel = .true.)
121
+ //"${k1}$ failed", thr = tol_dp , rel = .true.)
122
122
123
123
end do
124
124
end subroutine test_logfact_${t1[0]}$${k1}$
0 commit comments