double pow_ri(float *ap, int *bp) {
  double pow, x;
  int n;
  unsigned long u;
  pow = 1;
  x = *ap;
  n = *bp;
  if (n != 0) {
    if (n < 0) {
      n = -n;
      x = 1 / x;
    }
    for (u = n;;) {
      if (u & 01)
        pow *= x;
      if (u >>= 1)
        x *= x;
      else
        break;
    }
  }
  return (pow);
}


x86 assembly (clang, -O0):
<pow_ri>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  0xeec(%rip),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
je     11c3 <pow_ri+0xc3>
cmpl   $0x0,-0x24(%rbp)
jge    1161 <pow_ri+0x61>
xor    %eax,%eax
sub    -0x24(%rbp),%eax
mov    %eax,-0x24(%rbp)
movsd  0xea9(%rip),%xmm0
divsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movslq -0x24(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
and    $0x1,%rax
cmp    $0x0,%rax
je     118a <pow_ri+0x8a>
movsd  -0x20(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x30(%rbp),%rax
shr    $0x1,%rax
mov    %rax,-0x30(%rbp)
cmp    $0x0,%rax
je     11b4 <pow_ri+0xb4>
movsd  -0x20(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    11b9 <pow_ri+0xb9>
jmp    11be <pow_ri+0xbe>
jmp    1169 <pow_ri+0x69>
jmp    11c3 <pow_ri+0xc3>
movsd  -0x18(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<pow_ri>:
mov    (%rsi),%ecx
test   %ecx,%ecx
je     1133 <pow_ri+0x33>
movss  (%rdi),%xmm0
cvtss2sd %xmm0,%xmm1
jns    1122 <pow_ri+0x22>
neg    %ecx
movsd  0xee6(%rip),%xmm0
divsd  %xmm1,%xmm0
movapd %xmm0,%xmm1
mov    %ecx,%eax
test   $0x1,%al
je     113c <pow_ri+0x3c>
movapd %xmm1,%xmm0
cmp    $0x2,%ecx
jb     1149 <pow_ri+0x49>
jmp    114a <pow_ri+0x4a>
movsd  0xec5(%rip),%xmm0
ret
movsd  0xebc(%rip),%xmm0
cmp    $0x2,%ecx
jae    114a <pow_ri+0x4a>
ret
movsd  0xeae(%rip),%xmm2
jmp    1173 <pow_ri+0x73>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mulsd  %xmm3,%xmm0
mov    %rax,%rcx
shr    %rcx
cmp    $0x3,%rax
mov    %rcx,%rax
jbe    1149 <pow_ri+0x49>
mulsd  %xmm1,%xmm1
movapd %xmm2,%xmm3
test   $0x2,%al
je     1160 <pow_ri+0x60>
movapd %xmm1,%xmm3
jmp    1160 <pow_ri+0x60>

LLVM IR (-O0):
define double @pow_ri(float* %0, i32* %1) {
%3 = alloca float*
%4 = alloca i32*
%5 = alloca double
%6 = alloca double
%7 = alloca i32
%8 = alloca i64
store float* %0, float** %3
store i32* %1, i32** %4
store double 1.000000e+00, double* %5
%9 = load float*, float** %3
%10 = load float, float* %9
%11 = fpext float %10 to double
store double %11, double* %6
%12 = load i32*, i32** %4
%13 = load i32, i32* %12
store i32 %13, i32* %7
%14 = load i32, i32* %7
%15 = icmp ne i32 %14, 0
br i1 %15, label %16, label %46
16: ; preds = %2
%17 = load i32, i32* %7
%18 = icmp slt i32 %17, 0
br i1 %18, label %19, label %24
19: ; preds = %16
%20 = load i32, i32* %7
%21 = sub nsw i32 0, %20
store i32 %21, i32* %7
%22 = load double, double* %6
%23 = fdiv double 1.000000e+00, %22
store double %23, double* %6
br label %24
24: ; preds = %19, %16
%25 = load i32, i32* %7
%26 = sext i32 %25 to i64
store i64 %26, i64* %8
br label %27
27: ; preds = %44, %24
%28 = load i64, i64* %8
%29 = and i64 %28, 1
%30 = icmp ne i64 %29, 0
br i1 %30, label %31, label %35
31: ; preds = %27
%32 = load double, double* %6
%33 = load double, double* %5
%34 = fmul double %33, %32
store double %34, double* %5
br label %35
35: ; preds = %31, %27
%36 = load i64, i64* %8
%37 = lshr i64 %36, 1
store i64 %37, i64* %8
%38 = icmp ne i64 %37, 0
br i1 %38, label %39, label %43
39: ; preds = %35
%40 = load double, double* %6
%41 = load double, double* %6
%42 = fmul double %41, %40
store double %42, double* %6
br label %44
43: ; preds = %35
br label %45
44: ; preds = %39
br label %27
45: ; preds = %43
br label %46
46: ; preds = %45, %2
%47 = load double, double* %5
ret double %47
}

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