#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */
typedef int /*<<< orphan*/ __m128;

/* Variables and functions */
int /*<<< orphan*/ _mm_add_ps(int /*<<< orphan*/, int /*<<< orphan*/);
int /*<<< orphan*/ _mm_add_ss(int /*<<< orphan*/, int /*<<< orphan*/);
scalar_t__ _mm_cvtss_f32(int /*<<< orphan*/);
int /*<<< orphan*/ _mm_load_ps(float const *);
int /*<<< orphan*/ _mm_loadu_ps(float const *);
int /*<<< orphan*/ _mm_movehl_ps(int /*<<< orphan*/, int /*<<< orphan*/);
int /*<<< orphan*/ _mm_mul_ps(int /*<<< orphan*/, int /*<<< orphan*/);
int /*<<< orphan*/ _mm_setzero_ps();
int /*<<< orphan*/ _mm_shuffle_ps(int /*<<< orphan*/, int /*<<< orphan*/, int);
float sdot_sse2(const float *a, const float *b, int n) {
  float s = 0.0f;
  int aligned = ((((uintptr_t)(const void *)(a)) % (16)) == 0) &&
                ((((uintptr_t)(const void *)(b)) % (16)) == 0);
  if (!aligned && ((((uintptr_t)(const void *)(a)) % (16)) ==
                   (((uintptr_t)(const void *)(b)) % (16)))) {
    int k = 0;
    while (!aligned) {
      s += (*a) * (*b);
      n--;
      a++;
      b++;
      aligned = ((((uintptr_t)(const void *)(a)) % (16)) == 0) &&
                ((((uintptr_t)(const void *)(b)) % (16)) == 0);
      if (aligned || (++k > 2) || (n == 0))
        break;
    }
  }
  __m128 s4 = _mm_setzero_ps();
  if (aligned)
    for (int k = 0, nq = 4 * (n / 4); k < nq; k += 4)
      s4 = _mm_add_ps(s4, _mm_mul_ps(_mm_load_ps(a + k), _mm_load_ps(b + k)));
  else
    for (int k = 0, nq = 4 * (n / 4); k < nq; k += 4)
      s4 = _mm_add_ps(s4, _mm_mul_ps(_mm_loadu_ps(a + k), _mm_loadu_ps(b + k)));
  s4 = _mm_add_ps(s4, _mm_movehl_ps(s4, s4));
  s4 = _mm_add_ss(s4, _mm_shuffle_ps(s4, s4, 1));
  s += _mm_cvtss_f32(s4);
  for (int k = 4 * (n / 4); k < n; k++)
    s += a[k] * b[k];
  return s;
}


x86 assembly (clang, -O0):
<sdot_sse2>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0x8(%rbp),%rcx
and    $0xf,%rcx
xor    %eax,%eax
cmp    $0x0,%rcx
mov    %al,-0x39(%rbp)
jne    11d4 <sdot_sse2+0x44>
mov    -0x10(%rbp),%rax
and    $0xf,%rax
cmp    $0x0,%rax
sete   %al
mov    %al,-0x39(%rbp)
mov    -0x39(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x0,-0x1c(%rbp)
jne    12cb <sdot_sse2+0x13b>
mov    -0x8(%rbp),%rax
and    $0xf,%rax
mov    -0x10(%rbp),%rcx
and    $0xf,%rcx
cmp    %rcx,%rax
jne    12cb <sdot_sse2+0x13b>
movl   $0x0,-0x20(%rbp)
cmpl   $0x0,-0x1c(%rbp)
setne  %al
xor    $0xff,%al
test   $0x1,%al
jne    121f <sdot_sse2+0x8f>
jmp    12c6 <sdot_sse2+0x136>
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm2
movss  -0x18(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
mov    -0x8(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x8(%rbp)
mov    -0x10(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x10(%rbp)
mov    -0x8(%rbp),%rcx
and    $0xf,%rcx
xor    %eax,%eax
cmp    $0x0,%rcx
mov    %al,-0x3a(%rbp)
jne    128b <sdot_sse2+0xfb>
mov    -0x10(%rbp),%rax
and    $0xf,%rax
cmp    $0x0,%rax
sete   %al
mov    %al,-0x3a(%rbp)
mov    -0x3a(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x0,-0x1c(%rbp)
jne    12bc <sdot_sse2+0x12c>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
cmp    $0x2,%eax
jg     12bc <sdot_sse2+0x12c>
cmpl   $0x0,-0x14(%rbp)
jne    12c1 <sdot_sse2+0x131>
jmp    12c6 <sdot_sse2+0x136>
jmp    1209 <sdot_sse2+0x79>
jmp    12cb <sdot_sse2+0x13b>
mov    $0x0,%al
call   10b0 <_mm_setzero_ps@plt>
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x1c(%rbp)
je     135e <sdot_sse2+0x1ce>
movl   $0x0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
shl    $0x2,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    1359 <sdot_sse2+0x1c9>
mov    -0x24(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x8(%rbp),%rdi
movslq -0x28(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
call   1040 <_mm_load_ps@plt>
mov    %eax,-0x44(%rbp)
mov    -0x10(%rbp),%rdi
movslq -0x28(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
call   1040 <_mm_load_ps@plt>
mov    -0x44(%rbp),%edi
mov    %eax,%esi
call   1080 <_mm_mul_ps@plt>
mov    -0x40(%rbp),%edi
mov    %eax,%esi
call   1090 <_mm_add_ps@plt>
mov    %eax,-0x24(%rbp)
mov    -0x28(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x28(%rbp)
jmp    12f7 <sdot_sse2+0x167>
jmp    13dd <sdot_sse2+0x24d>
movl   $0x0,-0x30(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
shl    $0x2,%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    13d8 <sdot_sse2+0x248>
mov    -0x24(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x8(%rbp),%rdi
movslq -0x30(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
call   1060 <_mm_loadu_ps@plt>
mov    %eax,-0x4c(%rbp)
mov    -0x10(%rbp),%rdi
movslq -0x30(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
call   1060 <_mm_loadu_ps@plt>
mov    -0x4c(%rbp),%edi
mov    %eax,%esi
call   1080 <_mm_mul_ps@plt>
mov    -0x48(%rbp),%edi
mov    %eax,%esi
call   1090 <_mm_add_ps@plt>
mov    %eax,-0x24(%rbp)
mov    -0x30(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x30(%rbp)
jmp    1376 <sdot_sse2+0x1e6>
jmp    13dd <sdot_sse2+0x24d>
mov    -0x24(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x24(%rbp),%edi
mov    -0x24(%rbp),%esi
call   10a0 <_mm_movehl_ps@plt>
mov    -0x54(%rbp),%edi
mov    %eax,%esi
call   1090 <_mm_add_ps@plt>
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x24(%rbp),%edi
mov    -0x24(%rbp),%esi
mov    $0x1,%edx
call   1050 <_mm_shuffle_ps@plt>
mov    -0x50(%rbp),%edi
mov    %eax,%esi
call   1070 <_mm_add_ss@plt>
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%edi
call   1030 <_mm_cvtss_f32@plt>
cvtsi2ss %rax,%xmm0
addss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
shl    $0x2,%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    148c <sdot_sse2+0x2fc>
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x18(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1446 <sdot_sse2+0x2b6>
movss  -0x18(%rbp),%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sdot_sse2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %edx,%r14d
mov    %rsi,%rbp
mov    %rdi,%r15
mov    %rsi,%rbx
or     %rdi,%rbx
test   $0xf,%bl
je     1256 <sdot_sse2+0xc6>
mov    %ebp,%eax
xor    %r15d,%eax
and    $0xf,%eax
test   %rax,%rax
jne    1256 <sdot_sse2+0xc6>
lea    -0x1(%r14),%ebx
cmp    $0x2,%ebx
mov    $0x2,%eax
cmovb  %ebx,%eax
movss  (%r15),%xmm0
mulss  0x0(%rbp),%xmm0
xorps  %xmm1,%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x8(%rsp)
mov    %r15,%rdx
add    $0x4,%r15
mov    %rbp,%rsi
add    $0x4,%rbp
mov    %r15d,%ecx
or     %ebp,%ecx
test   $0xf,%cl
je     1247 <sdot_sse2+0xb7>
test   %eax,%eax
je     1331 <sdot_sse2+0x1a1>
mov    %rdx,%rcx
movss  0x4(%rdx),%xmm0
mov    %rsi,%rdi
mulss  0x4(%rsi),%xmm0
movss  0x8(%rsp),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x8(%rsp)
lea    -0x2(%r14),%ebx
lea    0x8(%rdx),%r15
lea    0x8(%rsi),%rbp
mov    %r15d,%ecx
or     %ebp,%ecx
test   $0xf,%cl
jne    12f0 <sdot_sse2+0x160>
xor    %eax,%eax
call   10b0 <_mm_setzero_ps@plt>
mov    %eax,%r13d
mov    %ebx,%r14d
jmp    1272 <sdot_sse2+0xe2>
xor    %eax,%eax
call   10b0 <_mm_setzero_ps@plt>
mov    %eax,%r13d
xorps  %xmm0,%xmm0
movss  %xmm0,0x8(%rsp)
test   $0xf,%bl
jne    133e <sdot_sse2+0x1ae>
lea    0x3(%r14),%r12d
test   %r14d,%r14d
cmovns %r14d,%r12d
and    $0xfffffffc,%r12d
cmp    $0x4,%r14d
jl     13bf <sdot_sse2+0x22f>
mov    %r14,0x10(%rsp)
mov    %r12d,0xc(%rsp)
movslq %r12d,%r14
xor    %r12d,%r12d
mov    %r15,0x20(%rsp)
mov    %r15,%rbx
mov    %rbp,0x18(%rsp)
nopl   0x0(%rax,%rax,1)
mov    %rbx,%rdi
call   1040 <_mm_load_ps@plt>
mov    %eax,%r15d
mov    %rbp,%rdi
call   1040 <_mm_load_ps@plt>
mov    %r15d,%edi
mov    %eax,%esi
call   1080 <_mm_mul_ps@plt>
mov    %r13d,%edi
mov    %eax,%esi
call   1090 <_mm_add_ps@plt>
mov    %eax,%r13d
add    $0x4,%r12
add    $0x10,%rbp
add    $0x10,%rbx
cmp    %r14,%r12
jl     12b0 <sdot_sse2+0x120>
jmp    13ab <sdot_sse2+0x21b>
cmp    $0x1,%eax
je     1331 <sdot_sse2+0x1a1>
mov    %rdx,%r15
movss  0x8(%rdx),%xmm0
mov    %rsi,%rbp
mulss  0x8(%rsi),%xmm0
movss  0x8(%rsp),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x8(%rsp)
add    $0xfffffffd,%r14d
add    $0xc,%r15
add    $0xc,%rbp
mov    %r15d,%eax
or     %ebp,%eax
mov    %r14d,%ebx
test   $0xf,%al
je     1247 <sdot_sse2+0xb7>
xor    %eax,%eax
call   10b0 <_mm_setzero_ps@plt>
mov    %eax,%r13d
mov    %ebx,%r14d
lea    0x3(%r14),%r12d
test   %r14d,%r14d
cmovns %r14d,%r12d
and    $0xfffffffc,%r12d
cmp    $0x4,%r14d
jl     13bf <sdot_sse2+0x22f>
mov    %r14,0x10(%rsp)
mov    %r12d,0xc(%rsp)
movslq %r12d,%r14
xor    %r12d,%r12d
mov    %r15,0x20(%rsp)
mov    %r15,%rbx
mov    %rbp,0x18(%rsp)
mov    %rbx,%rdi
call   1060 <_mm_loadu_ps@plt>
mov    %eax,%r15d
mov    %rbp,%rdi
call   1060 <_mm_loadu_ps@plt>
mov    %r15d,%edi
mov    %eax,%esi
call   1080 <_mm_mul_ps@plt>
mov    %r13d,%edi
mov    %eax,%esi
call   1090 <_mm_add_ps@plt>
mov    %eax,%r13d
add    $0x4,%r12
add    $0x10,%rbp
add    $0x10,%rbx
cmp    %r14,%r12
jl     1370 <sdot_sse2+0x1e0>
mov    0x20(%rsp),%r15
mov    0x18(%rsp),%rbp
mov    0x10(%rsp),%r14
mov    0xc(%rsp),%r12d
mov    %r13d,%edi
mov    %r13d,%esi
call   10a0 <_mm_movehl_ps@plt>
mov    %r13d,%edi
mov    %eax,%esi
call   1090 <_mm_add_ps@plt>
mov    %eax,%ebx
mov    %eax,%edi
mov    %eax,%esi
mov    $0x1,%edx
call   1050 <_mm_shuffle_ps@plt>
mov    %ebx,%edi
mov    %eax,%esi
call   1070 <_mm_add_ss@plt>
mov    %eax,%edi
call   1030 <_mm_cvtss_f32@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
addss  0x8(%rsp),%xmm0
cmp    %r14d,%r12d
jge    149c <sdot_sse2+0x30c>
movslq %r12d,%rax
movslq %r14d,%rcx
mov    %rax,%rdx
not    %rdx
add    %rcx,%rdx
mov    %rcx,%rsi
and    $0x3,%rsi
je     144a <sdot_sse2+0x2ba>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  (%r15,%rax,4),%xmm1
mulss  0x0(%rbp,%rax,4),%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rax
add    $0xffffffffffffffff,%rsi
jne    1430 <sdot_sse2+0x2a0>
cmp    $0x3,%rdx
jb     149c <sdot_sse2+0x30c>
movss  (%r15,%rax,4),%xmm1
movss  0x4(%r15,%rax,4),%xmm2
mulss  0x0(%rbp,%rax,4),%xmm1
mulss  0x4(%rbp,%rax,4),%xmm2
addss  %xmm0,%xmm1
movss  0x8(%r15,%rax,4),%xmm3
mulss  0x8(%rbp,%rax,4),%xmm3
addss  %xmm1,%xmm2
movss  0xc(%r15,%rax,4),%xmm0
mulss  0xc(%rbp,%rax,4),%xmm0
addss  %xmm2,%xmm3
addss  %xmm3,%xmm0
add    $0x4,%rax
cmp    %rax,%rcx
jne    1450 <sdot_sse2+0x2c0>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define float @sdot_sse2(float* %0, float* %1, i32 %2) {
%4 = alloca float*
%5 = alloca float*
%6 = alloca i32
%7 = alloca float
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
store float* %0, float** %4
store float* %1, float** %5
store i32 %2, i32* %6
store float 0.000000e+00, float* %7
%16 = load float*, float** %4
%17 = bitcast float* %16 to i8*
%18 = ptrtoint i8* %17 to i64
%19 = urem i64 %18, 16
%20 = icmp eq i64 %19, 0
br i1 %20, label %21, label %27
21: ; preds = %3
%22 = load float*, float** %5
%23 = bitcast float* %22 to i8*
%24 = ptrtoint i8* %23 to i64
%25 = urem i64 %24, 16
%26 = icmp eq i64 %25, 0
br label %27
27: ; preds = %21, %3
%28 = phi i1 [ false, %3 ], [ %26, %21 ]
%29 = zext i1 %28 to i32
store i32 %29, i32* %8
%30 = load i32, i32* %8
%31 = icmp ne i32 %30, 0
br i1 %31, label %86, label %32
32: ; preds = %27
%33 = load float*, float** %4
%34 = bitcast float* %33 to i8*
%35 = ptrtoint i8* %34 to i64
%36 = urem i64 %35, 16
%37 = load float*, float** %5
%38 = bitcast float* %37 to i8*
%39 = ptrtoint i8* %38 to i64
%40 = urem i64 %39, 16
%41 = icmp eq i64 %36, %40
br i1 %41, label %42, label %86
42: ; preds = %32
store i32 0, i32* %9
br label %43
43: ; preds = %84, %42
%44 = load i32, i32* %8
%45 = icmp ne i32 %44, 0
%46 = xor i1 %45, true
br i1 %46, label %47, label %85
47: ; preds = %43
%48 = load float*, float** %4
%49 = load float, float* %48
%50 = load float*, float** %5
%51 = load float, float* %50
%52 = load float, float* %7
%53 = call float @llvm.fmuladd.f32(float %49, float %51, float %52)
store float %53, float* %7
%54 = load i32, i32* %6
%55 = add nsw i32 %54, -1
store i32 %55, i32* %6
%56 = load float*, float** %4
%57 = getelementptr inbounds float, float* %56, i32 1
store float* %57, float** %4
%58 = load float*, float** %5
%59 = getelementptr inbounds float, float* %58, i32 1
store float* %59, float** %5
%60 = load float*, float** %4
%61 = bitcast float* %60 to i8*
%62 = ptrtoint i8* %61 to i64
%63 = urem i64 %62, 16
%64 = icmp eq i64 %63, 0
br i1 %64, label %65, label %71
65: ; preds = %47
%66 = load float*, float** %5
%67 = bitcast float* %66 to i8*
%68 = ptrtoint i8* %67 to i64
%69 = urem i64 %68, 16
%70 = icmp eq i64 %69, 0
br label %71
71: ; preds = %65, %47
%72 = phi i1 [ false, %47 ], [ %70, %65 ]
%73 = zext i1 %72 to i32
store i32 %73, i32* %8
%74 = load i32, i32* %8
%75 = icmp ne i32 %74, 0
br i1 %75, label %83, label %76
76: ; preds = %71
%77 = load i32, i32* %9
%78 = add nsw i32 %77, 1
store i32 %78, i32* %9
%79 = icmp sgt i32 %78, 2
br i1 %79, label %83, label %80
80: ; preds = %76
%81 = load i32, i32* %6
%82 = icmp eq i32 %81, 0
br i1 %82, label %83, label %84
83: ; preds = %80, %76, %71
br label %85
84: ; preds = %80
br label %43
85: ; preds = %83, %43
br label %86
86: ; preds = %85, %32, %27
%87 = call i32 (...) @_mm_setzero_ps()
store i32 %87, i32* %10
%88 = load i32, i32* %8
%89 = icmp ne i32 %88, 0
br i1 %89, label %90, label %116
90: ; preds = %86
store i32 0, i32* %11
%91 = load i32, i32* %6
%92 = sdiv i32 %91, 4
%93 = mul nsw i32 4, %92
store i32 %93, i32* %12
br label %94
94: ; preds = %112, %90
%95 = load i32, i32* %11
%96 = load i32, i32* %12
%97 = icmp slt i32 %95, %96
br i1 %97, label %98, label %115
98: ; preds = %94
%99 = load i32, i32* %10
%100 = load float*, float** %4
%101 = load i32, i32* %11
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds float, float* %100, i64 %102
%104 = call i32 @_mm_load_ps(float* %103)
%105 = load float*, float** %5
%106 = load i32, i32* %11
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds float, float* %105, i64 %107
%109 = call i32 @_mm_load_ps(float* %108)
%110 = call i32 @_mm_mul_ps(i32 %104, i32 %109)
%111 = call i32 @_mm_add_ps(i32 %99, i32 %110)
store i32 %111, i32* %10
br label %112
112: ; preds = %98
%113 = load i32, i32* %11
%114 = add nsw i32 %113, 4
store i32 %114, i32* %11
br label %94
115: ; preds = %94
br label %142
116: ; preds = %86
store i32 0, i32* %13
%117 = load i32, i32* %6
%118 = sdiv i32 %117, 4
%119 = mul nsw i32 4, %118
store i32 %119, i32* %14
br label %120
120: ; preds = %138, %116
%121 = load i32, i32* %13
%122 = load i32, i32* %14
%123 = icmp slt i32 %121, %122
br i1 %123, label %124, label %141
124: ; preds = %120
%125 = load i32, i32* %10
%126 = load float*, float** %4
%127 = load i32, i32* %13
%128 = sext i32 %127 to i64
%129 = getelementptr inbounds float, float* %126, i64 %128
%130 = call i32 @_mm_loadu_ps(float* %129)
%131 = load float*, float** %5
%132 = load i32, i32* %13
%133 = sext i32 %132 to i64
%134 = getelementptr inbounds float, float* %131, i64 %133
%135 = call i32 @_mm_loadu_ps(float* %134)
%136 = call i32 @_mm_mul_ps(i32 %130, i32 %135)
%137 = call i32 @_mm_add_ps(i32 %125, i32 %136)
store i32 %137, i32* %10
br label %138
138: ; preds = %124
%139 = load i32, i32* %13
%140 = add nsw i32 %139, 4
store i32 %140, i32* %13
br label %120
141: ; preds = %120
br label %142
142: ; preds = %141, %115
%143 = load i32, i32* %10
%144 = load i32, i32* %10
%145 = load i32, i32* %10
%146 = call i32 @_mm_movehl_ps(i32 %144, i32 %145)
%147 = call i32 @_mm_add_ps(i32 %143, i32 %146)
store i32 %147, i32* %10
%148 = load i32, i32* %10
%149 = load i32, i32* %10
%150 = load i32, i32* %10
%151 = call i32 @_mm_shuffle_ps(i32 %149, i32 %150, i32 1)
%152 = call i32 @_mm_add_ss(i32 %148, i32 %151)
store i32 %152, i32* %10
%153 = load i32, i32* %10
%154 = call i64 @_mm_cvtss_f32(i32 %153)
%155 = sitofp i64 %154 to float
%156 = load float, float* %7
%157 = fadd float %156, %155
store float %157, float* %7
%158 = load i32, i32* %6
%159 = sdiv i32 %158, 4
%160 = mul nsw i32 4, %159
store i32 %160, i32* %15
br label %161
161: ; preds = %178, %142
%162 = load i32, i32* %15
%163 = load i32, i32* %6
%164 = icmp slt i32 %162, %163
br i1 %164, label %165, label %181
165: ; preds = %161
%166 = load float*, float** %4
%167 = load i32, i32* %15
%168 = sext i32 %167 to i64
%169 = getelementptr inbounds float, float* %166, i64 %168
%170 = load float, float* %169
%171 = load float*, float** %5
%172 = load i32, i32* %15
%173 = sext i32 %172 to i64
%174 = getelementptr inbounds float, float* %171, i64 %173
%175 = load float, float* %174
%176 = load float, float* %7
%177 = call float @llvm.fmuladd.f32(float %170, float %175, float %176)
store float %177, float* %7
br label %178
178: ; preds = %165
%179 = load i32, i32* %15
%180 = add nsw i32 %179, 1
store i32 %180, i32* %15
br label %161
181: ; preds = %161
%182 = load float, float* %7
ret float %182
}

LLVM IR (-O2):
define float @sdot_sse2(float* %0, float* %1, i32 %2) {
%4 = ptrtoint float* %0 to i64
%5 = ptrtoint float* %1 to i64
%6 = or i64 %5, %4
%7 = and i64 %6, 15
%8 = icmp eq i64 %7, 0
%9 = xor i64 %5, %4
%10 = and i64 %9, 15
%11 = icmp ne i64 %10, 0
%12 = or i1 %8, %11
br i1 %12, label %73, label %13
13: ; preds = %3
%14 = add i32 %2, -1
%15 = call i32 @llvm.umin.i32(i32 %14, i32 2)
%16 = load float, float* %0
%17 = load float, float* %1
%18 = tail call float @llvm.fmuladd.f32(float %16, float %17, float 0.000000e+00)
%19 = add nsw i32 %2, -1
%20 = getelementptr inbounds float, float* %0, i64 1
%21 = getelementptr inbounds float, float* %1, i64 1
%22 = ptrtoint float* %20 to i64
%23 = and i64 %22, 15
%24 = icmp eq i64 %23, 0
%25 = ptrtoint float* %21 to i64
%26 = and i64 %25, 15
%27 = icmp eq i64 %26, 0
%28 = select i1 %24, i1 %27, i1 false
br i1 %28, label %67, label %29
29: ; preds = %13
%30 = icmp eq i32 %15, 0
br i1 %30, label %61, label %31
31: ; preds = %29
%32 = load float, float* %20
%33 = load float, float* %21
%34 = tail call float @llvm.fmuladd.f32(float %32, float %33, float %18)
%35 = add nsw i32 %2, -2
%36 = getelementptr inbounds float, float* %0, i64 2
%37 = getelementptr inbounds float, float* %1, i64 2
%38 = ptrtoint float* %36 to i64
%39 = and i64 %38, 15
%40 = icmp eq i64 %39, 0
%41 = ptrtoint float* %37 to i64
%42 = and i64 %41, 15
%43 = icmp eq i64 %42, 0
%44 = select i1 %40, i1 %43, i1 false
br i1 %44, label %67, label %45
45: ; preds = %31
%46 = icmp eq i32 %15, 1
br i1 %46, label %61, label %47
47: ; preds = %45
%48 = load float, float* %36
%49 = load float, float* %37
%50 = tail call float @llvm.fmuladd.f32(float %48, float %49, float %34)
%51 = add nsw i32 %2, -3
%52 = getelementptr inbounds float, float* %0, i64 3
%53 = getelementptr inbounds float, float* %1, i64 3
%54 = ptrtoint float* %52 to i64
%55 = and i64 %54, 15
%56 = icmp eq i64 %55, 0
%57 = ptrtoint float* %53 to i64
%58 = and i64 %57, 15
%59 = icmp eq i64 %58, 0
%60 = select i1 %56, i1 %59, i1 false
br i1 %60, label %67, label %61
61: ; preds = %47, %45, %29
%62 = phi float [ %18, %29 ], [ %34, %45 ], [ %50, %47 ]
%63 = phi i32 [ %19, %29 ], [ %35, %45 ], [ %51, %47 ]
%64 = phi float* [ %20, %29 ], [ %36, %45 ], [ %52, %47 ]
%65 = phi float* [ %21, %29 ], [ %37, %45 ], [ %53, %47 ]
%66 = tail call i32 (...) @_mm_setzero_ps()
br label %97
67: ; preds = %47, %31, %13
%68 = phi i32 [ %19, %13 ], [ %35, %31 ], [ %51, %47 ]
%69 = phi float [ %18, %13 ], [ %34, %31 ], [ %50, %47 ]
%70 = phi float* [ %21, %13 ], [ %37, %31 ], [ %53, %47 ]
%71 = phi float* [ %20, %13 ], [ %36, %31 ], [ %52, %47 ]
%72 = tail call i32 (...) @_mm_setzero_ps()
br label %75
73: ; preds = %3
%74 = tail call i32 (...) @_mm_setzero_ps()
br i1 %8, label %75, label %97
75: ; preds = %67, %73
%76 = phi i32 [ %72, %67 ], [ %74, %73 ]
%77 = phi float* [ %71, %67 ], [ %0, %73 ]
%78 = phi float* [ %70, %67 ], [ %1, %73 ]
%79 = phi float [ %69, %67 ], [ 0.000000e+00, %73 ]
%80 = phi i32 [ %68, %67 ], [ %2, %73 ]
%81 = sdiv i32 %80, 4
%82 = shl nsw i32 %81, 2
%83 = icmp sgt i32 %80, 3
br i1 %83, label %84, label %119
84: ; preds = %75
%85 = sext i32 %82 to i64
br label %86
86: ; preds = %84, %86
%87 = phi i64 [ 0, %84 ], [ %95, %86 ]
%88 = phi i32 [ %76, %84 ], [ %94, %86 ]
%89 = getelementptr inbounds float, float* %77, i64 %87
%90 = tail call i32 @_mm_load_ps(float* %89)
%91 = getelementptr inbounds float, float* %78, i64 %87
%92 = tail call i32 @_mm_load_ps(float* %91)
%93 = tail call i32 @_mm_mul_ps(i32 %90, i32 %92)
%94 = tail call i32 @_mm_add_ps(i32 %88, i32 %93)
%95 = add nuw nsw i64 %87, 4
%96 = icmp slt i64 %95, %85
br i1 %96, label %86, label %119
97: ; preds = %61, %73
%98 = phi i32 [ %66, %61 ], [ %74, %73 ]
%99 = phi float* [ %64, %61 ], [ %0, %73 ]
%100 = phi float* [ %65, %61 ], [ %1, %73 ]
%101 = phi float [ %62, %61 ], [ 0.000000e+00, %73 ]
%102 = phi i32 [ %63, %61 ], [ %2, %73 ]
%103 = sdiv i32 %102, 4
%104 = shl nsw i32 %103, 2
%105 = icmp sgt i32 %102, 3
br i1 %105, label %106, label %119
106: ; preds = %97
%107 = sext i32 %104 to i64
br label %108
108: ; preds = %106, %108
%109 = phi i64 [ 0, %106 ], [ %117, %108 ]
%110 = phi i32 [ %98, %106 ], [ %116, %108 ]
%111 = getelementptr inbounds float, float* %99, i64 %109
%112 = tail call i32 @_mm_loadu_ps(float* %111)
%113 = getelementptr inbounds float, float* %100, i64 %109
%114 = tail call i32 @_mm_loadu_ps(float* %113)
%115 = tail call i32 @_mm_mul_ps(i32 %112, i32 %114)
%116 = tail call i32 @_mm_add_ps(i32 %110, i32 %115)
%117 = add nuw nsw i64 %109, 4
%118 = icmp slt i64 %117, %107
br i1 %118, label %108, label %119
119: ; preds = %108, %86, %97, %75
%120 = phi i32 [ %104, %97 ], [ %82, %75 ], [ %82, %86 ], [ %104, %108 ]
%121 = phi float* [ %99, %97 ], [ %77, %75 ], [ %77, %86 ], [ %99, %108 ]
%122 = phi float* [ %100, %97 ], [ %78, %75 ], [ %78, %86 ], [ %100, %108 ]
%123 = phi float [ %101, %97 ], [ %79, %75 ], [ %79, %86 ], [ %101, %108 ]
%124 = phi i32 [ %102, %97 ], [ %80, %75 ], [ %80, %86 ], [ %102, %108 ]
%125 = phi i32 [ %98, %97 ], [ %76, %75 ], [ %94, %86 ], [ %116, %108 ]
%126 = tail call i32 @_mm_movehl_ps(i32 %125, i32 %125)
%127 = tail call i32 @_mm_add_ps(i32 %125, i32 %126)
%128 = tail call i32 @_mm_shuffle_ps(i32 %127, i32 %127, i32 1)
%129 = tail call i32 @_mm_add_ss(i32 %127, i32 %128)
%130 = tail call i64 @_mm_cvtss_f32(i32 %129)
%131 = sitofp i64 %130 to float
%132 = fadd float %123, %131
%133 = icmp slt i32 %120, %124
br i1 %133, label %134, label %158
134: ; preds = %119
%135 = sext i32 %120 to i64
%136 = sext i32 %124 to i64
%137 = xor i64 %135, -1
%138 = add nsw i64 %137, %136
%139 = and i64 %136, 3
%140 = icmp eq i64 %139, 0
br i1 %140, label %153, label %141
141: ; preds = %134, %141
%142 = phi i64 [ %150, %141 ], [ %135, %134 ]
%143 = phi float [ %149, %141 ], [ %132, %134 ]
%144 = phi i64 [ %151, %141 ], [ 0, %134 ]
%145 = getelementptr inbounds float, float* %121, i64 %142
%146 = load float, float* %145
%147 = getelementptr inbounds float, float* %122, i64 %142
%148 = load float, float* %147
%149 = tail call float @llvm.fmuladd.f32(float %146, float %148, float %143)
%150 = add nsw i64 %142, 1
%151 = add i64 %144, 1
%152 = icmp eq i64 %151, %139
br i1 %152, label %153, label %141
153: ; preds = %141, %134
%154 = phi float [ undef, %134 ], [ %149, %141 ]
%155 = phi i64 [ %135, %134 ], [ %150, %141 ]
%156 = phi float [ %132, %134 ], [ %149, %141 ]
%157 = icmp ult i64 %138, 3
br i1 %157, label %158, label %160
158: ; preds = %153, %160, %119
%159 = phi float [ %132, %119 ], [ %154, %153 ], [ %185, %160 ]
ret float %159
160: ; preds = %153, %160
%161 = phi i64 [ %186, %160 ], [ %155, %153 ]
%162 = phi float [ %185, %160 ], [ %156, %153 ]
%163 = getelementptr inbounds float, float* %121, i64 %161
%164 = load float, float* %163
%165 = getelementptr inbounds float, float* %122, i64 %161
%166 = load float, float* %165
%167 = tail call float @llvm.fmuladd.f32(float %164, float %166, float %162)
%168 = add nsw i64 %161, 1
%169 = getelementptr inbounds float, float* %121, i64 %168
%170 = load float, float* %169
%171 = getelementptr inbounds float, float* %122, i64 %168
%172 = load float, float* %171
%173 = tail call float @llvm.fmuladd.f32(float %170, float %172, float %167)
%174 = add nsw i64 %161, 2
%175 = getelementptr inbounds float, float* %121, i64 %174
%176 = load float, float* %175
%177 = getelementptr inbounds float, float* %122, i64 %174
%178 = load float, float* %177
%179 = tail call float @llvm.fmuladd.f32(float %176, float %178, float %173)
%180 = add nsw i64 %161, 3
%181 = getelementptr inbounds float, float* %121, i64 %180
%182 = load float, float* %181
%183 = getelementptr inbounds float, float* %122, i64 %180
%184 = load float, float* %183
%185 = tail call float @llvm.fmuladd.f32(float %182, float %184, float %179)
%186 = add nsw i64 %161, 4
%187 = icmp eq i64 %186, %136
br i1 %187, label %158, label %160
}