Go to the documentation of this file.
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_pow(Treal x, Treal y)
Treal template_blas_sqrt(Treal x)
float template_blas_sqrt< float >(float x)
Definition: template_blas_basicmath.cc:81
float template_blas_exp< float >(float x)
Definition: template_blas_basicmath.cc:107
double template_blas_sin< double >(double x)
Definition: template_blas_basicmath.cc:240
long double template_blas_log10< long double >(long double x)
Definition: template_blas_basicmath.cc:167
Treal template_blas_cos(Treal x)
double template_blas_log< double >(double x)
Definition: template_blas_basicmath.cc:136
Treal template_blas_exp(Treal x)
long double template_blas_cos< long double >(long double x)
Definition: template_blas_basicmath.cc:271
Treal template_blas_erfc(Treal x)
Treal template_blas_fabs(Treal x)
long double template_blas_sqrt< long double >(long double x)
Definition: template_blas_basicmath.cc:89
float template_blas_fabs< float >(float x)
Definition: template_blas_basicmath.cc:55
double template_blas_exp< double >(double x)
Definition: template_blas_basicmath.cc:110
float template_blas_erf< float >(float x)
Definition: template_blas_basicmath.cc:185
double template_blas_log10< double >(double x)
Definition: template_blas_basicmath.cc:162
double template_blas_fabs< double >(double x)
Definition: template_blas_basicmath.cc:58
float template_blas_erfc< float >(float x)
Definition: template_blas_basicmath.cc:211
Treal template_blas_log10(Treal x)
float template_blas_log< float >(float x)
Definition: template_blas_basicmath.cc:133
Treal template_blas_erf(Treal x)
double template_blas_erfc< double >(double x)
Definition: template_blas_basicmath.cc:214
double template_blas_pow< double >(double x, double y)
Definition: template_blas_basicmath.cc:292
long double template_blas_erfc< long double >(long double x)
Definition: template_blas_basicmath.cc:219
Treal template_blas_sin(Treal x)
float template_blas_cos< float >(float x)
Definition: template_blas_basicmath.cc:263
float template_blas_sin< float >(float x)
Definition: template_blas_basicmath.cc:237
Treal template_blas_log(Treal x)
double template_blas_erf< double >(double x)
Definition: template_blas_basicmath.cc:188
float template_blas_pow< float >(float x, float y)
Definition: template_blas_basicmath.cc:289
long double template_blas_erf< long double >(long double x)
Definition: template_blas_basicmath.cc:193
Treal template_blas_compute_pi_BBP(Treal dummy)
Definition: template_blas_basicmath.h:77
double template_blas_cos< double >(double x)
Definition: template_blas_basicmath.cc:266
long double template_blas_pow< long double >(long double x, long double y)
Definition: template_blas_basicmath.cc:297
long double template_blas_sin< long double >(long double x)
Definition: template_blas_basicmath.cc:245
long double template_blas_fabs< long double >(long double x)
Definition: template_blas_basicmath.cc:63
float template_blas_log10< float >(float x)
Definition: template_blas_basicmath.cc:159
long double template_blas_log< long double >(long double x)
Definition: template_blas_basicmath.cc:141
double template_blas_sqrt< double >(double x)
Definition: template_blas_basicmath.cc:84
long double template_blas_exp< long double >(long double x)
Definition: template_blas_basicmath.cc:115