Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F80910396
fftwmpi_engine.hh
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Tue, Sep 3, 05:18
Size
2 KB
Mime Type
text/x-c
Expires
Thu, Sep 5, 05:18 (2 d)
Engine
blob
Format
Raw Data
Handle
20470288
Attached To
rMUSPECTRE µSpectre
fftwmpi_engine.hh
View Options
/**
* @file fftwmpi_engine.hh
*
* @author Lars Pastewka <lars.pastewka@imtek.uni-freiburg.de>
*
* @date 06 Mar 2017
*
* @brief FFT engine using MPI-parallel FFTW
*
* Copyright © 2017 Till Junge
*
* µSpectre is free software; you can redistribute it and/or
* modify it under the terms of the GNU General Public License as
* published by the Free Software Foundation, either version 3, or (at
* your option) any later version.
*
* µSpectre is distributed in the hope that it will be useful, but
* WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with GNU Emacs; see the file COPYING. If not, write to the
* Free Software Foundation, Inc., 59 Temple Place - Suite 330,
* Boston, MA 02111-1307, USA.
*/
#ifndef FFTWMPI_ENGINE_H
#define FFTWMPI_ENGINE_H
#include "common/communicator.hh"
#include "fft/fft_engine_base.hh"
#include <fftw3-mpi.h>
namespace
muSpectre
{
/**
* implements the `muSpectre::Fft_Engine_Base` interface using the
* FFTW library
*/
template
<
Dim_t
DimS
,
Dim_t
DimM
>
class
FFTWMPIEngine
:
public
FFT_Engine_base
<
DimS
,
DimM
>
{
public
:
using
Parent
=
FFT_Engine_base
<
DimS
,
DimM
>
;
//!< base class
using
Ccoord
=
typename
Parent
::
Ccoord
;
//!< cell coordinates type
using
Rcoord
=
typename
Parent
::
Rcoord
;
//!< spatial coordinates type
//! field for Fourier transform of second-order tensor
using
Workspace_t
=
typename
Parent
::
Workspace_t
;
//! real-valued second-order tensor
using
Field_t
=
typename
Parent
::
Field_t
;
//! Default constructor
FFTWMPIEngine
()
=
delete
;
//! Constructor with system resolutions
FFTWMPIEngine
(
Ccoord
resolutions
,
Rcoord
lengths
,
Communicator
comm
=
Communicator
());
//! Copy constructor
FFTWMPIEngine
(
const
FFTWMPIEngine
&
other
)
=
delete
;
//! Move constructor
FFTWMPIEngine
(
FFTWMPIEngine
&&
other
)
=
default
;
//! Destructor
virtual
~
FFTWMPIEngine
()
noexcept
;
//! Copy assignment operator
FFTWMPIEngine
&
operator
=
(
const
FFTWMPIEngine
&
other
)
=
delete
;
//! Move assignment operator
FFTWMPIEngine
&
operator
=
(
FFTWMPIEngine
&&
other
)
=
default
;
// compute the plan, etc
virtual
void
initialise
(
FFT_PlanFlags
plan_flags
)
override
;
//! forward transform
virtual
Workspace_t
&
fft
(
Field_t
&
field
)
override
;
//! inverse transform
virtual
void
ifft
(
Field_t
&
field
)
const
override
;
protected
:
static
int
nb_engines
;
Communicator
comm
;
//!< communicator
fftw_plan
plan_fft
{};
//!< holds the plan for forward fourier transform
fftw_plan
plan_ifft
{};
//!< holds the plan for inverse fourier transform
bool
initialised
{
false
};
//!< to prevent double initialisation
private
:
};
}
// muSpectre
#endif
/* FFTWMPI_ENGINE_H */
Event Timeline
Log In to Comment