Rheolef
7.2
an efficient C++ finite element environment
Loading...
Searching...
No Matches
gauss_chebyschev.icc
Go to the documentation of this file.
1
21
#include <cmath>
22
#include <iterator>
23
template
<
class
Size,
class
OutputIterator1,
class
OutputIterator2>
24
void
gauss_chebyschev
(Size R, OutputIterator1 zeta, OutputIterator2 omega) {
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);
30
}
31
}
T
Expr1::float_type T
Definition
field_expr.h:230
gauss_chebyschev
void gauss_chebyschev(Size R, OutputIterator1 zeta, OutputIterator2 omega)
Definition
gauss_chebyschev.icc:24