@@ -100,8 +100,8 @@ subroutine test_fftfreq_2(error)
100
100
type (error_type), allocatable , intent (out ) :: error
101
101
102
102
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
105
105
106
106
integer , parameter :: n = 128
107
107
integer :: i
@@ -110,11 +110,11 @@ subroutine test_fftfreq_2(error)
110
110
111
111
do i = 1 , n
112
112
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
114
114
end do
115
115
116
116
xfou = fft(xvec)/ n
117
- xfou = imu* fftfreq(n)* xfou ! > derivative in Fourier space
117
+ xfou = imu* fftfreq(n)* xfou ! derivative in Fourier space
118
118
xvec = ifft(xfou)
119
119
call check(error, maxval (abs (xvec- xtrue)) < tol, &
120
120
" maxval(abs(xvec-xtrue)) < tol failed." )
@@ -125,8 +125,8 @@ subroutine test_fftfreq_3(error)
125
125
type (error_type), allocatable , intent (out ) :: error
126
126
127
127
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
130
130
131
131
integer , parameter :: n = 135
132
132
integer :: i
@@ -135,11 +135,11 @@ subroutine test_fftfreq_3(error)
135
135
136
136
do i = 1 , n
137
137
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
139
139
end do
140
140
141
141
xfou = fft(xvec)/ n
142
- xfou = imu* fftfreq(n)* xfou ! > derivative in Fourier space
142
+ xfou = imu* fftfreq(n)* xfou ! derivative in Fourier space
143
143
xvec = ifft(xfou)
144
144
call check(error, maxval (abs (xvec- xtrue)) < tol, &
145
145
" maxval(abs(xvec-xtrue)) < tol failed." )
0 commit comments