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.));
27 for (Size r = 0; r < R; r++) {
28 zeta [r] = cos((2*r+1)*pi/(2.*R));
29 omega[r] = pi/
T(1.*R);
void gauss_chebyschev(Size R, OutputIterator1 zeta, OutputIterator2 omega)