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);
int template_blas_scal(const integer *n, const Treal *da, Treal *dx, const integer *incx)
Definition: template_blas_scal.h:41
int integer
Definition: template_blas_common.h:38
#define minMACRO(a, b)
Definition: template_blas_common.h:44
int template_lapack_latrd(const char *uplo, const integer *n, const integer *nb, Treal *a, const integer *lda, Treal *e, Treal *tau, Treal *w, const integer *ldw)
Definition: template_lapack_latrd.h:40
int template_lapack_larfg(const integer *n, Treal *alpha, Treal *x, const integer *incx, Treal *tau)
Definition: template_lapack_larfg.h:41
int template_blas_gemv(const char *trans, const integer *m, const integer *n, const Treal *alpha, const Treal *a, const integer *lda, const Treal *x, const integer *incx, const Treal *beta, Treal *y, const integer *incy)
Definition: template_blas_gemv.h:41
int template_blas_axpy(const integer *n, const Treal *da, const Treal *dx, const integer *incx, Treal *dy, const integer *incy)
Definition: template_blas_axpy.h:41
int template_blas_symv(const char *uplo, const integer *n, const Treal *alpha, const Treal *a, const integer *lda, const Treal *x, const integer *incx, const Treal *beta, Treal *y, const integer *incy)
Definition: template_blas_symv.h:40
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:44
Treal template_blas_dot(const integer *n, const Treal *dx, const integer *incx, const Treal *dy, const integer *incy)
Definition: template_blas_dot.h:41