35 #ifndef TEMPLATE_LAPACK_LAEV2_HEADER 36 #define TEMPLATE_LAPACK_LAEV2_HEADER 41 Treal *rt1, Treal *rt2, Treal *cs1, Treal *sn1)
108 Treal acmn, acmx, ab, df, cs, ct, tb, sm, tn, rt, adf, acs;
128 }
else if (adf < ab) {
139 *rt1 = (sm - rt) * .5;
146 *rt2 = acmx / *rt1 * acmn - *b / *rt1 * *b;
147 }
else if (sm > 0.) {
148 *rt1 = (sm + rt) * .5;
155 *rt2 = acmx / *rt1 * acmn - *b / *rt1 * *b;
#define absMACRO(x)
Definition: template_blas_common.h:45
int template_lapack_laev2(Treal *a, Treal *b, Treal *c__, Treal *rt1, Treal *rt2, Treal *cs1, Treal *sn1)
Definition: template_lapack_laev2.h:40
int integer
Definition: template_blas_common.h:38
Treal template_blas_sqrt(Treal x)