double ran2(int *idum) {
  static int m = 714205, ia = 1366, ic = 150889;
  static double rm = 1.40015821787e-6;
  static int iff = 0;
  static int ir[97], iy;
  int j;
  double y;
  if (*idum < 0 || iff == 0) {
    iff = 1;
    *idum = (ic - *idum) % m;
    for (j = 0; j < 97; ++j) {
      *idum = (*idum * ia + ic) % m;
      ir[j] = *idum;
    }
    *idum = (*idum * ia + ic) % m;
    iy = *idum;
  }
  j = (iy * 97) / m;
  if (j > 96 || j < 0) {
    return (-1.0);
  }
  y = (double)ir[j];
  *idum = (*idum * ia + ic) % m;
  ir[j] = *idum;
  return ((double)iy * rm);
}


x86 assembly (clang, -O0):
<ran2>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmpl   $0x0,(%rax)
jl     1122 <ran2+0x22>
cmpl   $0x0,0x2f34(%rip)
jne    11c4 <ran2+0xc4>
movl   $0x1,0x2f24(%rip)
mov    0x2ef6(%rip),%eax
mov    -0x10(%rbp),%rcx
sub    (%rcx),%eax
cltd
idivl  0x2ee1(%rip)
mov    -0x10(%rbp),%rax
mov    %edx,(%rax)
movl   $0x0,-0x14(%rbp)
cmpl   $0x61,-0x14(%rbp)
jge    1198 <ran2+0x98>
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
imul   0x2ec1(%rip),%eax
add    0x2ebf(%rip),%eax
cltd
idivl  0x2eb0(%rip)
mov    -0x10(%rbp),%rax
mov    %edx,(%rax)
mov    -0x10(%rbp),%rax
mov    (%rax),%edx
movslq -0x14(%rbp),%rcx
lea    0x2ed9(%rip),%rax
mov    %edx,(%rax,%rcx,4)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    114c <ran2+0x4c>
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
imul   0x2e7f(%rip),%eax
add    0x2e7d(%rip),%eax
cltd
idivl  0x2e6e(%rip)
mov    -0x10(%rbp),%rax
mov    %edx,(%rax)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    %eax,0x3020(%rip)
imul   $0x61,0x3019(%rip),%eax
cltd
idivl  0x2e4e(%rip)
mov    %eax,-0x14(%rbp)
cmpl   $0x60,-0x14(%rbp)
jg     11e9 <ran2+0xe9>
cmpl   $0x0,-0x14(%rbp)
jge    11fb <ran2+0xfb>
movsd  0xe0f(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1259 <ran2+0x159>
movslq -0x14(%rbp),%rcx
lea    0x2e5a(%rip),%rax
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
imul   0x2e07(%rip),%eax
add    0x2e05(%rip),%eax
cltd
idivl  0x2df6(%rip)
mov    -0x10(%rbp),%rax
mov    %edx,(%rax)
mov    -0x10(%rbp),%rax
mov    (%rax),%edx
movslq -0x14(%rbp),%rcx
lea    0x2e1f(%rip),%rax
mov    %edx,(%rax,%rcx,4)
cvtsi2sdl 0x2f98(%rip),%xmm0
mulsd  0x2ddc(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<ran2>:
mov    (%rdi),%eax
test   %eax,%eax
js     111a <ran2+0x1a>
testb  $0x1,0x2f17(%rip)
je     111a <ran2+0x1a>
mov    0x2f13(%rip),%ecx
jmp    11dd <ran2+0xdd>
movb   $0x1,0x2f03(%rip)
mov    $0x24d69,%ecx
sub    %eax,%ecx
mov    $0x60,%eax
nopl   (%rax)
movslq %ecx,%rcx
imul   $0xffffffffbbed11ef,%rcx,%rdx
shr    $0x20,%rdx
add    %ecx,%edx
mov    %edx,%esi
shr    $0x1f,%esi
sar    $0x13,%edx
add    %esi,%edx
imul   $0xae5dd,%edx,%edx
sub    %edx,%ecx
imul   $0x556,%ecx,%edx
mov    %edx,%ecx
add    $0x24d69,%ecx
movslq %ecx,%rcx
imul   $0xffffffffbbed11ef,%rcx,%rsi
shr    $0x20,%rsi
add    %esi,%edx
add    $0x24d69,%edx
mov    %edx,%esi
shr    $0x1f,%esi
sar    $0x13,%edx
add    %esi,%edx
imul   $0xae5dd,%edx,%edx
sub    %edx,%ecx
test   %rax,%rax
je     119f <ran2+0x9f>
imul   $0x556,%ecx,%ecx
add    $0x24d69,%ecx
add    $0xfffffffffffffffe,%rax
jmp    1130 <ran2+0x30>
imul   $0x556,%ecx,%ecx
mov    %ecx,%eax
add    $0x24d69,%eax
cltq
imul   $0xffffffffbbed11ef,%rax,%rdx
shr    $0x20,%rdx
add    %edx,%ecx
add    $0x24d69,%ecx
mov    %ecx,%edx
shr    $0x1f,%edx
sar    $0x13,%ecx
add    %edx,%ecx
imul   $0xae5dd,%ecx,%ecx
sub    %ecx,%eax
mov    %eax,(%rdi)
mov    %eax,0x2e4d(%rip)
mov    %eax,%ecx
imul   $0x61,%ecx,%edx
cmp    $0xfff51a24,%edx
jge    11f1 <ran2+0xf1>
movsd  0xe10(%rip),%xmm0
ret
imul   $0x556,%eax,%eax
mov    %eax,%edx
add    $0x24d69,%edx
movslq %edx,%rdx
imul   $0xffffffffbbed11ef,%rdx,%rsi
shr    $0x20,%rsi
add    %esi,%eax
add    $0x24d69,%eax
mov    %eax,%esi
shr    $0x1f,%esi
sar    $0x13,%eax
add    %esi,%eax
imul   $0xae5dd,%eax,%eax
sub    %eax,%edx
mov    %edx,(%rdi)
cvtsi2sd %ecx,%xmm0
mulsd  0xdd4(%rip),%xmm0
ret

LLVM IR (-O0):
@ran2.m = internal global i32 714205
@ran2.ia = internal global i32 1366
@ran2.ic = internal global i32 150889
@ran2.rm = internal global double 0x3EB77DA23DC42526
@ran2.iff = internal global i32 0
@ran2.ir = internal global [97 x i32] zeroinitializer
@ran2.iy = internal global i32 0
define double @ran2(i32* %0) {
%2 = alloca double
%3 = alloca i32*
%4 = alloca i32
%5 = alloca double
store i32* %0, i32** %3
%6 = load i32*, i32** %3
%7 = load i32, i32* %6
%8 = icmp slt i32 %7, 0
br i1 %8, label %12, label %9
9: ; preds = %1
%10 = load i32, i32* @ran2.iff
%11 = icmp eq i32 %10, 0
br i1 %11, label %12, label %53
12: ; preds = %9, %1
store i32 1, i32* @ran2.iff
%13 = load i32, i32* @ran2.ic
%14 = load i32*, i32** %3
%15 = load i32, i32* %14
%16 = sub nsw i32 %13, %15
%17 = load i32, i32* @ran2.m
%18 = srem i32 %16, %17
%19 = load i32*, i32** %3
store i32 %18, i32* %19
store i32 0, i32* %4
br label %20
20: ; preds = %38, %12
%21 = load i32, i32* %4
%22 = icmp slt i32 %21, 97
br i1 %22, label %23, label %41
23: ; preds = %20
%24 = load i32*, i32** %3
%25 = load i32, i32* %24
%26 = load i32, i32* @ran2.ia
%27 = mul nsw i32 %25, %26
%28 = load i32, i32* @ran2.ic
%29 = add nsw i32 %27, %28
%30 = load i32, i32* @ran2.m
%31 = srem i32 %29, %30
%32 = load i32*, i32** %3
store i32 %31, i32* %32
%33 = load i32*, i32** %3
%34 = load i32, i32* %33
%35 = load i32, i32* %4
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds [97 x i32], [97 x i32]* @ran2.ir, i64 0, i64 %36
store i32 %34, i32* %37
br label %38
38: ; preds = %23
%39 = load i32, i32* %4
%40 = add nsw i32 %39, 1
store i32 %40, i32* %4
br label %20
41: ; preds = %20
%42 = load i32*, i32** %3
%43 = load i32, i32* %42
%44 = load i32, i32* @ran2.ia
%45 = mul nsw i32 %43, %44
%46 = load i32, i32* @ran2.ic
%47 = add nsw i32 %45, %46
%48 = load i32, i32* @ran2.m
%49 = srem i32 %47, %48
%50 = load i32*, i32** %3
store i32 %49, i32* %50
%51 = load i32*, i32** %3
%52 = load i32, i32* %51
store i32 %52, i32* @ran2.iy
br label %53
53: ; preds = %41, %9
%54 = load i32, i32* @ran2.iy
%55 = mul nsw i32 %54, 97
%56 = load i32, i32* @ran2.m
%57 = sdiv i32 %55, %56
store i32 %57, i32* %4
%58 = load i32, i32* %4
%59 = icmp sgt i32 %58, 96
br i1 %59, label %63, label %60
60: ; preds = %53
%61 = load i32, i32* %4
%62 = icmp slt i32 %61, 0
br i1 %62, label %63, label %64
63: ; preds = %60, %53
store double -1.000000e+00, double* %2
br label %88
64: ; preds = %60
%65 = load i32, i32* %4
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds [97 x i32], [97 x i32]* @ran2.ir, i64 0, i64 %66
%68 = load i32, i32* %67
%69 = sitofp i32 %68 to double
store double %69, double* %5
%70 = load i32*, i32** %3
%71 = load i32, i32* %70
%72 = load i32, i32* @ran2.ia
%73 = mul nsw i32 %71, %72
%74 = load i32, i32* @ran2.ic
%75 = add nsw i32 %73, %74
%76 = load i32, i32* @ran2.m
%77 = srem i32 %75, %76
%78 = load i32*, i32** %3
store i32 %77, i32* %78
%79 = load i32*, i32** %3
%80 = load i32, i32* %79
%81 = load i32, i32* %4
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds [97 x i32], [97 x i32]* @ran2.ir, i64 0, i64 %82
store i32 %80, i32* %83
%84 = load i32, i32* @ran2.iy
%85 = sitofp i32 %84 to double
%86 = load double, double* @ran2.rm
%87 = fmul double %85, %86
store double %87, double* %2
br label %88
88: ; preds = %64, %63
%89 = load double, double* %2
ret double %89
}

LLVM IR (-O2):
@ran2.iff = internal unnamed_addr global i1 false
@ran2.iy = internal unnamed_addr global i32 0
define double @ran2(i32* nocapture %0) {
%2 = load i32, i32* %0
%3 = icmp slt i32 %2, 0
%4 = load i1, i1* @ran2.iff
%5 = xor i1 %4, true
%6 = select i1 %3, i1 true, i1 %5
br i1 %6, label %9, label %7
7: ; preds = %1
%8 = load i32, i32* @ran2.iy
br label %27
9: ; preds = %1
store i1 true, i1* @ran2.iff
%10 = sub nsw i32 150889, %2
br label %11
11: ; preds = %19, %9
%12 = phi i64 [ 0, %9 ], [ %22, %19 ]
%13 = phi i32 [ %10, %9 ], [ %21, %19 ]
%14 = srem i32 %13, 714205
%15 = mul nsw i32 %14, 1366
%16 = add nsw i32 %15, 150889
%17 = srem i32 %16, 714205
%18 = icmp eq i64 %12, 96
br i1 %18, label %23, label %19
19: ; preds = %11
%20 = mul nsw i32 %17, 1366
%21 = add nsw i32 %20, 150889
%22 = add nuw nsw i64 %12, 2
br label %11
23: ; preds = %11
%24 = mul nsw i32 %17, 1366
%25 = add nsw i32 %24, 150889
%26 = srem i32 %25, 714205
store i32 %26, i32* %0
store i32 %26, i32* @ran2.iy
br label %27
27: ; preds = %7, %23
%28 = phi i32 [ %2, %7 ], [ %26, %23 ]
%29 = phi i32 [ %8, %7 ], [ %26, %23 ]
%30 = mul nsw i32 %29, 97
%31 = icmp slt i32 %30, -714204
br i1 %31, label %38, label %32
32: ; preds = %27
%33 = mul nsw i32 %28, 1366
%34 = add nsw i32 %33, 150889
%35 = srem i32 %34, 714205
store i32 %35, i32* %0
%36 = sitofp i32 %29 to double
%37 = fmul double %36, 0x3EB77DA23DC42526
br label %38
38: ; preds = %27, %32
%39 = phi double [ %37, %32 ], [ -1.000000e+00, %27 ]
ret double %39
}