23 template <
class Size,
class OutputIterator1,
class OutputIterator2>
25 typedef typename std::iterator_traits<OutputIterator1>::value_type
T;
26 static T pi = acos(
T(-1.));
28 omega [0] = pi/
T(2.*R-2);
29 for (Size r = 1; r < R-1; r++) {
30 zeta [r] = cos(r*pi/
T(R-1.));
31 omega[r] = pi/
T(R-1.);
34 omega [R-1] = pi/
T(2.*R-2);
void gauss_lobatto_chebyschev(Size R, OutputIterator1 zeta, OutputIterator2 omega)