37 #ifndef TEMPLATE_LAPACK_ORG2L_HEADER 38 #define TEMPLATE_LAPACK_ORG2L_HEADER 43 a,
const integer *lda,
const Treal *tau, Treal *work,
integer *info)
105 integer a_dim1, a_offset, i__1, i__2, i__3;
110 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1] 114 a_offset = 1 + a_dim1 * 1;
123 }
else if (*n < 0 || *n > *m) {
125 }
else if (*k < 0 || *k > *n) {
145 for (j = 1; j <= i__1; ++j) {
147 for (l = 1; l <= i__2; ++l) {
151 a_ref(*m - *n + j, j) = 1.;
156 for (i__ = 1; i__ <= i__1; ++i__) {
161 a_ref(*m - *n + ii, ii) = 1.;
165 a_offset], lda, &work[1]);
166 i__2 = *m - *n + ii - 1;
169 a_ref(*m - *n + ii, ii) = 1. - tau[i__];
174 for (l = *m - *n + ii + 1; l <= i__2; ++l) {
int template_blas_scal(const integer *n, const Treal *da, Treal *dx, const integer *incx)
Definition: template_blas_scal.h:43
int integer
Definition: template_blas_common.h:40
#define maxMACRO(a, b)
Definition: template_blas_common.h:45
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
int template_lapack_larf(const char *side, const integer *m, const integer *n, const Treal *v, const integer *incv, const Treal *tau, Treal *c__, const integer *ldc, Treal *work)
Definition: template_lapack_larf.h:42
int template_lapack_org2l(const integer *m, const integer *n, const integer *k, Treal *a, const integer *lda, const Treal *tau, Treal *work, integer *info)
Definition: template_lapack_org2l.h:42