double grib_power(long s, long n) {
  double divisor = 1.0;
  while (s < 0) {
    divisor /= n;
    s++;
  }
  while (s > 0) {
    divisor *= n;
    s--;
  }
  return divisor;
}


x86 assembly (clang, -O0):
<grib_power>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  0xeec(%rip),%xmm0
movsd  %xmm0,-0x18(%rbp)
cmpq   $0x0,-0x8(%rbp)
jge    1149 <grib_power+0x49>
cvtsi2sdq -0x10(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x8(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x8(%rbp)
jmp    1119 <grib_power+0x19>
jmp    114e <grib_power+0x4e>
cmpq   $0x0,-0x8(%rbp)
jle    117a <grib_power+0x7a>
cvtsi2sdq -0x10(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x8(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x8(%rbp)
jmp    114e <grib_power+0x4e>
movsd  -0x18(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<grib_power>:
test   %rdi,%rdi
js     117c <grib_power+0x7c>
je     11d7 <grib_power+0xd7>
cvtsi2sd %rsi,%xmm1
lea    -0x1(%rdi),%rax
mov    %rdi,%rcx
movsd  0xee1(%rip),%xmm0
and    $0x7,%rcx
je     1140 <grib_power+0x40>
xor    %edx,%edx
nopw   0x0(%rax,%rax,1)
mulsd  %xmm1,%xmm0
add    $0x1,%rdx
cmp    %rdx,%rcx
jne    1130 <grib_power+0x30>
sub    %rdx,%rdi
cmp    $0x7,%rax
jb     11d6 <grib_power+0xd6>
add    $0xffffffffffffffff,%rdi
xchg   %ax,%ax
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
add    $0xfffffffffffffff8,%rdi
cmp    $0xfffffffffffffffe,%rdi
jb     1150 <grib_power+0x50>
jmp    11d6 <grib_power+0xd6>
cvtsi2sd %rsi,%xmm1
mov    %edi,%eax
neg    %eax
and    $0x3,%rax
je     11e0 <grib_power+0xe0>
neg    %rax
movsd  0xe6a(%rip),%xmm0
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
divsd  %xmm1,%xmm0
add    $0xffffffffffffffff,%rcx
cmp    %rcx,%rax
jne    11a0 <grib_power+0xa0>
mov    %rdi,%rax
sub    %rcx,%rax
cmp    $0xfffffffffffffffc,%rdi
ja     11d6 <grib_power+0xd6>
nopl   0x0(%rax)
divsd  %xmm1,%xmm0
divsd  %xmm1,%xmm0
divsd  %xmm1,%xmm0
divsd  %xmm1,%xmm0
add    $0x4,%rax
jne    11c0 <grib_power+0xc0>
ret
movsd  0xe21(%rip),%xmm0
ret
movsd  0xe18(%rip),%xmm0
mov    %rdi,%rax
cmp    $0xfffffffffffffffc,%rdi
jbe    11c0 <grib_power+0xc0>
jmp    11d6 <grib_power+0xd6>

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

LLVM IR (-O2):
define double @grib_power(i64 %0, i64 %1) {
%3 = icmp slt i64 %0, 0
br i1 %3, label %4, label %22
4: ; preds = %2
%5 = sitofp i64 %1 to double
%6 = sub i64 0, %0
%7 = and i64 %6, 3
%8 = icmp eq i64 %7, 0
br i1 %8, label %17, label %9
9: ; preds = %4, %9
%10 = phi double [ %13, %9 ], [ 1.000000e+00, %4 ]
%11 = phi i64 [ %14, %9 ], [ %0, %4 ]
%12 = phi i64 [ %15, %9 ], [ 0, %4 ]
%13 = fdiv double %10, %5
%14 = add i64 %11, 1
%15 = add i64 %12, 1
%16 = icmp eq i64 %15, %7
br i1 %16, label %17, label %9
17: ; preds = %9, %4
%18 = phi double [ undef, %4 ], [ %13, %9 ]
%19 = phi double [ 1.000000e+00, %4 ], [ %13, %9 ]
%20 = phi i64 [ %0, %4 ], [ %14, %9 ]
%21 = icmp ugt i64 %0, -4
br i1 %21, label %65, label %42
22: ; preds = %2
%23 = icmp eq i64 %0, 0
br i1 %23, label %65, label %24
24: ; preds = %22
%25 = sitofp i64 %1 to double
%26 = add i64 %0, -1
%27 = and i64 %0, 7
%28 = icmp eq i64 %27, 0
br i1 %28, label %37, label %29
29: ; preds = %24, %29
%30 = phi double [ %33, %29 ], [ 1.000000e+00, %24 ]
%31 = phi i64 [ %34, %29 ], [ %0, %24 ]
%32 = phi i64 [ %35, %29 ], [ 0, %24 ]
%33 = fmul double %30, %25
%34 = add nsw i64 %31, -1
%35 = add i64 %32, 1
%36 = icmp eq i64 %35, %27
br i1 %36, label %37, label %29
37: ; preds = %29, %24
%38 = phi double [ 1.000000e+00, %24 ], [ %33, %29 ]
%39 = phi i64 [ %0, %24 ], [ %34, %29 ]
%40 = phi double [ undef, %24 ], [ %33, %29 ]
%41 = icmp ult i64 %26, 7
br i1 %41, label %65, label %51
42: ; preds = %17, %42
%43 = phi double [ %48, %42 ], [ %19, %17 ]
%44 = phi i64 [ %49, %42 ], [ %20, %17 ]
%45 = fdiv double %43, %5
%46 = fdiv double %45, %5
%47 = fdiv double %46, %5
%48 = fdiv double %47, %5
%49 = add i64 %44, 4
%50 = icmp eq i64 %49, 0
br i1 %50, label %65, label %42
51: ; preds = %37, %51
%52 = phi double [ %61, %51 ], [ %38, %37 ]
%53 = phi i64 [ %62, %51 ], [ %39, %37 ]
%54 = fmul double %52, %25
%55 = fmul double %54, %25
%56 = fmul double %55, %25
%57 = fmul double %56, %25
%58 = fmul double %57, %25
%59 = fmul double %58, %25
%60 = fmul double %59, %25
%61 = fmul double %60, %25
%62 = add nsw i64 %53, -8
%63 = add i64 %53, -9
%64 = icmp ult i64 %63, -2
br i1 %64, label %51, label %65
65: ; preds = %37, %51, %17, %42, %22
%66 = phi double [ 1.000000e+00, %22 ], [ %18, %17 ], [ %48, %42 ], [ %40, %37 ], [ %61, %51 ]
ret double %66
}