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) {
#define absMACRO(x)
Definition: template_blas_common.h:47
#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:42
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:42
int template_lapack_laev2(Treal *a, Treal *b, Treal *c__, Treal *rt1, Treal *rt2, Treal *cs1, Treal *sn1)
Definition: template_lapack_laev2.h:42
int integer
Definition: template_blas_common.h:40
Treal template_lapack_lapy2(Treal *x, Treal *y)
Definition: template_lapack_lapy2.h:42
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
int template_lapack_lasrt(const char *id, const integer *n, Treal *d__, integer *info)
Definition: template_lapack_lasrt.h:42
Treal template_lapack_d_sign(const Treal *a, const Treal *b)
Definition: template_lapack_lamch.h:48
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:42
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
int template_blas_swap(const integer *n, Treal *dx, const integer *incx, Treal *dy, const integer *incy)
Definition: template_blas_swap.h:42
int template_lapack_lartg(const Treal *f, const Treal *g, Treal *cs, Treal *sn, Treal *r__)
Definition: template_lapack_lartg.h:42
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:42
Treal template_lapack_lamch(const char *cmach, Treal dummyReal)
Definition: template_lapack_lamch.h:202
int template_lapack_lae2(const Treal *a, const Treal *b, const Treal *c__, Treal *rt1, Treal *rt2)
Definition: template_lapack_lae2.h:42
Treal template_lapack_lanst(const char *norm, const integer *n, const Treal *d__, const Treal *e)
Definition: template_lapack_lanst.h:42
Treal template_blas_sqrt(Treal x)
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:46