spherical_bessel_transformer.h

Go to the documentation of this file.

#ifndef SPHERICAL_BESSEL_TRANSFORMER_H_
#define SPHERICAL_BESSEL_TRANSFORMER_H_
 
#include <memory>
#include <fftw3.h>
 
namespace ModuleBase
{
 
class SphericalBesselTransformer
{
public:
    SphericalBesselTransformer(const bool cache_enabled = false);
    ~SphericalBesselTransformer() = default;
 
    SphericalBesselTransformer(SphericalBesselTransformer const&) = default;
    SphericalBesselTransformer(SphericalBesselTransformer &&) = default;
 
    SphericalBesselTransformer& operator=(const SphericalBesselTransformer&) = default;
    SphericalBesselTransformer& operator=(SphericalBesselTransformer&&) = default;
 
    void radrfft(
        const int l,
        const int ngrid,
        const double cutoff,
        const double* const in,
        double* const out,
        const int p = 0
    ) const;
 
 
    void direct(
        const int l,
        const int ngrid_in,
        const double* const grid_in,
        const double* const in,
        const int ngrid_out,
        const double* const grid_out,
        double* const out,
        const int p = 0
    ) const;
 
 
    size_t heap_usage() const;
 
    void clear();
 
    bool operator==(const SphericalBesselTransformer& rhs) const { return impl_ == rhs.impl_; }
 
 
private:
 
    class Impl; // forward declaration
    std::shared_ptr<Impl> impl_;
};
 
} // namespace ModuleBase
 
#endif