double power(double x, int n) {
  if (n == 0) {
    return 1;
  } else if (n == 1) {
    return x;
  } else if (n % 2 == 0) {
    return power(x * x, n / 2);
  } else {
    return x * power(x * x, (n - 1) / 2);
  }
}


x86 assembly (clang, -O0):
<power>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
movsd  %xmm0,-0x10(%rbp)
mov    %edi,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
jne    113c <power+0x2c>
movsd  0xece(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    11c9 <power+0xb9>
cmpl   $0x1,-0x14(%rbp)
jne    1155 <power+0x45>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    11c9 <power+0xb9>
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    118f <power+0x7f>
movsd  -0x10(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%edi
call   1030 <power@plt>
movsd  %xmm0,-0x8(%rbp)
jmp    11c9 <power+0xb9>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x10(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%edi
call   1030 <power@plt>
movaps %xmm0,%xmm1
movsd  -0x20(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<power>:
push   %rax
movapd %xmm0,%xmm1
test   %edi,%edi
je     1146 <power+0x36>
nopl   0x0(%rax)
cmp    $0x1,%edi
je     1172 <power+0x62>
movapd %xmm1,%xmm0
mulsd  %xmm1,%xmm0
test   $0x1,%dil
jne    1150 <power+0x40>
mov    %edi,%eax
shr    $0x1f,%eax
add    %edi,%eax
sar    %eax
movapd %xmm0,%xmm1
mov    %eax,%edi
test   %edi,%edi
jne    1120 <power+0x10>
movsd  0xeb2(%rip),%xmm1
jmp    1172 <power+0x62>
mov    %edi,%eax
add    $0xffffffff,%eax
shr    $0x1f,%eax
add    %eax,%edi
add    $0xffffffff,%edi
sar    %edi
movsd  %xmm1,(%rsp)
call   1030 <power@plt>
movsd  (%rsp),%xmm1
mulsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
pop    %rax
ret

LLVM IR (-O0):
define double @power(double %0, i32 %1) {
%3 = alloca double
%4 = alloca double
%5 = alloca i32
store double %0, double* %4
store i32 %1, i32* %5
%6 = load i32, i32* %5
%7 = icmp eq i32 %6, 0
br i1 %7, label %8, label %9
8: ; preds = %2
store double 1.000000e+00, double* %3
br label %35
9: ; preds = %2
%10 = load i32, i32* %5
%11 = icmp eq i32 %10, 1
br i1 %11, label %12, label %14
12: ; preds = %9
%13 = load double, double* %4
store double %13, double* %3
br label %35
14: ; preds = %9
%15 = load i32, i32* %5
%16 = srem i32 %15, 2
%17 = icmp eq i32 %16, 0
br i1 %17, label %18, label %25
18: ; preds = %14
%19 = load double, double* %4
%20 = load double, double* %4
%21 = fmul double %19, %20
%22 = load i32, i32* %5
%23 = sdiv i32 %22, 2
%24 = call double @power(double %21, i32 %23)
store double %24, double* %3
br label %35
25: ; preds = %14
%26 = load double, double* %4
%27 = load double, double* %4
%28 = load double, double* %4
%29 = fmul double %27, %28
%30 = load i32, i32* %5
%31 = sub nsw i32 %30, 1
%32 = sdiv i32 %31, 2
%33 = call double @power(double %29, i32 %32)
%34 = fmul double %26, %33
store double %34, double* %3
br label %35
35: ; preds = %25, %18, %12, %8
%36 = load double, double* %3
ret double %36
}

LLVM IR (-O2):
define double @power(double %0, i32 %1) {
br label %3
3: ; preds = %10, %2
%4 = phi double [ %0, %2 ], [ %9, %10 ]
%5 = phi i32 [ %1, %2 ], [ %11, %10 ]
switch i32 %5, label %6 [
i32 0, label %12
i32 1, label %13
]
6: ; preds = %3
%7 = and i32 %5, 1
%8 = icmp eq i32 %7, 0
%9 = fmul double %4, %4
br i1 %8, label %10, label %15
10: ; preds = %6
%11 = sdiv i32 %5, 2
br label %3
12: ; preds = %3
br label %13
13: ; preds = %3, %12, %15
%14 = phi double [ %19, %15 ], [ 1.000000e+00, %12 ], [ %4, %3 ]
ret double %14
15: ; preds = %6
%16 = add nsw i32 %5, -1
%17 = sdiv i32 %16, 2
%18 = tail call double @power(double %9, i32 %17)
%19 = fmul double %4, %18
br label %13
}