CubicBezierCurve.h

/***************************************************************************
 *   Copyright (C) 2004 by Jacques Gasselin                                *
 *   jacquesgasselin@hotmail.com                                           *
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU Library 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 Library General Public     *
 *   License along with this program; if not, write to the                 *
 *   Free Software Foundation, Inc.,                                       *
 *   59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.             *
 ***************************************************************************/
#ifndef GLCUBICBEZIERCURVE_H
#define GLCUBICBEZIERCURVE_H


#include "GLMatrix.h"

namespace mathglpp
{

template <typename T>
class CubicBezierCurve4
{
public:
    enum { D = 4 };
    typedef T value_type;
    typedef GLMatrix<value_type> matrix_type;

    CubicBezierCurve4(const value_type* cpoints)
    :basis(cpoints)
    {
        static value_type coeffs[D*D] =
            {  1, -0,  0, -0,
              -3,  3, -0,  0,
               3, -6,  3, -0,
              -1,  3, -3,  1 };

        basis *= GLMatrix<value_type>(coeffs);
    }

    CubicBezierCurve4(const value_type* c0, const value_type* c1, const value_type* c2, const value_type* c3)
    :basis(c0,c1,c2,c3)
    {
        static value_type coeffs[16] =
            {  1, -0,  0, -0,
              -3,  3, -0,  0,
               3, -6,  3, -0,
              -1,  3, -3,  1 };


        basis *= GLMatrix<value_type>(coeffs);
    }

    CubicBezierCurve4(const CubicBezierCurve4& bc)
    :basis(bc.basis){}

    ~CubicBezierCurve4(){}

    inline void interpolate(GLclampf interp, value_type* vertex) const
    {
        vertex[0] = 1;
        vertex[1] = value_type(interp);
        vertex[2] = value_type(interp * interp);
        vertex[3] = value_type(interp * interp * interp);
        basis.vdot4(vertex);
    }

private:
    GLMatrix<value_type> basis;
};


typedef CubicBezierCurve4<GLbyte> CubicBezierCurve4b;
typedef CubicBezierCurve4<GLshort> CubicBezierCurve4s;
typedef CubicBezierCurve4<GLint> CubicBezierCurve4i;
typedef CubicBezierCurve4<GLfloat> CubicBezierCurve4f;
typedef CubicBezierCurve4<GLdouble> CubicBezierCurve4d;

};

#endif

---