35 #ifndef TEMPLATE_LAPACK_SYTD2_HEADER
36 #define TEMPLATE_LAPACK_SYTD2_HEADER
42 lda, Treal *d__, Treal *e, Treal *tau,
integer *info)
161 integer a_dim1, a_offset, i__1, i__2, i__3;
167 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1]
171 a_offset = 1 + a_dim1 * 1;
203 for (i__ = *n - 1; i__ >= 1; --i__) {
210 e[i__] =
a_ref(i__, i__ + 1);
216 a_ref(i__, i__ + 1) = 1.;
221 1), &c__1, &c_b8, &tau[1], &c__1);
234 &c__1, &a[a_offset], lda);
236 a_ref(i__, i__ + 1) = e[i__];
238 d__[i__ + 1] =
a_ref(i__ + 1, i__ + 1);
242 d__[1] =
a_ref(1, 1);
248 for (i__ = 1; i__ <= i__1; ++i__) {
258 e[i__] =
a_ref(i__ + 1, i__);
264 a_ref(i__ + 1, i__) = 1.;
270 a_ref(i__ + 1, i__), &c__1, &c_b8, &tau[i__], &c__1);
276 i__ + 1, i__), &c__1);
286 i__], &c__1, &
a_ref(i__ + 1, i__ + 1), lda)
289 a_ref(i__ + 1, i__) = e[i__];
291 d__[i__] =
a_ref(i__, i__);
295 d__[*n] =
a_ref(*n, *n);