35 #ifndef TEMPLATE_BLAS_BASICMATH_HEADER
36 #define TEMPLATE_BLAS_BASICMATH_HEADER
65 Treal epsilon = std::numeric_limits<Treal>::epsilon();
66 Treal one_over_16 = (Treal)1 / (Treal)16;
67 Treal one_over_16_to_pow_k = 1;
73 (Treal)4 / (Treal)(8*k + 1) -
74 (Treal)2 / (Treal)(8*k + 4) -
75 (Treal)1 / (Treal)(8*k + 5) -
76 (Treal)1 / (Treal)(8*k + 6);
77 sum += one_over_16_to_pow_k * factor;
79 one_over_16_to_pow_k *= one_over_16;
81 while(one_over_16_to_pow_k > epsilon);