79 Treal epsilon = template_blas_get_machine_epsilon<Treal>();
80 Treal one_over_16 = (Treal)1 / (Treal)16;
81 Treal one_over_16_to_pow_k = 1;
87 (Treal)4 / (Treal)(8*k + 1) -
88 (Treal)2 / (Treal)(8*k + 4) -
89 (Treal)1 / (Treal)(8*k + 5) -
90 (Treal)1 / (Treal)(8*k + 6);
91 sum += one_over_16_to_pow_k * factor;
93 one_over_16_to_pow_k *= one_over_16;
95 while(one_over_16_to_pow_k > epsilon);