InterpolationCubic.h

Go to the documentation of this file.

/*=============================================================================
TexGen: Geometric textile modeller.
Copyright (C) 2006 Martin Sherburn
 
This program is free software; you can redistribute it and/or
modify it under the terms of the GNU General Public License
as published by the Free Software Foundation; either version 2
of the License, or (at your option) any later version.
 
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, write to the Free Software
Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
=============================================================================*/
 
#pragma once
#include "Interpolation.h"
 
namespace TexGen
{
    using namespace std;
 
    class CLASS_DECLSPEC CInterpolationCubic :  public CInterpolation
    {
    public:
        CInterpolationCubic(bool bPeriodic = true, bool bForceInPlaneTangent = false, bool bForceMasterNodeTangent = false);
        CInterpolationCubic(TiXmlElement &Element);
        ~CInterpolationCubic(void);
 
        CInterpolation* Copy() const { return new CInterpolationCubic(*this); }
        string GetType() const { return "CInterpolationCubic"; }
//      void PopulateTiXmlElement(TiXmlElement &Element, OUTPUT_TYPE OutputType) const;
 
        void Initialise(const vector<CNode> &MasterNodes) const;
 
        CSlaveNode GetNode(const vector<CNode> &MasterNodes, int iIndex, double t) const;
 
    protected:
        struct CUBICEQUATION
        {
            CUBICEQUATION(double a, double b, double c, double d)
            : m_a(a)
            , m_b(b)
            , m_c(c)
            , m_d(d)
            {}
 
            double Evaluate(double x) const
            {
                return (((m_d*x)+m_c)*x+m_b)*x+m_a;
            }
 
            double EvaluateDerivative(double x) const
            {
                return ((3*m_d*x)+2*m_c)*x+m_b;
            }
 
        protected:
            double m_a, m_b, m_c, m_d;
        };
 
        static void GetPeriodicCubicSplines(const vector<double> &Knots, vector<CUBICEQUATION> &Cubics);
        static void GetNaturalCubicSplines(const vector<double> &Knots, vector<CUBICEQUATION> &Cubics);
 
        mutable vector<CUBICEQUATION> m_XCubics;
        mutable vector<CUBICEQUATION> m_YCubics;
        mutable vector<CUBICEQUATION> m_ZCubics;
    };
 
};  // namespace TexGen