double ran2(int *idum) {
  int j;
  int k;
  static int idum2 = 123456789;
  static int iy = 0;
  static int iv[32];
  double tempran;
  if (*idum <= 0) {
    if (-(*idum) < 1) {
      *idum = 1;
    } else {
      *idum = -(*idum);
    }
    idum2 = (*idum);
    for (j = 32 + 7; j >= 0; j--) {
      k = (*idum) / 53668;
      *idum = 40014 * (*idum - k * 53668) - k * 12211;
      if (*idum < 0) {
        *idum += 2147483563;
      }
      if (j < 32) {
        iv[j] = *idum;
      }
    }
    iy = iv[0];
  }
  k = (*idum) / 53668;
  *idum = 40014 * (*idum - k * 53668) - k * 12211;
  if (*idum < 0) {
    *idum += 2147483563;
  }
  k = idum2 / 52774;
  idum2 = 40692 * (idum2 - k * 52774) - k * 3791;
  if (idum2 < 0) {
    idum2 += 2147483399;
  }
  j = iy / (1 + (2147483563 - 1) / 32);
  iy = iv[j] - idum2;
  iv[j] = *idum;
  if (iy < 1) {
    iy += (2147483563 - 1);
  }
  if ((tempran = (1.0 / 2147483563) * iy) > (1.0 - 1.2e-7)) {
    return (1.0 - 1.2e-7);
  } else
    return tempran;
}


x86 assembly (clang, -O0):
<ran2>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmpl   $0x0,(%rax)
jg     11ed <ran2+0xed>
mov    -0x10(%rbp),%rcx
xor    %eax,%eax
sub    (%rcx),%eax
cmp    $0x1,%eax
jge    1135 <ran2+0x35>
mov    -0x10(%rbp),%rax
movl   $0x1,(%rax)
jmp    1143 <ran2+0x43>
mov    -0x10(%rbp),%rax
xor    %ecx,%ecx
sub    (%rax),%ecx
mov    -0x10(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    %eax,0x2ed1(%rip)
movl   $0x27,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
jl     11e1 <ran2+0xe1>
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    $0xd1a4,%ecx
cltd
idiv   %ecx
mov    %eax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
imul   $0xd1a4,-0x18(%rbp),%ecx
sub    %ecx,%eax
imul   $0x9c4e,%eax,%ecx
imul   $0x2fb3,-0x18(%rbp),%eax
sub    %eax,%ecx
mov    -0x10(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x10(%rbp),%rax
cmpl   $0x0,(%rax)
jge    11b0 <ran2+0xb0>
mov    -0x10(%rbp),%rax
mov    (%rax),%ecx
add    $0x7fffffab,%ecx
mov    %ecx,(%rax)
cmpl   $0x20,-0x14(%rbp)
jge    11ce <ran2+0xce>
mov    -0x10(%rbp),%rax
mov    (%rax),%edx
movslq -0x14(%rbp),%rcx
lea    0x2e85(%rip),%rax
mov    %edx,(%rax,%rcx,4)
jmp    11d3 <ran2+0xd3>
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    1156 <ran2+0x56>
mov    0x2e69(%rip),%eax
mov    %eax,0x2e53(%rip)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    $0xd1a4,%ecx
cltd
idiv   %ecx
mov    %eax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
imul   $0xd1a4,-0x18(%rbp),%ecx
sub    %ecx,%eax
imul   $0x9c4e,%eax,%ecx
imul   $0x2fb3,-0x18(%rbp),%eax
sub    %eax,%ecx
mov    -0x10(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x10(%rbp),%rax
cmpl   $0x0,(%rax)
jge    123d <ran2+0x13d>
mov    -0x10(%rbp),%rax
mov    (%rax),%ecx
add    $0x7fffffab,%ecx
mov    %ecx,(%rax)
mov    0x2ddd(%rip),%eax
mov    $0xce26,%ecx
cltd
idiv   %ecx
mov    %eax,-0x18(%rbp)
mov    0x2dcc(%rip),%eax
imul   $0xce26,-0x18(%rbp),%ecx
sub    %ecx,%eax
imul   $0x9ef4,%eax,%eax
imul   $0xecf,-0x18(%rbp),%ecx
sub    %ecx,%eax
mov    %eax,0x2dae(%rip)
cmpl   $0x0,0x2da7(%rip)
jge    1290 <ran2+0x190>
mov    0x2d9b(%rip),%eax
add    $0x7fffff07,%eax
mov    %eax,0x2d90(%rip)
mov    0x2daa(%rip),%eax
mov    $0x3fffffe,%ecx
cltd
idiv   %ecx
mov    %eax,-0x14(%rbp)
movslq -0x14(%rbp),%rcx
lea    0x2da4(%rip),%rax
mov    (%rax,%rcx,4),%eax
sub    0x2d6b(%rip),%eax
mov    %eax,0x2d85(%rip)
mov    -0x10(%rbp),%rax
mov    (%rax),%edx
movslq -0x14(%rbp),%rcx
lea    0x2d84(%rip),%rax
mov    %edx,(%rax,%rcx,4)
cmpl   $0x1,0x2d6a(%rip)
jge    12ed <ran2+0x1ed>
mov    0x2d5e(%rip),%eax
add    $0x7fffffaa,%eax
mov    %eax,0x2d53(%rip)
cvtsi2sdl 0x2d4b(%rip),%xmm1
movsd  0xd0b(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  0xcf2(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    132a <ran2+0x22a>
movsd  0xce0(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1334 <ran2+0x234>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<ran2>:
mov    (%rdi),%edx
test   %edx,%edx
jle    11f0 <ran2+0xf0>
lea    0x2f2f(%rip),%r9
mov    0x2f08(%rip),%r8d
movslq (%r9),%rcx
movslq %edx,%rax
imul   $0x1389c755,%rax,%rdx
mov    %rdx,%rsi
shr    $0x3f,%rsi
shr    $0x20,%rdx
sar    $0xc,%edx
add    %esi,%edx
imul   $0xffff2e5c,%edx,%esi
add    %esi,%eax
imul   $0x9c4e,%eax,%eax
imul   $0xffffd04d,%edx,%esi
lea    (%rax,%rsi,1),%edx
add    $0x7fffffab,%edx
add    %esi,%eax
cmovns %eax,%edx
mov    %edx,(%rdi)
movslq %r8d,%rax
imul   $0x4f7a09cd,%rax,%rsi
mov    %rsi,%rdi
shr    $0x3f,%rdi
shr    $0x20,%rsi
sar    $0xe,%esi
add    %edi,%esi
imul   $0xffff31da,%esi,%edi
add    %edi,%eax
imul   $0x9ef4,%eax,%eax
imul   $0xfffff131,%esi,%esi
lea    (%rax,%rsi,1),%edi
add    $0x7fffff07,%edi
add    %esi,%eax
cmovns %eax,%edi
mov    %edi,0x2e85(%rip)
imul   $0xffffffff80000041,%rcx,%rax
shr    $0x20,%rax
add    %ecx,%eax
mov    %eax,%ecx
shr    $0x1f,%ecx
sar    $0x19,%eax
add    %ecx,%eax
movslq %eax,%rcx
lea    0x2e94(%rip),%rsi
mov    (%rsi,%rcx,4),%eax
mov    %edx,(%rsi,%rcx,4)
sub    %edi,%eax
mov    %eax,0x2e76(%rip)
jg     11d7 <ran2+0xd7>
add    $0x7fffffaa,%eax
mov    %eax,0x2e69(%rip)
cvtsi2sd %eax,%xmm1
mulsd  0xe1d(%rip),%xmm1
movsd  0xe1d(%rip),%xmm0
minsd  %xmm1,%xmm0
ret
neg    %edx
mov    $0x1,%r8d
cmovb  %edx,%r8d
mov    $0x27,%esi
lea    0x2e48(%rip),%r9
mov    %r8d,%edx
jmp    121a <ran2+0x11a>
nopl   (%rax)
add    $0xffffffffffffffff,%rsi
jae    1118 <ran2+0x18>
movslq %edx,%rdx
imul   $0x1389c755,%rdx,%rax
mov    %rax,%rcx
shr    $0x3f,%rcx
sar    $0x2c,%rax
add    %ecx,%eax
imul   $0xffff2e5c,%eax,%ecx
add    %ecx,%edx
imul   $0x9c4e,%edx,%ecx
imul   $0xffffd04d,%eax,%eax
lea    (%rcx,%rax,1),%edx
add    $0x7fffffab,%edx
add    %eax,%ecx
cmovns %ecx,%edx
cmp    $0x1f,%rsi
ja     1210 <ran2+0x110>
mov    %edx,(%r9,%rsi,4)
jmp    1210 <ran2+0x110>

LLVM IR (-O0):
@ran2.idum2 = internal global i32 123456789
@ran2.iy = internal global i32 0
@ran2.iv = internal global [32 x i32] zeroinitializer
define double @ran2(i32* %0) {
%2 = alloca double
%3 = alloca i32*
%4 = alloca i32
%5 = alloca i32
%6 = alloca double
store i32* %0, i32** %3
%7 = load i32*, i32** %3
%8 = load i32, i32* %7
%9 = icmp sle i32 %8, 0
br i1 %9, label %10, label %64
10: ; preds = %1
%11 = load i32*, i32** %3
%12 = load i32, i32* %11
%13 = sub nsw i32 0, %12
%14 = icmp slt i32 %13, 1
br i1 %14, label %15, label %17
15: ; preds = %10
%16 = load i32*, i32** %3
store i32 1, i32* %16
br label %22
17: ; preds = %10
%18 = load i32*, i32** %3
%19 = load i32, i32* %18
%20 = sub nsw i32 0, %19
%21 = load i32*, i32** %3
store i32 %20, i32* %21
br label %22
22: ; preds = %17, %15
%23 = load i32*, i32** %3
%24 = load i32, i32* %23
store i32 %24, i32* @ran2.idum2
store i32 39, i32* %4
br label %25
25: ; preds = %59, %22
%26 = load i32, i32* %4
%27 = icmp sge i32 %26, 0
br i1 %27, label %28, label %62
28: ; preds = %25
%29 = load i32*, i32** %3
%30 = load i32, i32* %29
%31 = sdiv i32 %30, 53668
store i32 %31, i32* %5
%32 = load i32*, i32** %3
%33 = load i32, i32* %32
%34 = load i32, i32* %5
%35 = mul nsw i32 %34, 53668
%36 = sub nsw i32 %33, %35
%37 = mul nsw i32 40014, %36
%38 = load i32, i32* %5
%39 = mul nsw i32 %38, 12211
%40 = sub nsw i32 %37, %39
%41 = load i32*, i32** %3
store i32 %40, i32* %41
%42 = load i32*, i32** %3
%43 = load i32, i32* %42
%44 = icmp slt i32 %43, 0
br i1 %44, label %45, label %49
45: ; preds = %28
%46 = load i32*, i32** %3
%47 = load i32, i32* %46
%48 = add nsw i32 %47, 2147483563
store i32 %48, i32* %46
br label %49
49: ; preds = %45, %28
%50 = load i32, i32* %4
%51 = icmp slt i32 %50, 32
br i1 %51, label %52, label %58
52: ; preds = %49
%53 = load i32*, i32** %3
%54 = load i32, i32* %53
%55 = load i32, i32* %4
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds [32 x i32], [32 x i32]* @ran2.iv, i64 0, i64 %56
store i32 %54, i32* %57
br label %58
58: ; preds = %52, %49
br label %59
59: ; preds = %58
%60 = load i32, i32* %4
%61 = add nsw i32 %60, -1
store i32 %61, i32* %4
br label %25
62: ; preds = %25
%63 = load i32, i32* getelementptr inbounds ([32 x i32], [32 x i32]* @ran2.iv, i64 0, i64 0)
store i32 %63, i32* @ran2.iy
br label %64
64: ; preds = %62, %1
%65 = load i32*, i32** %3
%66 = load i32, i32* %65
%67 = sdiv i32 %66, 53668
store i32 %67, i32* %5
%68 = load i32*, i32** %3
%69 = load i32, i32* %68
%70 = load i32, i32* %5
%71 = mul nsw i32 %70, 53668
%72 = sub nsw i32 %69, %71
%73 = mul nsw i32 40014, %72
%74 = load i32, i32* %5
%75 = mul nsw i32 %74, 12211
%76 = sub nsw i32 %73, %75
%77 = load i32*, i32** %3
store i32 %76, i32* %77
%78 = load i32*, i32** %3
%79 = load i32, i32* %78
%80 = icmp slt i32 %79, 0
br i1 %80, label %81, label %85
81: ; preds = %64
%82 = load i32*, i32** %3
%83 = load i32, i32* %82
%84 = add nsw i32 %83, 2147483563
store i32 %84, i32* %82
br label %85
85: ; preds = %81, %64
%86 = load i32, i32* @ran2.idum2
%87 = sdiv i32 %86, 52774
store i32 %87, i32* %5
%88 = load i32, i32* @ran2.idum2
%89 = load i32, i32* %5
%90 = mul nsw i32 %89, 52774
%91 = sub nsw i32 %88, %90
%92 = mul nsw i32 40692, %91
%93 = load i32, i32* %5
%94 = mul nsw i32 %93, 3791
%95 = sub nsw i32 %92, %94
store i32 %95, i32* @ran2.idum2
%96 = load i32, i32* @ran2.idum2
%97 = icmp slt i32 %96, 0
br i1 %97, label %98, label %101
98: ; preds = %85
%99 = load i32, i32* @ran2.idum2
%100 = add nsw i32 %99, 2147483399
store i32 %100, i32* @ran2.idum2
br label %101
101: ; preds = %98, %85
%102 = load i32, i32* @ran2.iy
%103 = sdiv i32 %102, 67108862
store i32 %103, i32* %4
%104 = load i32, i32* %4
%105 = sext i32 %104 to i64
%106 = getelementptr inbounds [32 x i32], [32 x i32]* @ran2.iv, i64 0, i64 %105
%107 = load i32, i32* %106
%108 = load i32, i32* @ran2.idum2
%109 = sub nsw i32 %107, %108
store i32 %109, i32* @ran2.iy
%110 = load i32*, i32** %3
%111 = load i32, i32* %110
%112 = load i32, i32* %4
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds [32 x i32], [32 x i32]* @ran2.iv, i64 0, i64 %113
store i32 %111, i32* %114
%115 = load i32, i32* @ran2.iy
%116 = icmp slt i32 %115, 1
br i1 %116, label %117, label %120
117: ; preds = %101
%118 = load i32, i32* @ran2.iy
%119 = add nsw i32 %118, 2147483562
store i32 %119, i32* @ran2.iy
br label %120
120: ; preds = %117, %101
%121 = load i32, i32* @ran2.iy
%122 = sitofp i32 %121 to double
%123 = fmul double 0x3E0000000AA00007, %122
store double %123, double* %6
%124 = fcmp ogt double %123, 0x3FEFFFFFBF935359
br i1 %124, label %125, label %126
125: ; preds = %120
store double 0x3FEFFFFFBF935359, double* %2
br label %128
126: ; preds = %120
%127 = load double, double* %6
store double %127, double* %2
br label %128
128: ; preds = %126, %125
%129 = load double, double* %2
ret double %129
}

LLVM IR (-O2):
@ran2.idum2 = internal unnamed_addr global i32 123456789
@ran2.iy = internal unnamed_addr global i32 0
@ran2.iv = internal unnamed_addr global [32 x i32] zeroinitializer
define double @ran2(i32* nocapture %0) {
%2 = load i32, i32* %0
%3 = icmp slt i32 %2, 1
br i1 %3, label %6, label %4
4: ; preds = %1
%5 = load i32, i32* @ran2.idum2
br label %28
6: ; preds = %1
%7 = icmp eq i32 %2, 0
%8 = sub nsw i32 0, %2
%9 = select i1 %7, i32 1, i32 %8
br label %10
10: ; preds = %6, %25
%11 = phi i64 [ 39, %6 ], [ %26, %25 ]
%12 = phi i32 [ %9, %6 ], [ %21, %25 ]
%13 = sdiv i32 %12, 53668
%14 = mul nsw i32 %13, -53668
%15 = add i32 %14, %12
%16 = mul nsw i32 %15, 40014
%17 = mul nsw i32 %13, -12211
%18 = add i32 %16, %17
%19 = icmp slt i32 %18, 0
%20 = add nsw i32 %18, 2147483563
%21 = select i1 %19, i32 %20, i32 %18
%22 = icmp ult i64 %11, 32
br i1 %22, label %23, label %25
23: ; preds = %10
%24 = getelementptr inbounds [32 x i32], [32 x i32]* @ran2.iv, i64 0, i64 %11
store i32 %21, i32* %24
br label %25
25: ; preds = %10, %23
%26 = add nsw i64 %11, -1
%27 = icmp eq i64 %11, 0
br i1 %27, label %28, label %10
28: ; preds = %25, %4
%29 = phi i32* [ @ran2.iy, %4 ], [ getelementptr inbounds ([32 x i32], [32 x i32]* @ran2.iv, i64 0, i64 0), %25 ]
%30 = phi i32 [ %5, %4 ], [ %9, %25 ]
%31 = phi i32 [ %2, %4 ], [ %21, %25 ]
%32 = load i32, i32* %29
%33 = sdiv i32 %31, 53668
%34 = mul nsw i32 %33, -53668
%35 = add nsw i32 %34, %31
%36 = mul nsw i32 %35, 40014
%37 = mul nsw i32 %33, -12211
%38 = add i32 %36, %37
%39 = icmp slt i32 %38, 0
%40 = add nsw i32 %38, 2147483563
%41 = select i1 %39, i32 %40, i32 %38
store i32 %41, i32* %0
%42 = sdiv i32 %30, 52774
%43 = mul nsw i32 %42, -52774
%44 = add i32 %43, %30
%45 = mul nsw i32 %44, 40692
%46 = mul nsw i32 %42, -3791
%47 = add i32 %45, %46
%48 = icmp slt i32 %47, 0
%49 = add nsw i32 %47, 2147483399
%50 = select i1 %48, i32 %49, i32 %47
store i32 %50, i32* @ran2.idum2
%51 = sdiv i32 %32, 67108862
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds [32 x i32], [32 x i32]* @ran2.iv, i64 0, i64 %52
%54 = load i32, i32* %53
%55 = sub nsw i32 %54, %50
store i32 %55, i32* @ran2.iy
store i32 %41, i32* %53
%56 = icmp slt i32 %55, 1
br i1 %56, label %57, label %59
57: ; preds = %28
%58 = add nsw i32 %55, 2147483562
store i32 %58, i32* @ran2.iy
br label %59
59: ; preds = %57, %28
%60 = phi i32 [ %58, %57 ], [ %55, %28 ]
%61 = sitofp i32 %60 to double
%62 = fmul double %61, 0x3E0000000AA00007
%63 = fcmp ogt double %62, 0x3FEFFFFFBF935359
%64 = select i1 %63, double 0x3FEFFFFFBF935359, double %62
ret double %64
}