37 #ifndef TEMPLATE_LAPACK_SYGST_HEADER
38 #define TEMPLATE_LAPACK_SYGST_HEADER
124 integer a_dim1, a_offset, b_dim1, b_offset, i__1, i__2, i__3;
130 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1]
131 #define b_ref(a_1,a_2) b[(a_2)*b_dim1 + a_1]
135 a_offset = 1 + a_dim1 * 1;
138 b_offset = 1 + b_dim1 * 1;
144 if (*itype < 1 || *itype > 3) {
172 if (nb <= 1 || nb >= *n) {
188 for (k = 1; i__2 < 0 ? k >= i__1 : k <= i__1; k += i__2) {
198 i__3 = *n - k - kb + 1;
202 i__3 = *n - k - kb + 1;
206 i__3 = *n - k - kb + 1;
208 k, k + kb), lda, &
b_ref(k, k + kb), ldb, &
209 c_b14, &
a_ref(k + kb, k + kb), lda);
210 i__3 = *n - k - kb + 1;
214 i__3 = *n - k - kb + 1;
216 &i__3, &c_b14, &
b_ref(k + kb, k + kb), ldb, &
217 a_ref(k, k + kb), lda);
227 for (k = 1; i__1 < 0 ? k >= i__2 : k <= i__2; k += i__1) {
237 i__3 = *n - k - kb + 1;
239 &kb, &c_b14, &
b_ref(k, k), ldb, &
a_ref(k + kb,
241 i__3 = *n - k - kb + 1;
243 , lda, &
b_ref(k + kb, k), ldb, &c_b14, &
a_ref(
245 i__3 = *n - k - kb + 1;
247 a_ref(k + kb, k), lda, &
b_ref(k + kb, k), ldb,
248 &c_b14, &
a_ref(k + kb, k + kb), lda);
249 i__3 = *n - k - kb + 1;
251 , lda, &
b_ref(k + kb, k), ldb, &c_b14, &
a_ref(
253 i__3 = *n - k - kb + 1;
255 i__3, &kb, &c_b14, &
b_ref(k + kb, k + kb),
256 ldb, &
a_ref(k + kb, k), lda);
268 for (k = 1; i__2 < 0 ? k >= i__1 : k <= i__1; k += i__2) {
277 kb, &c_b14, &b[b_offset], ldb, &
a_ref(1, k), lda);
280 lda, &
b_ref(1, k), ldb, &c_b14, &
a_ref(1, k), lda);
283 1, k), lda, &
b_ref(1, k), ldb, &c_b14, &a[
287 lda, &
b_ref(1, k), ldb, &c_b14, &
a_ref(1, k), lda);
290 &c_b14, &
b_ref(k, k), ldb, &
a_ref(1, k), lda);
301 for (k = 1; i__1 < 0 ? k >= i__2 : k <= i__2; k += i__1) {
310 i__3, &c_b14, &b[b_offset], ldb, &
a_ref(k, 1),
314 lda, &
b_ref(k, 1), ldb, &c_b14, &
a_ref(k, 1), lda);
317 1), lda, &
b_ref(k, 1), ldb, &c_b14, &a[a_offset],
321 lda, &
b_ref(k, 1), ldb, &c_b14, &
a_ref(k, 1), lda);
324 &c_b14, &
b_ref(k, k), ldb, &
a_ref(k, 1), lda);