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) {
#define absMACRO(x)
Definition: template_blas_common.h:45
#define z___ref(a_1, a_2)
int template_lapack_lascl(const char *type__, const integer *kl, const integer *ku, const Treal *cfrom, const Treal *cto, const integer *m, const integer *n, Treal *a, const integer *lda, integer *info)
Definition: template_lapack_lascl.h:40
int template_lapack_laset(const char *uplo, const integer *m, const integer *n, const Treal *alpha, const Treal *beta, Treal *a, const integer *lda)
Definition: template_lapack_laset.h:40
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_lapack_lapy2(Treal *x, Treal *y)
Definition: template_lapack_lapy2.h:40
#define maxMACRO(a, b)
Definition: template_blas_common.h:43
int template_lapack_lasrt(const char *id, const integer *n, Treal *d__, integer *info)
Definition: template_lapack_lasrt.h:40
Treal template_lapack_d_sign(const Treal *a, const Treal *b)
Definition: template_lapack_lamch.h:46
int template_lapack_steqr(const char *compz, const integer *n, Treal *d__, Treal *e, Treal *z__, const integer *ldz, Treal *work, integer *info)
Definition: template_lapack_steqr.h:40
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:144
int template_blas_swap(const integer *n, Treal *dx, const integer *incx, Treal *dy, const integer *incy)
Definition: template_blas_swap.h:40
int template_lapack_lartg(const Treal *f, const Treal *g, Treal *cs, Treal *sn, Treal *r__)
Definition: template_lapack_lartg.h:40
int template_lapack_lasr(const char *side, const char *pivot, const char *direct, const integer *m, const integer *n, const Treal *c__, const Treal *s, Treal *a, const integer *lda)
Definition: template_lapack_lasr.h:40
Treal template_lapack_lamch(const char *cmach, Treal dummyReal)
Definition: template_lapack_lamch.h:199
int template_lapack_lae2(const Treal *a, const Treal *b, const Treal *c__, Treal *rt1, Treal *rt2)
Definition: template_lapack_lae2.h:40
Treal template_lapack_lanst(const char *norm, const integer *n, const Treal *d__, const Treal *e)
Definition: template_lapack_lanst.h:40
Treal template_blas_sqrt(Treal x)
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:44