CRNLinearInterpolation.h

Go to the documentation of this file.

/* Copyright 2012-2016 CoReNum, ENS-Lyon
 * 
 * This file is part of libcrn.
 * 
 * libcrn is free software: you can redistribute it and/or modify
 * it under the terms of the GNU Lesser General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 * 
 * libcrn 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 Lesser General Public License for more details.
 * 
 * You should have received a copy of the GNU Lesser General Public License
 * along with libcrn.  If not, see <http://www.gnu.org/licenses/>.
 * 
 * file: CRNLinearInterpolation.h
 * \author Yann LEYDIER
 */

#ifndef CRNLinearInterpolation_HEADER
#define CRNLinearInterpolation_HEADER

#include <CRNMath/CRNCartesian2DFunction.h>
#include <CRNGeometry/CRNPoint2DDouble.h>
#include <CRNException.h>

namespace crn
{
    class LinearInterpolation: public crn::Cartesian2DFunction
    {
        public:
            template<typename Iter> LinearInterpolation(Iter beg, Iter en):
                data(std::distance(beg, en))
            {
                if (data.size() < 2)
                    throw crn::ExceptionDimension("There must be at least 2 samples to interpolate anything.");
                size_t cnt = 0;
                for (Iter tmp = beg; tmp != en; ++tmp)
                    data[cnt++] = makePoint(*tmp);
                std::sort(data.begin(), data.end(), Point2DDouble::Comparer(Direction::LEFT));
            }
            LinearInterpolation(const LinearInterpolation&) = default;
            LinearInterpolation(LinearInterpolation&&) = default;
            virtual ~LinearInterpolation() override {}
            LinearInterpolation& operator=(const LinearInterpolation&) = default;
            LinearInterpolation& operator=(LinearInterpolation&&) = default;

            virtual double operator[](double x) const override;
            using Cartesian2DFunction::operator[];

            const std::vector<crn::Point2DDouble>& GetData() const noexcept { return data; }

            void Deserialize(xml::Element &el);
            xml::Element Serialize(xml::Element &parent) const;

        private:
            template<typename T> Point2DDouble makePoint(const crn::Point2D<T> &p)
            { return crn::Point2DDouble(double(p.X), double(p.Y)); }
            template<typename T> Point2DDouble makePoint(const crn::Point2D<T> *p)
            { return crn::Point2DDouble(double(p->X), double(p->Y)); }
            template<typename T, typename Y> Point2DDouble makePoint(const std::pair<T, Y> &p)
            { return crn::Point2DDouble(double(p.first), double(p.second)); }
            
            std::vector<crn::Point2DDouble> data; 
        CRN_DECLARE_CLASS_CONSTRUCTOR(LinearInterpolation)          
        CRN_SERIALIZATION_CONSTRUCTOR(LinearInterpolation)          
    };
    template<> struct IsSerializable<LinearInterpolation> : public std::true_type {};
    template<> struct IsClonable<LinearInterpolation> : public std::true_type {};

    CRN_ALIAS_SMART_PTR(LinearInterpolation)
}

#endif