35 #ifndef TEMPLATE_LAPACK_SYGST_HEADER
36 #define TEMPLATE_LAPACK_SYGST_HEADER
122 integer a_dim1, a_offset, b_dim1, b_offset, i__1, i__2, i__3;
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) {
170 if (nb <= 1 || nb >= *n) {
186 for (k = 1; i__2 < 0 ? k >= i__1 : k <= i__1; k += i__2) {
196 i__3 = *n - k - kb + 1;
200 i__3 = *n - k - kb + 1;
204 i__3 = *n - k - kb + 1;
206 k, k + kb), lda, &
b_ref(k, k + kb), ldb, &
207 c_b14, &
a_ref(k + kb, k + kb), lda);
208 i__3 = *n - k - kb + 1;
212 i__3 = *n - k - kb + 1;
214 &i__3, &c_b14, &
b_ref(k + kb, k + kb), ldb, &
215 a_ref(k, k + kb), lda);
225 for (k = 1; i__1 < 0 ? k >= i__2 : k <= i__2; k += i__1) {
235 i__3 = *n - k - kb + 1;
237 &kb, &c_b14, &
b_ref(k, k), ldb, &
a_ref(k + kb,
239 i__3 = *n - k - kb + 1;
241 , lda, &
b_ref(k + kb, k), ldb, &c_b14, &
a_ref(
243 i__3 = *n - k - kb + 1;
245 a_ref(k + kb, k), lda, &
b_ref(k + kb, k), ldb,
246 &c_b14, &
a_ref(k + kb, k + kb), lda);
247 i__3 = *n - k - kb + 1;
249 , lda, &
b_ref(k + kb, k), ldb, &c_b14, &
a_ref(
251 i__3 = *n - k - kb + 1;
253 i__3, &kb, &c_b14, &
b_ref(k + kb, k + kb),
254 ldb, &
a_ref(k + kb, k), lda);
266 for (k = 1; i__2 < 0 ? k >= i__1 : k <= i__1; k += i__2) {
275 kb, &c_b14, &b[b_offset], ldb, &
a_ref(1, k), lda);
278 lda, &
b_ref(1, k), ldb, &c_b14, &
a_ref(1, k), lda);
281 1, k), lda, &
b_ref(1, k), ldb, &c_b14, &a[
285 lda, &
b_ref(1, k), ldb, &c_b14, &
a_ref(1, k), lda);
288 &c_b14, &
b_ref(k, k), ldb, &
a_ref(1, k), lda);
299 for (k = 1; i__1 < 0 ? k >= i__2 : k <= i__2; k += i__1) {
308 i__3, &c_b14, &b[b_offset], ldb, &
a_ref(k, 1),
312 lda, &
b_ref(k, 1), ldb, &c_b14, &
a_ref(k, 1), lda);
315 1), lda, &
b_ref(k, 1), ldb, &c_b14, &a[a_offset],
319 lda, &
b_ref(k, 1), ldb, &c_b14, &
a_ref(k, 1), lda);
322 &c_b14, &
b_ref(k, k), ldb, &
a_ref(k, 1), lda);