1 #include "rheolef/compiler.h"
22 #include "rheolef/gauss_chebyschev.h"
24 template <
class Iterator1,
class Iterator2,
class Size>
25 typename std::iterator_traits<Iterator1>::value_type
27 Size R_e, Iterator1 zeta_e, Iterator2 omega_e,
29 typedef typename std::iterator_traits<Iterator1>::value_type
T;
31 for (Size k = 0; k <=
order; k++) {
32 T exact = 0, value = 0;
33 for (Size r = 0; r < R_e; r++)
34 exact +=
pow (zeta_e[r], k)*omega_e[r];
35 for (Size r = 0; r < R; r++)
36 value +=
pow (zeta[r], k)*omega[r];
37 error = std::max(error, fabs (exact - value));