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);
int template_blas_scal(const integer *n, const Treal *da, Treal *dx, const integer *incx)
Definition: template_blas_scal.h:43
int integer
Definition: template_blas_common.h:40
#define minMACRO(a, b)
Definition: template_blas_common.h:46
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:42
int template_lapack_larfg(const integer *n, Treal *alpha, Treal *x, const integer *incx, Treal *tau)
Definition: template_lapack_larfg.h:43
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:43
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:43
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:42
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:46
Treal template_blas_dot(const integer *n, const Treal *dx, const integer *incx, const Treal *dy, const integer *incy)
Definition: template_blas_dot.h:43