double pdoublexpdouble(double *ai, double *y, unsigned long n) {
  double sum = 0;
  for (n -= 1; n > 0; n--)
    sum += *(ai + n) * *(y + n);
  return sum;
}


x86 assembly (clang, -O0):
<pdoublexpdouble>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x18(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x18(%rbp)
cmpq   $0x0,-0x18(%rbp)
jbe    116c <pdoublexpdouble+0x6c>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x20(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x18(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x18(%rbp)
jmp    1124 <pdoublexpdouble+0x24>
movsd  -0x20(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<pdoublexpdouble>:
mov    %rdx,%rax
add    $0xffffffffffffffff,%rax
je     1185 <pdoublexpdouble+0x85>
add    $0xfffffffffffffffe,%rdx
mov    %rax,%rcx
xorpd  %xmm0,%xmm0
and    $0x3,%rcx
je     1138 <pdoublexpdouble+0x38>
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rax,8),%xmm1
mulsd  (%rsi,%rax,8),%xmm1
addsd  %xmm1,%xmm0
add    $0xffffffffffffffff,%rax
add    $0xffffffffffffffff,%rcx
jne    1120 <pdoublexpdouble+0x20>
cmp    $0x3,%rdx
jb     1184 <pdoublexpdouble+0x84>
xchg   %ax,%ax
movsd  (%rdi,%rax,8),%xmm1
mulsd  (%rsi,%rax,8),%xmm1
addsd  %xmm0,%xmm1
movsd  -0x8(%rdi,%rax,8),%xmm2
mulsd  -0x8(%rsi,%rax,8),%xmm2
addsd  %xmm1,%xmm2
movsd  -0x18(%rdi,%rax,8),%xmm0
movsd  -0x10(%rdi,%rax,8),%xmm1
mulsd  -0x10(%rsi,%rax,8),%xmm1
addsd  %xmm2,%xmm1
mulsd  -0x18(%rsi,%rax,8),%xmm0
addsd  %xmm1,%xmm0
add    $0xfffffffffffffffc,%rax
jne    1140 <pdoublexpdouble+0x40>
ret
xorps  %xmm0,%xmm0
ret

LLVM IR (-O0):
define double @pdoublexpdouble(double* %0, double* %1, i64 %2) {
%4 = alloca double*
%5 = alloca double*
%6 = alloca i64
%7 = alloca double
store double* %0, double** %4
store double* %1, double** %5
store i64 %2, i64* %6
store double 0.000000e+00, double* %7
%8 = load i64, i64* %6
%9 = sub i64 %8, 1
store i64 %9, i64* %6
br label %10
10: ; preds = %24, %3
%11 = load i64, i64* %6
%12 = icmp ugt i64 %11, 0
br i1 %12, label %13, label %27
13: ; preds = %10
%14 = load double*, double** %4
%15 = load i64, i64* %6
%16 = getelementptr inbounds double, double* %14, i64 %15
%17 = load double, double* %16
%18 = load double*, double** %5
%19 = load i64, i64* %6
%20 = getelementptr inbounds double, double* %18, i64 %19
%21 = load double, double* %20
%22 = load double, double* %7
%23 = call double @llvm.fmuladd.f64(double %17, double %21, double %22)
store double %23, double* %7
br label %24
24: ; preds = %13
%25 = load i64, i64* %6
%26 = add i64 %25, -1
store i64 %26, i64* %6
br label %10
27: ; preds = %10
%28 = load double, double* %7
ret double %28
}

LLVM IR (-O2):
define double @pdoublexpdouble(double* nocapture %0, double* nocapture %1, i64 %2) {
%4 = add i64 %2, -1
%5 = icmp eq i64 %4, 0
br i1 %5, label %55, label %6
6: ; preds = %3
%7 = add i64 %2, -2
%8 = and i64 %4, 3
%9 = icmp eq i64 %8, 0
br i1 %9, label %22, label %10
10: ; preds = %6, %10
%11 = phi i64 [ %19, %10 ], [ %4, %6 ]
%12 = phi double [ %18, %10 ], [ 0.000000e+00, %6 ]
%13 = phi i64 [ %20, %10 ], [ 0, %6 ]
%14 = getelementptr inbounds double, double* %0, i64 %11
%15 = load double, double* %14
%16 = getelementptr inbounds double, double* %1, i64 %11
%17 = load double, double* %16
%18 = tail call double @llvm.fmuladd.f64(double %15, double %17, double %12)
%19 = add i64 %11, -1
%20 = add i64 %13, 1
%21 = icmp eq i64 %20, %8
br i1 %21, label %22, label %10
22: ; preds = %10, %6
%23 = phi double [ undef, %6 ], [ %18, %10 ]
%24 = phi i64 [ %4, %6 ], [ %19, %10 ]
%25 = phi double [ 0.000000e+00, %6 ], [ %18, %10 ]
%26 = icmp ult i64 %7, 3
br i1 %26, label %55, label %27
27: ; preds = %22, %27
%28 = phi i64 [ %53, %27 ], [ %24, %22 ]
%29 = phi double [ %52, %27 ], [ %25, %22 ]
%30 = getelementptr inbounds double, double* %0, i64 %28
%31 = load double, double* %30
%32 = getelementptr inbounds double, double* %1, i64 %28
%33 = load double, double* %32
%34 = tail call double @llvm.fmuladd.f64(double %31, double %33, double %29)
%35 = add i64 %28, -1
%36 = getelementptr inbounds double, double* %0, i64 %35
%37 = load double, double* %36
%38 = getelementptr inbounds double, double* %1, i64 %35
%39 = load double, double* %38
%40 = tail call double @llvm.fmuladd.f64(double %37, double %39, double %34)
%41 = add i64 %28, -2
%42 = getelementptr inbounds double, double* %0, i64 %41
%43 = load double, double* %42
%44 = getelementptr inbounds double, double* %1, i64 %41
%45 = load double, double* %44
%46 = tail call double @llvm.fmuladd.f64(double %43, double %45, double %40)
%47 = add i64 %28, -3
%48 = getelementptr inbounds double, double* %0, i64 %47
%49 = load double, double* %48
%50 = getelementptr inbounds double, double* %1, i64 %47
%51 = load double, double* %50
%52 = tail call double @llvm.fmuladd.f64(double %49, double %51, double %46)
%53 = add i64 %28, -4
%54 = icmp eq i64 %53, 0
br i1 %54, label %55, label %27
55: ; preds = %22, %27, %3
%56 = phi double [ 0.000000e+00, %3 ], [ %23, %22 ], [ %52, %27 ]
ret double %56
}