spline_function.hpp

Go to the documentation of this file.

/*****************************************************************************
** Ifdefs
*****************************************************************************/

#ifndef ECL_GEOMETRY_SPLINE_FUNCTION_HPP_
#define ECL_GEOMETRY_SPLINE_FUNCTION_HPP_

/*****************************************************************************
** Includes
*****************************************************************************/

#include <ecl/config/macros.hpp>
#include <ecl/containers/array.hpp>
#include <ecl/exceptions/standard_exception.hpp>

/*****************************************************************************
** Namespaces
*****************************************************************************/

namespace ecl {

/*****************************************************************************
** Interface [GenericSplineFunction]
*****************************************************************************/
class ECL_PUBLIC GenericSplineFunction {
    public:
                virtual ~GenericSplineFunction() {};
        virtual double operator()(const double &x) const = 0;
        virtual double derivative(const double &x) const = 0;
        virtual double dderivative(const double &x) const = 0;

        const Array<double,2>& domain() { return time_domain; }

    protected:
        Array<double,2> time_domain;
};
/*****************************************************************************
** Interface [SplineFunction]
*****************************************************************************/
template <typename Function>
class ECL_PUBLIC SplineFunction : public GenericSplineFunction {
    public:
        SplineFunction (const double& time_begin, const double& time_end, const Function& f) : function(f) {
            time_domain << time_begin, time_end;
        };
        virtual ~SplineFunction() {};

        double operator()(const double &x) const ecl_assert_throw_decl(StandardException) {
            ecl_assert_throw( ( x >= time_domain[0] ) && x <= time_domain[1], StandardException(LOC,OutOfRangeError ) );
            return function(x);
        }

        double derivative(const double &x) const ecl_assert_throw_decl(StandardException) {
            ecl_assert_throw( ( x >= time_domain[0] ) && x <= time_domain[1], StandardException(LOC,OutOfRangeError ) );
            return function.derivative(x);
        }

        double dderivative(const double &x) const ecl_assert_throw_decl(StandardException) {
            ecl_assert_throw( ( x >= time_domain[0] ) && x <= time_domain[1], StandardException(LOC,OutOfRangeError ) );
            return function.dderivative(x);
        }

    private:
        Function function;
};

}; // namespace ecl


#endif /* SPLINE_FUNCTION_HPP_ */