37 #ifndef TEMPLATE_BLAS_BASICMATH_HEADER 38 #define TEMPLATE_BLAS_BASICMATH_HEADER 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);
Treal template_blas_cos(Treal x)
Treal template_blas_pow(Treal x, Treal y)
Treal template_blas_erfc(Treal x)
Treal template_blas_fabs(Treal x)
Treal template_blas_exp(Treal x)
Treal template_blas_log(Treal x)
Treal template_blas_sin(Treal x)
Treal template_blas_log10(Treal x)
Treal template_blas_erf(Treal x)
Treal template_blas_compute_pi_BBP(Treal dummy)
Definition: template_blas_basicmath.h:77
Treal template_blas_sqrt(Treal x)