lib/polyfit.h

Go to the documentation of this file.

//
// canola - canon canola 1614p emulator
// Copyright (C) 2012 Peter Miller
//
// This program is free software; you can redistribute it and/or modify
// it under the terms of the GNU General Public License, version 3, as
// published by the Free Software Foundation.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// General Public License for more details.
//
// You should have received a copy of the GNU General Public License along
// with this program. If not, see <http://www.gnu.org/licenses/>.
//

#ifndef LIB_POLYFIT_H
#define LIB_POLYFIT_H

class polyfit
{
public:
    virtual ~polyfit();

    polyfit(int order = 1);

    void enter(double x, double y);

    enum solve_t
    {
        solve_ok,
        solve_singular,
        solve_insufficient
    };

    solve_t solve(void);

    double coefficient(int n);

    double evaluate(double x);

    void print(void) const;

    bool enough_points(void) const { return matrix[0][0] >= rows; }

private:
    unsigned int rows;

    unsigned int cols;

    typedef long double data_t;

    data_t *coefficients;

    data_t **matrix;

    static data_t epsilon;

    data_t **mat_alloc();

    void mat_dump(data_t **data);

    void mat_free(data_t **data);

    polyfit();

    polyfit(const polyfit &);

    polyfit &operator=(const polyfit &);
};

#endif // LIB_POLYFIT_H