Go to the documentation of this file.
37 #ifndef TEMPLATE_LAPACK_ORG2R_HEADER
38 #define TEMPLATE_LAPACK_ORG2R_HEADER
43 a,
const integer *lda,
const Treal *tau, Treal *work,
integer *info)
105 integer a_dim1, a_offset, i__1, i__2;
109 #define a_ref(a_1,a_2) a[(a_2)*a_dim1 + a_1]
113 a_offset = 1 + a_dim1 * 1;
122 }
else if (*n < 0 || *n > *m) {
124 }
else if (*k < 0 || *k > *n) {
144 for (j = *k + 1; j <= i__1; ++j) {
146 for (l = 1; l <= i__2; ++l) {
154 for (i__ = *k; i__ >= 1; --i__) {
159 a_ref(i__, i__) = 1.;
163 a_ref(i__, i__ + 1), lda, &work[1]);
170 a_ref(i__, i__) = 1. - tau[i__];
175 for (l = 1; l <= i__1; ++l) {
int template_blas_scal(const integer *n, const Treal *da, Treal *dx, const integer *incx)
Definition: template_blas_scal.h:43
int template_lapack_org2r(const integer *m, const integer *n, const integer *k, Treal *a, const integer *lda, const Treal *tau, Treal *work, integer *info)
Definition: template_lapack_org2r.h:42
int template_blas_erbla(const char *srname, integer *info)
Definition: template_blas_common.cc:146
int integer
Definition: template_blas_common.h:40
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
#define maxMACRO(a, b)
Definition: template_blas_common.h:45