37 #ifndef TEMPLATE_LAPACK_STEQR_HEADER
38 #define TEMPLATE_LAPACK_STEQR_HEADER
43 Treal *e, Treal *z__,
const integer *ldz, Treal *work,
126 integer z_dim1, z_offset, i__1, i__2;
148 #define z___ref(a_1,a_2) z__[(a_2)*z_dim1 + a_1]
154 z_offset = 1 + z_dim1 * 1;
174 }
else if (*ldz < 1 || (icompz > 0 && *ldz <
maxMACRO(1,*n) ) ) {
203 safmax = 1. / safmin;
233 for (m = l1; m <= i__1; ++m) {
234 tst = (d__1 = e[m],
absMACRO(d__1));
266 if (anorm > ssfmax) {
274 }
else if (anorm < ssfmin) {
301 for (m = l; m <= i__1; ++m) {
303 d__2 = (d__1 = e[m],
absMACRO(d__1));
305 if (tst <= eps2 * (d__1 = d__[m],
absMACRO(d__1)) * (d__2 = d__[m
331 work[*n - 1 + l] = s;
347 if (jtot == nmaxit) {
354 g = (d__[l + 1] - p) / (e[l] * 2.);
366 for (i__ = mm1; i__ >= i__1; --i__) {
373 g = d__[i__ + 1] - p;
374 r__ = (d__[i__] - g) * s + c__ * 2. * b;
376 d__[i__ + 1] = g + p;
383 work[*n - 1 + i__] = -s;
422 for (m = l; m >= i__1; --m) {
424 d__2 = (d__1 = e[m - 1],
absMACRO(d__1));
426 if (tst <= eps2 * (d__1 = d__[m],
absMACRO(d__1)) * (d__2 = d__[m
453 work[*n - 1 + m] = s;
469 if (jtot == nmaxit) {
476 g = (d__[l - 1] - p) / (e[l - 1] * 2.);
488 for (i__ = m; i__ <= i__1; ++i__) {
496 r__ = (d__[i__ + 1] - g) * s + c__ * 2. * b;
505 work[*n - 1 + i__] = s;
540 i__1 = lendsv - lsv + 1;
546 }
else if (iscale == 2) {
547 i__1 = lendsv - lsv + 1;
562 for (i__ = 1; i__ <= i__1; ++i__) {
584 for (ii = 2; ii <= i__1; ++ii) {
589 for (j = ii; j <= i__2; ++j) {