37 #ifndef TEMPLATE_LAPACK_SYTRD_HEADER 38 #define TEMPLATE_LAPACK_SYTRD_HEADER 44 lda, Treal *d__, Treal *e, Treal *tau, Treal *
179 integer a_dim1, a_offset, i__1, i__2, i__3;
188 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1] 192 a_offset = 1 + a_dim1 * 1;
204 lquery = *lwork == -1;
211 }
else if (*lwork < 1 && ! lquery) {
222 work[1] = (Treal) lwkopt;
242 if (nb > 1 && nb < *n) {
264 i__1 = *lwork / ldwork;
284 kk = *n - (*n - nx + nb - 1) / nb * nb;
287 for (i__ = *n - nb + 1; i__2 < 0 ? i__ >= i__1 : i__ <= i__1; i__ +=
303 lda, &work[1], &ldwork, &c_b23, &a[a_offset], lda);
309 for (j = i__; j <= i__3; ++j) {
310 a_ref(j - 1, j) = e[j - 1];
311 d__[j] =
a_ref(j, j);
326 for (i__ = 1; i__1 < 0 ? i__ >= i__2 : i__ <= i__2; i__ += i__1) {
334 i__], &work[1], &ldwork);
339 i__3 = *n - i__ - nb + 1;
341 i__), lda, &work[nb + 1], &ldwork, &c_b23, &
a_ref(i__ +
348 for (j = i__; j <= i__3; ++j) {
349 a_ref(j + 1, j) = e[j];
350 d__[j] =
a_ref(j, j);
363 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:43
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:43
int integer
Definition: template_blas_common.h:40
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:281
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
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:42
bool logical
Definition: template_blas_common.h:41
int ftnlen
Definition: template_blas_common.h:42
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:42
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:46