35 #ifndef TEMPLATE_LAPACK_LAE2_HEADER 36 #define TEMPLATE_LAPACK_LAE2_HEADER 41 Treal *rt1, Treal *rt2)
97 Treal acmn, acmx, ab, df, tb, sm, rt, adf;
116 }
else if (adf < ab) {
127 *rt1 = (sm - rt) * .5;
133 *rt2 = acmx / *rt1 * acmn - *b / *rt1 * *b;
134 }
else if (sm > 0.) {
135 *rt1 = (sm + rt) * .5;
141 *rt2 = acmx / *rt1 * acmn - *b / *rt1 * *b;
#define absMACRO(x)
Definition: template_blas_common.h:45
int template_lapack_lae2(const Treal *a, const Treal *b, const Treal *c__, Treal *rt1, Treal *rt2)
Definition: template_lapack_lae2.h:40
Treal template_blas_sqrt(Treal x)