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( ( x >= time_domain[0] ) && x <= time_domain[1], StandardException(LOC,OutOfRangeError ) );
            return function(x);
        }
 
        double derivative(const double &x) const {
            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( ( x >= time_domain[0] ) && x <= time_domain[1], StandardException(LOC,OutOfRangeError ) );
            return function.dderivative(x);
        }
 
    private:
        Function function;
};
 
}; // namespace ecl
 
 
#endif /* SPLINE_FUNCTION_HPP_ */