PUApproximant_fwd.hpp

Go to the documentation of this file.

#ifndef MEDUSA_BITS_INTERPOLANTS_PUAPPROXIMANT_FWD_HPP_
#define MEDUSA_BITS_INTERPOLANTS_PUAPPROXIMANT_FWD_HPP_
 
#include <Eigen/Core>
#include <medusa/bits/domains/DomainDiscretization_fwd.hpp>
#include <medusa/bits/types/Range_fwd.hpp>
 
namespace mm {
 
template <typename vec_t>
class PUApproximant {
    template <typename scalar_t>
    using eigen_vt = Eigen::Matrix<scalar_t, Eigen::Dynamic, 1>;  
    using d_scalar_t = typename vec_t::scalar_t;                  
 
  public:
    static typename vec_t::scalar_t weight(typename vec_t::scalar_t r) {
        assert_msg(r >= 0 && r <= 1,
                   "Effective radius must be larger than 0 and smaller than 1, got %g", r);
        return ipow<3>(1 - r) * (6 * r * r + 3 * r + 1);
    }
 
    template <typename scalar_t, typename engine_t>
    static eigen_vt<scalar_t> evaluate(const DomainDiscretization<vec_t>& domain,
                                       const eigen_vt<scalar_t>& values,
                                       const Range<vec_t>& query_points, d_scalar_t radius_factor,
                                       const engine_t& engine);
 
    template <typename scalar_t, typename engine_t>
    static eigen_vt<scalar_t> evaluate(const DomainDiscretization<vec_t>& domain,
                                       const eigen_vt<scalar_t>& values,
                                       const KDTree<vec_t>& query_points_tree,
                                       d_scalar_t radius_factor, const engine_t& engine);
};
 
}  // namespace mm
 
#endif  // MEDUSA_BITS_INTERPOLANTS_PUAPPROXIMANT_FWD_HPP_