---

riesz.h

Go to the documentation of this file.

#ifndef __RIESZ__
#define __RIESZ__

#include "base_2D.h"

namespace Function_2D{

  class Riesz_1:public Filter_interface_2D{

  public:
    inline Riesz_1(const real_type center_x = 0,
                   const real_type center_y = 0,
                   const real_type sigma_x = 1,
                   const real_type sigma_y = 1,
                   const real_type orientation_angle = 0);

    virtual inline real_type base_value_proxy(const real_type x,
                                              const real_type y) const;

  };

  class Riesz_2:public Filter_interface_2D{

  public:
    inline Riesz_2(const real_type center_x = 0,
                   const real_type center_y = 0,
                   const real_type sigma_x = 1,
                   const real_type sigma_y = 1,
                   const real_type orientation_angle = 0);

    virtual inline real_type base_value_proxy(const real_type x,
                                              const real_type y) const;

  };


/*
  ###########################################################
  ###########################################################
  ###########################################################
  ##############                               ##############
  ##############  I M P L E M E N T A T I O N  ##############
  ##############                               ##############
  ###########################################################
  ###########################################################
  ###########################################################
*/

  Riesz_1::Riesz_1(const real_type center_x,
                   const real_type center_y,
                   const real_type sigma_x,
                   const real_type sigma_y,
                   const real_type orientation_angle):
    Filter_interface_2D(center_x,center_y,
                        sigma_x,sigma_y,
                        1,1,
                        orientation_angle,
                        1){}

  Riesz_1::real_type Riesz_1::base_value_proxy(const real_type x,
                                               const real_type y) const{

    const real_type r = x*x+y*y;
    
    return (r>0) ? x/pow(r,1.5) : 0;
  }

    
  Riesz_2::Riesz_2(const real_type center_x,
                   const real_type center_y,
                   const real_type sigma_x,
                   const real_type sigma_y,
                   const real_type orientation_angle):
    Filter_interface_2D(center_x,center_y,
                        sigma_x,sigma_y,
                        1,1,
                        orientation_angle,
                        1){}

  Riesz_2::real_type Riesz_2::base_value_proxy(const real_type x,
                                               const real_type y) const{

    const real_type r = x*x+y*y;

    return (r>0) ? y/pow(r,1.5) : 0;
  }
}

#endif

---