helmholtz_3d_hf_level_data.h

Go to the documentation of this file.

 
#ifndef FMM_HELMHOLTZ_3D_HF_LEVEL_DATA_H_INCLUDED
#define FMM_HELMHOLTZ_3D_HF_LEVEL_DATA_H_INCLUDED
 
#include "unit_sphere.h"
#include "unit_sphere_interpolator.h"
 
#include <Eigen/Dense>
 
#include <vector>
 
namespace NiHu
{
namespace fmm
{
 
class helmholtz_3d_hf_level_data
{
public:
    typedef Eigen::Matrix<std::complex<double>, Eigen::Dynamic, 1> cvector_t;
 
    helmholtz_3d_hf_level_data()
        : m_interp_ups(1)
        , m_interp_dns(1)
    {
    }
 
    void set_expansion_length(size_t L);
 
    size_t get_expansion_length() const;
 
    unit_sphere const &get_unit_sphere() const;
 
    void set_num_threads(size_t n);
 
    cvector_t const &interp_up(cvector_t const &x) const;
 
    cvector_t const &interp_down(cvector_t const &x) const;
 
    void set_interp_up(interpolator const &interp);
 
    void set_interp_dn(interpolator const &interp);
 
private:
    size_t m_expansion_length; // L
    unit_sphere m_unit_sphere;  // the unit phere on this level
    std::vector<interpolator> m_interp_ups; // upward interpolators fall all threads
    std::vector<interpolator> m_interp_dns; // downward interpolators for all threads
};
 
} // end of namespace fmm
} // end of namespace NiHu
 
#endif /* FMM_HELMHOLTZ_3D_HF_LEVEL_DATA_H_INCLUDED */