@@ -100,8 +100,8 @@ subroutine test_fftfreq_2(error)
100100 type (error_type), allocatable , intent (out ) :: error
101101
102102 real (rk), parameter :: tol = 1.0e-12_rk
103- real (rk), parameter :: twopi = 8 * atan (1.0_rk ) ! > 2*pi
104- complex (rk), parameter :: imu = (0 ,1 ) ! > imaginary unit
103+ real (rk), parameter :: twopi = 8 * atan (1.0_rk ) ! 2*pi
104+ complex (rk), parameter :: imu = (0 ,1 ) ! imaginary unit
105105
106106 integer , parameter :: n = 128
107107 integer :: i
@@ -110,11 +110,11 @@ subroutine test_fftfreq_2(error)
110110
111111 do i = 1 , n
112112 xvec(i) = cos (twopi* (i-1 )/ n)
113- xtrue(i) = - sin (twopi* (i-1 )/ n) ! > derivative in physical space
113+ xtrue(i) = - sin (twopi* (i-1 )/ n) ! derivative in physical space
114114 end do
115115
116116 xfou = fft(xvec)/ n
117- xfou = imu* fftfreq(n)* xfou ! > derivative in Fourier space
117+ xfou = imu* fftfreq(n)* xfou ! derivative in Fourier space
118118 xvec = ifft(xfou)
119119 call check(error, maxval (abs (xvec- xtrue)) < tol, &
120120 " maxval(abs(xvec-xtrue)) < tol failed." )
@@ -125,8 +125,8 @@ subroutine test_fftfreq_3(error)
125125 type (error_type), allocatable , intent (out ) :: error
126126
127127 real (rk), parameter :: tol = 1.0e-12_rk
128- real (rk), parameter :: twopi = 8 * atan (1.0_rk ) ! > 2*pi
129- complex (rk), parameter :: imu = (0 ,1 ) ! > imaginary unit
128+ real (rk), parameter :: twopi = 8 * atan (1.0_rk ) ! 2*pi
129+ complex (rk), parameter :: imu = (0 ,1 ) ! imaginary unit
130130
131131 integer , parameter :: n = 135
132132 integer :: i
@@ -135,11 +135,11 @@ subroutine test_fftfreq_3(error)
135135
136136 do i = 1 , n
137137 xvec(i) = cos (twopi* (i-1 )/ n)
138- xtrue(i) = - sin (twopi* (i-1 )/ n) ! > derivative in physical space
138+ xtrue(i) = - sin (twopi* (i-1 )/ n) ! derivative in physical space
139139 end do
140140
141141 xfou = fft(xvec)/ n
142- xfou = imu* fftfreq(n)* xfou ! > derivative in Fourier space
142+ xfou = imu* fftfreq(n)* xfou ! derivative in Fourier space
143143 xvec = ifft(xfou)
144144 call check(error, maxval (abs (xvec- xtrue)) < tol, &
145145 " maxval(abs(xvec-xtrue)) < tol failed." )
0 commit comments