unit_sphere.h

Go to the documentation of this file.

 
#ifndef NIHU_UNIT_SPHERE_H_INCLUDED
#define NIHU_UNIT_SPHERE_H_INCLUDED
 
#include "Eigen/Dense"
 
namespace NiHu
{
namespace fmm
{
 
class unit_sphere
{
    typedef Eigen::Matrix<std::complex<double>, Eigen::Dynamic, 1> cvector_t;
    typedef Eigen::Matrix<double, Eigen::Dynamic, 1> dvector_t;
    typedef Eigen::Matrix<std::complex<double>, Eigen::Dynamic, Eigen::Dynamic> cmatrix_t;
    typedef Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic> dmatrix_t;
 
    size_t m_P; // quadrature density
    dvector_t theta, phi, w, wtheta;
    Eigen::Matrix<double, 3, Eigen::Dynamic> s;
 
public:
    unit_sphere()
    {
    }
 
    unit_sphere(size_t P);
 
    size_t get_P(void) const
    {
        return m_P;
    }
 
    dvector_t const &get_w(void) const
    {
        return w;
    }
 
    dvector_t const &get_theta(void) const
    {
        return this->theta;
    }
 
    dvector_t const &get_wtheta(void) const
    {
        return this->wtheta;
    }
 
    dvector_t const &get_phi(void) const
    {
        return phi;
    }
 
    Eigen::Matrix<double, 3, Eigen::Dynamic> const &get_s(void) const
    {
        return s;
    }
 
    cmatrix_t spht(cvector_t const &f, int L) const;
 
    cvector_t ispht(cmatrix_t const &Flm) const;
};
 
} // namespace fmm
} // namespace NiHu
 
#endif /* NIHU_UNIT_SPHERE_H_INCLUDED */