double ran3(long *idum) {
  int j;
  long k;
  static long iy = 0;
  static long iv[32];
  float temp;
  if (*idum <= 0 || !iy) {
    if (-(*idum) < 1)
      *idum = 1;
    else
      *idum = -(*idum);
    for (j = 32 + 7; j >= 0; j--) {
      k = (*idum) / 127773;
      *idum = 16807 * (*idum - k * 127773) - 2836 * k;
      if (*idum < 0)
        *idum += 2147483647;
      if (j < 32)
        iv[j] = *idum;
    }
    iy = iv[0];
  }
  k = (*idum) / 127773;
  *idum = 16807 * (*idum - k * 127773) - 2836 * k;
  if (*idum < 0)
    *idum += 2147483647;
  j = iy / (1 + (2147483563 - 1) / 32);
  iy = iv[j];
  iv[j] = *idum;
  if ((temp = (1.0 / 2147483647) * iy) > (1.0 - 1.2e-7))
    return (1.0 - 1.2e-7);
  else
    return temp;
}


x86 assembly (clang, -O0):
<ran3>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmpq   $0x0,(%rax)
jle    1124 <ran3+0x24>
cmpq   $0x0,0x2f12(%rip)
jne    1208 <ran3+0x108>
mov    -0x10(%rbp),%rcx
xor    %eax,%eax
sub    (%rcx),%rax
cmp    $0x1,%rax
jge    1147 <ran3+0x47>
mov    -0x10(%rbp),%rax
movq   $0x1,(%rax)
jmp    1157 <ran3+0x57>
mov    -0x10(%rbp),%rax
xor    %ecx,%ecx
sub    (%rax),%rcx
mov    -0x10(%rbp),%rax
mov    %rcx,(%rax)
movl   $0x27,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
jl     11fa <ran3+0xfa>
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    $0x1f31d,%ecx
cqto
idiv   %rcx
mov    %rax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
imul   $0x1f31d,-0x20(%rbp),%rcx
sub    %rcx,%rax
imul   $0x41a7,%rax,%rcx
imul   $0xb14,-0x20(%rbp),%rax
sub    %rax,%rcx
mov    -0x10(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x10(%rbp),%rax
cmpq   $0x0,(%rax)
jge    11c7 <ran3+0xc7>
mov    -0x10(%rbp),%rax
mov    (%rax),%rcx
add    $0x7fffffff,%rcx
mov    %rcx,(%rax)
cmpl   $0x20,-0x14(%rbp)
jge    11e7 <ran3+0xe7>
mov    -0x10(%rbp),%rax
mov    (%rax),%rdx
movslq -0x14(%rbp),%rcx
lea    0x2e5d(%rip),%rax
mov    %rdx,(%rax,%rcx,8)
jmp    11ec <ran3+0xec>
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    115e <ran3+0x5e>
mov    0x2e3f(%rip),%rax
mov    %rax,0x2e28(%rip)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    $0x1f31d,%ecx
cqto
idiv   %rcx
mov    %rax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
imul   $0x1f31d,-0x20(%rbp),%rcx
sub    %rcx,%rax
imul   $0x41a7,%rax,%rcx
imul   $0xb14,-0x20(%rbp),%rax
sub    %rax,%rcx
mov    -0x10(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x10(%rbp),%rax
cmpq   $0x0,(%rax)
jge    1267 <ran3+0x167>
mov    -0x10(%rbp),%rax
mov    (%rax),%rcx
add    $0x7fffffff,%rcx
mov    %rcx,(%rax)
mov    0x2dc2(%rip),%rax
mov    $0x3fffffe,%ecx
cqto
idiv   %rcx
mov    %eax,-0x14(%rbp)
movslq -0x14(%rbp),%rcx
lea    0x2dba(%rip),%rax
mov    (%rax,%rcx,8),%rax
mov    %rax,0x2d9f(%rip)
mov    -0x10(%rbp),%rax
mov    (%rax),%rdx
movslq -0x14(%rbp),%rcx
lea    0x2d9d(%rip),%rax
mov    %rdx,(%rax,%rcx,8)
cvtsi2sdq 0x2d80(%rip),%xmm1
movsd  0xd50(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
cvtss2sd %xmm0,%xmm0
movsd  0xd2f(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    12ed <ran3+0x1ed>
movsd  0xd1d(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    12fb <ran3+0x1fb>
movss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<ran3>:
mov    (%rdi),%rcx
movabs $0x41a705af1fe3fb79,%r9
05 a7 41
test   %rcx,%rcx
jle    111e <ran3+0x1e>
mov    0x2f17(%rip),%r8
test   %r8,%r8
jne    1197 <ran3+0x97>
mov    %rcx,%rax
neg    %rax
test   %rcx,%rcx
mov    $0x1,%ecx
cmovs  %rax,%rcx
mov    $0x27,%esi
lea    0x2f04(%rip),%r8
jmp    1146 <ran3+0x46>
xchg   %ax,%ax
add    $0xffffffffffffffff,%rsi
jae    1190 <ran3+0x90>
mov    %rcx,%rax
imul   %r9
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0xf,%rdx
add    %rax,%rdx
imul   $0xfffffffffffe0ce3,%rdx,%rax
add    %rcx,%rax
imul   $0x41a7,%rax,%rax
imul   $0xfffffffffffff4ec,%rdx,%rdx
lea    (%rax,%rdx,1),%rcx
add    $0x7fffffff,%rcx
add    %rdx,%rax
cmovns %rax,%rcx
cmp    $0x1f,%rsi
ja     1140 <ran3+0x40>
mov    %rcx,(%r8,%rsi,8)
jmp    1140 <ran3+0x40>
mov    0x2ea9(%rip),%r8
mov    %rcx,%rax
imul   %r9
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0xf,%rdx
add    %rax,%rdx
imul   $0xfffffffffffe0ce3,%rdx,%rax
add    %rcx,%rax
imul   $0x41a7,%rax,%rax
imul   $0xfffffffffffff4ec,%rdx,%rcx
lea    (%rax,%rcx,1),%rsi
add    $0x7fffffff,%rsi
add    %rcx,%rax
cmovns %rax,%rsi
mov    %rsi,(%rdi)
movabs $0x8000004000002001,%rcx
mov    %r8,%rax
imul   %rcx
add    %r8,%rdx
mov    %rdx,%rax
shr    $0x3f,%rax
shr    $0x19,%rdx
add    %eax,%edx
movslq %edx,%rax
lea    0x2e3e(%rip),%rcx
mov    (%rcx,%rax,8),%rdx
mov    %rdx,0x2e23(%rip)
cvtsi2sd %rdx,%xmm0
mulsd  0xde6(%rip),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    %rsi,(%rcx,%rax,8)
cvtss2sd %xmm0,%xmm1
movsd  0xdda(%rip),%xmm0
minsd  %xmm1,%xmm0
ret

LLVM IR (-O0):
@ran3.iy = internal global i64 0
@ran3.iv = internal global [32 x i64] zeroinitializer
define double @ran3(i64* %0) {
%2 = alloca double
%3 = alloca i64*
%4 = alloca i32
%5 = alloca i64
%6 = alloca float
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 %13, label %10
10: ; preds = %1
%11 = load i64, i64* @ran3.iy
%12 = icmp ne i64 %11, 0
br i1 %12, label %65, label %13
13: ; preds = %10, %1
%14 = load i64*, i64** %3
%15 = load i64, i64* %14
%16 = sub nsw i64 0, %15
%17 = icmp slt i64 %16, 1
br i1 %17, label %18, label %20
18: ; preds = %13
%19 = load i64*, i64** %3
store i64 1, i64* %19
br label %25
20: ; preds = %13
%21 = load i64*, i64** %3
%22 = load i64, i64* %21
%23 = sub nsw i64 0, %22
%24 = load i64*, i64** %3
store i64 %23, i64* %24
br label %25
25: ; preds = %20, %18
store i32 39, i32* %4
br label %26
26: ; preds = %60, %25
%27 = load i32, i32* %4
%28 = icmp sge i32 %27, 0
br i1 %28, label %29, label %63
29: ; preds = %26
%30 = load i64*, i64** %3
%31 = load i64, i64* %30
%32 = sdiv i64 %31, 127773
store i64 %32, i64* %5
%33 = load i64*, i64** %3
%34 = load i64, i64* %33
%35 = load i64, i64* %5
%36 = mul nsw i64 %35, 127773
%37 = sub nsw i64 %34, %36
%38 = mul nsw i64 16807, %37
%39 = load i64, i64* %5
%40 = mul nsw i64 2836, %39
%41 = sub nsw i64 %38, %40
%42 = load i64*, i64** %3
store i64 %41, i64* %42
%43 = load i64*, i64** %3
%44 = load i64, i64* %43
%45 = icmp slt i64 %44, 0
br i1 %45, label %46, label %50
46: ; preds = %29
%47 = load i64*, i64** %3
%48 = load i64, i64* %47
%49 = add nsw i64 %48, 2147483647
store i64 %49, i64* %47
br label %50
50: ; preds = %46, %29
%51 = load i32, i32* %4
%52 = icmp slt i32 %51, 32
br i1 %52, label %53, label %59
53: ; preds = %50
%54 = load i64*, i64** %3
%55 = load i64, i64* %54
%56 = load i32, i32* %4
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds [32 x i64], [32 x i64]* @ran3.iv, i64 0, i64 %57
store i64 %55, i64* %58
br label %59
59: ; preds = %53, %50
br label %60
60: ; preds = %59
%61 = load i32, i32* %4
%62 = add nsw i32 %61, -1
store i32 %62, i32* %4
br label %26
63: ; preds = %26
%64 = load i64, i64* getelementptr inbounds ([32 x i64], [32 x i64]* @ran3.iv, i64 0, i64 0)
store i64 %64, i64* @ran3.iy
br label %65
65: ; preds = %63, %10
%66 = load i64*, i64** %3
%67 = load i64, i64* %66
%68 = sdiv i64 %67, 127773
store i64 %68, i64* %5
%69 = load i64*, i64** %3
%70 = load i64, i64* %69
%71 = load i64, i64* %5
%72 = mul nsw i64 %71, 127773
%73 = sub nsw i64 %70, %72
%74 = mul nsw i64 16807, %73
%75 = load i64, i64* %5
%76 = mul nsw i64 2836, %75
%77 = sub nsw i64 %74, %76
%78 = load i64*, i64** %3
store i64 %77, i64* %78
%79 = load i64*, i64** %3
%80 = load i64, i64* %79
%81 = icmp slt i64 %80, 0
br i1 %81, label %82, label %86
82: ; preds = %65
%83 = load i64*, i64** %3
%84 = load i64, i64* %83
%85 = add nsw i64 %84, 2147483647
store i64 %85, i64* %83
br label %86
86: ; preds = %82, %65
%87 = load i64, i64* @ran3.iy
%88 = sdiv i64 %87, 67108862
%89 = trunc i64 %88 to i32
store i32 %89, i32* %4
%90 = load i32, i32* %4
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds [32 x i64], [32 x i64]* @ran3.iv, i64 0, i64 %91
%93 = load i64, i64* %92
store i64 %93, i64* @ran3.iy
%94 = load i64*, i64** %3
%95 = load i64, i64* %94
%96 = load i32, i32* %4
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds [32 x i64], [32 x i64]* @ran3.iv, i64 0, i64 %97
store i64 %95, i64* %98
%99 = load i64, i64* @ran3.iy
%100 = sitofp i64 %99 to double
%101 = fmul double 0x3E00000000200000, %100
%102 = fptrunc double %101 to float
store float %102, float* %6
%103 = fpext float %102 to double
%104 = fcmp ogt double %103, 0x3FEFFFFFBF935359
br i1 %104, label %105, label %106
105: ; preds = %86
store double 0x3FEFFFFFBF935359, double* %2
br label %109
106: ; preds = %86
%107 = load float, float* %6
%108 = fpext float %107 to double
store double %108, double* %2
br label %109
109: ; preds = %106, %105
%110 = load double, double* %2
ret double %110
}

LLVM IR (-O2):
@ran3.iy = internal unnamed_addr global i64 0
@ran3.iv = internal unnamed_addr global [32 x i64] zeroinitializer
define double @ran3(i64* nocapture %0) {
%2 = load i64, i64* %0
%3 = icmp sgt i64 %2, 0
%4 = load i64, i64* @ran3.iy
%5 = icmp ne i64 %4, 0
%6 = select i1 %3, i1 %5, i1 false
br i1 %6, label %31, label %7
7: ; preds = %1
%8 = icmp sgt i64 %2, -1
%9 = sub nsw i64 0, %2
%10 = select i1 %8, i64 1, i64 %9
br label %11
11: ; preds = %7, %26
%12 = phi i64 [ 39, %7 ], [ %27, %26 ]
%13 = phi i64 [ %10, %7 ], [ %22, %26 ]
%14 = sdiv i64 %13, 127773
%15 = mul nsw i64 %14, -127773
%16 = add i64 %15, %13
%17 = mul nsw i64 %16, 16807
%18 = mul nsw i64 %14, -2836
%19 = add i64 %17, %18
%20 = icmp slt i64 %19, 0
%21 = add nsw i64 %19, 2147483647
%22 = select i1 %20, i64 %21, i64 %19
%23 = icmp ult i64 %12, 32
br i1 %23, label %24, label %26
24: ; preds = %11
%25 = getelementptr inbounds [32 x i64], [32 x i64]* @ran3.iv, i64 0, i64 %12
store i64 %22, i64* %25
br label %26
26: ; preds = %11, %24
%27 = add nsw i64 %12, -1
%28 = icmp eq i64 %12, 0
br i1 %28, label %29, label %11
29: ; preds = %26
%30 = load i64, i64* getelementptr inbounds ([32 x i64], [32 x i64]* @ran3.iv, i64 0, i64 0)
br label %31
31: ; preds = %1, %29
%32 = phi i64 [ %4, %1 ], [ %30, %29 ]
%33 = phi i64 [ %2, %1 ], [ %22, %29 ]
%34 = sdiv i64 %33, 127773
%35 = mul nsw i64 %34, -127773
%36 = add i64 %35, %33
%37 = mul nsw i64 %36, 16807
%38 = mul nsw i64 %34, -2836
%39 = add i64 %37, %38
%40 = icmp slt i64 %39, 0
%41 = add nsw i64 %39, 2147483647
%42 = select i1 %40, i64 %41, i64 %39
store i64 %42, i64* %0
%43 = sdiv i64 %32, 67108862
%44 = shl i64 %43, 32
%45 = ashr exact i64 %44, 32
%46 = getelementptr inbounds [32 x i64], [32 x i64]* @ran3.iv, i64 0, i64 %45
%47 = load i64, i64* %46
store i64 %47, i64* @ran3.iy
store i64 %42, i64* %46
%48 = sitofp i64 %47 to double
%49 = fmul double %48, 0x3E00000000200000
%50 = fptrunc double %49 to float
%51 = fpext float %50 to double
%52 = fcmp ogt double %51, 0x3FEFFFFFBF935359
%53 = select i1 %52, double 0x3FEFFFFFBF935359, double %51
ret double %53
}