35 #ifndef TEMPLATE_LAPACK_SYGS2_HEADER
36 #define TEMPLATE_LAPACK_SYGS2_HEADER
119 integer a_dim1, a_offset, b_dim1, b_offset, i__1, i__2;
126 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1]
127 #define b_ref(a_1,a_2) b[(a_2)*b_dim1 + a_1]
131 a_offset = 1 + a_dim1 * 1;
134 b_offset = 1 + b_dim1 * 1;
140 if (*itype < 1 || *itype > 3) {
163 for (k = 1; k <= i__1; ++k) {
183 lda, &
b_ref(k, k + 1), ldb,
184 &
a_ref(k + 1, k + 1), lda);
190 k + 1), ldb, &
a_ref(k, k + 1), lda);
199 for (k = 1; k <= i__1; ++k) {
219 &c__1, &
b_ref(k + 1, k),
220 &c__1, &
a_ref(k + 1, k + 1),
228 + 1, k + 1), ldb, &
a_ref(k + 1, k), &c__1);
239 for (k = 1; k <= i__1; ++k) {
247 ldb, &
a_ref(1, k), &c__1);
253 &c__1, &a[a_offset], lda);
260 a_ref(k, k) = akk * (d__1 * d__1);
268 for (k = 1; k <= i__1; ++k) {
276 ldb, &
a_ref(k, 1), lda);
282 ldb, &a[a_offset], lda);
289 a_ref(k, k) = akk * (d__1 * d__1);