Polyharmonic_fwd.hpp

Go to the documentation of this file.

#ifndef MEDUSA_BITS_APPROXIMATIONS_POLYHARMONIC_FWD_HPP_
#define MEDUSA_BITS_APPROXIMATIONS_POLYHARMONIC_FWD_HPP_
 
#include <medusa/Config.hpp>
#include <medusa/bits/utils/numutils.hpp>
#include <cmath>
#include "Operators_fwd.hpp"
 
namespace mm {
 
template <typename scal_t, int k = -1>
class Polyharmonic {
  public:
    static_assert((k > 0 && k % 2 == 1) || k == -1, "k must be odd or -1 (Dynamic)");
    typedef scal_t scalar_t;  
 
  private:
    const int order_;  
 
  public:
    Polyharmonic();
    Polyharmonic(int order);
 
    int order() const { return order_; }
 
    inline scalar_t operator()(scalar_t r2, int derivative) const;
 
    template <int dimension>
    inline scalar_t operator()(scalar_t r2, Lap<dimension> lap) const;
 
    inline scalar_t operator()(scalar_t r2) const { return ipow(std::sqrt(r2), order_); }
 
    template <typename S, int K>
    friend std::ostream& operator<<(std::ostream& os, const Polyharmonic<S, K>& m);
};
 
}  // namespace mm
 
#endif  // MEDUSA_BITS_APPROXIMATIONS_POLYHARMONIC_FWD_HPP_