1 #include "gmx_lapack.h"
4 F77_FUNC(sorm2r,SORM2R)(const char *side,
17 int a_dim1, a_offset, c_dim1, c_offset, i__1, i__2;
19 int i__, i1, i2, i3, ic, jc, mi, ni, nq;
26 a_offset = 1 + a_dim1;
30 c_offset = 1 + c_dim1;
34 left = (*side=='L' || *side=='l');
35 notran = (*trans=='N' || *trans=='n');
45 if (*m <= 0 || *n <= 0 || *k <= 0) {
49 if ((left && !notran) || (!left && notran)) {
69 for (i__ = i1; i__2 < 0 ? i__ >= i__1 : i__ <= i__1; i__ += i__2) {
81 aii = a[i__ + i__ * a_dim1];
82 a[i__ + i__ * a_dim1] = 1.;
83 F77_FUNC(slarf,SLARF)(side, &mi, &ni, &a[i__ + i__ * a_dim1], &c__1, &tau[i__], &c__[
84 ic + jc * c_dim1], ldc, &work[1]);
85 a[i__ + i__ * a_dim1] = aii;