aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--test.f9028
1 files changed, 25 insertions, 3 deletions
diff --git a/test.f90 b/test.f90
index abc2e58..b5beb0b 100644
--- a/test.f90
+++ b/test.f90
@@ -1,10 +1,32 @@
PROGRAM test
use, intrinsic :: iso_c_binding
- implicit none
+ IMPLICIT none
- include 'fftw3.f03'
+ INCLUDE 'fftw3.f03'
- write(*,*) "ok"
+ integer, parameter :: npoints = 1024
+ real(C_DOUBLE), parameter :: start = 0, end = 1
+ integer :: i
+ real(C_DOUBLE) :: t, &
+ delta = (end - start) / real(npoints - 1, C_DOUBLE)
+ DO i = 0, npoints - 1
+
+ t = i * delta + start
+
+ write(*,*) t, " ", f1(t)
+
+ END DO
+
+CONTAINS
+
+ real(C_DOUBLE) FUNCTION f1(t)
+ real(C_DOUBLE), intent(in) :: t
+ real(C_DOUBLE), parameter :: pi = acos(-1.0)
+
+ f1 = sin(2 * pi * t * 200) + 2 * sin(2 * pi * t * 400)
+
+ END FUNCTION f1
+
END PROGRAM test