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;