37 #ifndef TEMPLATE_LAPACK_LATRD_HEADER
38 #define TEMPLATE_LAPACK_LATRD_HEADER
43 a,
const integer *lda, Treal *e, Treal *tau, Treal *w,
186 integer a_dim1, a_offset, w_dim1, w_offset, i__1, i__2, i__3;
191 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1]
192 #define w_ref(a_1,a_2) w[(a_2)*w_dim1 + a_1]
196 a_offset = 1 + a_dim1 * 1;
201 w_offset = 1 + w_dim1 * 1;
214 for (i__ = *n; i__ >= i__1; --i__) {
222 lda, &
w_ref(i__, iw + 1), ldw, &c_b6, &
a_ref(1, i__),
226 ldw, &
a_ref(i__, i__ + 1), lda, &c_b6, &
a_ref(1, i__),
237 e[i__ - 1] =
a_ref(i__ - 1, i__);
238 a_ref(i__ - 1, i__) = 1.;
244 i__), &c__1, &c_b16, &
w_ref(1, iw), &c__1);
249 , ldw, &
a_ref(1, i__), &c__1, &c_b16, &
w_ref(i__
254 + 1), lda, &
w_ref(i__ + 1, iw), &c__1, &c_b6, &
255 w_ref(1, iw), &c__1);
259 1), lda, &
a_ref(1, i__), &c__1, &c_b16, &
w_ref(
260 i__ + 1, iw), &c__1);
264 + 1), ldw, &
w_ref(i__ + 1, iw), &c__1, &c_b6, &
265 w_ref(1, iw), &c__1);
271 c__1, &
a_ref(1, i__), &c__1);
284 for (i__ = 1; i__ <= i__1; ++i__) {
291 w_ref(i__, 1), ldw, &c_b6, &
a_ref(i__, i__), &c__1);
295 a_ref(i__, 1), lda, &c_b6, &
a_ref(i__, i__), &c__1);
306 e[i__] =
a_ref(i__ + 1, i__);
307 a_ref(i__ + 1, i__) = 1.;
313 a_ref(i__ + 1, i__), &c__1, &c_b16, &
w_ref(i__ + 1,
318 ldw, &
a_ref(i__ + 1, i__), &c__1, &c_b16, &
w_ref(1,
323 , lda, &
w_ref(1, i__), &c__1, &c_b6, &
w_ref(i__ + 1,
328 lda, &
a_ref(i__ + 1, i__), &c__1, &c_b16, &
w_ref(1,
333 , ldw, &
w_ref(1, i__), &c__1, &c_b6, &
w_ref(i__ + 1,
339 c__1, &
a_ref(i__ + 1, i__), &c__1);