37 #ifndef TEMPLATE_LAPACK_LAE2_HEADER 38 #define TEMPLATE_LAPACK_LAE2_HEADER 43 Treal *rt1, Treal *rt2)
99 Treal acmn, acmx, ab, df, tb, sm, rt, adf;
118 }
else if (adf < ab) {
129 *rt1 = (sm - rt) * .5;
135 *rt2 = acmx / *rt1 * acmn - *b / *rt1 * *b;
136 }
else if (sm > 0.) {
137 *rt1 = (sm + rt) * .5;
143 *rt2 = acmx / *rt1 * acmn - *b / *rt1 * *b;
#define absMACRO(x)
Definition: template_blas_common.h:47
int template_lapack_lae2(const Treal *a, const Treal *b, const Treal *c__, Treal *rt1, Treal *rt2)
Definition: template_lapack_lae2.h:42
Treal template_blas_sqrt(Treal x)