aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--src/fourier.f9014
1 files changed, 12 insertions, 2 deletions
diff --git a/src/fourier.f90 b/src/fourier.f90
index be0a16f..b1d2b84 100644
--- a/src/fourier.f90
+++ b/src/fourier.f90
@@ -87,10 +87,20 @@ PROGRAM test
DO i = 1, ncomplex
write(*,*) i - 1, " ", abs(arr_complex(i))
-
END DO
ELSE
- STOP "not implemented"
+ DO i = 1, ncomplex
+
+ IF (abs(arr_complex(i)) < 50) arr_complex(i) = 0
+ END DO
+
+ call fftw_execute_dft_c2r(plan_b, arr_complex, arr_real)
+
+ DO i = 1, npoints
+
+ write(*,*) (i - 1) * delta + start, " ", arr_real(i)
+ END DO
+
END IF
call fftw_free(p_real)