From 94e4f02728e102d69a1a2378869c22519751fe47 Mon Sep 17 00:00:00 2001 From: Wojtek Kosior Date: Sat, 1 Jun 2019 12:08:17 +0200 Subject: move: text.f90 → src/f1.f90 MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- src/f1.f90 | 32 ++++++++++++++++++++++++++++++++ 1 file changed, 32 insertions(+) create mode 100644 src/f1.f90 (limited to 'src') diff --git a/src/f1.f90 b/src/f1.f90 new file mode 100644 index 0000000..38f360f --- /dev/null +++ b/src/f1.f90 @@ -0,0 +1,32 @@ +PROGRAM test + + use, intrinsic :: iso_c_binding + IMPLICIT none + + INCLUDE 'fftw3.f03' + + 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 + + PURE 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