37 #ifndef TEMPLATE_LAPACK_GGHRD_HEADER
38 #define TEMPLATE_LAPACK_GGHRD_HEADER
163 integer a_dim1, a_offset, b_dim1, b_offset, q_dim1, q_offset, z_dim1,
164 z_offset, i__1, i__2, i__3;
172 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1]
173 #define b_ref(a_1,a_2) b[(a_2)*b_dim1 + a_1]
174 #define q_ref(a_1,a_2) q[(a_2)*q_dim1 + a_1]
175 #define z___ref(a_1,a_2) z__[(a_2)*z_dim1 + a_1]
179 a_offset = 1 + a_dim1 * 1;
182 b_offset = 1 + b_dim1 * 1;
185 q_offset = 1 + q_dim1 * 1;
188 z_offset = 1 + z_dim1 * 1;
227 }
else if (icompz <= 0) {
231 }
else if (*ilo < 1) {
233 }
else if (*ihi > *n || *ihi < *ilo - 1) {
239 }
else if ( ( ilq && *ldq < *n ) || *ldq < 1) {
241 }
else if ( ( ilz && *ldz < *n ) || *ldz < 1) {
268 for (jcol = 1; jcol <= i__1; ++jcol) {
270 for (jrow = jcol + 1; jrow <= i__2; ++jrow) {
271 b_ref(jrow, jcol) = 0.;
280 for (jcol = *ilo; jcol <= i__1; ++jcol) {
283 for (jrow = *ihi; jrow >= i__2; --jrow) {
287 temp =
a_ref(jrow - 1, jcol);
290 a_ref(jrow, jcol) = 0.;
294 i__3 = *n + 2 - jrow;
304 temp =
b_ref(jrow, jrow);
307 b_ref(jrow, jrow - 1) = 0.;