35 #ifndef TEMPLATE_LAPACK_STEQR_HEADER
36 #define TEMPLATE_LAPACK_STEQR_HEADER
41 Treal *e, Treal *z__,
const integer *ldz, Treal *work,
124 integer z_dim1, z_offset, i__1, i__2;
146 #define z___ref(a_1,a_2) z__[(a_2)*z_dim1 + a_1]
152 z_offset = 1 + z_dim1 * 1;
172 }
else if (*ldz < 1 || (icompz > 0 && *ldz <
maxMACRO(1,*n) ) ) {
201 safmax = 1. / safmin;
231 for (m = l1; m <= i__1; ++m) {
232 tst = (d__1 = e[m],
absMACRO(d__1));
264 if (anorm > ssfmax) {
272 }
else if (anorm < ssfmin) {
299 for (m = l; m <= i__1; ++m) {
301 d__2 = (d__1 = e[m],
absMACRO(d__1));
303 if (tst <= eps2 * (d__1 = d__[m],
absMACRO(d__1)) * (d__2 = d__[m
329 work[*n - 1 + l] = s;
345 if (jtot == nmaxit) {
352 g = (d__[l + 1] - p) / (e[l] * 2.);
364 for (i__ = mm1; i__ >= i__1; --i__) {
371 g = d__[i__ + 1] - p;
372 r__ = (d__[i__] - g) * s + c__ * 2. * b;
374 d__[i__ + 1] = g + p;
381 work[*n - 1 + i__] = -s;
420 for (m = l; m >= i__1; --m) {
422 d__2 = (d__1 = e[m - 1],
absMACRO(d__1));
424 if (tst <= eps2 * (d__1 = d__[m],
absMACRO(d__1)) * (d__2 = d__[m
451 work[*n - 1 + m] = s;
467 if (jtot == nmaxit) {
474 g = (d__[l - 1] - p) / (e[l - 1] * 2.);
486 for (i__ = m; i__ <= i__1; ++i__) {
494 r__ = (d__[i__ + 1] - g) * s + c__ * 2. * b;
503 work[*n - 1 + i__] = s;
538 i__1 = lendsv - lsv + 1;
544 }
else if (iscale == 2) {
545 i__1 = lendsv - lsv + 1;
560 for (i__ = 1; i__ <= i__1; ++i__) {
582 for (ii = 2; ii <= i__1; ++ii) {
587 for (j = ii; j <= i__2; ++j) {