13 #ifndef eli_geom_curve_torsion_hpp
14 #define eli_geom_curve_torsion_hpp
23 template<
typename curve__>
24 void torsion(
typename curve__::data_type &tau,
const curve__ &c,
const typename curve__::data_type &t)
30 assert((t>=0) && (t<=1));
32 typename curve__::point_type xp(c.fp(t)), xpp(c.fpp(t)), xppp(c.fppp(t)), tmp1;
33 typename curve__::data_type tmp2;
37 tau=tmp1.dot(xppp)/(tmp2*tmp2);
void torsion(typename curve__::data_type &tau, const curve__ &c, const typename curve__::data_type &t)
Definition: torsion.hpp:24