35 #ifndef TEMPLATE_LAPACK_SYTRD_HEADER 36 #define TEMPLATE_LAPACK_SYTRD_HEADER 42 lda, Treal *d__, Treal *e, Treal *tau, Treal *
177 integer a_dim1, a_offset, i__1, i__2, i__3;
186 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1] 190 a_offset = 1 + a_dim1 * 1;
202 lquery = *lwork == -1;
209 }
else if (*lwork < 1 && ! lquery) {
220 work[1] = (Treal) lwkopt;
240 if (nb > 1 && nb < *n) {
262 i__1 = *lwork / ldwork;
282 kk = *n - (*n - nx + nb - 1) / nb * nb;
285 for (i__ = *n - nb + 1; i__2 < 0 ? i__ >= i__1 : i__ <= i__1; i__ +=
301 lda, &work[1], &ldwork, &c_b23, &a[a_offset], lda);
307 for (j = i__; j <= i__3; ++j) {
308 a_ref(j - 1, j) = e[j - 1];
309 d__[j] =
a_ref(j, j);
324 for (i__ = 1; i__1 < 0 ? i__ >= i__2 : i__ <= i__2; i__ += i__1) {
332 i__], &work[1], &ldwork);
337 i__3 = *n - i__ - nb + 1;
339 i__), lda, &work[nb + 1], &ldwork, &c_b23, &
a_ref(i__ +
346 for (j = i__; j <= i__3; ++j) {
347 a_ref(j + 1, j) = e[j];
348 d__[j] =
a_ref(j, j);
361 work[1] = (Treal) lwkopt;
int template_lapack_sytd2(const char *uplo, const integer *n, Treal *a, const integer *lda, Treal *d__, Treal *e, Treal *tau, integer *info)
Definition: template_lapack_sytd2.h:41
int template_lapack_sytrd(const char *uplo, const integer *n, Treal *a, const integer *lda, Treal *d__, Treal *e, Treal *tau, Treal *work, const integer *lwork, integer *info)
Definition: template_lapack_sytrd.h:41
int integer
Definition: template_blas_common.h:38
integer template_lapack_ilaenv(const integer *ispec, const char *name__, const char *opts, const integer *n1, const integer *n2, const integer *n3, const integer *n4, ftnlen name_len, ftnlen opts_len)
Definition: template_lapack_common.cc:279
#define maxMACRO(a, b)
Definition: template_blas_common.h:43
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:144
int template_lapack_latrd(const char *uplo, const integer *n, const integer *nb, Treal *a, const integer *lda, Treal *e, Treal *tau, Treal *w, const integer *ldw)
Definition: template_lapack_latrd.h:40
bool logical
Definition: template_blas_common.h:39
int ftnlen
Definition: template_blas_common.h:40
int template_blas_syr2k(const char *uplo, const char *trans, const integer *n, const integer *k, const Treal *alpha, const Treal *a, const integer *lda, const Treal *b, const integer *ldb, const Treal *beta, Treal *c__, const integer *ldc)
Definition: template_blas_syr2k.h:40
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:44