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


x86 assembly (clang, -O0):
<power>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
mov    %edi,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
jne    1128 <power+0x28>
movsd  0xee2(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    118f <power+0x8f>
cmpl   $0x0,-0x14(%rbp)
jge    114c <power+0x4c>
movsd  0xec6(%rip),%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
xor    %eax,%eax
sub    -0x14(%rbp),%eax
mov    %eax,-0x14(%rbp)
jmp    1151 <power+0x51>
movsd  0xea7(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
cmpl   $0x0,-0x14(%rbp)
jle    1185 <power+0x85>
movsd  -0x10(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    115e <power+0x5e>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<power>:
test   %edi,%edi
je     1141 <power+0x41>
jns    1118 <power+0x18>
movsd  0xef2(%rip),%xmm1
divsd  %xmm0,%xmm1
neg    %edi
movapd %xmm1,%xmm0
add    $0x1,%edi
movsd  0xedd(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mulsd  %xmm0,%xmm1
add    $0xffffffff,%edi
cmp    $0x1,%edi
ja     1130 <power+0x30>
movapd %xmm1,%xmm0
ret
movsd  0xeb7(%rip),%xmm0
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 eq i32 %7, 0
br i1 %8, label %9, label %10
9: ; preds = %2
store double 1.000000e+00, double* %3
br label %31
10: ; preds = %2
%11 = load i32, i32* %5
%12 = icmp slt i32 %11, 0
br i1 %12, label %13, label %18
13: ; preds = %10
%14 = load double, double* %4
%15 = fdiv double 1.000000e+00, %14
store double %15, double* %4
%16 = load i32, i32* %5
%17 = sub nsw i32 0, %16
store i32 %17, i32* %5
br label %18
18: ; preds = %13, %10
br label %19
19: ; preds = %18
store double 1.000000e+00, double* %6
br label %20
20: ; preds = %23, %19
%21 = load i32, i32* %5
%22 = icmp sgt i32 %21, 0
br i1 %22, label %23, label %29
23: ; preds = %20
%24 = load double, double* %4
%25 = load double, double* %6
%26 = fmul double %25, %24
store double %26, double* %6
%27 = load i32, i32* %5
%28 = add nsw i32 %27, -1
store i32 %28, i32* %5
br label %20
29: ; preds = %20
%30 = load double, double* %6
store double %30, double* %3
br label %31
31: ; preds = %29, %9
%32 = load double, double* %3
ret double %32
}

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