double power(double base, int exponent) {
  if (exponent < 0)
    return power(1.0 / base, -exponent);
  else {
    double answer;
    answer = 1.0;
    while (exponent > 0) {
      answer = answer * base;
      exponent = exponent - 1;
    }
    return answer;
  }
}


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)
jge    114b <power+0x3b>
movsd  0xece(%rip),%xmm0
divsd  -0x10(%rbp),%xmm0
xor    %edi,%edi
sub    -0x14(%rbp),%edi
call   1030 <power@plt>
movsd  %xmm0,-0x8(%rbp)
jmp    1189 <power+0x79>
movsd  0xead(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
cmpl   $0x0,-0x14(%rbp)
jle    117f <power+0x6f>
movsd  -0x20(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1158 <power+0x48>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<power>:
test   %edi,%edi
js     110c <power+0xc>
movapd %xmm0,%xmm1
mov    %edi,%eax
jmp    1134 <power+0x34>
movsd  0xeec(%rip),%xmm2
mov    %edi,%eax
cs nopw 0x0(%rax,%rax,1)
movapd %xmm2,%xmm1
divsd  %xmm0,%xmm1
neg    %eax
movapd %xmm1,%xmm0
test   %edi,%edi
mov    %eax,%edi
jg     1120 <power+0x20>
test   %eax,%eax
jle    1164 <power+0x64>
lea    -0x1(%rax),%ecx
mov    %eax,%edx
and    $0x7,%edx
je     116d <power+0x6d>
xor    %esi,%esi
movsd  0xeb4(%rip),%xmm0
nopl   0x0(%rax)
mulsd  %xmm1,%xmm0
add    $0x1,%esi
cmp    %esi,%edx
jne    1150 <power+0x50>
sub    %esi,%eax
cmp    $0x7,%ecx
jae    117a <power+0x7a>
jmp    11a8 <power+0xa8>
movsd  0xe94(%rip),%xmm0
ret
movsd  0xe8b(%rip),%xmm0
cmp    $0x7,%ecx
jb     11a8 <power+0xa8>
add    $0xffffffff,%eax
nopl   (%rax)
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
add    $0xfffffff8,%eax
cmp    $0xfffffffe,%eax
jb     1180 <power+0x80>
ret

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

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