double ran2(long *idum) {
  long j;
  long k;
  static long idum2 = 123456789;
  static long iy = 0;
  static long iv[32];
  double temp;
  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)
    idum += 2147483399;
  j = iy / (1 + (2147483563 - 1) / 32);
  iy = iv[j] - idum2;
  iv[j] = *idum;
  if (iy < 1)
    iy += (2147483563 - 1);
  if ((temp = (1.0 / 2147483563) * iy) > (1.0 - 1.2e-7))
    return (1.0 - 1.2e-7);
  else
    return temp;
}


x86 assembly (clang, -O0):
<ran2>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmpq   $0x0,(%rax)
jg     120e <ran2+0x10e>
mov    -0x10(%rbp),%rcx
xor    %eax,%eax
sub    (%rcx),%rax
cmp    $0x1,%rax
jge    1139 <ran2+0x39>
mov    -0x10(%rbp),%rax
movq   $0x1,(%rax)
jmp    1149 <ran2+0x49>
mov    -0x10(%rbp),%rax
xor    %ecx,%ecx
sub    (%rax),%rcx
mov    -0x10(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    %rax,0x2ec9(%rip)
movq   $0x27,-0x18(%rbp)
cmpq   $0x0,-0x18(%rbp)
jl     1200 <ran2+0x100>
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    $0xd1a4,%ecx
cqto
idiv   %rcx
mov    %rax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
imul   $0xd1a4,-0x20(%rbp),%rcx
sub    %rcx,%rax
imul   $0x9c4e,%rax,%rcx
imul   $0x2fb3,-0x20(%rbp),%rax
sub    %rax,%rcx
mov    -0x10(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x10(%rbp),%rax
cmpq   $0x0,(%rax)
jge    11c9 <ran2+0xc9>
mov    -0x10(%rbp),%rax
mov    (%rax),%rcx
add    $0x7fffffab,%rcx
mov    %rcx,(%rax)
cmpq   $0x20,-0x18(%rbp)
jge    11ea <ran2+0xea>
mov    -0x10(%rbp),%rax
mov    (%rax),%rdx
mov    -0x18(%rbp),%rcx
lea    0x2e6a(%rip),%rax
mov    %rdx,(%rax,%rcx,8)
jmp    11ef <ran2+0xef>
mov    -0x18(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x18(%rbp)
jmp    115f <ran2+0x5f>
mov    0x2e49(%rip),%rax
mov    %rax,0x2e32(%rip)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    $0xd1a4,%ecx
cqto
idiv   %rcx
mov    %rax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
imul   $0xd1a4,-0x20(%rbp),%rcx
sub    %rcx,%rax
imul   $0x9c4e,%rax,%rcx
imul   $0x2fb3,-0x20(%rbp),%rax
sub    %rax,%rcx
mov    -0x10(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x10(%rbp),%rax
cmpq   $0x0,(%rax)
jge    126d <ran2+0x16d>
mov    -0x10(%rbp),%rax
mov    (%rax),%rcx
add    $0x7fffffab,%rcx
mov    %rcx,(%rax)
mov    0x2dac(%rip),%rax
mov    $0xce26,%ecx
cqto
idiv   %rcx
mov    %rax,-0x20(%rbp)
mov    0x2d97(%rip),%rax
imul   $0xce26,-0x20(%rbp),%rcx
sub    %rcx,%rax
imul   $0x9ef4,%rax,%rax
imul   $0xecf,-0x20(%rbp),%rcx
sub    %rcx,%rax
mov    %rax,0x2d73(%rip)
cmpq   $0x0,0x2d6b(%rip)
jge    12cd <ran2+0x1cd>
movabs $0x3fffff838,%rax
add    -0x10(%rbp),%rax
mov    %rax,-0x10(%rbp)
mov    0x2d6c(%rip),%rax
mov    $0x3fffffe,%ecx
cqto
idiv   %rcx
mov    %rax,-0x18(%rbp)
mov    -0x18(%rbp),%rcx
lea    0x2d63(%rip),%rax
mov    (%rax,%rcx,8),%rax
sub    0x2d28(%rip),%rax
mov    %rax,0x2d41(%rip)
mov    -0x10(%rbp),%rax
mov    (%rax),%rdx
mov    -0x18(%rbp),%rcx
lea    0x2d3f(%rip),%rax
mov    %rdx,(%rax,%rcx,8)
cmpq   $0x1,0x2d23(%rip)
01
jge    1337 <ran2+0x237>
mov    0x2d16(%rip),%rax
add    $0x7fffffaa,%rax
mov    %rax,0x2d09(%rip)
cvtsi2sdq 0x2d00(%rip),%xmm1
movsd  0xcc0(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xca7(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1375 <ran2+0x275>
movsd  0xc95(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    137f <ran2+0x27f>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<ran2>:
mov    (%rdi),%rsi
movabs $0x9c4e3aa71ae25487,%r10
3a 4e 9c
test   %rsi,%rsi
jle    1226 <ran2+0x126>
lea    0x2f23(%rip),%r9
mov    0x2efc(%rip),%r8
mov    (%r9),%r9
mov    %rsi,%rax
imul   %r10
add    %rsi,%rdx
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0xf,%rdx
add    %rax,%rdx
imul   $0xffffffffffff2e5c,%rdx,%rax
add    %rsi,%rax
imul   $0x9c4e,%rax,%rax
imul   $0xffffffffffffd04d,%rdx,%rcx
lea    (%rax,%rcx,1),%rdx
add    $0x7fffffab,%rdx
add    %rcx,%rax
cmovns %rax,%rdx
mov    %rdx,(%rdi)
movabs $0x4f7a09cc7bd54605,%rcx
09 7a 4f
mov    %r8,%rax
imul   %rcx
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0xe,%rdx
add    %rax,%rdx
imul   $0xffffffffffff31da,%rdx,%rax
add    %r8,%rax
imul   $0x9ef4,%rax,%rax
imul   $0xfffffffffffff131,%rdx,%rcx
add    %rax,%rcx
mov    %rcx,0x2e75(%rip)
mov    %rcx,%rax
sar    $0x3f,%rax
movabs $0x3fffff838,%rsi
and    %rax,%rsi
movabs $0x8000004000002001,%rdx
mov    %r9,%rax
imul   %rdx
add    %r9,%rdx
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0x19,%rdx
add    %rax,%rdx
lea    0x2e69(%rip),%r8
mov    (%r8,%rdx,8),%rax
sub    %rcx,%rax
mov    %rax,0x2e4b(%rip)
mov    (%rdi,%rsi,1),%rcx
mov    %rcx,(%r8,%rdx,8)
jg     120c <ran2+0x10c>
add    $0x7fffffaa,%rax
mov    %rax,0x2e34(%rip)
cvtsi2sd %rax,%xmm1
mulsd  0xde7(%rip),%xmm1
movsd  0xde7(%rip),%xmm0
minsd  %xmm1,%xmm0
ret
neg    %rsi
mov    $0x1,%r8d
cmovb  %rsi,%r8
mov    $0x27,%ecx
lea    0x2e11(%rip),%r9
mov    %r8,%rsi
jmp    125a <ran2+0x15a>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0xffffffffffffffff,%rcx
jae    1124 <ran2+0x24>
mov    %rsi,%rax
imul   %r10
add    %rsi,%rdx
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0xf,%rdx
add    %rax,%rdx
imul   $0xffffffffffff2e5c,%rdx,%rax
add    %rsi,%rax
imul   $0x9c4e,%rax,%rax
imul   $0xffffffffffffd04d,%rdx,%rdx
lea    (%rax,%rdx,1),%rsi
add    $0x7fffffab,%rsi
add    %rdx,%rax
cmovns %rax,%rsi
cmp    $0x1f,%rcx
ja     1250 <ran2+0x150>
mov    %rsi,(%r9,%rcx,8)
jmp    1250 <ran2+0x150>

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

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