geometry/ruled_surface.hh

#ifndef VTESTBED_GEOMETRY_RULED_SURFACE_HH
#define VTESTBED_GEOMETRY_RULED_SURFACE_HH

#include <algorithm>
#include <vector>

namespace vtb {

    namespace geometry {

        template <class T, class Curve>
        class Ruled_surface {

        public:
            using curve_type = Curve;
            using value_type = typename curve_type::value_type;
            using size_type = size_t;

        private:
            curve_type _curve1, _curve2;

        public:

            inline explicit
            Ruled_surface(const curve_type& c, const curve_type& d):
            _curve1(c), _curve2(d) {}

            inline value_type
            operator()(T u, T v) const {
                const auto& c = this->_curve1;
                const auto& d = this->_curve2;
                return (T{1}-v)*c(u) + v*d(u);
            }

            inline size_type
            max_size() const {
                return std::max(this->_curve1.size(), this->_curve2.size());
            }

        };

    }

}

#endif // vim:filetype=cpp