double ran2d(idum)
long *idum;
{
  static long iy, ir[98];
  static int iff = 0;
  int j;
  void nrerror();
  if (*idum < 0 || iff == 0) {
    iff = 1;
    if ((*idum = (150889 - (*idum)) % 714025) < 0)
      *idum = -(*idum);
    for (j = 1; j <= 97; j++) {
      *idum = (1366 * (*idum) + 150889) % 714025;
      ir[j] = (*idum);
    }
    *idum = (1366 * (*idum) + 150889) % 714025;
    iy = (*idum);
  }
  j = 1 + 97.0 * iy / 714025;
  if (j > 97 || j < 1)
    nrerror("RAN2: This cannot happen.");
  iy = ir[j];
  *idum = (1366 * (*idum) + 150889) % 714025;
  ir[j] = (*idum);
  return (double)iy / 714025;
}


x86 assembly (clang, -O0):
<ran2d>:
push   %rbp
mov    %rsp,%rbp
sub    $0x10,%rsp
mov    %rdi,-0x8(%rbp)
mov    -0x8(%rbp),%rax
cmpq   $0x0,(%rax)
jl     1137 <ran2d+0x27>
cmpl   $0x0,0x322f(%rip)
jne    11ff <ran2d+0xef>
movl   $0x1,0x321f(%rip)
mov    -0x8(%rbp),%rcx
mov    $0x24d69,%eax
sub    (%rcx),%rax
mov    $0xae529,%ecx
cqto
idiv   %rcx
mov    -0x8(%rbp),%rax
mov    %rdx,(%rax)
cmp    $0x0,%rdx
jge    1178 <ran2d+0x68>
mov    -0x8(%rbp),%rax
xor    %ecx,%ecx
sub    (%rax),%rcx
mov    -0x8(%rbp),%rax
mov    %rcx,(%rax)
movl   $0x1,-0xc(%rbp)
cmpl   $0x61,-0xc(%rbp)
jg     11cf <ran2d+0xbf>
mov    -0x8(%rbp),%rax
imul   $0x556,(%rax),%rax
add    $0x24d69,%rax
mov    $0xae529,%ecx
cqto
idiv   %rcx
mov    -0x8(%rbp),%rax
mov    %rdx,(%rax)
mov    -0x8(%rbp),%rax
mov    (%rax),%rdx
movslq -0xc(%rbp),%rcx
lea    0x2e93(%rip),%rax
mov    %rdx,(%rax,%rcx,8)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    117f <ran2d+0x6f>
mov    -0x8(%rbp),%rax
imul   $0x556,(%rax),%rax
add    $0x24d69,%rax
mov    $0xae529,%ecx
cqto
idiv   %rcx
mov    -0x8(%rbp),%rax
mov    %rdx,(%rax)
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
mov    %rax,0x2e41(%rip)
cvtsi2sdq 0x2e38(%rip),%xmm0
movsd  0xe00(%rip),%xmm1
mulsd  %xmm0,%xmm1
movsd  0xdec(%rip),%xmm0
divsd  %xmm0,%xmm1
movsd  0xdd8(%rip),%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0xc(%rbp)
cmpl   $0x61,-0xc(%rbp)
jg     1247 <ran2d+0x137>
cmpl   $0x1,-0xc(%rbp)
jge    1255 <ran2d+0x145>
lea    0xdca(%rip),%rdi
mov    $0x0,%al
call   1030 <nrerror@plt>
movslq -0xc(%rbp),%rcx
lea    0x2df0(%rip),%rax
mov    (%rax,%rcx,8),%rax
mov    %rax,0x2dd5(%rip)
mov    -0x8(%rbp),%rax
imul   $0x556,(%rax),%rax
add    $0x24d69,%rax
mov    $0xae529,%ecx
cqto
idiv   %rcx
mov    -0x8(%rbp),%rax
mov    %rdx,(%rax)
mov    -0x8(%rbp),%rax
mov    (%rax),%rdx
movslq -0xc(%rbp),%rcx
lea    0x2db1(%rip),%rax
mov    %rdx,(%rax,%rcx,8)
cvtsi2sdq 0x2d94(%rip),%xmm0
movsd  0xd54(%rip),%xmm1
divsd  %xmm1,%xmm0
add    $0x10,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ran2d>:
push   %r15
push   %r14
push   %rbx
mov    %rdi,%r14
mov    (%rdi),%rax
movabs $0x5dfc995681872319,%rbx
99 fc 5d
test   %rax,%rax
js     113f <ran2d+0x2f>
testb  $0x1,0x322f(%rip)
je     113f <ran2d+0x2f>
mov    0x2f06(%rip),%rcx
jmp    1235 <ran2d+0x125>
movb   $0x1,0x321a(%rip)
mov    $0x24d69,%esi
sub    %rax,%rsi
mov    %rsi,%rax
imul   %rbx
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0x12,%rdx
add    %rax,%rdx
imul   $0xae529,%rdx,%rax
sub    %rax,%rsi
mov    %rsi,%rcx
neg    %rcx
cmovs  %rsi,%rcx
mov    $0x2,%esi
lea    0x2ece(%rip),%rdi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
imul   $0x556,%rcx,%rcx
add    $0x24d69,%rcx
mov    %rcx,%rax
imul   %rbx
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0x12,%rdx
add    %rax,%rdx
imul   $0xae529,%rdx,%rax
sub    %rax,%rcx
mov    %rcx,-0x8(%rdi,%rsi,8)
cmp    $0x62,%rsi
je     11fd <ran2d+0xed>
imul   $0x556,%rcx,%rcx
add    $0x24d69,%rcx
mov    %rcx,%rax
imul   %rbx
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0x12,%rdx
add    %rax,%rdx
imul   $0xae529,%rdx,%rax
sub    %rax,%rcx
mov    %rcx,(%rdi,%rsi,8)
add    $0x2,%rsi
jmp    1190 <ran2d+0x80>
imul   $0x556,%ecx,%eax
add    $0x24d69,%eax
cltq
imul   $0x5dfc9957,%rax,%rcx
mov    %rcx,%rdx
shr    $0x3f,%rdx
sar    $0x32,%rcx
add    %edx,%ecx
imul   $0xae529,%ecx,%ecx
sub    %ecx,%eax
cltq
mov    %rax,(%r14)
mov    %rax,0x2e0e(%rip)
mov    %rax,%rcx
cvtsi2sd %rcx,%xmm0
mulsd  0xdbe(%rip),%xmm0
divsd  0xdbe(%rip),%xmm0
addsd  0xdbe(%rip),%xmm0
cvttsd2si %xmm0,%r15d
lea    -0x62(%r15),%ecx
cmp    $0xffffff9e,%ecx
ja     1271 <ran2d+0x161>
lea    0xdb1(%rip),%rdi
xor    %eax,%eax
call   1030 <nrerror@plt>
mov    (%r14),%rax
movslq %r15d,%rsi
lea    0x2dd5(%rip),%r8
mov    (%r8,%rsi,8),%rdi
imul   $0x556,%rax,%rcx
add    $0x24d69,%rcx
mov    %rcx,%rax
imul   %rbx
mov    %rdi,0x2da6(%rip)
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0x12,%rdx
add    %rax,%rdx
imul   $0xae529,%rdx,%rax
sub    %rax,%rcx
mov    %rcx,(%r14)
mov    %rcx,(%r8,%rsi,8)
xorps  %xmm0,%xmm0
cvtsi2sd %rdi,%xmm0
divsd  0xd3f(%rip),%xmm0
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
@ran2d.iy = internal global i64 0
@ran2d.ir = internal global [98 x i64] zeroinitializer
@ran2d.iff = internal global i32 0
@.str = [26 x i8] c"RAN2: This cannot happen.\00"
define double @ran2d(i64* %0) {
%2 = alloca i64*
%3 = alloca i32
store i64* %0, i64** %2
%4 = load i64*, i64** %2
%5 = load i64, i64* %4
%6 = icmp slt i64 %5, 0
br i1 %6, label %10, label %7
7: ; preds = %1
%8 = load i32, i32* @ran2d.iff
%9 = icmp eq i32 %8, 0
br i1 %9, label %10, label %50
10: ; preds = %7, %1
store i32 1, i32* @ran2d.iff
%11 = load i64*, i64** %2
%12 = load i64, i64* %11
%13 = sub nsw i64 150889, %12
%14 = srem i64 %13, 714025
%15 = load i64*, i64** %2
store i64 %14, i64* %15
%16 = icmp slt i64 %14, 0
br i1 %16, label %17, label %22
17: ; preds = %10
%18 = load i64*, i64** %2
%19 = load i64, i64* %18
%20 = sub nsw i64 0, %19
%21 = load i64*, i64** %2
store i64 %20, i64* %21
br label %22
22: ; preds = %17, %10
store i32 1, i32* %3
br label %23
23: ; preds = %38, %22
%24 = load i32, i32* %3
%25 = icmp sle i32 %24, 97
br i1 %25, label %26, label %41
26: ; preds = %23
%27 = load i64*, i64** %2
%28 = load i64, i64* %27
%29 = mul nsw i64 1366, %28
%30 = add nsw i64 %29, 150889
%31 = srem i64 %30, 714025
%32 = load i64*, i64** %2
store i64 %31, i64* %32
%33 = load i64*, i64** %2
%34 = load i64, i64* %33
%35 = load i32, i32* %3
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds [98 x i64], [98 x i64]* @ran2d.ir, i64 0, i64 %36
store i64 %34, i64* %37
br label %38
38: ; preds = %26
%39 = load i32, i32* %3
%40 = add nsw i32 %39, 1
store i32 %40, i32* %3
br label %23
41: ; preds = %23
%42 = load i64*, i64** %2
%43 = load i64, i64* %42
%44 = mul nsw i64 1366, %43
%45 = add nsw i64 %44, 150889
%46 = srem i64 %45, 714025
%47 = load i64*, i64** %2
store i64 %46, i64* %47
%48 = load i64*, i64** %2
%49 = load i64, i64* %48
store i64 %49, i64* @ran2d.iy
br label %50
50: ; preds = %41, %7
%51 = load i64, i64* @ran2d.iy
%52 = sitofp i64 %51 to double
%53 = fmul double 9.700000e+01, %52
%54 = fdiv double %53, 7.140250e+05
%55 = fadd double 1.000000e+00, %54
%56 = fptosi double %55 to i32
store i32 %56, i32* %3
%57 = load i32, i32* %3
%58 = icmp sgt i32 %57, 97
br i1 %58, label %62, label %59
59: ; preds = %50
%60 = load i32, i32* %3
%61 = icmp slt i32 %60, 1
br i1 %61, label %62, label %63
62: ; preds = %59, %50
br label %63
63: ; preds = %62, %59
%64 = load i32, i32* %3
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds [98 x i64], [98 x i64]* @ran2d.ir, i64 0, i64 %65
%67 = load i64, i64* %66
store i64 %67, i64* @ran2d.iy
%68 = load i64*, i64** %2
%69 = load i64, i64* %68
%70 = mul nsw i64 1366, %69
%71 = add nsw i64 %70, 150889
%72 = srem i64 %71, 714025
%73 = load i64*, i64** %2
store i64 %72, i64* %73
%74 = load i64*, i64** %2
%75 = load i64, i64* %74
%76 = load i32, i32* %3
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds [98 x i64], [98 x i64]* @ran2d.ir, i64 0, i64 %77
store i64 %75, i64* %78
%79 = load i64, i64* @ran2d.iy
%80 = sitofp i64 %79 to double
%81 = fdiv double %80, 7.140250e+05
ret double %81
}

LLVM IR (-O2):
@ran2d.iy = internal unnamed_addr global i64 0
@ran2d.ir = internal unnamed_addr global [98 x i64] zeroinitializer
@ran2d.iff = internal unnamed_addr global i1 false
@.str = [26 x i8] c"RAN2: This cannot happen.\00"
define double @ran2d(i64* nocapture %0) {
%2 = load i64, i64* %0
%3 = icmp slt i64 %2, 0
%4 = load i1, i1* @ran2d.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 i64, i64* @ran2d.iy
br label %34
9: ; preds = %1
store i1 true, i1* @ran2d.iff
%10 = sub nsw i64 150889, %2
%11 = srem i64 %10, 714025
%12 = tail call i64 @llvm.abs.i64(i64 %11, i1 true)
br label %13
13: ; preds = %22, %9
%14 = phi i64 [ 1, %9 ], [ %27, %22 ]
%15 = phi i64 [ %12, %9 ], [ %25, %22 ]
%16 = mul nsw i64 %15, 1366
%17 = add nsw i64 %16, 150889
%18 = srem i64 %17, 714025
%19 = getelementptr inbounds [98 x i64], [98 x i64]* @ran2d.ir, i64 0, i64 %14
store i64 %18, i64* %19
%20 = add nuw nsw i64 %14, 1
%21 = icmp eq i64 %20, 98
br i1 %21, label %28, label %22
22: ; preds = %13
%23 = mul nsw i64 %18, 1366
%24 = add nsw i64 %23, 150889
%25 = srem i64 %24, 714025
%26 = getelementptr inbounds [98 x i64], [98 x i64]* @ran2d.ir, i64 0, i64 %20
store i64 %25, i64* %26
%27 = add nuw nsw i64 %14, 2
br label %13
28: ; preds = %13
%29 = trunc i64 %18 to i32
%30 = mul i32 %29, 1366
%31 = add i32 %30, 150889
%32 = srem i32 %31, 714025
%33 = sext i32 %32 to i64
store i64 %33, i64* %0
store i64 %33, i64* @ran2d.iy
br label %34
34: ; preds = %7, %28
%35 = phi i64 [ %2, %7 ], [ %33, %28 ]
%36 = phi i64 [ %8, %7 ], [ %33, %28 ]
%37 = sitofp i64 %36 to double
%38 = fmul double %37, 9.700000e+01
%39 = fdiv double %38, 7.140250e+05
%40 = fadd double %39, 1.000000e+00
%41 = fptosi double %40 to i32
%42 = add i32 %41, -98
%43 = icmp ult i32 %42, -97
br i1 %43, label %44, label %46
44: ; preds = %34
tail call void (i8*, ...) bitcast (void (...)* @nrerror to void (i8*, ...)*)(i8* getelementptr inbounds ([26 x i8], [26 x i8]* @.str, i64 0, i64 0))
%45 = load i64, i64* %0
br label %46
46: ; preds = %34, %44
%47 = phi i64 [ %35, %34 ], [ %45, %44 ]
%48 = sext i32 %41 to i64
%49 = getelementptr inbounds [98 x i64], [98 x i64]* @ran2d.ir, i64 0, i64 %48
%50 = load i64, i64* %49
store i64 %50, i64* @ran2d.iy
%51 = mul nsw i64 %47, 1366
%52 = add nsw i64 %51, 150889
%53 = srem i64 %52, 714025
store i64 %53, i64* %0
store i64 %53, i64* %49
%54 = sitofp i64 %50 to double
%55 = fdiv double %54, 7.140250e+05
ret double %55
}