src/include/PointerSystem.h

/*
 * PointerSystem.h
 *
 *  Created on: Jan 22, 2013
 *      Author: romain
 */

#ifndef POINTERSYSTEM_H_
#define POINTERSYSTEM_H_

#include "DynamicalSystem.h"
#include "utils.h"
#include "Box.h"
#include "stdexcept"

namespace nltool {


class PointerSystem: public DynamicalSystem {
public:

        PointerSystem();

        PointerSystem(int dim, Number (*f)(const int x, const stdVector &pt, void* data),
                        Number (*df)(const int x, const stdVector &pt, const int dx, void* data),
                        Number (*ddf)(const int x, const stdVector &pt, const int dx,
                                        const int dy, void* data), void* data, bool discretized = false);

        virtual ~PointerSystem();

        Vector eval(const Vector& x);

        Vector evalCurvature(const Vector& x, int dx, int dy);

        Number getMaxCurvature(const unsigned i, const unsigned dx,
                        const unsigned dy, const Box& domain);

        static double CurvatureFunc(const std::vector<double> &x,
                        std::vector<double> &grad, void* f_data);
        // variable used by CurvatureFunc
        unsigned curvatureFunctionX;
        unsigned curvatureFunctionDx;
        unsigned curvatureFunctionDy;
        void* dataPtr;
private:
        Number (*fPtr)(const int x, const stdVector &pt, void* data); //function studied f
        Number (*dfPtr)(const int x, const stdVector &pt, const int dx, void* data); // df/dx
        Number (*ddfPtr)(const int x, const stdVector &pt, const int dx,
                        const int dy, void* data); // d^2f/dxdy


};

}  // namespace nltool


#endif /* POINTERSYSTEM_H_ */