double power(double r, int p) {
  double tr;
  int neg;
  if (r == 0.0)
    return ((p == 0) ? 1.0 : 0.0);
  if ((neg = (p < 0)))
    p = -p;
  tr = 1.0;
  while (p > 0) {
    if (p & 1)
      tr *= r;
    p >>= 1;
    if (p)
      r *= r;
  }
  return ((neg ? 1.0 / tr : tr));
}


x86 assembly (clang, -O0):
<power>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
mov    %edi,-0x14(%rbp)
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    115e <power+0x5e>
jp     115e <power+0x5e>
mov    -0x14(%rbp),%eax
movsd  0xed1(%rip),%xmm0
xorps  %xmm1,%xmm1
movsd  %xmm1,-0x38(%rbp)
cmp    $0x0,%eax
movsd  %xmm0,-0x30(%rbp)
je     114f <power+0x4f>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1216 <power+0x116>
cmpl   $0x0,-0x14(%rbp)
setl   %al
mov    %al,%cl
and    $0x1,%cl
movzbl %cl,%ecx
mov    %ecx,-0x24(%rbp)
test   $0x1,%al
jne    117d <power+0x7d>
jmp    1185 <power+0x85>
xor    %eax,%eax
sub    -0x14(%rbp),%eax
mov    %eax,-0x14(%rbp)
movsd  0xe73(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
cmpl   $0x0,-0x14(%rbp)
jle    11e1 <power+0xe1>
mov    -0x14(%rbp),%eax
and    $0x1,%eax
cmp    $0x0,%eax
je     11ba <power+0xba>
movsd  -0x10(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
je     11dc <power+0xdc>
movsd  -0x10(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
jmp    1192 <power+0x92>
cmpl   $0x0,-0x24(%rbp)
je     1202 <power+0x102>
movsd  0xe0d(%rip),%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    120c <power+0x10c>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<power>:
movapd %xmm0,%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jne    111d <power+0x1d>
jp     111d <power+0x1d>
test   %edi,%edi
jne    1177 <power+0x77>
movsd  0xee4(%rip),%xmm0
ret
mov    %edi,%eax
neg    %eax
cmovs  %edi,%eax
test   %eax,%eax
je     116b <power+0x6b>
movsd  0xed0(%rip),%xmm2
movapd %xmm2,%xmm0
jmp    114c <power+0x4c>
cs nopw 0x0(%rax,%rax,1)
mulsd  %xmm3,%xmm0
mulsd  %xmm4,%xmm1
mov    %ecx,%eax
jbe    1173 <power+0x73>
movapd %xmm2,%xmm3
test   $0x1,%al
je     1158 <power+0x58>
movapd %xmm1,%xmm3
mov    %eax,%ecx
shr    %ecx
movapd %xmm2,%xmm4
cmp    $0x1,%eax
je     1140 <power+0x40>
movapd %xmm1,%xmm4
jmp    1140 <power+0x40>
movsd  0xe8d(%rip),%xmm0
test   %edi,%edi
js     1178 <power+0x78>
ret
movsd  0xe80(%rip),%xmm1
divsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @power(double %0, i32 %1) {
%3 = alloca double
%4 = alloca double
%5 = alloca i32
%6 = alloca double
%7 = alloca i32
store double %0, double* %4
store i32 %1, i32* %5
%8 = load double, double* %4
%9 = fcmp oeq double %8, 0.000000e+00
br i1 %9, label %10, label %15
10: ; preds = %2
%11 = load i32, i32* %5
%12 = icmp eq i32 %11, 0
%13 = zext i1 %12 to i64
%14 = select i1 %12, double 1.000000e+00, double 0.000000e+00
store double %14, double* %3
br label %54
15: ; preds = %2
%16 = load i32, i32* %5
%17 = icmp slt i32 %16, 0
%18 = zext i1 %17 to i32
store i32 %18, i32* %7
br i1 %17, label %19, label %22
19: ; preds = %15
%20 = load i32, i32* %5
%21 = sub nsw i32 0, %20
store i32 %21, i32* %5
br label %22
22: ; preds = %19, %15
store double 1.000000e+00, double* %6
br label %23
23: ; preds = %43, %22
%24 = load i32, i32* %5
%25 = icmp sgt i32 %24, 0
br i1 %25, label %26, label %44
26: ; preds = %23
%27 = load i32, i32* %5
%28 = and i32 %27, 1
%29 = icmp ne i32 %28, 0
br i1 %29, label %30, label %34
30: ; preds = %26
%31 = load double, double* %4
%32 = load double, double* %6
%33 = fmul double %32, %31
store double %33, double* %6
br label %34
34: ; preds = %30, %26
%35 = load i32, i32* %5
%36 = ashr i32 %35, 1
store i32 %36, i32* %5
%37 = load i32, i32* %5
%38 = icmp ne i32 %37, 0
br i1 %38, label %39, label %43
39: ; preds = %34
%40 = load double, double* %4
%41 = load double, double* %4
%42 = fmul double %41, %40
store double %42, double* %4
br label %43
43: ; preds = %39, %34
br label %23
44: ; preds = %23
%45 = load i32, i32* %7
%46 = icmp ne i32 %45, 0
br i1 %46, label %47, label %50
47: ; preds = %44
%48 = load double, double* %6
%49 = fdiv double 1.000000e+00, %48
br label %52
50: ; preds = %44
%51 = load double, double* %6
br label %52
52: ; preds = %50, %47
%53 = phi double [ %49, %47 ], [ %51, %50 ]
store double %53, double* %3
br label %54
54: ; preds = %52, %10
%55 = load double, double* %3
ret double %55
}

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