2 #include "gmx_lapack.h"
5 F77_FUNC(slarnv,SLARNV)(int *idist,
20 for (iv = 1; iv <= i__1; iv += 64) {
21 i__2 = 64, i__3 = *n - iv + 1;
22 il = (i__2<i__3) ? i__2 : i__3;
29 F77_FUNC(slaruv,SLARUV)(&iseed[1], &il2, u);
34 for (i__ = 1; i__ <= i__2; ++i__) {
35 x[iv + i__ - 1] = u[i__ - 1];
37 } else if (*idist == 2) {
40 for (i__ = 1; i__ <= i__2; ++i__) {
41 x[iv + i__ - 1] = u[i__ - 1] * 2. - 1.;
43 } else if (*idist == 3) {
46 for (i__ = 1; i__ <= i__2; ++i__) {
47 x[iv + i__ - 1] = sqrt(log(u[(i__ << 1) - 2]) * -2.) *
48 cos(u[(i__ << 1) - 1] *
49 (float)6.2831853071795864769252867663);