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.;
int template_lapack_laset(const char *uplo, const integer *m, const integer *n, const Treal *alpha, const Treal *beta, Treal *a, const integer *lda)
Definition: template_lapack_laset.h:42
int integer
Definition: template_blas_common.h:40
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
int template_lapack_gghrd(const char *compq, const char *compz, const integer *n, const integer *ilo, const integer *ihi, Treal *a, const integer *lda, Treal *b, const integer *ldb, Treal *q, const integer *ldq, Treal *z__, const integer *ldz, integer *info)
Definition: template_lapack_gghrd.h:42
int template_blas_rot(const integer *n, Treal *dx, const integer *incx, Treal *dy, const integer *incy, const Treal *c__, const Treal *s)
Definition: template_blas_rot.h:42
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
int template_lapack_lartg(const Treal *f, const Treal *g, Treal *cs, Treal *sn, Treal *r__)
Definition: template_lapack_lartg.h:42
#define z___ref(a_1, a_2)
bool logical
Definition: template_blas_common.h:41
#define TRUE_
Definition: template_lapack_common.h:42
#define FALSE_
Definition: template_lapack_common.h:43
logical template_blas_lsame(const char *ca, const char *cb)
Definition: template_blas_common.cc:46