Homec4science
Diffusion tamaas (petsc) Browse

tamaas/src/mpi/fft_mpi_engine.hh
8facf8e0cda6petsc

Last Change
Raw File

fft_mpi_engine.hh
View Options

/**
* @file
* @section LICENSE
*
* Copyright (©) 2016-2020 EPFL (École Polytechnique Fédérale de Lausanne),
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* This program is free software: you can redistribute it and/or modify
* it under the terms of the GNU Affero General Public License as published
* by the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* This program 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 Affero General Public License for more details.
*
* You should have received a copy of the GNU Affero General Public License
* along with this program. If not, see <https://www.gnu.org/licenses/>.
*
*/
/* -------------------------------------------------------------------------- */
#ifndef FFT_MPI_ENGINE_HH
#define FFT_MPI_ENGINE_HH
/* -------------------------------------------------------------------------- */
#include "fft_engine.hh"
#include "fftw_interface.hh"
#include "grid.hh"
#include "grid_hermitian.hh"
#include <map>
/* -------------------------------------------------------------------------- */
namespace tamaas {
class FFTMPIEngine : public FFTEngine {
public:
void forward(Grid<Real, 1>& /*real*/, GridHermitian<Real, 1>& /*spectral*/) {
TAMAAS_EXCEPTION("FFTW/MPI does not support 1D transforms");
}
void backward(Grid<Real, 1>& /*real*/, GridHermitian<Real, 1>& /*spectral*/) {
TAMAAS_EXCEPTION("FFTW/MPI does not support 1D transforms");
}
/// FFTW/MPI forward (r2c) transform
void forward(Grid<Real, 2>& real, GridHermitian<Real, 2>& spectral);
/// FFTW/MPI backward (c2r) transform
void backward(Grid<Real, 2>& real, GridHermitian<Real, 2>& spectral);
protected:
/// Make a transform signature from a pair of grids
static key_t make_key(Grid<Real, 2>& real, GridHermitian<Real, 2>& spectral);
/// Return the plans pair for a given transform signature
plan_t& getPlans(key_t key);
/// Get FFTW local sizes from an hermitian grid
static auto local_size(const key_t& key);
protected:
/// Buffer for real data because of FFTW/MPI layout
std::map<key_t, Grid<Real, 2>> workspaces;
};
} // namespace tamaas
#endif // FFT_MPI_ENGINE_HH
c4science · Help