double power(double x, int n) {
  double y;
  if (n == 0)
    return 1;
  if (x == 0) {
    if (n < 0) {
      return 1e200;
    }
    return 0;
  }
  if (n < 0) {
    x = 1. / x;
    n = -n;
  }
  y = 1.;
  while (n > 0) {
    if (n & 1) {
      y *= x;
    }
    n /= 2;
    x *= x;
  }
  return y;
}


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    11ee <power+0xee>
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1169 <power+0x69>
jp     1169 <power+0x69>
cmpl   $0x0,-0x14(%rbp)
jge    115c <power+0x5c>
movsd  0xeb6(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    11ee <power+0xee>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    11ee <power+0xee>
cmpl   $0x0,-0x14(%rbp)
jge    118d <power+0x8d>
movsd  0xe85(%rip),%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
xor    %eax,%eax
sub    -0x14(%rbp),%eax
mov    %eax,-0x14(%rbp)
movsd  0xe6b(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
cmpl   $0x0,-0x14(%rbp)
jle    11e4 <power+0xe4>
mov    -0x14(%rbp),%eax
and    $0x1,%eax
cmp    $0x0,%eax
je     11c2 <power+0xc2>
movsd  -0x10(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x14(%rbp)
movsd  -0x10(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
jmp    119a <power+0x9a>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<power>:
test   %edi,%edi
je     1119 <power+0x19>
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1122 <power+0x22>
jp     1122 <power+0x22>
test   %edi,%edi
js     1173 <power+0x73>
movapd %xmm1,%xmm0
ret
movsd  0xedf(%rip),%xmm0
ret
test   %edi,%edi
jns    1138 <power+0x38>
movsd  0xed2(%rip),%xmm1
divsd  %xmm0,%xmm1
neg    %edi
movapd %xmm1,%xmm0
movsd  0xec0(%rip),%xmm2
movapd %xmm2,%xmm1
jmp    1163 <power+0x63>
cs nopw 0x0(%rax,%rax,1)
mulsd  %xmm3,%xmm1
mov    %edi,%eax
shr    %eax
mulsd  %xmm0,%xmm0
cmp    $0x1,%edi
mov    %eax,%edi
jbe    1114 <power+0x14>
movapd %xmm2,%xmm3
test   $0x1,%dil
je     1150 <power+0x50>
movapd %xmm0,%xmm3
jmp    1150 <power+0x50>
movsd  0xe8d(%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 %46
10: ; preds = %2
%11 = load double, double* %4
%12 = fcmp oeq double %11, 0.000000e+00
br i1 %12, label %13, label %18
13: ; preds = %10
%14 = load i32, i32* %5
%15 = icmp slt i32 %14, 0
br i1 %15, label %16, label %17
16: ; preds = %13
store double 0x6974E718D7D7625A, double* %3
br label %46
17: ; preds = %13
store double 0.000000e+00, double* %3
br label %46
18: ; preds = %10
%19 = load i32, i32* %5
%20 = icmp slt i32 %19, 0
br i1 %20, label %21, label %26
21: ; preds = %18
%22 = load double, double* %4
%23 = fdiv double 1.000000e+00, %22
store double %23, double* %4
%24 = load i32, i32* %5
%25 = sub nsw i32 0, %24
store i32 %25, i32* %5
br label %26
26: ; preds = %21, %18
store double 1.000000e+00, double* %6
br label %27
27: ; preds = %38, %26
%28 = load i32, i32* %5
%29 = icmp sgt i32 %28, 0
br i1 %29, label %30, label %44
30: ; preds = %27
%31 = load i32, i32* %5
%32 = and i32 %31, 1
%33 = icmp ne i32 %32, 0
br i1 %33, label %34, label %38
34: ; preds = %30
%35 = load double, double* %4
%36 = load double, double* %6
%37 = fmul double %36, %35
store double %37, double* %6
br label %38
38: ; preds = %34, %30
%39 = load i32, i32* %5
%40 = sdiv i32 %39, 2
store i32 %40, i32* %5
%41 = load double, double* %4
%42 = load double, double* %4
%43 = fmul double %42, %41
store double %43, double* %4
br label %27
44: ; preds = %27
%45 = load double, double* %6
store double %45, double* %3
br label %46
46: ; preds = %44, %17, %16, %9
%47 = load double, double* %3
ret double %47
}

LLVM IR (-O2):
define double @power(double %0, i32 %1) {
%3 = icmp eq i32 %1, 0
br i1 %3, label %27, label %4
4: ; preds = %2
%5 = fcmp oeq double %0, 0.000000e+00
%6 = icmp slt i32 %1, 0
br i1 %5, label %7, label %9
7: ; preds = %4
%8 = select i1 %6, double 0x6974E718D7D7625A, double 0.000000e+00
br label %27
9: ; preds = %4
br i1 %6, label %10, label %13
10: ; preds = %9
%11 = fdiv double 1.000000e+00, %0
%12 = sub nsw i32 0, %1
br label %13
13: ; preds = %9, %10
%14 = phi i32 [ %12, %10 ], [ %1, %9 ]
%15 = phi double [ %11, %10 ], [ %0, %9 ]
br label %16
16: ; preds = %13, %16
%17 = phi double [ %23, %16 ], [ 1.000000e+00, %13 ]
%18 = phi i32 [ %24, %16 ], [ %14, %13 ]
%19 = phi double [ %25, %16 ], [ %15, %13 ]
%20 = and i32 %18, 1
%21 = icmp eq i32 %20, 0
%22 = select i1 %21, double 1.000000e+00, double %19
%23 = fmul double %17, %22
%24 = lshr i32 %18, 1
%25 = fmul double %19, %19
%26 = icmp ult i32 %18, 2
br i1 %26, label %27, label %16
27: ; preds = %16, %7, %2
%28 = phi double [ 1.000000e+00, %2 ], [ %8, %7 ], [ %23, %16 ]
ret double %28
}