@@ -73,8 +73,8 @@ subroutine test_eye
73
73
msg= " sum(rye - diag([(1.0_sp,i=1,6)])) < sptol failed." ,warn= warn)
74
74
75
75
cye = eye(7 )
76
- call check(abs (trace(cye) - cmplx (7.0_sp ,0.0_sp )) < sptol, &
77
- msg= " abs(trace(cye) - cmplx(7.0_sp,0.0_sp)) < sptol failed." ,warn= warn)
76
+ call check(abs (trace(cye) - cmplx (7.0_sp ,0.0_sp ,kind = sp )) < sptol, &
77
+ msg= " abs(trace(cye) - cmplx(7.0_sp,0.0_sp,kind=sp )) < sptol failed." ,warn= warn)
78
78
end subroutine
79
79
80
80
subroutine test_diag_rsp
@@ -153,7 +153,7 @@ subroutine test_diag_rqp
153
153
subroutine test_diag_csp
154
154
integer , parameter :: n = 3
155
155
complex (sp) :: v(n), a(n,n), b(n,n)
156
- complex (sp), parameter :: i_ = cmplx (0 ,1 )
156
+ complex (sp), parameter :: i_ = cmplx (0 ,1 ,kind = sp )
157
157
integer :: i,j
158
158
write (* ,* ) " test_diag_csp"
159
159
a = diag([(i,i= 1 ,n)]) + diag([(i_,i= 1 ,n)])
@@ -170,7 +170,7 @@ subroutine test_diag_csp
170
170
subroutine test_diag_cdp
171
171
integer , parameter :: n = 3
172
172
complex (dp) :: v(n), a(n,n), b(n,n)
173
- complex (dp), parameter :: i_ = cmplx (0 ,1 )
173
+ complex (dp), parameter :: i_ = cmplx (0 ,1 ,kind = dp )
174
174
integer :: i,j
175
175
write (* ,* ) " test_diag_cdp"
176
176
a = diag([i_],- 2 ) + diag([i_],2 )
@@ -181,7 +181,7 @@ subroutine test_diag_cdp
181
181
subroutine test_diag_cqp
182
182
integer , parameter :: n = 3
183
183
complex (qp) :: v(n), a(n,n), b(n,n)
184
- complex (qp), parameter :: i_ = cmplx (0 ,1 )
184
+ complex (qp), parameter :: i_ = cmplx (0 ,1 ,kind = qp )
185
185
integer :: i,j
186
186
write (* ,* ) " test_diag_cqp"
187
187
a = diag([i_,i_],- 1 ) + diag([i_,i_],1 )
@@ -333,7 +333,7 @@ subroutine test_trace_csp
333
333
integer , parameter :: n = 5
334
334
real (sp) :: re(n,n), im(n,n)
335
335
complex (sp) :: a(n,n), b(n,n)
336
- complex (sp), parameter :: i_ = cmplx (0 ,1 )
336
+ complex (sp), parameter :: i_ = cmplx (0 ,1 ,kind = sp )
337
337
write (* ,* ) " test_trace_csp"
338
338
339
339
call random_number (re)
@@ -352,12 +352,12 @@ subroutine test_trace_csp
352
352
subroutine test_trace_cdp
353
353
integer , parameter :: n = 3
354
354
complex (dp) :: a(n,n), ans
355
- complex (dp), parameter :: i_ = cmplx (0 ,1 )
355
+ complex (dp), parameter :: i_ = cmplx (0 ,1 ,kind = dp )
356
356
integer :: j
357
357
write (* ,* ) " test_trace_cdp"
358
358
359
359
a = reshape ([(j + (n** 2 - (j-1 ))* i_,j= 1 ,n** 2 )],[n,n])
360
- ans = cmplx (15 ,15 ) ! (1 + 5 + 9) + (9 + 5 + 1)i
360
+ ans = cmplx (15 ,15 ,kind = dp ) ! (1 + 5 + 9) + (9 + 5 + 1)i
361
361
362
362
call check(abs (trace(a) - ans) < dptol, &
363
363
msg= " abs(trace(a) - ans) < dptol failed." ,warn= warn)
@@ -366,7 +366,7 @@ subroutine test_trace_cdp
366
366
subroutine test_trace_cqp
367
367
integer , parameter :: n = 3
368
368
complex (qp) :: a(n,n)
369
- complex (qp), parameter :: i_ = cmplx (0 ,1 )
369
+ complex (qp), parameter :: i_ = cmplx (0 ,1 ,kind = qp )
370
370
write (* ,* ) " test_trace_cqp"
371
371
a = 3 * eye(n) + 4 * eye(n)* i_ ! pythagorean triple
372
372
call check(abs (trace(a)) - 3 * 5.0_qp < qptol, &
0 commit comments