double ddot_vec(const double *X, const double *Y, int incX, int incY, int N) {
  register int i;
  register double c0, c1;
  c0 = 0.0;
  c1 = 0.0;
  for (i = 0; i < N - 1; i += 2) {
    c0 += X[0] * Y[0];
    X += incX;
    Y += incY;
    c1 += X[0] * Y[0];
    X += incX;
    Y += incY;
  }
  if (i == N) {
    return c0 + c1;
  }
  c0 += X[0] * Y[0];
  return c0 + c1;
}


x86 assembly (clang, -O0):
<ddot_vec>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %ecx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    11e4 <ddot_vec+0xe4>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x1c(%rbp),%ecx
mov    -0x10(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x10(%rbp)
mov    -0x20(%rbp),%ecx
mov    -0x18(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x1c(%rbp),%ecx
mov    -0x10(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x10(%rbp)
mov    -0x20(%rbp),%ecx
mov    -0x18(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x28(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x28(%rbp)
jmp    112d <ddot_vec+0x2d>
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jne    1204 <ddot_vec+0x104>
movsd  -0x30(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1235 <ddot_vec+0x135>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<ddot_vec>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
cmp    $0x2,%r8d
jl     113c <ddot_vec+0x3c>
movslq %edx,%r9
movslq %ecx,%r10
lea    -0x2(%r8),%eax
mov    %eax,%r15d
shr    %r15d
add    $0x1,%r15d
cmp    $0x2,%eax
jae    1147 <ddot_vec+0x47>
xorpd  %xmm1,%xmm1
test   $0x1,%r15b
jne    11e9 <ddot_vec+0xe9>
jmp    1211 <ddot_vec+0x111>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
jmp    121f <ddot_vec+0x11f>
mov    %r15d,%r12d
and    $0xfffffffe,%r12d
mov    %r9,%r11
shl    $0x5,%r11
mov    %r10,%r14
shl    $0x5,%r14
lea    0x0(,%r9,8),%rbx
lea    0x0(,%r10,8),%rbp
xorpd  %xmm1,%xmm1
mov    %rdi,%rcx
mov    %rsi,%rax
lea    (%rdi,%rbx,1),%rdx
lea    (%rsi,%rbp,1),%r13
movsd  (%rdi),%xmm0
movhpd (%rdi,%rbx,1),%xmm0
movsd  (%rsi),%xmm2
movhpd (%rsi,%rbp,1),%xmm2
mulpd  %xmm0,%xmm2
addpd  %xmm1,%xmm2
lea    (%rdx,%rbx,1),%rdi
mov    %r13,%rsi
add    %rbp,%rsi
movsd  (%rbx,%rdx,1),%xmm0
movhpd (%rbx,%rdi,1),%xmm0
add    %rbx,%rdi
movsd  0x0(%rbp,%r13,1),%xmm3
movhpd 0x0(%rbp,%rsi,1),%xmm3
add    %rbp,%rsi
movapd %xmm0,%xmm1
mulpd  %xmm3,%xmm1
addpd  %xmm2,%xmm1
add    %rbx,%rdi
add    %rbp,%rsi
add    $0xfffffffe,%r12d
jne    1170 <ddot_vec+0x70>
add    %r11,%rcx
add    %r14,%rax
mov    %rcx,%rdi
mov    %rax,%rsi
test   $0x1,%r15b
je     1211 <ddot_vec+0x111>
lea    (%rdi,%r9,8),%rax
lea    (%rsi,%r10,8),%rcx
movsd  (%rdi),%xmm0
movhpd (%rdi,%r9,8),%xmm0
movsd  (%rsi),%xmm3
movhpd (%rsi,%r10,8),%xmm3
lea    (%rax,%r9,8),%rdi
lea    (%rcx,%r10,8),%rsi
movapd %xmm1,%xmm2
mulpd  %xmm3,%xmm0
addpd  %xmm2,%xmm0
mov    %r8d,%eax
and    $0xfffffffe,%eax
movapd %xmm0,%xmm1
cmp    %r8d,%eax
je     1234 <ddot_vec+0x134>
movsd  (%rdi),%xmm1
mulsd  (%rsi),%xmm1
addsd  %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @ddot_vec(double* %0, double* %1, i32 %2, i32 %3, i32 %4) {
%6 = alloca double
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
%14 = alloca double
store double* %0, double** %7
store double* %1, double** %8
store i32 %2, i32* %9
store i32 %3, i32* %10
store i32 %4, i32* %11
store double 0.000000e+00, double* %13
store double 0.000000e+00, double* %14
store i32 0, i32* %12
br label %15
15: ; preds = %53, %5
%16 = load i32, i32* %12
%17 = load i32, i32* %11
%18 = sub nsw i32 %17, 1
%19 = icmp slt i32 %16, %18
br i1 %19, label %20, label %56
20: ; preds = %15
%21 = load double*, double** %7
%22 = getelementptr inbounds double, double* %21, i64 0
%23 = load double, double* %22
%24 = load double*, double** %8
%25 = getelementptr inbounds double, double* %24, i64 0
%26 = load double, double* %25
%27 = load double, double* %13
%28 = call double @llvm.fmuladd.f64(double %23, double %26, double %27)
store double %28, double* %13
%29 = load i32, i32* %9
%30 = load double*, double** %7
%31 = sext i32 %29 to i64
%32 = getelementptr inbounds double, double* %30, i64 %31
store double* %32, double** %7
%33 = load i32, i32* %10
%34 = load double*, double** %8
%35 = sext i32 %33 to i64
%36 = getelementptr inbounds double, double* %34, i64 %35
store double* %36, double** %8
%37 = load double*, double** %7
%38 = getelementptr inbounds double, double* %37, i64 0
%39 = load double, double* %38
%40 = load double*, double** %8
%41 = getelementptr inbounds double, double* %40, i64 0
%42 = load double, double* %41
%43 = load double, double* %14
%44 = call double @llvm.fmuladd.f64(double %39, double %42, double %43)
store double %44, double* %14
%45 = load i32, i32* %9
%46 = load double*, double** %7
%47 = sext i32 %45 to i64
%48 = getelementptr inbounds double, double* %46, i64 %47
store double* %48, double** %7
%49 = load i32, i32* %10
%50 = load double*, double** %8
%51 = sext i32 %49 to i64
%52 = getelementptr inbounds double, double* %50, i64 %51
store double* %52, double** %8
br label %53
53: ; preds = %20
%54 = load i32, i32* %12
%55 = add nsw i32 %54, 2
store i32 %55, i32* %12
br label %15
56: ; preds = %15
%57 = load i32, i32* %12
%58 = load i32, i32* %11
%59 = icmp eq i32 %57, %58
br i1 %59, label %60, label %64
60: ; preds = %56
%61 = load double, double* %13
%62 = load double, double* %14
%63 = fadd double %61, %62
store double %63, double* %6
br label %76
64: ; preds = %56
%65 = load double*, double** %7
%66 = getelementptr inbounds double, double* %65, i64 0
%67 = load double, double* %66
%68 = load double*, double** %8
%69 = getelementptr inbounds double, double* %68, i64 0
%70 = load double, double* %69
%71 = load double, double* %13
%72 = call double @llvm.fmuladd.f64(double %67, double %70, double %71)
store double %72, double* %13
%73 = load double, double* %13
%74 = load double, double* %14
%75 = fadd double %73, %74
store double %75, double* %6
br label %76
76: ; preds = %64, %60
%77 = load double, double* %6
ret double %77
}

LLVM IR (-O2):
define double @ddot_vec(double* nocapture %0, double* nocapture %1, i32 %2, i32 %3, i32 %4) {
%6 = icmp sgt i32 %4, 1
br i1 %6, label %7, label %81
7: ; preds = %5
%8 = sext i32 %2 to i64
%9 = sext i32 %3 to i64
%10 = add i32 %4, -2
%11 = lshr i32 %10, 1
%12 = add nuw i32 %11, 1
%13 = and i32 %12, 1
%14 = icmp ult i32 %10, 2
br i1 %14, label %50, label %15
15: ; preds = %7
%16 = and i32 %12, -2
br label %17
17: ; preds = %17, %15
%18 = phi double* [ %0, %15 ], [ %46, %17 ]
%19 = phi double* [ %1, %15 ], [ %47, %17 ]
%20 = phi <2 x double> [ zeroinitializer, %15 ], [ %45, %17 ]
%21 = phi i32 [ 0, %15 ], [ %48, %17 ]
%22 = load double, double* %18
%23 = load double, double* %19
%24 = getelementptr inbounds double, double* %18, i64 %8
%25 = getelementptr inbounds double, double* %19, i64 %9
%26 = load double, double* %24
%27 = load double, double* %25
%28 = insertelement <2 x double> poison, double %22, i64 0
%29 = insertelement <2 x double> %28, double %26, i64 1
%30 = insertelement <2 x double> poison, double %23, i64 0
%31 = insertelement <2 x double> %30, double %27, i64 1
%32 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %29, <2 x double> %31, <2 x double> %20)
%33 = getelementptr inbounds double, double* %24, i64 %8
%34 = getelementptr inbounds double, double* %25, i64 %9
%35 = load double, double* %33
%36 = load double, double* %34
%37 = getelementptr inbounds double, double* %33, i64 %8
%38 = getelementptr inbounds double, double* %34, i64 %9
%39 = load double, double* %37
%40 = load double, double* %38
%41 = insertelement <2 x double> poison, double %35, i64 0
%42 = insertelement <2 x double> %41, double %39, i64 1
%43 = insertelement <2 x double> poison, double %36, i64 0
%44 = insertelement <2 x double> %43, double %40, i64 1
%45 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %42, <2 x double> %44, <2 x double> %32)
%46 = getelementptr inbounds double, double* %37, i64 %8
%47 = getelementptr inbounds double, double* %38, i64 %9
%48 = add nuw nsw i32 %21, 2
%49 = icmp eq i32 %48, %16
br i1 %49, label %50, label %17
50: ; preds = %17, %7
%51 = phi <2 x double> [ undef, %7 ], [ %32, %17 ]
%52 = phi <2 x double> [ undef, %7 ], [ %42, %17 ]
%53 = phi <2 x double> [ undef, %7 ], [ %44, %17 ]
%54 = phi double* [ undef, %7 ], [ %46, %17 ]
%55 = phi double* [ undef, %7 ], [ %47, %17 ]
%56 = phi double* [ %0, %7 ], [ %46, %17 ]
%57 = phi double* [ %1, %7 ], [ %47, %17 ]
%58 = phi <2 x double> [ zeroinitializer, %7 ], [ %45, %17 ]
%59 = icmp eq i32 %13, 0
br i1 %59, label %73, label %60
60: ; preds = %50
%61 = load double, double* %56
%62 = load double, double* %57
%63 = getelementptr inbounds double, double* %56, i64 %8
%64 = getelementptr inbounds double, double* %57, i64 %9
%65 = load double, double* %63
%66 = load double, double* %64
%67 = insertelement <2 x double> poison, double %61, i64 0
%68 = insertelement <2 x double> %67, double %65, i64 1
%69 = insertelement <2 x double> poison, double %62, i64 0
%70 = insertelement <2 x double> %69, double %66, i64 1
%71 = getelementptr inbounds double, double* %63, i64 %8
%72 = getelementptr inbounds double, double* %64, i64 %9
br label %73
73: ; preds = %50, %60
%74 = phi <2 x double> [ %51, %50 ], [ %58, %60 ]
%75 = phi <2 x double> [ %52, %50 ], [ %68, %60 ]
%76 = phi <2 x double> [ %53, %50 ], [ %70, %60 ]
%77 = phi double* [ %54, %50 ], [ %71, %60 ]
%78 = phi double* [ %55, %50 ], [ %72, %60 ]
%79 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %75, <2 x double> %76, <2 x double> %74)
%80 = and i32 %4, -2
br label %81
81: ; preds = %73, %5
%82 = phi double* [ %1, %5 ], [ %78, %73 ]
%83 = phi double* [ %0, %5 ], [ %77, %73 ]
%84 = phi i32 [ 0, %5 ], [ %80, %73 ]
%85 = phi <2 x double> [ zeroinitializer, %5 ], [ %79, %73 ]
%86 = icmp eq i32 %84, %4
%87 = extractelement <2 x double> %85, i64 0
br i1 %86, label %92, label %88
88: ; preds = %81
%89 = load double, double* %83
%90 = load double, double* %82
%91 = tail call double @llvm.fmuladd.f64(double %89, double %90, double %87)
br label %92
92: ; preds = %81, %88
%93 = phi double [ %91, %88 ], [ %87, %81 ]
%94 = extractelement <2 x double> %85, i64 1
%95 = fadd double %94, %93
ret double %95
}