#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*/ __m128d;

/* Variables and functions */
int /*<<< orphan*/ _mm_add_pd(int /*<<< orphan*/, int /*<<< orphan*/);
int /*<<< orphan*/ _mm_cvtps_pd(int /*<<< orphan*/);
scalar_t__ _mm_cvtsd_f64(int /*<<< orphan*/);
int /*<<< orphan*/ _mm_load_ps(float const *);
int /*<<< orphan*/ _mm_loadu_ps(float const *);
int /*<<< orphan*/ _mm_mul_ps(int /*<<< orphan*/, int /*<<< orphan*/);
int /*<<< orphan*/ _mm_setzero_pd();
int /*<<< orphan*/ _mm_shuffle_pd(int /*<<< orphan*/, int /*<<< orphan*/, int);
double dsdot_sse2(const float *a, const float *b, int n) {
  double s = 0.0;
  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;
    }
  }
  __m128d s2 = _mm_setzero_pd();
  __m128d s2u = _mm_setzero_pd();
  if (aligned) {
    for (int k = 0, nq = 4 * (n / 4); k < nq; k += 4) {
      s2 = _mm_add_pd(
          s2, _mm_cvtps_pd(_mm_mul_ps(_mm_load_ps(a + k), _mm_load_ps(b + k))));
      s2u = _mm_add_pd(s2u, _mm_cvtps_pd(_mm_mul_ps(_mm_loadu_ps(a + k + 2),
                                                    _mm_loadu_ps(b + k + 2))));
    }
  } else {
    for (int k = 0, nq = 4 * (n / 4); k < nq; k += 4) {
      s2 = _mm_add_pd(s2, _mm_cvtps_pd(_mm_mul_ps(_mm_loadu_ps(a + k),
                                                  _mm_loadu_ps(b + k))));
      s2u = _mm_add_pd(s2u, _mm_cvtps_pd(_mm_mul_ps(_mm_loadu_ps(a + k + 2),
                                                    _mm_loadu_ps(b + k + 2))));
    }
  }
  s2 = _mm_add_pd(s2, s2u);
  s2 = _mm_add_pd(s2, _mm_shuffle_pd(s2, s2, 1));
  s += _mm_cvtsd_f64(s2);
  for (int k = 4 * (n / 4); k < n; k++)
    s += a[k] * b[k];
  return s;
}


x86 assembly (clang, -O0):
<dsdot_sse2>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rcx
and    $0xf,%rcx
xor    %eax,%eax
cmp    $0x0,%rcx
mov    %al,-0x45(%rbp)
jne    11c4 <dsdot_sse2+0x44>
mov    -0x10(%rbp),%rax
and    $0xf,%rax
cmp    $0x0,%rax
sete   %al
mov    %al,-0x45(%rbp)
mov    -0x45(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jne    12b7 <dsdot_sse2+0x137>
mov    -0x8(%rbp),%rax
and    $0xf,%rax
mov    -0x10(%rbp),%rcx
and    $0xf,%rcx
cmp    %rcx,%rax
jne    12b7 <dsdot_sse2+0x137>
movl   $0x0,-0x28(%rbp)
cmpl   $0x0,-0x24(%rbp)
setne  %al
xor    $0xff,%al
test   $0x1,%al
jne    120f <dsdot_sse2+0x8f>
jmp    12b2 <dsdot_sse2+0x132>
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x10(%rbp),%rax
mulss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%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,-0x46(%rbp)
jne    1277 <dsdot_sse2+0xf7>
mov    -0x10(%rbp),%rax
and    $0xf,%rax
cmp    $0x0,%rax
sete   %al
mov    %al,-0x46(%rbp)
mov    -0x46(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jne    12a8 <dsdot_sse2+0x128>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
cmp    $0x2,%eax
jg     12a8 <dsdot_sse2+0x128>
cmpl   $0x0,-0x14(%rbp)
jne    12ad <dsdot_sse2+0x12d>
jmp    12b2 <dsdot_sse2+0x132>
jmp    11f9 <dsdot_sse2+0x79>
jmp    12b7 <dsdot_sse2+0x137>
mov    $0x0,%al
call   10a0 <_mm_setzero_pd@plt>
mov    %eax,-0x2c(%rbp)
mov    $0x0,%al
call   10a0 <_mm_setzero_pd@plt>
mov    %eax,-0x30(%rbp)
cmpl   $0x0,-0x24(%rbp)
je     13b2 <dsdot_sse2+0x232>
movl   $0x0,-0x34(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
shl    $0x2,%eax
mov    %eax,-0x38(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    13ad <dsdot_sse2+0x22d>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x8(%rbp),%rdi
movslq -0x34(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
call   1030 <_mm_load_ps@plt>
mov    %eax,-0x58(%rbp)
mov    -0x10(%rbp),%rdi
movslq -0x34(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
call   1030 <_mm_load_ps@plt>
mov    -0x58(%rbp),%edi
mov    %eax,%esi
call   1080 <_mm_mul_ps@plt>
mov    %eax,%edi
call   1090 <_mm_cvtps_pd@plt>
mov    -0x54(%rbp),%edi
mov    %eax,%esi
call   1050 <_mm_add_pd@plt>
mov    %eax,-0x2c(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x8(%rbp),%rdi
movslq -0x34(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
add    $0x8,%rdi
call   1060 <_mm_loadu_ps@plt>
mov    %eax,-0x50(%rbp)
mov    -0x10(%rbp),%rdi
movslq -0x34(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
add    $0x8,%rdi
call   1060 <_mm_loadu_ps@plt>
mov    -0x50(%rbp),%edi
mov    %eax,%esi
call   1080 <_mm_mul_ps@plt>
mov    %eax,%edi
call   1090 <_mm_cvtps_pd@plt>
mov    -0x4c(%rbp),%edi
mov    %eax,%esi
call   1050 <_mm_add_pd@plt>
mov    %eax,-0x30(%rbp)
mov    -0x34(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x34(%rbp)
jmp    12ed <dsdot_sse2+0x16d>
jmp    148f <dsdot_sse2+0x30f>
movl   $0x0,-0x3c(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
shl    $0x2,%eax
mov    %eax,-0x40(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    148a <dsdot_sse2+0x30a>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x8(%rbp),%rdi
movslq -0x3c(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
call   1060 <_mm_loadu_ps@plt>
mov    %eax,-0x68(%rbp)
mov    -0x10(%rbp),%rdi
movslq -0x3c(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
call   1060 <_mm_loadu_ps@plt>
mov    -0x68(%rbp),%edi
mov    %eax,%esi
call   1080 <_mm_mul_ps@plt>
mov    %eax,%edi
call   1090 <_mm_cvtps_pd@plt>
mov    -0x64(%rbp),%edi
mov    %eax,%esi
call   1050 <_mm_add_pd@plt>
mov    %eax,-0x2c(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x8(%rbp),%rdi
movslq -0x3c(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
add    $0x8,%rdi
call   1060 <_mm_loadu_ps@plt>
mov    %eax,-0x60(%rbp)
mov    -0x10(%rbp),%rdi
movslq -0x3c(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
add    $0x8,%rdi
call   1060 <_mm_loadu_ps@plt>
mov    -0x60(%rbp),%edi
mov    %eax,%esi
call   1080 <_mm_mul_ps@plt>
mov    %eax,%edi
call   1090 <_mm_cvtps_pd@plt>
mov    -0x5c(%rbp),%edi
mov    %eax,%esi
call   1050 <_mm_add_pd@plt>
mov    %eax,-0x30(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x3c(%rbp)
jmp    13ca <dsdot_sse2+0x24a>
jmp    148f <dsdot_sse2+0x30f>
mov    -0x2c(%rbp),%edi
mov    -0x30(%rbp),%esi
call   1050 <_mm_add_pd@plt>
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x6c(%rbp)
mov    -0x2c(%rbp),%edi
mov    -0x2c(%rbp),%esi
mov    $0x1,%edx
call   1070 <_mm_shuffle_pd@plt>
mov    -0x6c(%rbp),%edi
mov    %eax,%esi
call   1050 <_mm_add_pd@plt>
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%edi
call   1040 <_mm_cvtsd_f64@plt>
cvtsi2sd %rax,%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
shl    $0x2,%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    152a <dsdot_sse2+0x3aa>
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    14e8 <dsdot_sse2+0x368>
movsd  -0x20(%rbp),%xmm0
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dsdot_sse2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %edx,%ebp
mov    %rsi,%r12
mov    %rdi,%r14
mov    %r12d,%eax
or     %r14d,%eax
xor    %ebx,%ebx
test   $0xf,%al
sete   %bl
je     11fa <dsdot_sse2+0x7a>
mov    %r12d,%eax
xor    %r14d,%eax
and    $0xf,%eax
test   %rax,%rax
jne    11fa <dsdot_sse2+0x7a>
lea    -0x1(%rbp),%ecx
cmp    $0x2,%ecx
mov    $0x2,%eax
cmovb  %ecx,%eax
movss  (%r14),%xmm0
mulss  (%r12),%xmm0
cvtss2sd %xmm0,%xmm0
xorpd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
lea    0x4(%r14),%rdx
lea    0x4(%r12),%rsi
mov    %edx,%edi
or     %esi,%edi
mov    $0x1,%ebx
test   $0xf,%dil
jne    1208 <dsdot_sse2+0x88>
movsd  %xmm1,0x8(%rsp)
mov    %ecx,%ebp
jmp    124b <dsdot_sse2+0xcb>
xorps  %xmm0,%xmm0
movsd  %xmm0,0x8(%rsp)
jmp    129d <dsdot_sse2+0x11d>
test   %eax,%eax
je     1241 <dsdot_sse2+0xc1>
movss  0x4(%r14),%xmm0
mulss  0x4(%r12),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  %xmm0,%xmm1
lea    -0x2(%rbp),%esi
lea    0x8(%r14),%rcx
lea    0x8(%r12),%rdx
mov    %ecx,%edi
or     %edx,%edi
test   $0xf,%dil
jne    1253 <dsdot_sse2+0xd3>
movsd  %xmm1,0x8(%rsp)
mov    %esi,%ebp
jmp    1262 <dsdot_sse2+0xe2>
movsd  %xmm1,0x8(%rsp)
mov    %ecx,%ebp
xor    %ebx,%ebx
mov    %rsi,%r12
mov    %rdx,%r14
jmp    129d <dsdot_sse2+0x11d>
cmp    $0x1,%eax
jne    126a <dsdot_sse2+0xea>
movsd  %xmm1,0x8(%rsp)
mov    %esi,%ebp
xor    %ebx,%ebx
mov    %rdx,%r12
mov    %rcx,%r14
jmp    129d <dsdot_sse2+0x11d>
movss  0x8(%r14),%xmm0
mulss  0x8(%r12),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
add    $0xfffffffd,%ebp
add    $0xc,%r14
add    $0xc,%r12
mov    %r14d,%eax
or     %r12d,%eax
xor    %ebx,%ebx
test   $0xf,%al
sete   %bl
xor    %eax,%eax
call   10a0 <_mm_setzero_pd@plt>
mov    %eax,%r13d
xor    %eax,%eax
call   10a0 <_mm_setzero_pd@plt>
mov    %eax,%r15d
lea    0x3(%rbp),%eax
test   %ebp,%ebp
cmovns %ebp,%eax
and    $0xfffffffc,%eax
test   %ebx,%ebx
mov    %r12,0x28(%rsp)
mov    %r14,0x20(%rsp)
mov    %rbp,0x30(%rsp)
mov    %eax,0x1c(%rsp)
je     1378 <dsdot_sse2+0x1f8>
cmp    $0x4,%ebp
jl     1417 <dsdot_sse2+0x297>
cltq
mov    %rax,0x10(%rsp)
mov    %r12,%rax
xor    %r12d,%r12d
mov    %r14,%rbp
mov    %rax,%r14
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %rbp,%rdi
call   1030 <_mm_load_ps@plt>
mov    %eax,%ebx
mov    %r14,%rdi
call   1030 <_mm_load_ps@plt>
mov    %ebx,%edi
mov    %eax,%esi
call   1080 <_mm_mul_ps@plt>
mov    %eax,%edi
call   1090 <_mm_cvtps_pd@plt>
mov    %r13d,%edi
mov    %eax,%esi
call   1050 <_mm_add_pd@plt>
mov    %eax,%r13d
lea    0x8(%rbp),%rdi
call   1060 <_mm_loadu_ps@plt>
mov    %eax,%ebx
lea    0x8(%r14),%rdi
call   1060 <_mm_loadu_ps@plt>
mov    %ebx,%edi
mov    %eax,%esi
call   1080 <_mm_mul_ps@plt>
mov    %eax,%edi
call   1090 <_mm_cvtps_pd@plt>
mov    %r15d,%edi
mov    %eax,%esi
call   1050 <_mm_add_pd@plt>
mov    %eax,%r15d
add    $0x4,%r12
add    $0x10,%r14
add    $0x10,%rbp
cmp    0x10(%rsp),%r12
jl     1300 <dsdot_sse2+0x180>
jmp    1417 <dsdot_sse2+0x297>
cmp    $0x4,%ebp
jl     1417 <dsdot_sse2+0x297>
cltq
mov    %rax,0x10(%rsp)
mov    %r12,%rbx
xor    %r12d,%r12d
mov    %r14,%rbp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rbp,%rdi
call   1060 <_mm_loadu_ps@plt>
mov    %eax,%r14d
mov    %rbx,%rdi
call   1060 <_mm_loadu_ps@plt>
mov    %r14d,%edi
mov    %eax,%esi
call   1080 <_mm_mul_ps@plt>
mov    %eax,%edi
call   1090 <_mm_cvtps_pd@plt>
mov    %r13d,%edi
mov    %eax,%esi
call   1050 <_mm_add_pd@plt>
mov    %eax,%r13d
lea    0x8(%rbp),%rdi
call   1060 <_mm_loadu_ps@plt>
mov    %eax,%r14d
lea    0x8(%rbx),%rdi
call   1060 <_mm_loadu_ps@plt>
mov    %r14d,%edi
mov    %eax,%esi
call   1080 <_mm_mul_ps@plt>
mov    %eax,%edi
call   1090 <_mm_cvtps_pd@plt>
mov    %r15d,%edi
mov    %eax,%esi
call   1050 <_mm_add_pd@plt>
mov    %eax,%r15d
add    $0x4,%r12
add    $0x10,%rbx
add    $0x10,%rbp
cmp    0x10(%rsp),%r12
jl     13a0 <dsdot_sse2+0x220>
mov    %r13d,%edi
mov    %r15d,%esi
call   1050 <_mm_add_pd@plt>
mov    %eax,%ebx
mov    %eax,%edi
mov    %eax,%esi
mov    $0x1,%edx
call   1070 <_mm_shuffle_pd@plt>
mov    %ebx,%edi
mov    %eax,%esi
call   1050 <_mm_add_pd@plt>
mov    %eax,%edi
call   1040 <_mm_cvtsd_f64@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %rax,%xmm1
movsd  0x8(%rsp),%xmm0
addsd  %xmm1,%xmm0
mov    0x30(%rsp),%rcx
mov    0x1c(%rsp),%eax
cmp    %ecx,%eax
jge    151e <dsdot_sse2+0x39e>
cltq
movslq %ecx,%rcx
mov    %rax,%rdx
not    %rdx
add    %rcx,%rdx
mov    %rcx,%rsi
and    $0x3,%rsi
je     152d <dsdot_sse2+0x3ad>
mov    0x28(%rsp),%rdi
mov    0x20(%rsp),%rbp
nopw   0x0(%rax,%rax,1)
movss  0x0(%rbp,%rax,4),%xmm1
mulss  (%rdi,%rax,4),%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rax
add    $0xffffffffffffffff,%rsi
jne    1490 <dsdot_sse2+0x310>
cmp    $0x3,%rdx
jb     151e <dsdot_sse2+0x39e>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  0x0(%rbp,%rax,4),%xmm2
movss  0x4(%rbp,%rax,4),%xmm1
mulss  (%rdi,%rax,4),%xmm2
cvtss2sd %xmm2,%xmm2
addsd  %xmm0,%xmm2
mulss  0x4(%rdi,%rax,4),%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm2,%xmm1
movss  0x8(%rbp,%rax,4),%xmm0
mulss  0x8(%rdi,%rax,4),%xmm0
xorps  %xmm2,%xmm2
cvtss2sd %xmm0,%xmm2
addsd  %xmm1,%xmm2
movss  0xc(%rbp,%rax,4),%xmm1
mulss  0xc(%rdi,%rax,4),%xmm1
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
addsd  %xmm2,%xmm0
add    $0x4,%rax
cmp    %rax,%rcx
jne    14c0 <dsdot_sse2+0x340>
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    0x28(%rsp),%rdi
mov    0x20(%rsp),%rbp
cmp    $0x3,%rdx
jae    14c0 <dsdot_sse2+0x340>
jmp    151e <dsdot_sse2+0x39e>

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

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