gauss_chebyschev.icc
Go to the documentation of this file.
#include <cmath>
#include <iterator>
template <class Size, class OutputIterator1, class OutputIterator2>
void gauss_chebyschev (Size R, OutputIterator1 zeta, OutputIterator2 omega) {
typedef typename std::iterator_traits<OutputIterator1>::value_type T;
static T pi = acos(T(-1.));
for (Size r = 0; r < R; r++) {
zeta [r] = cos((2*r+1)*pi/(2.*R));
omega[r] = pi/T(1.*R);
}
}
gauss_chebyschev
void gauss_chebyschev(Size R, OutputIterator1 zeta, OutputIterator2 omega)
Definition: gauss_chebyschev.icc:24
T
Expr1::float_type T
Definition: field_expr.h:218