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