4 #include "gmx_lapack.h"
7 F77_FUNC(dorg2r,DORG2R)(int *m,
16 int a_dim1, a_offset, i__1, i__2;
23 a_offset = 1 + a_dim1;
35 for (j = *k + 1; j <= i__1; ++j) {
37 for (l = 1; l <= i__2; ++l) {
38 a[l + j * a_dim1] = 0.0;
40 a[j + j * a_dim1] = 1.0;
42 for (i__ = *k; i__ >= 1; --i__) {
44 a[i__ + i__ * a_dim1] = 1.0;
47 F77_FUNC(dlarf,DLARF)("L", &i__1, &i__2, &a[i__ + i__ * a_dim1], &c__1,
48 &tau[i__], &a[i__ + (i__ + 1) * a_dim1], lda, &work[1]);
53 F77_FUNC(dscal,DSCAL)(&i__1, &r__1, &a[i__ + 1 + i__ * a_dim1], &c__1);
55 a[i__ + i__ * a_dim1] = 1.0 - tau[i__];
57 for (l = 1; l <= i__1; ++l) {
58 a[l + i__ * a_dim1] = 0.0;