37 #ifndef TEMPLATE_LAPACK_LARFT_HEADER
38 #define TEMPLATE_LAPACK_LARFT_HEADER
43 k, Treal *v,
const integer *ldv,
const Treal *tau, Treal *t,
151 integer t_dim1, t_offset, v_dim1, v_offset, i__1, i__2, i__3;
156 #define t_ref(a_1,a_2) t[(a_2)*t_dim1 + a_1]
157 #define v_ref(a_1,a_2) v[(a_2)*v_dim1 + a_1]
161 v_offset = 1 + v_dim1 * 1;
165 t_offset = 1 + t_dim1 * 1;
175 for (i__ = 1; i__ <= i__1; ++i__) {
176 if (tau[i__] == 0.) {
181 for (j = 1; j <= i__2; ++j) {
189 vii =
v_ref(i__, i__);
190 v_ref(i__, i__) = 1.;
199 ldv, &
v_ref(i__, i__), &c__1, &c_b8, &
t_ref(1,
209 , ldv, &
v_ref(i__, i__), ldv, &c_b8, &
t_ref(1,
212 v_ref(i__, i__) = vii;
218 t_offset], ldt, &
t_ref(1, i__), &c__1);
219 t_ref(i__, i__) = tau[i__];
224 for (i__ = *k; i__ >= 1; --i__) {
225 if (tau[i__] == 0.) {
230 for (j = i__; j <= i__1; ++j) {
240 vii =
v_ref(*n - *k + i__, i__);
241 v_ref(*n - *k + i__, i__) = 1.;
246 i__1 = *n - *k + i__;
250 i__ + 1), ldv, &
v_ref(1, i__), &c__1, &c_b8, &
251 t_ref(i__ + 1, i__), &c__1);
252 v_ref(*n - *k + i__, i__) = vii;
254 vii =
v_ref(i__, *n - *k + i__);
255 v_ref(i__, *n - *k + i__) = 1.;
261 i__2 = *n - *k + i__;
264 i__ + 1, 1), ldv, &
v_ref(i__, 1), ldv, &c_b8,
265 &
t_ref(i__ + 1, i__), &c__1);
266 v_ref(i__, *n - *k + i__) = vii;
273 i__ + 1, i__ + 1), ldt, &
t_ref(i__ + 1, i__), &
276 t_ref(i__, i__) = tau[i__];