double umath_double_pow(double x, int n) {
  double result = 1.0;
  int neg = 0;
  if (n < 0) {
    n = -n;
    neg = 1;
  }
  while (n-- != 0) {
    result *= x;
  }
  if (neg != 0) {
    result = 1.0 / result;
  }
  return result;
}


x86 assembly (clang, -O0):
<umath_double_pow>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %edi,-0xc(%rbp)
movsd  0xeec(%rip),%xmm0
movsd  %xmm0,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
cmpl   $0x0,-0xc(%rbp)
jge    1139 <umath_double_pow+0x39>
xor    %eax,%eax
sub    -0xc(%rbp),%eax
mov    %eax,-0xc(%rbp)
movl   $0x1,-0x1c(%rbp)
jmp    113e <umath_double_pow+0x3e>
mov    -0xc(%rbp),%eax
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0xc(%rbp)
cmp    $0x0,%eax
je     1166 <umath_double_pow+0x66>
movsd  -0x8(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
jmp    113e <umath_double_pow+0x3e>
cmpl   $0x0,-0x1c(%rbp)
je     1182 <umath_double_pow+0x82>
movsd  0xe88(%rip),%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<umath_double_pow>:
mov    %edi,%eax
neg    %eax
cmovs  %edi,%eax
test   %eax,%eax
je     1134 <umath_double_pow+0x34>
lea    -0x1(%rax),%ecx
mov    %eax,%edx
and    $0x7,%edx
je     113e <umath_double_pow+0x3e>
xor    %esi,%esi
movsd  0xee1(%rip),%xmm1
nop
mulsd  %xmm0,%xmm1
add    $0x1,%esi
cmp    %esi,%edx
jne    1120 <umath_double_pow+0x20>
sub    %esi,%eax
cmp    $0x7,%ecx
jae    1150 <umath_double_pow+0x50>
jmp    1175 <umath_double_pow+0x75>
movsd  0xec4(%rip),%xmm1
jmp    1175 <umath_double_pow+0x75>
movsd  0xeba(%rip),%xmm1
cmp    $0x7,%ecx
jb     1175 <umath_double_pow+0x75>
nopl   0x0(%rax,%rax,1)
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
add    $0xfffffff8,%eax
jne    1150 <umath_double_pow+0x50>
test   %edi,%edi
js     117e <umath_double_pow+0x7e>
movapd %xmm1,%xmm0
ret
movsd  0xe7a(%rip),%xmm0
divsd  %xmm1,%xmm0
movapd %xmm0,%xmm1
movapd %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @umath_double_pow(double %0, i32 %1) {
%3 = alloca double
%4 = alloca i32
%5 = alloca double
%6 = alloca i32
store double %0, double* %3
store i32 %1, i32* %4
store double 1.000000e+00, double* %5
store i32 0, i32* %6
%7 = load i32, i32* %4
%8 = icmp slt i32 %7, 0
br i1 %8, label %9, label %12
9: ; preds = %2
%10 = load i32, i32* %4
%11 = sub nsw i32 0, %10
store i32 %11, i32* %4
store i32 1, i32* %6
br label %12
12: ; preds = %9, %2
br label %13
13: ; preds = %17, %12
%14 = load i32, i32* %4
%15 = add nsw i32 %14, -1
store i32 %15, i32* %4
%16 = icmp ne i32 %14, 0
br i1 %16, label %17, label %21
17: ; preds = %13
%18 = load double, double* %3
%19 = load double, double* %5
%20 = fmul double %19, %18
store double %20, double* %5
br label %13
21: ; preds = %13
%22 = load i32, i32* %6
%23 = icmp ne i32 %22, 0
br i1 %23, label %24, label %27
24: ; preds = %21
%25 = load double, double* %5
%26 = fdiv double 1.000000e+00, %25
store double %26, double* %5
br label %27
27: ; preds = %24, %21
%28 = load double, double* %5
ret double %28
}

LLVM IR (-O2):
define double @umath_double_pow(double %0, i32 %1) {
%3 = icmp slt i32 %1, 0
%4 = sub nsw i32 0, %1
%5 = select i1 %3, i32 %4, i32 %1
%6 = icmp eq i32 %5, 0
br i1 %6, label %37, label %7
7: ; preds = %2
%8 = add nsw i32 %5, -1
%9 = and i32 %5, 7
%10 = icmp eq i32 %9, 0
br i1 %10, label %19, label %11
11: ; preds = %7, %11
%12 = phi double [ %16, %11 ], [ 1.000000e+00, %7 ]
%13 = phi i32 [ %15, %11 ], [ %5, %7 ]
%14 = phi i32 [ %17, %11 ], [ 0, %7 ]
%15 = add nsw i32 %13, -1
%16 = fmul double %12, %0
%17 = add i32 %14, 1
%18 = icmp eq i32 %17, %9
br i1 %18, label %19, label %11
19: ; preds = %11, %7
%20 = phi double [ undef, %7 ], [ %16, %11 ]
%21 = phi double [ 1.000000e+00, %7 ], [ %16, %11 ]
%22 = phi i32 [ %5, %7 ], [ %15, %11 ]
%23 = icmp ult i32 %8, 7
br i1 %23, label %37, label %24
24: ; preds = %19, %24
%25 = phi double [ %35, %24 ], [ %21, %19 ]
%26 = phi i32 [ %34, %24 ], [ %22, %19 ]
%27 = fmul double %25, %0
%28 = fmul double %27, %0
%29 = fmul double %28, %0
%30 = fmul double %29, %0
%31 = fmul double %30, %0
%32 = fmul double %31, %0
%33 = fmul double %32, %0
%34 = add nsw i32 %26, -8
%35 = fmul double %33, %0
%36 = icmp eq i32 %34, 0
br i1 %36, label %37, label %24
37: ; preds = %19, %24, %2
%38 = phi double [ 1.000000e+00, %2 ], [ %20, %19 ], [ %35, %24 ]
%39 = fdiv double 1.000000e+00, %38
%40 = select i1 %3, double %39, double %38
ret double %40
}