1
1
#:include "common.fypp"
2
- #:set R_KINDS_TYPES = [KT for KT in REAL_KINDS_TYPES if KT[0] in ["sp","dp"]]
3
- #:set C_KINDS_TYPES = [KT for KT in CMPLX_KINDS_TYPES if KT[0] in ["sp","dp"] ]
4
- #:set CI_KINDS_TYPES = INT_KINDS_TYPES + C_KINDS_TYPES
2
+ #:set CI_KINDS_TYPES = INT_KINDS_TYPES + CMPLX_KINDS_TYPES
3
+ #:set IDX_CMPLX_KINDS_TYPES = [(i, CMPLX_KINDS[i], CMPLX_TYPES[i], CMPLX_INIT[i]) for i in range(len(CMPLX_KINDS)) ]
4
+ #:set IDX_REAL_KINDS_TYPES = [(i, REAL_KINDS[i], REAL_TYPES[i], REAL_INIT[i]) for i in range(len(REAL_KINDS))]
5
5
module test_specialfunctions_gamma
6
6
use testdrive, only : new_unittest, unittest_type, error_type, check
7
- use stdlib_kinds, only: sp, dp, int8, int16, int32, int64
7
+ use stdlib_kinds, only: sp, dp, xdp, qp, int8, int16, int32, int64
8
8
use stdlib_specialfunctions_gamma, only: gamma, log_gamma, log_factorial, &
9
9
lower_incomplete_gamma, &
10
10
upper_incomplete_gamma, &
@@ -18,7 +18,7 @@ module test_specialfunctions_gamma
18
18
19
19
public :: collect_specialfunctions_gamma
20
20
21
- #:for k1, t1 in R_KINDS_TYPES
21
+ #:for k1, t1 in REAL_KINDS_TYPES
22
22
${t1}$, parameter :: tol_${k1}$ = 1000 * epsilon(1.0_${k1}$)
23
23
#:endfor
24
24
@@ -35,15 +35,15 @@ contains
35
35
test_logfact_${t1[0]}$${k1}$) &
36
36
#:endfor
37
37
38
- #:for k1, t1 in CI_KINDS_TYPES
38
+ #:for k1, t1 in CI_KINDS_TYPES[:-1]
39
39
, new_unittest("gamma_${t1[0]}$${k1}$", &
40
40
test_gamma_${t1[0]}$${k1}$) &
41
41
, new_unittest("log_gamma_${t1[0]}$${k1}$", &
42
42
test_loggamma_${t1[0]}$${k1}$) &
43
43
#:endfor
44
44
45
45
#:for k1, t1 in INT_KINDS_TYPES
46
- #:for k2, t2 in R_KINDS_TYPES
46
+ #:for k2, t2 in REAL_KINDS_TYPES[:-1]
47
47
, new_unittest("lower_incomplete_gamma_${t1[0]}$${k1}$${k2}$", &
48
48
test_lincgamma_${t1[0]}$${k1}$${k2}$) &
49
49
, new_unittest("log_lower_incomplete_gamma_${t1[0]}$${k1}$${k2}$", &
@@ -59,7 +59,7 @@ contains
59
59
#:endfor
60
60
#:endfor
61
61
62
- #:for k1, t1 in R_KINDS_TYPES
62
+ #:for k1, t1 in REAL_KINDS_TYPES[:-1]
63
63
, new_unittest("lower_incomplete_gamma_${t1[0]}$${k1}$", &
64
64
test_lincgamma_${t1[0]}$${k1}$) &
65
65
, new_unittest("log_lower_incomplete_gamma_${t1[0]}$${k1}$", &
@@ -124,7 +124,7 @@ contains
124
124
125
125
126
126
127
- #:for k1, t1 in CI_KINDS_TYPES
127
+ #:for k1, t1 in CI_KINDS_TYPES[:-1]
128
128
129
129
subroutine test_gamma_${t1[0]}$${k1}$(error)
130
130
type(error_type), allocatable, intent(out) :: error
@@ -262,7 +262,7 @@ contains
262
262
263
263
264
264
#:for k1, t1 in INT_KINDS_TYPES
265
- #:for k2, t2 in R_KINDS_TYPES
265
+ #:for k2, t2 in REAL_KINDS_TYPES[:-1]
266
266
267
267
subroutine test_lincgamma_${t1[0]}$${k1}$${k2}$(error)
268
268
type(error_type), allocatable, intent(out) :: error
@@ -411,7 +411,7 @@ contains
411
411
412
412
413
413
414
- #:for k1, t1 in R_KINDS_TYPES
414
+ #:for k1, t1 in REAL_KINDS_TYPES[:-1]
415
415
416
416
subroutine test_lincgamma_${t1[0]}$${k1}$(error)
417
417
type(error_type), allocatable, intent(out) :: error
0 commit comments