s2_poly_wave_approx.hpp

/* s2_poly_wave_approx.hpp */
/* $Id: s2_poly_wave_approx.hpp 21022 2011-07-11 11:45:03Z urij $ */
#ifndef __s2_poly_wave_approx_hpp
#define __s2_poly_wave_approx_hpp

#include <s2proc/wave_field.hpp>
#include <s2proc/s2_corr.hpp>
#include <mth/least_squares_2d.hpp>

#include <utility>

class S2PolyWaveApprox
{
public:
  
  struct Params
  {
    double twin_start; 
    double twin_len;   
    double dt;         
    size_t pdeg;       
    size_t tbnum;      
    void info() const
    {
      printf("S2PolyWaveApprox::Params: win(%g..%g), dt=%g, pdeg=%zu, tbnum=%zu\n",
         twin_start, twin_start+twin_len, dt, pdeg, tbnum);
    }
  };
  
  typedef std::pair<double, double> Node; 
  
  S2PolyWaveApprox(const S2PolyWaveApprox::Params &params) : m_params(params) {};
  
  virtual ~S2PolyWaveApprox() {};
  
  LeastSquares2D::PolyCoeffs operator() (const WaveField &field) const {return calc_poly(field);}
  
  LeastSquares2D::PolyCoeffs calc_poly(const WaveField &field) const;
  
  LeastSquares2D::PolyCoeffs calc_poly(const std::vector<S2PolyWaveApprox::Node>& nodes) const;
  
  virtual std::vector<S2PolyWaveApprox::Node> calc_nodes(const WaveField &field) const;
  
  void setParams(const S2PolyWaveApprox::Params &params) {m_params = params;};
  
  const S2PolyWaveApprox::Params& getParams() const {return m_params;}
  
protected:
  
  Params m_params;
  
private:
  
};

#endif /* s2_poly_wave_approx.hpp */