From 21a2abb375380d745d648ba868ef0a0de4d3cbbe Mon Sep 17 00:00:00 2001 From: Wojtek Kosior Date: Sat, 1 Jun 2019 12:03:01 +0200 Subject: compute sinuses signal --- test.f90 | 28 +++++++++++++++++++++++++--- 1 file 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 -- cgit v1.2.3