37 #ifndef TEMPLATE_LAPACK_SYTD2_HEADER 38 #define TEMPLATE_LAPACK_SYTD2_HEADER 44 lda, Treal *d__, Treal *e, Treal *tau,
integer *info)
163 integer a_dim1, a_offset, i__1, i__2, i__3;
169 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1] 173 a_offset = 1 + a_dim1 * 1;
205 for (i__ = *n - 1; i__ >= 1; --i__) {
212 e[i__] =
a_ref(i__, i__ + 1);
218 a_ref(i__, i__ + 1) = 1.;
223 1), &c__1, &c_b8, &tau[1], &c__1);
236 &c__1, &a[a_offset], lda);
238 a_ref(i__, i__ + 1) = e[i__];
240 d__[i__ + 1] =
a_ref(i__ + 1, i__ + 1);
244 d__[1] =
a_ref(1, 1);
250 for (i__ = 1; i__ <= i__1; ++i__) {
260 e[i__] =
a_ref(i__ + 1, i__);
266 a_ref(i__ + 1, i__) = 1.;
272 a_ref(i__ + 1, i__), &c__1, &c_b8, &tau[i__], &c__1);
278 i__ + 1, i__), &c__1);
288 i__], &c__1, &
a_ref(i__ + 1, i__ + 1), lda)
291 a_ref(i__ + 1, i__) = e[i__];
293 d__[i__] =
a_ref(i__, i__);
297 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:43
int integer
Definition: template_blas_common.h:40
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:42
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
#define minMACRO(a, b)
Definition: template_blas_common.h:46
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
int template_lapack_larfg(const integer *n, Treal *alpha, Treal *x, const integer *incx, Treal *tau)
Definition: template_lapack_larfg.h:43
bool logical
Definition: template_blas_common.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: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