1 #include "rheolef/compiler.h"
23 template <
class Iterator1,
class Iterator2,
class Size>
24 typename std::iterator_traits<Iterator1>::value_type
26 typedef typename std::iterator_traits<Iterator1>::value_type
T;
28 for (Size k = 0; k <=
order; k++) {
30 for (Size r = 0; r < R; r++) {
31 value +=
pow (zeta[r], k)*omega[r];
33 T exact = (k % 2 == 1) ?
T(0.) : 2/
T(k+1.);
34 error = std::max (error, fabs (exact - value));
36 typedef typename std::iterator_traits<Iterator1>::value_type
T;
std::iterator_traits< Iterator1 >::value_type gauss_legendre_check(Size R, Iterator1 zeta, Iterator2 omega, Size order)
space_mult_list< T, M > pow(const space_basic< T, M > &X, size_t n)