@@ -168,7 +168,7 @@ contains
168168 do k = 1 , sys_size
169169 do j = 0 , m
170170 do l = 1 , cmplx_size
171- data_fltr_cmplx_gpu(l + j* cmplx_size + (k - 1 )* cmplx_size* x_size) = (0_dp , 0_dp )
171+ data_fltr_cmplx_gpu(l + j* cmplx_size + (k - 1 )* cmplx_size* x_size) = (0._dp , 0._dp )
172172 end do
173173 end do
174174 end do
@@ -224,7 +224,7 @@ contains
224224 do k = 1 , sys_size
225225 do j = 0 , m
226226 do l = 1 , cmplx_size
227- data_fltr_cmplx_gpu(l + j* cmplx_size + (k - 1 )* cmplx_size* x_size) = (0_dp , 0_dp )
227+ data_fltr_cmplx_gpu(l + j* cmplx_size + (k - 1 )* cmplx_size* x_size) = (0._dp , 0._dp )
228228 end do
229229 end do
230230 end do
@@ -243,7 +243,7 @@ contains
243243 $:GPU_MKL_DISPATCH()
244244 ierr = DftiComputeForward(fwd_plan_mkl, data_real_gpu, data_cmplx_gpu)
245245
246- Nfq = min (floor(2_dp * real (i, dp)* pi), cmplx_size)
246+ Nfq = min (floor(2._dp * real (i, dp)* pi), cmplx_size)
247247 $:GPU_UPDATE(device= ' [Nfq]' )
248248
249249 $:GPU_PARALLEL_LOOP(collapse= 3 )
@@ -277,7 +277,7 @@ contains
277277 do k = 1 , sys_size
278278 do j = 0 , m
279279 do l = 1 , cmplx_size
280- data_fltr_cmplx_gpu(l + j* cmplx_size + (k - 1 )* cmplx_size* x_size) = (0_dp , 0_dp )
280+ data_fltr_cmplx_gpu(l + j* cmplx_size + (k - 1 )* cmplx_size* x_size) = (0._dp , 0._dp )
281281 end do
282282 end do
283283 end do
@@ -344,7 +344,7 @@ contains
344344 do k = 1 , sys_size
345345 do j = 0 , m
346346 do l = 1 , cmplx_size
347- data_fltr_cmplx_gpu(l + j* cmplx_size + (k - 1 )* cmplx_size* x_size) = (0_dp , 0_dp )
347+ data_fltr_cmplx_gpu(l + j* cmplx_size + (k - 1 )* cmplx_size* x_size) = (0._dp , 0._dp )
348348 end do
349349 end do
350350 end do
@@ -369,7 +369,7 @@ contains
369369#endif
370370 #:endcall GPU_HOST_DATA
371371
372- Nfq = min (floor(2_dp * real (i, dp)* pi), cmplx_size)
372+ Nfq = min (floor(2._dp * real (i, dp)* pi), cmplx_size)
373373 $:GPU_UPDATE(device= ' [Nfq]' )
374374
375375 $:GPU_PARALLEL_LOOP(collapse= 3 )
@@ -408,7 +408,7 @@ contains
408408 Nfq = 3
409409 do j = 0 , m
410410 do k = 1 , sys_size
411- data_fltr_cmplx(:) = (0_dp , 0_dp )
411+ data_fltr_cmplx(:) = (0._dp , 0._dp )
412412 data_real(1 :p + 1 ) = q_cons_vf(k)%sf(j, 0 ,0 :p)
413413 call fftw_execute_dft_r2c(fwd_plan, data_real, data_cmplx)
414414 data_fltr_cmplx(1 :Nfq) = data_cmplx(1 :Nfq)
@@ -420,10 +420,10 @@ contains
420420
421421 ! Apply Fourier filter to additional rings
422422 do i = 1 , fourier_rings
423- Nfq = min (floor(2_dp * real (i, dp)* pi), cmplx_size)
423+ Nfq = min (floor(2._dp * real (i, dp)* pi), cmplx_size)
424424 do j = 0 , m
425425 do k = 1 , sys_size
426- data_fltr_cmplx(:) = (0_dp , 0_dp )
426+ data_fltr_cmplx(:) = (0._dp , 0._dp )
427427 data_real(1 :p + 1 ) = q_cons_vf(k)%sf(j, i,0 :p)
428428 call fftw_execute_dft_r2c(fwd_plan, data_real, data_cmplx)
429429 data_fltr_cmplx(1 :Nfq) = data_cmplx(1 :Nfq)
0 commit comments