least_squares_2d.hpp

// $Id: least_squares_2d.hpp,v 1.2 2009/08/03 15:58:36 urij Exp $

#ifndef __least_squares_2d_hpp
#define __least_squares_2d_hpp

#include <vector>
#include <cstring>

namespace LeastSquares2D
{
  using std::vector;
  
  struct PolyCoeff 
  {
    PolyCoeff() : c(0.), fixed(false) {};
    PolyCoeff(double i_c, bool i_fixed = false) {c = i_c; fixed = i_fixed;}
    
    double c;     
    bool fixed;   
  };
  
  typedef std::vector<PolyCoeff> PolyCoeffs; 
  
  bool getPoly(const vector<double> &x, 
           const vector<double> &y, 
           const vector<double> &w, 
           PolyCoeffs &cfs);
  
  bool getPoly(const vector<double> &x, 
           const vector<double> &y, 
           PolyCoeffs &cfs);
  
  double calcPoly(const PolyCoeffs &coeffs, double x);
    
  double ls_m_elem(size_t i, size_t j, const vector<double> &x, const vector<double> &w);
  
  double ls_b_elem(size_t i, const vector<double> &x, const vector<double> &y, const vector<double> &w);
};

#endif //__least_squares_2d_hpp