void _inner_vec_ddot(double *y0, int incY, const double *a0, const double *x0,
                     int incX, double alpha, int nC) {
  register int i;
  register double ytmp;
  ytmp = 0.0;
  for (i = 0; i < nC - 3; i += 4) {
    ytmp += a0[0] * x0[0];
    ytmp += a0[1] * x0[incX];
    ytmp += a0[2] * x0[2 * incX];
    ytmp += a0[3] * x0[3 * incX];
    x0 += incX << 2;
    a0 += 4;
  }
  if (i == nC)
    goto update;
  if (i < nC - 1) {
    ytmp += a0[0] * x0[0];
    ytmp += a0[1] * x0[incX];
    x0 += incX << 1;
    a0 += 2;
    i += 2;
  }
  if (i < nC) {
    ytmp += a0[0] * x0[0];
    x0 += incX;
    i++;
  }
update:
  y0[0] += ytmp * alpha;
}


x86 assembly (clang, -O0):
<_inner_vec_ddot>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
movsd  %xmm0,-0x30(%rbp)
mov    %r9d,-0x34(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x34(%rbp),%ecx
sub    $0x3,%ecx
cmp    %ecx,%eax
jge    1214 <_inner_vec_ddot+0x114>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x20(%rbp),%rax
imul   $0x3,-0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x24(%rbp),%ecx
shl    $0x2,%ecx
mov    -0x20(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x18(%rbp),%rax
add    $0x20,%rax
mov    %rax,-0x18(%rbp)
mov    -0x38(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x38(%rbp)
jmp    112f <_inner_vec_ddot+0x2f>
mov    -0x38(%rbp),%eax
cmp    -0x34(%rbp),%eax
jne    1225 <_inner_vec_ddot+0x125>
jmp    12fe <_inner_vec_ddot+0x1fe>
mov    -0x38(%rbp),%eax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12ad <_inner_vec_ddot+0x1ad>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x24(%rbp),%ecx
shl    $0x1,%ecx
mov    -0x20(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x18(%rbp),%rax
add    $0x10,%rax
mov    %rax,-0x18(%rbp)
mov    -0x38(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    12f9 <_inner_vec_ddot+0x1f9>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x24(%rbp),%ecx
mov    -0x20(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    12fe <_inner_vec_ddot+0x1fe>
movsd  -0x40(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<_inner_vec_ddot>:
push   %rbp
push   %rbx
cmp    $0x4,%r9d
jl     1183 <_inner_vec_ddot+0x83>
lea    -0x3(%r9),%r10d
movslq %r8d,%rsi
lea    (%rsi,%rsi,1),%eax
movslq %eax,%r11
lea    (%rsi,%rsi,2),%eax
cltq
lea    0x0(,%rsi,4),%ebx
movslq %ebx,%rbx
shl    $0x3,%rbx
xor    %ebp,%ebp
xorpd  %xmm1,%xmm1
xchg   %ax,%ax
movsd  (%rdx),%xmm2
movsd  0x8(%rdx),%xmm3
mulsd  (%rcx),%xmm2
addsd  %xmm1,%xmm2
mulsd  (%rcx,%rsi,8),%xmm3
addsd  %xmm2,%xmm3
movsd  0x10(%rdx),%xmm2
mulsd  (%rcx,%r11,8),%xmm2
addsd  %xmm3,%xmm2
movsd  0x18(%rdx),%xmm1
mulsd  (%rcx,%rax,8),%xmm1
addsd  %xmm2,%xmm1
add    $0x20,%rdx
add    $0x4,%ebp
add    %rbx,%rcx
cmp    %r10d,%ebp
jl     1130 <_inner_vec_ddot+0x30>
mov    %r9d,%eax
and    $0xfffffffc,%eax
cmp    %r9d,%eax
jne    118e <_inner_vec_ddot+0x8e>
jmp    11d5 <_inner_vec_ddot+0xd5>
xorpd  %xmm1,%xmm1
xor    %eax,%eax
cmp    %r9d,%eax
je     11d5 <_inner_vec_ddot+0xd5>
lea    -0x1(%r9),%esi
cmp    %esi,%eax
jge    11c4 <_inner_vec_ddot+0xc4>
movsd  (%rdx),%xmm2
mulsd  (%rcx),%xmm2
addsd  %xmm1,%xmm2
movsd  0x8(%rdx),%xmm1
movslq %r8d,%rsi
mulsd  (%rcx,%rsi,8),%xmm1
addsd  %xmm2,%xmm1
add    %r8d,%r8d
movslq %r8d,%rsi
lea    (%rcx,%rsi,8),%rcx
add    $0x10,%rdx
or     $0x2,%eax
cmp    %r9d,%eax
jge    11d5 <_inner_vec_ddot+0xd5>
movsd  (%rdx),%xmm2
mulsd  (%rcx),%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdi),%xmm1
movsd  %xmm1,(%rdi)
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
define void @_inner_vec_ddot(double* %0, i32 %1, double* %2, double* %3, i32 %4, double %5, i32 %6) {
%8 = alloca double*
%9 = alloca i32
%10 = alloca double*
%11 = alloca double*
%12 = alloca i32
%13 = alloca double
%14 = alloca i32
%15 = alloca i32
%16 = alloca double
store double* %0, double** %8
store i32 %1, i32* %9
store double* %2, double** %10
store double* %3, double** %11
store i32 %4, i32* %12
store double %5, double* %13
store i32 %6, i32* %14
store double 0.000000e+00, double* %16
store i32 0, i32* %15
br label %17
17: ; preds = %70, %7
%18 = load i32, i32* %15
%19 = load i32, i32* %14
%20 = sub nsw i32 %19, 3
%21 = icmp slt i32 %18, %20
br i1 %21, label %22, label %73
22: ; preds = %17
%23 = load double*, double** %10
%24 = getelementptr inbounds double, double* %23, i64 0
%25 = load double, double* %24
%26 = load double*, double** %11
%27 = getelementptr inbounds double, double* %26, i64 0
%28 = load double, double* %27
%29 = load double, double* %16
%30 = call double @llvm.fmuladd.f64(double %25, double %28, double %29)
store double %30, double* %16
%31 = load double*, double** %10
%32 = getelementptr inbounds double, double* %31, i64 1
%33 = load double, double* %32
%34 = load double*, double** %11
%35 = load i32, i32* %12
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %34, i64 %36
%38 = load double, double* %37
%39 = load double, double* %16
%40 = call double @llvm.fmuladd.f64(double %33, double %38, double %39)
store double %40, double* %16
%41 = load double*, double** %10
%42 = getelementptr inbounds double, double* %41, i64 2
%43 = load double, double* %42
%44 = load double*, double** %11
%45 = load i32, i32* %12
%46 = mul nsw i32 2, %45
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %44, i64 %47
%49 = load double, double* %48
%50 = load double, double* %16
%51 = call double @llvm.fmuladd.f64(double %43, double %49, double %50)
store double %51, double* %16
%52 = load double*, double** %10
%53 = getelementptr inbounds double, double* %52, i64 3
%54 = load double, double* %53
%55 = load double*, double** %11
%56 = load i32, i32* %12
%57 = mul nsw i32 3, %56
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %55, i64 %58
%60 = load double, double* %59
%61 = load double, double* %16
%62 = call double @llvm.fmuladd.f64(double %54, double %60, double %61)
store double %62, double* %16
%63 = load i32, i32* %12
%64 = shl i32 %63, 2
%65 = load double*, double** %11
%66 = sext i32 %64 to i64
%67 = getelementptr inbounds double, double* %65, i64 %66
store double* %67, double** %11
%68 = load double*, double** %10
%69 = getelementptr inbounds double, double* %68, i64 4
store double* %69, double** %10
br label %70
70: ; preds = %22
%71 = load i32, i32* %15
%72 = add nsw i32 %71, 4
store i32 %72, i32* %15
br label %17
73: ; preds = %17
%74 = load i32, i32* %15
%75 = load i32, i32* %14
%76 = icmp eq i32 %74, %75
br i1 %76, label %77, label %78
77: ; preds = %73
br label %131
78: ; preds = %73
%79 = load i32, i32* %15
%80 = load i32, i32* %14
%81 = sub nsw i32 %80, 1
%82 = icmp slt i32 %79, %81
br i1 %82, label %83, label %111
83: ; preds = %78
%84 = load double*, double** %10
%85 = getelementptr inbounds double, double* %84, i64 0
%86 = load double, double* %85
%87 = load double*, double** %11
%88 = getelementptr inbounds double, double* %87, i64 0
%89 = load double, double* %88
%90 = load double, double* %16
%91 = call double @llvm.fmuladd.f64(double %86, double %89, double %90)
store double %91, double* %16
%92 = load double*, double** %10
%93 = getelementptr inbounds double, double* %92, i64 1
%94 = load double, double* %93
%95 = load double*, double** %11
%96 = load i32, i32* %12
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds double, double* %95, i64 %97
%99 = load double, double* %98
%100 = load double, double* %16
%101 = call double @llvm.fmuladd.f64(double %94, double %99, double %100)
store double %101, double* %16
%102 = load i32, i32* %12
%103 = shl i32 %102, 1
%104 = load double*, double** %11
%105 = sext i32 %103 to i64
%106 = getelementptr inbounds double, double* %104, i64 %105
store double* %106, double** %11
%107 = load double*, double** %10
%108 = getelementptr inbounds double, double* %107, i64 2
store double* %108, double** %10
%109 = load i32, i32* %15
%110 = add nsw i32 %109, 2
store i32 %110, i32* %15
br label %111
111: ; preds = %83, %78
%112 = load i32, i32* %15
%113 = load i32, i32* %14
%114 = icmp slt i32 %112, %113
br i1 %114, label %115, label %130
115: ; preds = %111
%116 = load double*, double** %10
%117 = getelementptr inbounds double, double* %116, i64 0
%118 = load double, double* %117
%119 = load double*, double** %11
%120 = getelementptr inbounds double, double* %119, i64 0
%121 = load double, double* %120
%122 = load double, double* %16
%123 = call double @llvm.fmuladd.f64(double %118, double %121, double %122)
store double %123, double* %16
%124 = load i32, i32* %12
%125 = load double*, double** %11
%126 = sext i32 %124 to i64
%127 = getelementptr inbounds double, double* %125, i64 %126
store double* %127, double** %11
%128 = load i32, i32* %15
%129 = add nsw i32 %128, 1
store i32 %129, i32* %15
br label %130
130: ; preds = %115, %111
br label %131
131: ; preds = %130, %77
%132 = load double, double* %16
%133 = load double, double* %13
%134 = load double*, double** %8
%135 = getelementptr inbounds double, double* %134, i64 0
%136 = load double, double* %135
%137 = call double @llvm.fmuladd.f64(double %132, double %133, double %136)
store double %137, double* %135
ret void
}

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