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 ++++++++++++++++++++++++++++++++ test.f90 | 32 -------------------------------- 2 files changed, 32 insertions(+), 32 deletions(-) create mode 100644 src/f1.f90 delete mode 100644 test.f90 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 diff --git a/test.f90 b/test.f90 deleted file mode 100644 index 38f360f..0000000 --- a/test.f90 +++ /dev/null @@ -1,32 +0,0 @@ -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