float dot(float *a, unsigned int a_length, float *history, float *b,
          unsigned int b_last_index) {
  float dotprod = 0.0f;
  unsigned int i = 0;
  if (a_length > b_last_index) {
    for (; i < (a_length - b_last_index - 1); i++)
      dotprod += a[i] * history[b_last_index + i];
  }
  for (; i < a_length; i++)
    dotprod += a[i] * b[(i + 1UL + b_last_index) - a_length];
  return dotprod;
}


x86 assembly (clang, -O0):
<dot>:
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)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x24(%rbp),%eax
jbe    118d <dot+0x8d>
jmp    1137 <dot+0x37>
mov    -0x2c(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    -0x24(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jae    1188 <dot+0x88>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x2c(%rbp),%ecx
mov    %ecx,%ecx
movss  (%rax,%rcx,4),%xmm2
movss  -0x28(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1137 <dot+0x37>
jmp    118d <dot+0x8d>
jmp    1192 <dot+0x92>
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    11e6 <dot+0xe6>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%rcx
mov    -0x24(%rbp),%edx
add    %rdx,%rcx
mov    -0xc(%rbp),%edx
sub    %rdx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x28(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1192 <dot+0x92>
movss  -0x28(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<dot>:
push   %r15
push   %r14
push   %rbx
xor    %ebx,%ebx
xorps  %xmm0,%xmm0
mov    %esi,%r9d
sub    %r8d,%r9d
jbe    1203 <dot+0x103>
mov    %r8d,%eax
not    %eax
add    %esi,%eax
je     1203 <dot+0x103>
mov    %eax,%r11d
lea    -0x1(%r11),%rax
mov    %r11d,%r10d
and    $0x3,%r10d
cmp    $0x3,%rax
jae    1142 <dot+0x42>
xorps  %xmm0,%xmm0
xor    %r15d,%r15d
jmp    11c5 <dot+0xc5>
and    $0xfffffffc,%r11d
neg    %r11
xorps  %xmm0,%xmm0
mov    $0x3,%r15d
mov    %r8d,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
mov    %ebx,%r14d
movss  -0xc(%rdi,%r15,4),%xmm1
mulss  (%rdx,%r14,4),%xmm1
movss  -0x8(%rdi,%r15,4),%xmm2
lea    0x1(%rbx),%eax
mulss  (%rdx,%rax,4),%xmm2
addss  %xmm0,%xmm1
addss  %xmm1,%xmm2
movss  -0x4(%rdi,%r15,4),%xmm1
lea    0x2(%rbx),%eax
mulss  (%rdx,%rax,4),%xmm1
addss  %xmm2,%xmm1
movss  (%rdi,%r15,4),%xmm0
lea    0x3(%rbx),%eax
mulss  (%rdx,%rax,4),%xmm0
addss  %xmm1,%xmm0
add    $0x4,%ebx
lea    (%r11,%r15,1),%rax
add    $0x4,%rax
add    $0x4,%r15
cmp    $0x3,%rax
jne    1160 <dot+0x60>
add    $0xfffffffffffffffd,%r15
add    $0xfffffffe,%r9d
test   %r10,%r10
je     11fc <dot+0xfc>
lea    (%r8,%r15,1),%r11d
lea    (%rdi,%r15,4),%r14
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
movss  (%r14,%rbx,4),%xmm1
lea    (%r11,%rbx,1),%eax
mulss  (%rdx,%rax,4),%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rbx
cmp    %rbx,%r10
jne    11e0 <dot+0xe0>
add    $0x1,%r9d
mov    %r9d,%ebx
cmp    %esi,%ebx
jae    12a7 <dot+0x1a7>
mov    %r8d,%r8d
mov    %esi,%eax
sub    %rax,%r8
mov    %ebx,%edx
sub    %ebx,%esi
mov    %rdx,%r9
not    %r9
add    %rax,%r9
and    $0x3,%rsi
je     1248 <dot+0x148>
lea    0x1(%r8),%rbx
lea    (%rcx,%rbx,4),%rbx
xchg   %ax,%ax
movss  (%rdi,%rdx,4),%xmm1
mulss  (%rbx,%rdx,4),%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rdx
add    $0xffffffffffffffff,%rsi
jne    1230 <dot+0x130>
cmp    $0x3,%r9
jb     12a7 <dot+0x1a7>
lea    (%rcx,%r8,4),%rcx
add    $0x10,%rcx
cs nopw 0x0(%rax,%rax,1)
movss  (%rdi,%rdx,4),%xmm1
movss  0x4(%rdi,%rdx,4),%xmm2
mulss  -0xc(%rcx,%rdx,4),%xmm1
mulss  -0x8(%rcx,%rdx,4),%xmm2
addss  %xmm0,%xmm1
movss  0x8(%rdi,%rdx,4),%xmm3
mulss  -0x4(%rcx,%rdx,4),%xmm3
addss  %xmm1,%xmm2
movss  0xc(%rdi,%rdx,4),%xmm0
mulss  (%rcx,%rdx,4),%xmm0
addss  %xmm2,%xmm3
addss  %xmm3,%xmm0
add    $0x4,%rdx
cmp    %rdx,%rax
jne    1260 <dot+0x160>
pop    %rbx
pop    %r14
pop    %r15
ret

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

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