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;