35 #ifndef TEMPLATE_LAPACK_SYTD2_HEADER
36 #define TEMPLATE_LAPACK_SYTD2_HEADER
42 lda, Treal *d__, Treal *e, Treal *tau,
integer *info)
161 integer a_dim1, a_offset, i__1, i__2, i__3;
167 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1]
171 a_offset = 1 + a_dim1 * 1;
203 for (i__ = *n - 1; i__ >= 1; --i__) {
210 e[i__] =
a_ref(i__, i__ + 1);
216 a_ref(i__, i__ + 1) = 1.;
221 1), &c__1, &c_b8, &tau[1], &c__1);
234 &c__1, &a[a_offset], lda);
236 a_ref(i__, i__ + 1) = e[i__];
238 d__[i__ + 1] =
a_ref(i__ + 1, i__ + 1);
242 d__[1] =
a_ref(1, 1);
248 for (i__ = 1; i__ <= i__1; ++i__) {
258 e[i__] =
a_ref(i__ + 1, i__);
264 a_ref(i__ + 1, i__) = 1.;
270 a_ref(i__ + 1, i__), &c__1, &c_b8, &tau[i__], &c__1);
276 i__ + 1, i__), &c__1);
286 i__], &c__1, &
a_ref(i__ + 1, i__ + 1), lda)
289 a_ref(i__ + 1, i__) = e[i__];
291 d__[i__] =
a_ref(i__, i__);
295 d__[*n] =
a_ref(*n, *n);
int template_lapack_sytd2(const char *uplo, const integer *n, Treal *a, const integer *lda, Treal *d__, Treal *e, Treal *tau, integer *info)
Definition: template_lapack_sytd2.h:41
int integer
Definition: template_blas_common.h:38
int template_blas_syr2(const char *uplo, const integer *n, const Treal *alpha, const Treal *x, const integer *incx, const Treal *y, const integer *incy, Treal *a, const integer *lda)
Definition: template_blas_syr2.h:40
#define maxMACRO(a, b)
Definition: template_blas_common.h:43
#define minMACRO(a, b)
Definition: template_blas_common.h:44
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:144
int template_lapack_larfg(const integer *n, Treal *alpha, Treal *x, const integer *incx, Treal *tau)
Definition: template_lapack_larfg.h:41
bool logical
Definition: template_blas_common.h:39
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