fft_plan.F

Go to the documentation of this file.

!--------------------------------------------------------------------------------------------------!
!   CP2K: A general program to perform molecular dynamics simulations                              !
!   Copyright 2000-2024 CP2K developers group <https://cp2k.org>                                   !
!                                                                                                  !
!   SPDX-License-Identifier: GPL-2.0-or-later                                                      !
!--------------------------------------------------------------------------------------------------!
 
! **************************************************************************************************
!> \brief Type to store data about a (1D or 3D) FFT, including FFTW plan
!> \par History
!>      IAB 09-Jan-2009 : initial version
!>                        (c) The Numerical Algorithms Group (NAG) Ltd, 2009 on behalf of the HECToR project
!>      IAB 09-Oct-2009 : Added additional fields needed when using OpenMP
!>                        (c) The Numerical Algorithms Group (NAG) Ltd, 2009 on behalf of the HECToR project
!> \author JGH
! **************************************************************************************************
 
MODULE fft_plan
   USE iso_c_binding,                   ONLY: c_null_ptr,&
                                              c_ptr
 
   IMPLICIT NONE
   PRIVATE
 
   PUBLIC :: fft_plan_type
 
   TYPE fft_plan_type
 
      INTEGER                             :: fft_type = -1
      INTEGER                             :: fsign = 0
      LOGICAL                             :: trans = .false., fft_in_place = .false., valid = .false., separated_plans = .false.
      INTEGER                             :: n = -1, m = -1
      INTEGER, DIMENSION(3)               :: n_3d = -1
 
!   Handle for the FFTW plan
      TYPE(c_ptr)                         :: fftw_plan = c_null_ptr
 
!   Plan for the remaining rows for 1D FFT when number of threads does not divide the number of rows exactly
!$    TYPE(C_PTR)                         :: alt_fftw_plan = C_NULL_PTR
!$    LOGICAL                             :: need_alt_plan = .FALSE.
!$    INTEGER                             :: num_threads_needed = -1, num_rows = -1, alt_num_rows = -1
 
!   Individual plans (used by hand-optimised 3D FFT)
      TYPE(c_ptr)                     :: fftw_plan_nx = c_null_ptr, fftw_plan_ny = c_null_ptr, fftw_plan_nz = c_null_ptr
!   Plans for the remaining rows (when the number of threads does not divide the number of rows exactly)
      TYPE(c_ptr)                     :: fftw_plan_nx_r = c_null_ptr, fftw_plan_ny_r = c_null_ptr, fftw_plan_nz_r = c_null_ptr
 
   END TYPE fft_plan_type
 
END MODULE fft_plan