double DoubleFactorial(int n) {
  double r = 1;
  for (int i = 1; i <= n; ++i) {
    r *= (2 * i - 1);
  }
  return r;
}


x86 assembly (clang, -O0):
<DoubleFactorial>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movsd  0xef1(%rip),%xmm0
movsd  %xmm0,-0x10(%rbp)
movl   $0x1,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     114c <DoubleFactorial+0x4c>
mov    -0x14(%rbp),%eax
shl    $0x1,%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    111b <DoubleFactorial+0x1b>
movsd  -0x10(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<DoubleFactorial>:
test   %edi,%edi
jle    1120 <DoubleFactorial+0x20>
lea    -0x1(%rdi),%ecx
mov    %edi,%eax
and    $0x3,%eax
cmp    $0x3,%ecx
jae    1129 <DoubleFactorial+0x29>
movsd  0xee7(%rip),%xmm0
mov    $0x1,%esi
jmp    118d <DoubleFactorial+0x8d>
movsd  0xed8(%rip),%xmm0
ret
and    $0xfffffffc,%edi
neg    %edi
movsd  0xeca(%rip),%xmm0
mov    $0x1,%edx
mov    $0x7,%ecx
lea    -0x6(%rcx),%esi
xorps  %xmm1,%xmm1
cvtsi2sd %esi,%xmm1
mulsd  %xmm0,%xmm1
lea    -0x4(%rcx),%esi
xorps  %xmm0,%xmm0
cvtsi2sd %esi,%xmm0
mulsd  %xmm1,%xmm0
lea    -0x2(%rcx),%esi
xorps  %xmm1,%xmm1
cvtsi2sd %esi,%xmm1
mulsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
mulsd  %xmm1,%xmm0
add    $0x8,%ecx
mov    %edx,%esi
add    $0x4,%esi
lea    (%rdi,%rdx,1),%r8d
add    $0x4,%r8d
mov    %esi,%edx
cmp    $0x1,%r8d
jne    1140 <DoubleFactorial+0x40>
test   %eax,%eax
je     11b3 <DoubleFactorial+0xb3>
lea    (%rsi,%rsi,1),%ecx
add    $0xffffffff,%ecx
nopw   0x0(%rax,%rax,1)
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
mulsd  %xmm1,%xmm0
add    $0x2,%ecx
add    $0xffffffff,%eax
jne    11a0 <DoubleFactorial+0xa0>
ret

LLVM IR (-O0):
define double @DoubleFactorial(i32 %0) {
%2 = alloca i32
%3 = alloca double
%4 = alloca i32
store i32 %0, i32* %2
store double 1.000000e+00, double* %3
store i32 1, i32* %4
br label %5
5: ; preds = %16, %1
%6 = load i32, i32* %4
%7 = load i32, i32* %2
%8 = icmp sle i32 %6, %7
br i1 %8, label %9, label %19
9: ; preds = %5
%10 = load i32, i32* %4
%11 = mul nsw i32 2, %10
%12 = sub nsw i32 %11, 1
%13 = sitofp i32 %12 to double
%14 = load double, double* %3
%15 = fmul double %14, %13
store double %15, double* %3
br label %16
16: ; preds = %9
%17 = load i32, i32* %4
%18 = add nsw i32 %17, 1
store i32 %18, i32* %4
br label %5
19: ; preds = %5
%20 = load double, double* %3
ret double %20
}

LLVM IR (-O2):
define double @DoubleFactorial(i32 %0) {
%2 = icmp slt i32 %0, 1
br i1 %2, label %25, label %3
3: ; preds = %1
%4 = add i32 %0, -1
%5 = and i32 %0, 3
%6 = icmp ult i32 %4, 3
br i1 %6, label %9, label %7
7: ; preds = %3
%8 = and i32 %0, -4
br label %27
9: ; preds = %27, %3
%10 = phi double [ undef, %3 ], [ %46, %27 ]
%11 = phi i32 [ 1, %3 ], [ %47, %27 ]
%12 = phi double [ 1.000000e+00, %3 ], [ %46, %27 ]
%13 = icmp eq i32 %5, 0
br i1 %13, label %25, label %14
14: ; preds = %9, %14
%15 = phi i32 [ %22, %14 ], [ %11, %9 ]
%16 = phi double [ %21, %14 ], [ %12, %9 ]
%17 = phi i32 [ %23, %14 ], [ 0, %9 ]
%18 = shl nuw nsw i32 %15, 1
%19 = add nsw i32 %18, -1
%20 = sitofp i32 %19 to double
%21 = fmul double %16, %20
%22 = add nuw i32 %15, 1
%23 = add i32 %17, 1
%24 = icmp eq i32 %23, %5
br i1 %24, label %25, label %14
25: ; preds = %9, %14, %1
%26 = phi double [ 1.000000e+00, %1 ], [ %10, %9 ], [ %21, %14 ]
ret double %26
27: ; preds = %27, %7
%28 = phi i32 [ 1, %7 ], [ %47, %27 ]
%29 = phi double [ 1.000000e+00, %7 ], [ %46, %27 ]
%30 = phi i32 [ 0, %7 ], [ %48, %27 ]
%31 = shl nuw nsw i32 %28, 1
%32 = add nsw i32 %31, -1
%33 = sitofp i32 %32 to double
%34 = fmul double %29, %33
%35 = shl i32 %28, 1
%36 = or i32 %35, 1
%37 = sitofp i32 %36 to double
%38 = fmul double %34, %37
%39 = shl i32 %28, 1
%40 = add i32 %39, 3
%41 = sitofp i32 %40 to double
%42 = fmul double %38, %41
%43 = shl i32 %28, 1
%44 = add i32 %43, 5
%45 = sitofp i32 %44 to double
%46 = fmul double %42, %45
%47 = add nuw i32 %28, 4
%48 = add i32 %30, 4
%49 = icmp eq i32 %48, %8
br i1 %49, label %9, label %27
}