1 #include "gmx_lapack.h"
4 F77_FUNC(sorm2l,SORM2L)(const char *side,
17 int a_dim1, a_offset, c_dim1, c_offset, i__1, i__2;
20 int i__, i1, i2, i3, mi, ni, nq;
26 a_offset = 1 + a_dim1;
30 c_offset = 1 + c_dim1;
36 left = (*side=='L' || *side=='l');
37 notran = (*trans=='N' || *trans=='n');
49 if (*m == 0 || *n == 0 || *k == 0) {
53 if ((left && notran) || (! left && ! notran)) {
71 for (i__ = i1; i__2 < 0 ? i__ >= i__1 : i__ <= i__1; i__ += i__2) {
80 aii = a[nq - *k + i__ + i__ * a_dim1];
81 a[nq - *k + i__ + i__ * a_dim1] = 1.;
82 F77_FUNC(slarf,SLARF)(side, &mi, &ni, &a[i__ * a_dim1 + 1], &c__1, &tau[i__], &c__[
83 c_offset], ldc, &work[1]);
84 a[nq - *k + i__ + i__ * a_dim1] = aii;