double powi(b, x)
double b;
long x;
{
  double temp;
  long i;
  if (b == 0.0)
    return ((double)0.0);
  if (x == 0L)
    return ((double)1.0);
  if (x > 0L) {
    temp = b;
    for (i = x - 1; i > 0L; i--)
      temp *= b;
    return temp;
  }
  if (x < 0L) {
    temp = 1.0 / b;
    for (i = x + 1; i < 0L; i++)
      temp *= (1.0 / b);
    return temp;
  }
  return 0;
}


x86 assembly (clang, -O0):
<powi>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
mov    %rdi,-0x18(%rbp)
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1132 <powi+0x32>
jp     1132 <powi+0x32>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    121d <powi+0x11d>
cmpq   $0x0,-0x18(%rbp)
jne    114f <powi+0x4f>
movsd  0xebb(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    121d <powi+0x11d>
cmpq   $0x0,-0x18(%rbp)
jle    11aa <powi+0xaa>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x18(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x28(%rbp)
cmpq   $0x0,-0x28(%rbp)
jle    119b <powi+0x9b>
movsd  -0x10(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x28(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x28(%rbp)
jmp    1170 <powi+0x70>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    121d <powi+0x11d>
cmpq   $0x0,-0x18(%rbp)
jge    1215 <powi+0x115>
movsd  0xe43(%rip),%xmm0
divsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x18(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
cmpq   $0x0,-0x28(%rbp)
jge    1206 <powi+0x106>
movsd  0xe1a(%rip),%xmm0
divsd  -0x10(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    11d3 <powi+0xd3>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    121d <powi+0x11d>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<powi>:
movapd %xmm0,%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jne    1111 <powi+0x11>
jp     1111 <powi+0x11>
ret
test   %rdi,%rdi
je     1123 <powi+0x23>
jle    112c <powi+0x2c>
cmp    $0x1,%rdi
jne    1175 <powi+0x75>
movapd %xmm1,%xmm0
ret
movsd  0xed5(%rip),%xmm0
ret
movsd  0xecc(%rip),%xmm2
divsd  %xmm1,%xmm2
cmp    $0xffffffffffffffff,%rdi
je     1190 <powi+0x90>
mov    %edi,%eax
and    $0x7,%eax
cmp    $0x7,%rax
je     1195 <powi+0x95>
xor    $0x7,%rax
neg    %rax
xor    %ecx,%ecx
movapd %xmm2,%xmm0
cs nopw 0x0(%rax,%rax,1)
mulsd  %xmm2,%xmm0
add    $0xffffffffffffffff,%rcx
cmp    %rcx,%rax
jne    1160 <powi+0x60>
mov    %rdi,%rax
sub    %rcx,%rax
jmp    119c <powi+0x9c>
add    $0x1,%rdi
movapd %xmm1,%xmm0
nopl   (%rax)
mulsd  %xmm1,%xmm0
add    $0xffffffffffffffff,%rdi
cmp    $0x2,%rdi
ja     1180 <powi+0x80>
jmp    1110 <powi+0x10>
movapd %xmm2,%xmm0
ret
mov    %rdi,%rax
movapd %xmm2,%xmm0
add    $0x8,%rdi
cmp    $0x7,%rdi
jb     1110 <powi+0x10>
add    $0x1,%rax
xchg   %ax,%ax
mulsd  %xmm2,%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm2,%xmm0
add    $0x8,%rax
jne    11b0 <powi+0xb0>
jmp    1110 <powi+0x10>

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

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