35 #ifndef TEMPLATE_LAPACK_LARFT_HEADER
36 #define TEMPLATE_LAPACK_LARFT_HEADER
41 k, Treal *v,
const integer *ldv,
const Treal *tau, Treal *t,
149 integer t_dim1, t_offset, v_dim1, v_offset, i__1, i__2, i__3;
154 #define t_ref(a_1,a_2) t[(a_2)*t_dim1 + a_1]
155 #define v_ref(a_1,a_2) v[(a_2)*v_dim1 + a_1]
159 v_offset = 1 + v_dim1 * 1;
163 t_offset = 1 + t_dim1 * 1;
173 for (i__ = 1; i__ <= i__1; ++i__) {
174 if (tau[i__] == 0.) {
179 for (j = 1; j <= i__2; ++j) {
187 vii =
v_ref(i__, i__);
188 v_ref(i__, i__) = 1.;
197 ldv, &
v_ref(i__, i__), &c__1, &c_b8, &
t_ref(1,
207 , ldv, &
v_ref(i__, i__), ldv, &c_b8, &
t_ref(1,
210 v_ref(i__, i__) = vii;
216 t_offset], ldt, &
t_ref(1, i__), &c__1);
217 t_ref(i__, i__) = tau[i__];
222 for (i__ = *k; i__ >= 1; --i__) {
223 if (tau[i__] == 0.) {
228 for (j = i__; j <= i__1; ++j) {
238 vii =
v_ref(*n - *k + i__, i__);
239 v_ref(*n - *k + i__, i__) = 1.;
244 i__1 = *n - *k + i__;
248 i__ + 1), ldv, &
v_ref(1, i__), &c__1, &c_b8, &
249 t_ref(i__ + 1, i__), &c__1);
250 v_ref(*n - *k + i__, i__) = vii;
252 vii =
v_ref(i__, *n - *k + i__);
253 v_ref(i__, *n - *k + i__) = 1.;
259 i__2 = *n - *k + i__;
262 i__ + 1, 1), ldv, &
v_ref(i__, 1), ldv, &c_b8,
263 &
t_ref(i__ + 1, i__), &c__1);
264 v_ref(i__, *n - *k + i__) = vii;
271 i__ + 1, i__ + 1), ldt, &
t_ref(i__ + 1, i__), &
274 t_ref(i__, i__) = tau[i__];