double cfrac_tan(double x, unsigned int n) {
  double k = n * 2 - 1;
  double result = 1 / (k / x);
  for (unsigned int i = n - 1; i > 0; i--) {
    k -= 2;
    result = 1 / (k / x - result);
  }
  return result;
}


x86 assembly (clang, -O0):
<cfrac_tan>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %edi,-0xc(%rbp)
mov    -0xc(%rbp),%ecx
mov    %ecx,%eax
lea    -0x1(%rax,%rax,1),%eax
cvtsi2sd %rax,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm1
divsd  -0x8(%rbp),%xmm1
movsd  0xecf(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jbe    1191 <cfrac_tan+0x91>
movsd  -0x18(%rbp),%xmm0
movsd  0xeae(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm1
divsd  -0x8(%rbp),%xmm1
subsd  -0x20(%rbp),%xmm1
movsd  0xe86(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1143 <cfrac_tan+0x43>
movsd  -0x20(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<cfrac_tan>:
lea    (%rdi,%rdi,1),%eax
add    $0xffffffff,%eax
cvtsi2sd %rax,%xmm2
movapd %xmm2,%xmm4
divsd  %xmm0,%xmm4
mov    %edi,%eax
movsd  0xee3(%rip),%xmm3
movapd %xmm3,%xmm1
divsd  %xmm4,%xmm1
add    $0xffffffff,%eax
je     1195 <cfrac_tan+0x95>
test   $0x1,%al
je     114d <cfrac_tan+0x4d>
addsd  0xed2(%rip),%xmm2
movapd %xmm2,%xmm4
divsd  %xmm0,%xmm4
subsd  %xmm1,%xmm4
lea    -0x2(%rdi),%eax
movapd %xmm3,%xmm1
divsd  %xmm4,%xmm1
cmp    $0x2,%edi
je     1195 <cfrac_tan+0x95>
movsd  0xeae(%rip),%xmm4
nopw   0x0(%rax,%rax,1)
addsd  %xmm4,%xmm2
movapd %xmm2,%xmm5
divsd  %xmm0,%xmm5
subsd  %xmm1,%xmm5
movapd %xmm3,%xmm1
divsd  %xmm5,%xmm1
addsd  %xmm4,%xmm2
movapd %xmm2,%xmm5
divsd  %xmm0,%xmm5
subsd  %xmm1,%xmm5
movapd %xmm3,%xmm1
divsd  %xmm5,%xmm1
add    $0xfffffffe,%eax
jne    1160 <cfrac_tan+0x60>
movapd %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @cfrac_tan(double %0, i32 %1) {
%3 = alloca double
%4 = alloca i32
%5 = alloca double
%6 = alloca double
%7 = alloca i32
store double %0, double* %3
store i32 %1, i32* %4
%8 = load i32, i32* %4
%9 = mul i32 %8, 2
%10 = sub i32 %9, 1
%11 = uitofp i32 %10 to double
store double %11, double* %5
%12 = load double, double* %5
%13 = load double, double* %3
%14 = fdiv double %12, %13
%15 = fdiv double 1.000000e+00, %14
store double %15, double* %6
%16 = load i32, i32* %4
%17 = sub i32 %16, 1
store i32 %17, i32* %7
br label %18
18: ; preds = %30, %2
%19 = load i32, i32* %7
%20 = icmp ugt i32 %19, 0
br i1 %20, label %21, label %33
21: ; preds = %18
%22 = load double, double* %5
%23 = fsub double %22, 2.000000e+00
store double %23, double* %5
%24 = load double, double* %5
%25 = load double, double* %3
%26 = fdiv double %24, %25
%27 = load double, double* %6
%28 = fsub double %26, %27
%29 = fdiv double 1.000000e+00, %28
store double %29, double* %6
br label %30
30: ; preds = %21
%31 = load i32, i32* %7
%32 = add i32 %31, -1
store i32 %32, i32* %7
br label %18
33: ; preds = %18
%34 = load double, double* %6
ret double %34
}

LLVM IR (-O2):
define double @cfrac_tan(double %0, i32 %1) {
%3 = shl i32 %1, 1
%4 = add i32 %3, -1
%5 = uitofp i32 %4 to double
%6 = fdiv double %5, %0
%7 = add i32 %1, -1
%8 = fdiv double 1.000000e+00, %6
%9 = icmp eq i32 %7, 0
br i1 %9, label %25, label %10
10: ; preds = %2
%11 = and i32 %7, 1
%12 = icmp eq i32 %11, 0
br i1 %12, label %19, label %13
13: ; preds = %10
%14 = fadd double %5, -2.000000e+00
%15 = fdiv double %14, %0
%16 = fsub double %15, %8
%17 = add i32 %1, -2
%18 = fdiv double 1.000000e+00, %16
br label %19
19: ; preds = %13, %10
%20 = phi double [ undef, %10 ], [ %18, %13 ]
%21 = phi double [ %8, %10 ], [ %18, %13 ]
%22 = phi i32 [ %7, %10 ], [ %17, %13 ]
%23 = phi double [ %5, %10 ], [ %14, %13 ]
%24 = icmp eq i32 %1, 2
br i1 %24, label %25, label %27
25: ; preds = %19, %27, %2
%26 = phi double [ %8, %2 ], [ %20, %19 ], [ %39, %27 ]
ret double %26
27: ; preds = %19, %27
%28 = phi double [ %39, %27 ], [ %21, %19 ]
%29 = phi i32 [ %38, %27 ], [ %22, %19 ]
%30 = phi double [ %35, %27 ], [ %23, %19 ]
%31 = fadd double %30, -2.000000e+00
%32 = fdiv double %31, %0
%33 = fsub double %32, %28
%34 = fdiv double 1.000000e+00, %33
%35 = fadd double %31, -2.000000e+00
%36 = fdiv double %35, %0
%37 = fsub double %36, %34
%38 = add i32 %29, -2
%39 = fdiv double 1.000000e+00, %37
%40 = icmp eq i32 %38, 0
br i1 %40, label %25, label %27
}