types.h

Go to the documentation of this file.

#ifndef SIRIUS_FFTW_TYPES_H_
#define SIRIUS_FFTW_TYPES_H_

#include <memory>

#include <fftw3.h>

#include "sirius/utils/log.h"

namespace sirius {
namespace fftw {

namespace detail {

struct ComplexDeleter {
    void operator()(::fftw_complex* complex) {
        if (complex == nullptr) {
            return;
        }
        ::fftw_free(complex);
    }
};

struct RealDeleter {
    void operator()(double* real) { ::fftw_free(real); }
};

}  // namespace detail

using ComplexUPtr = std::unique_ptr<::fftw_complex[], detail::ComplexDeleter>;

#if (!defined(__GNUC__) && __cplusplus <= 201402L) || \
      (defined(__GNUC__) && __GNUC__ < 7 && __cplusplus <= 201402L)

// C++14: no shared_ptr array syntax, classic definition
using ComplexSPtr = std::shared_ptr<::fftw_complex>;

#else

// C++17: std::shared_ptr array syntax

// compiling with GCC7 breaks the build if C++14 syntax is used
// definition of std::shared_ptr<T>::element_type has changed:
// <7: typedef _Tp element_type;
// >=7: using element_type = typename remove_extent<_Tp>::type;
// With GCC <7:
//   std::shared_ptr<double[2]>::element_type <=> double[2]
// With GCC >=7:
//   std::shared_ptr<double[2]>::element_type <=> double
//   std::shared_ptr<double[][2]>::element_type <=> double[2]

using ComplexSPtr = std::shared_ptr<::fftw_complex[]>;

#endif  // (!defined(__GNUC__) && __cplusplus <= 201402L) ||
        //   (defined(__GNUC__) && __GNUC__ < 7 && __cplusplus <= 201402L)

using RealUPtr = std::unique_ptr<double[], detail::RealDeleter>;

}  // namespace fftw
}  // namespace sirius

#endif  // SIRIUS_FFTW_TYPES_H_