37 #ifndef TEMPLATE_LAPACK_SYGS2_HEADER 38 #define TEMPLATE_LAPACK_SYGS2_HEADER 121 integer a_dim1, a_offset, b_dim1, b_offset, i__1, i__2;
128 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1] 129 #define b_ref(a_1,a_2) b[(a_2)*b_dim1 + a_1] 133 a_offset = 1 + a_dim1 * 1;
136 b_offset = 1 + b_dim1 * 1;
142 if (*itype < 1 || *itype > 3) {
165 for (k = 1; k <= i__1; ++k) {
185 lda, &
b_ref(k, k + 1), ldb,
186 &
a_ref(k + 1, k + 1), lda);
192 k + 1), ldb, &
a_ref(k, k + 1), lda);
201 for (k = 1; k <= i__1; ++k) {
221 &c__1, &
b_ref(k + 1, k),
222 &c__1, &
a_ref(k + 1, k + 1),
230 + 1, k + 1), ldb, &
a_ref(k + 1, k), &c__1);
241 for (k = 1; k <= i__1; ++k) {
249 ldb, &
a_ref(1, k), &c__1);
255 &c__1, &a[a_offset], lda);
262 a_ref(k, k) = akk * (d__1 * d__1);
270 for (k = 1; k <= i__1; ++k) {
278 ldb, &
a_ref(k, 1), lda);
284 ldb, &a[a_offset], lda);
291 a_ref(k, k) = akk * (d__1 * d__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
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
int template_blas_trmv(const char *uplo, const char *trans, const char *diag, const integer *n, const Treal *a, const integer *lda, Treal *x, const integer *incx)
Definition: template_blas_trmv.h:42
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
int template_blas_trsv(const char *uplo, const char *trans, const char *diag, const integer *n, const Treal *a, const integer *lda, Treal *x, const integer *incx)
Definition: template_blas_trsv.h:42
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_lapack_sygs2(const integer *itype, const char *uplo, const integer *n, Treal *a, const integer *lda, Treal *b, const integer *ldb, integer *info)
Definition: template_lapack_sygs2.h:43
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:46