double TIMESTE(pow, a)
double a;
int pow;
{
  int p;
  register double n, t, r, s;
  p = pow < 0 ? -pow : pow;
  n = 5.0;
  t = 2.0;
  if ((p & 1) == 0) {
    r = 1.0;
    s = 1.0;
  } else {
    r = 5.0;
    s = 2.0;
  }
  p >>= 1;
  while (p != 0) {
    n *= n;
    t *= t;
    if ((p & 1) != 0) {
      r *= n;
      s *= t;
    }
    p >>= 1;
  }
  if (pow < 0) {
    return ((a / r) / s);
  } else {
    return ((a * r) * s);
  }
}


x86 assembly (clang, -O0):
<TIMESTE>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
cmpl   $0x0,-0xc(%rbp)
jge    1123 <TIMESTE+0x23>
xor    %eax,%eax
sub    -0xc(%rbp),%eax
mov    %eax,-0x44(%rbp)
jmp    1129 <TIMESTE+0x29>
mov    -0xc(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x1c(%rbp)
movsd  0xed1(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xebc(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x1c(%rbp),%eax
and    $0x1,%eax
cmp    $0x0,%eax
jne    1177 <TIMESTE+0x77>
movsd  0xeb0(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xea3(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    1191 <TIMESTE+0x91>
movsd  0xe89(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe74(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x1c(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x0,-0x1c(%rbp)
je     11fd <TIMESTE+0xfd>
movsd  -0x28(%rbp),%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x1c(%rbp),%eax
and    $0x1,%eax
cmp    $0x0,%eax
je     11ef <TIMESTE+0xef>
movsd  -0x28(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x1c(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    119a <TIMESTE+0x9a>
cmpl   $0x0,-0xc(%rbp)
jge    1220 <TIMESTE+0x120>
movsd  -0x18(%rbp),%xmm0
divsd  -0x38(%rbp),%xmm0
divsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1234 <TIMESTE+0x134>
movsd  -0x18(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<TIMESTE>:
mov    %edi,%eax
neg    %eax
cmovs  %edi,%eax
xor    %ecx,%ecx
test   $0x1,%al
sete   %cl
lea    0xeeb(%rip),%rdx
movsd  (%rdx,%rcx,8),%xmm1
lea    0xeef(%rip),%rdx
movsd  (%rdx,%rcx,8),%xmm2
cmp    $0x2,%eax
jae    113c <TIMESTE+0x3c>
test   %edi,%edi
js     1181 <TIMESTE+0x81>
mulsd  %xmm0,%xmm1
mulsd  %xmm2,%xmm1
movapd %xmm1,%xmm0
ret
movapd 0xedc(%rip),%xmm3
movsd  0xee4(%rip),%xmm4
mov    %eax,%ecx
jmp    115f <TIMESTE+0x5f>
mulsd  %xmm6,%xmm2
mulsd  %xmm5,%xmm1
cmp    $0x3,%eax
mov    %ecx,%eax
jbe    112b <TIMESTE+0x2b>
shr    %ecx
mulpd  %xmm3,%xmm3
movapd %xmm4,%xmm5
test   $0x2,%al
je     1175 <TIMESTE+0x75>
movapd %xmm3,%xmm5
unpckhpd %xmm3,%xmm5
movapd %xmm4,%xmm6
je     1150 <TIMESTE+0x50>
movapd %xmm3,%xmm6
jmp    1150 <TIMESTE+0x50>
divsd  %xmm1,%xmm0
divsd  %xmm2,%xmm0
ret

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

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