35 #ifndef TEMPLATE_LAPACK_GGHRD_HEADER
36 #define TEMPLATE_LAPACK_GGHRD_HEADER
161 integer a_dim1, a_offset, b_dim1, b_offset, q_dim1, q_offset, z_dim1,
162 z_offset, i__1, i__2, i__3;
170 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1]
171 #define b_ref(a_1,a_2) b[(a_2)*b_dim1 + a_1]
172 #define q_ref(a_1,a_2) q[(a_2)*q_dim1 + a_1]
173 #define z___ref(a_1,a_2) z__[(a_2)*z_dim1 + a_1]
177 a_offset = 1 + a_dim1 * 1;
180 b_offset = 1 + b_dim1 * 1;
183 q_offset = 1 + q_dim1 * 1;
186 z_offset = 1 + z_dim1 * 1;
225 }
else if (icompz <= 0) {
229 }
else if (*ilo < 1) {
231 }
else if (*ihi > *n || *ihi < *ilo - 1) {
237 }
else if ( ( ilq && *ldq < *n ) || *ldq < 1) {
239 }
else if ( ( ilz && *ldz < *n ) || *ldz < 1) {
266 for (jcol = 1; jcol <= i__1; ++jcol) {
268 for (jrow = jcol + 1; jrow <= i__2; ++jrow) {
269 b_ref(jrow, jcol) = 0.;
278 for (jcol = *ilo; jcol <= i__1; ++jcol) {
281 for (jrow = *ihi; jrow >= i__2; --jrow) {
285 temp =
a_ref(jrow - 1, jcol);
288 a_ref(jrow, jcol) = 0.;
292 i__3 = *n + 2 - jrow;
302 temp =
b_ref(jrow, jrow);
305 b_ref(jrow, jrow - 1) = 0.;