source: trunk/lib/statistics/Polynomial.h @ 389

// $Id: Polynomial.h 389 2005-08-15 11:37:07Z peter $

#ifndef _theplu_statistics_regression_polynomial_
#define _theplu_statistics_regression_polynomial_

#include <c++_tools/statistics/OneDimensional.h>
#include <c++_tools/statistics/MultiDimensional.h>

#include <c++_tools/gslapi/vector.h>

#include <gsl/gsl_multifit.h>


namespace theplu {
namespace statistics {
namespace regression {

  ///
  /// @todo document
  ///
  class Polynomial : public OneDimensional
  {
  public:

    ///
    ///
    ///
    inline Polynomial(size_t power)
      : OneDimensional(), power_(power) {}

    ///
    /// @brief Destructor
    ///
    inline ~Polynomial(void) {};

    ///
    ///
    ///
    void fit(const gslapi::vector& x, const gslapi::vector& y);

    ///
    ///
    ///
    gslapi::vector fit_parameters(void) { return md_.fit_parameters(); }

    inline void fit(const gslapi::vector& x, const gslapi::vector& y,
                    const gslapi::vector& w) { assert(0); }

    inline void predict(const double x, double& y, double& y_err,
                         const double w=1) { assert(0); }

    inline std::ostream& print(std::ostream& s) const
      { assert(0); return s; }

    inline std::ostream& print_header(std::ostream& s) const
      { assert(0); return s; }


  private:
    MultiDimensional md_;
    size_t power_;

  };


}}} // of namespaces regression, statisitcs and thep

#endif