Loading...
Searching...
No Matches
gauss_radau_chebyschev.icc
Go to the documentation of this file.
#include <cmath>
#include <iterator>
template <class Size, class OutputIterator1, class OutputIterator2>
void gauss_radau_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-1);
for (Size r = 1; r < R; r++) {
zeta [r] = cos(2*r*pi/T(2.*R-1));
omega[r] = 2*pi/T(2.*R-1);
}
}
T
Expr1::float_type T
Definition field_expr.h:230
gauss_radau_chebyschev
void gauss_radau_chebyschev(Size R, OutputIterator1 zeta, OutputIterator2 omega)
Definition gauss_radau_chebyschev.icc:24