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;