bezier_curve.hh

#ifndef VTESTBED_GEOMETRY_BEZIER_CURVE_HH
#define VTESTBED_GEOMETRY_BEZIER_CURVE_HH

#include <vtestbed/base/blitz.hh>

namespace vtb {

    namespace geometry {

        template <class T, int N>
        class Bezier_curve {

        public:
            using value_type = Vertex<T,N>;
            using reference = value_type&;
            using const_reference = const value_type&;
            using container_type = std::vector<value_type>;
            using size_type = size_t;

        private:
            const container_type& _values;

        public:

            inline explicit
            Bezier_curve(const container_type& values):
            _values(values) {}

            inline value_type
            operator()(T t) const {
                T t1 = T{1} - t;
                auto p = this->_values;
                int n = p.size()-1;
                for (int j=1; j<=n; ++j) {
                    for (int i=0; i<=n-j; ++i) {
                        p[i] = p[i]*t1 + p[i+1]*t;
                    }
                }
                return p[0];
            }

        };

    }

}

#endif // vim:filetype=cpp