35 #ifndef TEMPLATE_LAPACK_LATRD_HEADER
36 #define TEMPLATE_LAPACK_LATRD_HEADER
41 a,
const integer *lda, Treal *e, Treal *tau, Treal *w,
184 integer a_dim1, a_offset, w_dim1, w_offset, i__1, i__2, i__3;
189 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1]
190 #define w_ref(a_1,a_2) w[(a_2)*w_dim1 + a_1]
194 a_offset = 1 + a_dim1 * 1;
199 w_offset = 1 + w_dim1 * 1;
212 for (i__ = *n; i__ >= i__1; --i__) {
220 lda, &
w_ref(i__, iw + 1), ldw, &c_b6, &
a_ref(1, i__),
224 ldw, &
a_ref(i__, i__ + 1), lda, &c_b6, &
a_ref(1, i__),
235 e[i__ - 1] =
a_ref(i__ - 1, i__);
236 a_ref(i__ - 1, i__) = 1.;
242 i__), &c__1, &c_b16, &
w_ref(1, iw), &c__1);
247 , ldw, &
a_ref(1, i__), &c__1, &c_b16, &
w_ref(i__
252 + 1), lda, &
w_ref(i__ + 1, iw), &c__1, &c_b6, &
253 w_ref(1, iw), &c__1);
257 1), lda, &
a_ref(1, i__), &c__1, &c_b16, &
w_ref(
258 i__ + 1, iw), &c__1);
262 + 1), ldw, &
w_ref(i__ + 1, iw), &c__1, &c_b6, &
263 w_ref(1, iw), &c__1);
269 c__1, &
a_ref(1, i__), &c__1);
282 for (i__ = 1; i__ <= i__1; ++i__) {
289 w_ref(i__, 1), ldw, &c_b6, &
a_ref(i__, i__), &c__1);
293 a_ref(i__, 1), lda, &c_b6, &
a_ref(i__, i__), &c__1);
304 e[i__] =
a_ref(i__ + 1, i__);
305 a_ref(i__ + 1, i__) = 1.;
311 a_ref(i__ + 1, i__), &c__1, &c_b16, &
w_ref(i__ + 1,
316 ldw, &
a_ref(i__ + 1, i__), &c__1, &c_b16, &
w_ref(1,
321 , lda, &
w_ref(1, i__), &c__1, &c_b6, &
w_ref(i__ + 1,
326 lda, &
a_ref(i__ + 1, i__), &c__1, &c_b16, &
w_ref(1,
331 , ldw, &
w_ref(1, i__), &c__1, &c_b6, &
w_ref(i__ + 1,
337 c__1, &
a_ref(i__ + 1, i__), &c__1);