double dChoose(long lNum, long lDenom) {
  long lStart;
  double dResult;
  long i, j;
  if (lNum < lDenom)
    return 0.0;
  if (lNum < 0.0l || lDenom < 0.0l)
    return 0.0;
  if (lDenom > lNum - lDenom) {
    lStart = lDenom;
    lDenom = lNum - lDenom;
  } else {
    lStart = lNum - lDenom;
  }
  dResult = 1.0;
  i = lNum;
  j = lDenom;
  while (i > lStart || j > 1) {
    if (i > lStart) {
      if (j > 1)
        dResult *= ((double)i-- / j--);
      else
        dResult *= i--;
    } else if (j > 1)
      dResult /= j--;
  }
  return (dResult);
}


x86 assembly (clang, -O0):
<dChoose>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x18(%rbp),%rax
jge    1127 <dChoose+0x27>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    12a9 <dChoose+0x1a9>
mov    -0x10(%rbp),%rax
mov    %rax,-0x40(%rbp)
fildll -0x40(%rbp)
fldz
fucomip %st(1),%st
fstp   %st(0)
ja     115f <dChoose+0x5f>
jmp    1143 <dChoose+0x43>
mov    -0x18(%rbp),%rax
mov    %rax,-0x48(%rbp)
fildll -0x48(%rbp)
fldz
fucomip %st(1),%st
fstp   %st(0)
jbe    116c <dChoose+0x6c>
jmp    115f <dChoose+0x5f>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    12a9 <dChoose+0x1a9>
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
sub    -0x18(%rbp),%rcx
cmp    %rcx,%rax
jle    119a <dChoose+0x9a>
mov    -0x18(%rbp),%rax
mov    %rax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
sub    -0x18(%rbp),%rax
mov    %rax,-0x18(%rbp)
jmp    11a6 <dChoose+0xa6>
mov    -0x10(%rbp),%rax
sub    -0x18(%rbp),%rax
mov    %rax,-0x20(%rbp)
movsd  0xe52(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x30(%rbp),%rcx
mov    $0x1,%al
cmp    -0x20(%rbp),%rcx
mov    %al,-0x49(%rbp)
jg     11e1 <dChoose+0xe1>
cmpq   $0x1,-0x38(%rbp)
setg   %al
mov    %al,-0x49(%rbp)
mov    -0x49(%rbp),%al
test   $0x1,%al
jne    11f1 <dChoose+0xf1>
jmp    129f <dChoose+0x19f>
mov    -0x30(%rbp),%rax
cmp    -0x20(%rbp),%rax
jle    1268 <dChoose+0x168>
cmpq   $0x1,-0x38(%rbp)
jle    1245 <dChoose+0x145>
mov    -0x30(%rbp),%rax
mov    %rax,%rcx
add    $0xffffffffffffffff,%rcx
mov    %rcx,-0x30(%rbp)
cvtsi2sd %rax,%xmm0
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0xffffffffffffffff,%rcx
mov    %rcx,-0x38(%rbp)
cvtsi2sd %rax,%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    1263 <dChoose+0x163>
mov    -0x30(%rbp),%rax
mov    %rax,%rcx
add    $0xffffffffffffffff,%rcx
mov    %rcx,-0x30(%rbp)
cvtsi2sd %rax,%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    129a <dChoose+0x19a>
cmpq   $0x1,-0x38(%rbp)
jle    1295 <dChoose+0x195>
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0xffffffffffffffff,%rcx
mov    %rcx,-0x38(%rbp)
cvtsi2sd %rax,%xmm1
movsd  -0x28(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    129a <dChoose+0x19a>
jmp    11c3 <dChoose+0xc3>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<dChoose>:
mov    %rdi,%r8
xorpd  %xmm0,%xmm0
sub    %rsi,%r8
jl     11aa <dChoose+0xaa>
mov    %rsi,%rax
or     %rdi,%rax
js     11aa <dChoose+0xaa>
cmp    %rsi,%r8
mov    %rsi,%rcx
cmovl  %r8,%rcx
cmovl  %rsi,%r8
cmp    $0x2,%rcx
setge  %dl
cmp    %rdi,%r8
setl   %al
jl     1148 <dChoose+0x48>
cmp    $0x2,%rcx
jge    1148 <dChoose+0x48>
movsd  0xeb9(%rip),%xmm0
ret
movsd  0xeb0(%rip),%xmm0
test   $0x1,%al
je     1180 <dChoose+0x80>
xorps  %xmm1,%xmm1
cvtsi2sd %rdi,%xmm1
add    $0xffffffffffffffff,%rdi
test   $0x1,%dl
je     1175 <dChoose+0x75>
xorps  %xmm2,%xmm2
cvtsi2sd %rcx,%xmm2
add    $0xffffffffffffffff,%rcx
divsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
jmp    1195 <dChoose+0x95>
nopl   0x0(%rax,%rax,1)
test   $0x1,%dl
je     1195 <dChoose+0x95>
xorps  %xmm1,%xmm1
cvtsi2sd %rcx,%xmm1
add    $0xffffffffffffffff,%rcx
divsd  %xmm1,%xmm0
cmp    $0x2,%rcx
setge  %dl
cmp    %r8,%rdi
setg   %al
jg     1150 <dChoose+0x50>
cmp    $0x1,%rcx
jg     1150 <dChoose+0x50>
ret

LLVM IR (-O0):
define double @dChoose(i64 %0, i64 %1) {
%3 = alloca double
%4 = alloca i64
%5 = alloca i64
%6 = alloca i64
%7 = alloca double
%8 = alloca i64
%9 = alloca i64
store i64 %0, i64* %4
store i64 %1, i64* %5
%10 = load i64, i64* %4
%11 = load i64, i64* %5
%12 = icmp slt i64 %10, %11
br i1 %12, label %13, label %14
13: ; preds = %2
store double 0.000000e+00, double* %3
br label %87
14: ; preds = %2
%15 = load i64, i64* %4
%16 = sitofp i64 %15 to x86_fp80
%17 = fcmp olt x86_fp80 %16, 0xK00000000000000000000
br i1 %17, label %22, label %18
18: ; preds = %14
%19 = load i64, i64* %5
%20 = sitofp i64 %19 to x86_fp80
%21 = fcmp olt x86_fp80 %20, 0xK00000000000000000000
br i1 %21, label %22, label %23
22: ; preds = %18, %14
store double 0.000000e+00, double* %3
br label %87
23: ; preds = %18
%24 = load i64, i64* %5
%25 = load i64, i64* %4
%26 = load i64, i64* %5
%27 = sub nsw i64 %25, %26
%28 = icmp sgt i64 %24, %27
br i1 %28, label %29, label %34
29: ; preds = %23
%30 = load i64, i64* %5
store i64 %30, i64* %6
%31 = load i64, i64* %4
%32 = load i64, i64* %5
%33 = sub nsw i64 %31, %32
store i64 %33, i64* %5
br label %38
34: ; preds = %23
%35 = load i64, i64* %4
%36 = load i64, i64* %5
%37 = sub nsw i64 %35, %36
store i64 %37, i64* %6
br label %38
38: ; preds = %34, %29
store double 1.000000e+00, double* %7
%39 = load i64, i64* %4
store i64 %39, i64* %8
%40 = load i64, i64* %5
store i64 %40, i64* %9
br label %41
41: ; preds = %84, %38
%42 = load i64, i64* %8
%43 = load i64, i64* %6
%44 = icmp sgt i64 %42, %43
br i1 %44, label %48, label %45
45: ; preds = %41
%46 = load i64, i64* %9
%47 = icmp sgt i64 %46, 1
br label %48
48: ; preds = %45, %41
%49 = phi i1 [ true, %41 ], [ %47, %45 ]
br i1 %49, label %50, label %85
50: ; preds = %48
%51 = load i64, i64* %8
%52 = load i64, i64* %6
%53 = icmp sgt i64 %51, %52
br i1 %53, label %54, label %74
54: ; preds = %50
%55 = load i64, i64* %9
%56 = icmp sgt i64 %55, 1
br i1 %56, label %57, label %67
57: ; preds = %54
%58 = load i64, i64* %8
%59 = add nsw i64 %58, -1
store i64 %59, i64* %8
%60 = sitofp i64 %58 to double
%61 = load i64, i64* %9
%62 = add nsw i64 %61, -1
store i64 %62, i64* %9
%63 = sitofp i64 %61 to double
%64 = fdiv double %60, %63
%65 = load double, double* %7
%66 = fmul double %65, %64
store double %66, double* %7
br label %73
67: ; preds = %54
%68 = load i64, i64* %8
%69 = add nsw i64 %68, -1
store i64 %69, i64* %8
%70 = sitofp i64 %68 to double
%71 = load double, double* %7
%72 = fmul double %71, %70
store double %72, double* %7
br label %73
73: ; preds = %67, %57
br label %84
74: ; preds = %50
%75 = load i64, i64* %9
%76 = icmp sgt i64 %75, 1
br i1 %76, label %77, label %83
77: ; preds = %74
%78 = load i64, i64* %9
%79 = add nsw i64 %78, -1
store i64 %79, i64* %9
%80 = sitofp i64 %78 to double
%81 = load double, double* %7
%82 = fdiv double %81, %80
store double %82, double* %7
br label %83
83: ; preds = %77, %74
br label %84
84: ; preds = %83, %73
br label %41
85: ; preds = %48
%86 = load double, double* %7
store double %86, double* %3
br label %87
87: ; preds = %85, %22, %13
%88 = load double, double* %3
ret double %88
}

LLVM IR (-O2):
define double @dChoose(i64 %0, i64 %1) {
%3 = icmp sge i64 %0, %1
%4 = or i64 %1, %0
%5 = icmp sgt i64 %4, -1
%6 = and i1 %3, %5
br i1 %6, label %7, label %43
7: ; preds = %2
%8 = sub nsw i64 %0, %1
%9 = icmp slt i64 %8, %1
%10 = select i1 %9, i64 %8, i64 %1
%11 = select i1 %9, i64 %1, i64 %8
%12 = icmp slt i64 %11, %0
%13 = icmp sgt i64 %10, 1
%14 = select i1 %12, i1 true, i1 %13
br i1 %14, label %15, label %43
15: ; preds = %7, %36
%16 = phi i1 [ %41, %36 ], [ %13, %7 ]
%17 = phi i1 [ %40, %36 ], [ %12, %7 ]
%18 = phi i64 [ %39, %36 ], [ %10, %7 ]
%19 = phi i64 [ %38, %36 ], [ %0, %7 ]
%20 = phi double [ %37, %36 ], [ 1.000000e+00, %7 ]
br i1 %17, label %21, label %31
21: ; preds = %15
%22 = add nsw i64 %19, -1
%23 = sitofp i64 %19 to double
br i1 %16, label %24, label %29
24: ; preds = %21
%25 = add nsw i64 %18, -1
%26 = sitofp i64 %18 to double
%27 = fdiv double %23, %26
%28 = fmul double %20, %27
br label %36
29: ; preds = %21
%30 = fmul double %20, %23
br label %36
31: ; preds = %15
br i1 %16, label %32, label %36
32: ; preds = %31
%33 = add nsw i64 %18, -1
%34 = sitofp i64 %18 to double
%35 = fdiv double %20, %34
br label %36
36: ; preds = %31, %32, %24, %29
%37 = phi double [ %28, %24 ], [ %30, %29 ], [ %35, %32 ], [ %20, %31 ]
%38 = phi i64 [ %22, %24 ], [ %22, %29 ], [ %19, %32 ], [ %19, %31 ]
%39 = phi i64 [ %25, %24 ], [ %18, %29 ], [ %33, %32 ], [ %18, %31 ]
%40 = icmp sgt i64 %38, %11
%41 = icmp sgt i64 %39, 1
%42 = select i1 %40, i1 true, i1 %41
br i1 %42, label %15, label %43
43: ; preds = %36, %7, %2
%44 = phi double [ 0.000000e+00, %2 ], [ 1.000000e+00, %7 ], [ %37, %36 ]
ret double %44
}