37 #ifndef TEMPLATE_LAPACK_LAEV2_HEADER 38 #define TEMPLATE_LAPACK_LAEV2_HEADER 43 Treal *rt1, Treal *rt2, Treal *cs1, Treal *sn1)
110 Treal acmn, acmx, ab, df, cs, ct, tb, sm, tn, rt, adf, acs;
130 }
else if (adf < ab) {
141 *rt1 = (sm - rt) * .5;
148 *rt2 = acmx / *rt1 * acmn - *b / *rt1 * *b;
149 }
else if (sm > 0.) {
150 *rt1 = (sm + rt) * .5;
157 *rt2 = acmx / *rt1 * acmn - *b / *rt1 * *b;
#define absMACRO(x)
Definition: template_blas_common.h:47
int template_lapack_laev2(Treal *a, Treal *b, Treal *c__, Treal *rt1, Treal *rt2, Treal *cs1, Treal *sn1)
Definition: template_lapack_laev2.h:42
int integer
Definition: template_blas_common.h:40
Treal template_blas_sqrt(Treal x)