2 #include "gmx_lapack.h"
6 F77_FUNC(dlae2,DLAE2)(double *a,
13 double ab, df, tb, sm, rt, adf, acmn, acmx;
21 if (fabs(*a) > fabs(*c__)) {
30 rt = adf * sqrt(d__1 * d__1 + 1.);
31 } else if (adf < ab) {
33 rt = ab * sqrt(d__1 * d__1 + 1.);
39 *rt1 = (sm - rt) * .5;
40 *rt2 = acmx / *rt1 * acmn - *b / *rt1 * *b;
42 *rt1 = (sm + rt) * .5;
43 *rt2 = acmx / *rt1 * acmn - *b / *rt1 * *b;