Gaussian_fwd.hpp

Go to the documentation of this file.

#ifndef MEDUSA_BITS_APPROXIMATIONS_GAUSSIAN_FWD_HPP_
#define MEDUSA_BITS_APPROXIMATIONS_GAUSSIAN_FWD_HPP_
 
#include <medusa/Config.hpp>
#include <iosfwd>
#include <cmath>
#include "Operators_fwd.hpp"
namespace mm {
 
template <class scal_t>
class Gaussian {
  public:
    typedef scal_t scalar_t;  
 
  private:
    scalar_t shape_;  
 
  public:
    Gaussian(scalar_t shape = 1.0);
 
    scalar_t shape() const { return shape_; }
 
    void setShape(scalar_t shape) { shape_ = shape; }
 
    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 std::exp(-r2/shape_/shape_); }
 
    template <typename V>
    friend std::ostream& operator<<(std::ostream& os, const Gaussian<V>& m);
};
 
}  // namespace mm
 
#endif  // MEDUSA_BITS_APPROXIMATIONS_GAUSSIAN_FWD_HPP_