#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*/);
scalar_t__ _mm_cvtsd_f64(int /*<<< orphan*/);
int /*<<< orphan*/ _mm_load_pd(double const *);
int /*<<< orphan*/ _mm_loadu_pd(double const *);
int /*<<< orphan*/ _mm_mul_pd(int /*<<< orphan*/, int /*<<< orphan*/);
int /*<<< orphan*/ _mm_setzero_pd();
int /*<<< orphan*/ _mm_shuffle_pd(int /*<<< orphan*/, int /*<<< orphan*/, int);
double ddot_sse2(const double *a, const double *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)))) {
    s += (*a) * (*b);
    n--;
    a++;
    b++;
    aligned = ((((uintptr_t)(const void *)(a)) % (16)) == 0) &&
              ((((uintptr_t)(const void *)(b)) % (16)) == 0);
  }
  __m128d s2 = _mm_setzero_pd();
  if (aligned)
    for (int k = 0, nq = 2 * (n / 2); k < nq; k += 2)
      s2 = _mm_add_pd(s2, _mm_mul_pd(_mm_load_pd(a + k), _mm_load_pd(b + k)));
  else
    for (int k = 0, nq = 2 * (n / 2); k < nq; k += 2)
      s2 = _mm_add_pd(s2, _mm_mul_pd(_mm_loadu_pd(a + k), _mm_loadu_pd(b + k)));
  s2 = _mm_add_pd(s2, _mm_shuffle_pd(s2, s2, 1));
  s += _mm_cvtsd_f64(s2);
  for (int k = 2 * (n / 2); k < n; k++)
    s += a[k] * b[k];
  return s;
}


x86 assembly (clang, -O0):
<ddot_sse2>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%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,-0x3d(%rbp)
jne    11b4 <ddot_sse2+0x44>
mov    -0x10(%rbp),%rax
and    $0xf,%rax
cmp    $0x0,%rax
sete   %al
mov    %al,-0x3d(%rbp)
mov    -0x3d(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jne    1259 <ddot_sse2+0xe9>
mov    -0x8(%rbp),%rax
and    $0xf,%rax
mov    -0x10(%rbp),%rcx
and    $0xf,%rcx
cmp    %rcx,%rax
jne    1259 <ddot_sse2+0xe9>
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm2
movsd  -0x20(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
mov    -0x8(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x8(%rbp)
mov    -0x10(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x10(%rbp)
mov    -0x8(%rbp),%rcx
and    $0xf,%rcx
xor    %eax,%eax
cmp    $0x0,%rcx
mov    %al,-0x3e(%rbp)
jne    124e <ddot_sse2+0xde>
mov    -0x10(%rbp),%rax
and    $0xf,%rax
cmp    $0x0,%rax
sete   %al
mov    %al,-0x3e(%rbp)
mov    -0x3e(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x24(%rbp)
mov    $0x0,%al
call   1070 <_mm_setzero_pd@plt>
mov    %eax,-0x28(%rbp)
cmpl   $0x0,-0x24(%rbp)
je     12ec <ddot_sse2+0x17c>
movl   $0x0,-0x2c(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
shl    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jge    12e7 <ddot_sse2+0x177>
mov    -0x28(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x8(%rbp),%rdi
movslq -0x2c(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
call   1060 <_mm_load_pd@plt>
mov    %eax,-0x48(%rbp)
mov    -0x10(%rbp),%rdi
movslq -0x2c(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
call   1060 <_mm_load_pd@plt>
mov    -0x48(%rbp),%edi
mov    %eax,%esi
call   1080 <_mm_mul_pd@plt>
mov    -0x44(%rbp),%edi
mov    %eax,%esi
call   1040 <_mm_add_pd@plt>
mov    %eax,-0x28(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x2c(%rbp)
jmp    1285 <ddot_sse2+0x115>
jmp    136b <ddot_sse2+0x1fb>
movl   $0x0,-0x34(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
shl    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    1366 <ddot_sse2+0x1f6>
mov    -0x28(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x8(%rbp),%rdi
movslq -0x34(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
call   1090 <_mm_loadu_pd@plt>
mov    %eax,-0x50(%rbp)
mov    -0x10(%rbp),%rdi
movslq -0x34(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
call   1090 <_mm_loadu_pd@plt>
mov    -0x50(%rbp),%edi
mov    %eax,%esi
call   1080 <_mm_mul_pd@plt>
mov    -0x4c(%rbp),%edi
mov    %eax,%esi
call   1040 <_mm_add_pd@plt>
mov    %eax,-0x28(%rbp)
mov    -0x34(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x34(%rbp)
jmp    1304 <ddot_sse2+0x194>
jmp    136b <ddot_sse2+0x1fb>
mov    -0x28(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x28(%rbp),%edi
mov    -0x28(%rbp),%esi
mov    $0x1,%edx
call   1050 <_mm_shuffle_pd@plt>
mov    -0x54(%rbp),%edi
mov    %eax,%esi
call   1040 <_mm_add_pd@plt>
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%edi
call   1030 <_mm_cvtsd_f64@plt>
cvtsi2sd %rax,%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
shl    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    13fc <ddot_sse2+0x28c>
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x20(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    13b6 <ddot_sse2+0x246>
movsd  -0x20(%rbp),%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ddot_sse2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %edx,%r12d
mov    %rsi,%rbp
mov    %rdi,%rbx
mov    %rsi,%r14
or     %rdi,%r14
test   $0xf,%r14b
je     1262 <ddot_sse2+0xf2>
mov    %ebp,%eax
xor    %ebx,%eax
and    $0xf,%eax
test   %rax,%rax
jne    1262 <ddot_sse2+0xf2>
movsd  (%rbx),%xmm0
mulsd  0x0(%rbp),%xmm0
xorpd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x20(%rsp)
add    $0xffffffff,%r12d
add    $0x8,%rbx
add    $0x8,%rbp
mov    %ebx,%r14d
or     %ebp,%r14d
xor    %eax,%eax
call   1070 <_mm_setzero_pd@plt>
mov    %eax,%r13d
test   $0xf,%r14b
je     1280 <ddot_sse2+0x110>
mov    %r12d,%r14d
shr    $0x1f,%r14d
add    %r12d,%r14d
and    $0xfffffffe,%r14d
cmp    $0x2,%r12d
jl     12fc <ddot_sse2+0x18c>
mov    %r14d,0x4(%rsp)
movslq %r14d,%rax
mov    %rax,0x18(%rsp)
xor    %r15d,%r15d
mov    %rbx,0x8(%rsp)
mov    %rbp,0x10(%rsp)
cs nopw 0x0(%rax,%rax,1)
mov    %rbx,%rdi
call   1090 <_mm_loadu_pd@plt>
mov    %eax,%r14d
mov    %rbp,%rdi
call   1090 <_mm_loadu_pd@plt>
mov    %r14d,%edi
mov    %eax,%esi
call   1080 <_mm_mul_pd@plt>
mov    %r13d,%edi
mov    %eax,%esi
call   1040 <_mm_add_pd@plt>
mov    %eax,%r13d
add    $0x2,%r15
add    $0x10,%rbp
add    $0x10,%rbx
cmp    0x18(%rsp),%r15
jl     1220 <ddot_sse2+0xb0>
jmp    12ed <ddot_sse2+0x17d>
xor    %eax,%eax
call   1070 <_mm_setzero_pd@plt>
mov    %eax,%r13d
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x20(%rsp)
test   $0xf,%r14b
jne    11e4 <ddot_sse2+0x74>
mov    %r12d,%r14d
shr    $0x1f,%r14d
add    %r12d,%r14d
and    $0xfffffffe,%r14d
cmp    $0x2,%r12d
jl     12fc <ddot_sse2+0x18c>
mov    %r14d,0x4(%rsp)
movslq %r14d,%rax
mov    %rax,0x18(%rsp)
xor    %r15d,%r15d
mov    %rbx,0x8(%rsp)
mov    %rbp,0x10(%rsp)
xchg   %ax,%ax
mov    %rbx,%rdi
call   1060 <_mm_load_pd@plt>
mov    %eax,%r14d
mov    %rbp,%rdi
call   1060 <_mm_load_pd@plt>
mov    %r14d,%edi
mov    %eax,%esi
call   1080 <_mm_mul_pd@plt>
mov    %r13d,%edi
mov    %eax,%esi
call   1040 <_mm_add_pd@plt>
mov    %eax,%r13d
add    $0x2,%r15
add    $0x10,%rbp
add    $0x10,%rbx
cmp    0x18(%rsp),%r15
jl     12b0 <ddot_sse2+0x140>
mov    0x10(%rsp),%rbp
mov    0x8(%rsp),%rbx
mov    0x4(%rsp),%r14d
mov    %r13d,%edi
mov    %r13d,%esi
mov    $0x1,%edx
call   1050 <_mm_shuffle_pd@plt>
mov    %r13d,%edi
mov    %eax,%esi
call   1040 <_mm_add_pd@plt>
mov    %eax,%edi
call   1030 <_mm_cvtsd_f64@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
addsd  0x20(%rsp),%xmm0
cmp    %r12d,%r14d
jge    13b8 <ddot_sse2+0x248>
movslq %r14d,%rax
movslq %r12d,%rcx
sub    %r14d,%r12d
mov    %rax,%rdx
not    %rdx
add    %rcx,%rdx
and    $0x3,%r12
je     1369 <ddot_sse2+0x1f9>
nopl   0x0(%rax)
movsd  (%rbx,%rax,8),%xmm1
mulsd  0x0(%rbp,%rax,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rax
add    $0xffffffffffffffff,%r12
jne    1350 <ddot_sse2+0x1e0>
cmp    $0x3,%rdx
jb     13b8 <ddot_sse2+0x248>
nop
movsd  (%rbx,%rax,8),%xmm1
movsd  0x8(%rbx,%rax,8),%xmm2
mulsd  0x0(%rbp,%rax,8),%xmm1
mulsd  0x8(%rbp,%rax,8),%xmm2
addsd  %xmm0,%xmm1
movsd  0x10(%rbx,%rax,8),%xmm3
mulsd  0x10(%rbp,%rax,8),%xmm3
addsd  %xmm1,%xmm2
movsd  0x18(%rbx,%rax,8),%xmm0
mulsd  0x18(%rbp,%rax,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%rax
cmp    %rax,%rcx
jne    1370 <ddot_sse2+0x200>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define double @ddot_sse2(double* %0, double* %1, i32 %2) {
%4 = ptrtoint double* %0 to i64
%5 = ptrtoint double* %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 %28, label %13
13: ; preds = %3
%14 = load double, double* %0
%15 = load double, double* %1
%16 = tail call double @llvm.fmuladd.f64(double %14, double %15, double 0.000000e+00)
%17 = add nsw i32 %2, -1
%18 = getelementptr inbounds double, double* %0, i64 1
%19 = getelementptr inbounds double, double* %1, i64 1
%20 = ptrtoint double* %18 to i64
%21 = and i64 %20, 15
%22 = icmp eq i64 %21, 0
%23 = ptrtoint double* %19 to i64
%24 = and i64 %23, 15
%25 = icmp eq i64 %24, 0
%26 = select i1 %22, i1 %25, i1 false
%27 = tail call i32 (...) @_mm_setzero_pd()
br i1 %26, label %30, label %52
28: ; preds = %3
%29 = tail call i32 (...) @_mm_setzero_pd()
br i1 %8, label %30, label %52
30: ; preds = %13, %28
%31 = phi i32 [ %27, %13 ], [ %29, %28 ]
%32 = phi double* [ %18, %13 ], [ %0, %28 ]
%33 = phi double* [ %19, %13 ], [ %1, %28 ]
%34 = phi double [ %16, %13 ], [ 0.000000e+00, %28 ]
%35 = phi i32 [ %17, %13 ], [ %2, %28 ]
%36 = sdiv i32 %35, 2
%37 = shl nsw i32 %36, 1
%38 = icmp sgt i32 %35, 1
br i1 %38, label %39, label %74
39: ; preds = %30
%40 = sext i32 %37 to i64
br label %41
41: ; preds = %39, %41
%42 = phi i64 [ 0, %39 ], [ %50, %41 ]
%43 = phi i32 [ %31, %39 ], [ %49, %41 ]
%44 = getelementptr inbounds double, double* %32, i64 %42
%45 = tail call i32 @_mm_load_pd(double* %44)
%46 = getelementptr inbounds double, double* %33, i64 %42
%47 = tail call i32 @_mm_load_pd(double* %46)
%48 = tail call i32 @_mm_mul_pd(i32 %45, i32 %47)
%49 = tail call i32 @_mm_add_pd(i32 %43, i32 %48)
%50 = add nuw nsw i64 %42, 2
%51 = icmp slt i64 %50, %40
br i1 %51, label %41, label %74
52: ; preds = %13, %28
%53 = phi i32 [ %27, %13 ], [ %29, %28 ]
%54 = phi double* [ %18, %13 ], [ %0, %28 ]
%55 = phi double* [ %19, %13 ], [ %1, %28 ]
%56 = phi double [ %16, %13 ], [ 0.000000e+00, %28 ]
%57 = phi i32 [ %17, %13 ], [ %2, %28 ]
%58 = sdiv i32 %57, 2
%59 = shl nsw i32 %58, 1
%60 = icmp sgt i32 %57, 1
br i1 %60, label %61, label %74
61: ; preds = %52
%62 = sext i32 %59 to i64
br label %63
63: ; preds = %61, %63
%64 = phi i64 [ 0, %61 ], [ %72, %63 ]
%65 = phi i32 [ %53, %61 ], [ %71, %63 ]
%66 = getelementptr inbounds double, double* %54, i64 %64
%67 = tail call i32 @_mm_loadu_pd(double* %66)
%68 = getelementptr inbounds double, double* %55, i64 %64
%69 = tail call i32 @_mm_loadu_pd(double* %68)
%70 = tail call i32 @_mm_mul_pd(i32 %67, i32 %69)
%71 = tail call i32 @_mm_add_pd(i32 %65, i32 %70)
%72 = add nuw nsw i64 %64, 2
%73 = icmp slt i64 %72, %62
br i1 %73, label %63, label %74
74: ; preds = %63, %41, %52, %30
%75 = phi i32 [ %59, %52 ], [ %37, %30 ], [ %37, %41 ], [ %59, %63 ]
%76 = phi double* [ %54, %52 ], [ %32, %30 ], [ %32, %41 ], [ %54, %63 ]
%77 = phi double* [ %55, %52 ], [ %33, %30 ], [ %33, %41 ], [ %55, %63 ]
%78 = phi double [ %56, %52 ], [ %34, %30 ], [ %34, %41 ], [ %56, %63 ]
%79 = phi i32 [ %57, %52 ], [ %35, %30 ], [ %35, %41 ], [ %57, %63 ]
%80 = phi i32 [ %53, %52 ], [ %31, %30 ], [ %49, %41 ], [ %71, %63 ]
%81 = tail call i32 @_mm_shuffle_pd(i32 %80, i32 %80, i32 1)
%82 = tail call i32 @_mm_add_pd(i32 %80, i32 %81)
%83 = tail call i64 @_mm_cvtsd_f64(i32 %82)
%84 = sitofp i64 %83 to double
%85 = fadd double %78, %84
%86 = icmp slt i32 %75, %79
br i1 %86, label %87, label %112
87: ; preds = %74
%88 = sext i32 %75 to i64
%89 = sext i32 %79 to i64
%90 = sub nsw i64 %89, %88
%91 = xor i64 %88, -1
%92 = add nsw i64 %91, %89
%93 = and i64 %90, 3
%94 = icmp eq i64 %93, 0
br i1 %94, label %107, label %95
95: ; preds = %87, %95
%96 = phi i64 [ %104, %95 ], [ %88, %87 ]
%97 = phi double [ %103, %95 ], [ %85, %87 ]
%98 = phi i64 [ %105, %95 ], [ 0, %87 ]
%99 = getelementptr inbounds double, double* %76, i64 %96
%100 = load double, double* %99
%101 = getelementptr inbounds double, double* %77, i64 %96
%102 = load double, double* %101
%103 = tail call double @llvm.fmuladd.f64(double %100, double %102, double %97)
%104 = add nsw i64 %96, 1
%105 = add i64 %98, 1
%106 = icmp eq i64 %105, %93
br i1 %106, label %107, label %95
107: ; preds = %95, %87
%108 = phi double [ undef, %87 ], [ %103, %95 ]
%109 = phi i64 [ %88, %87 ], [ %104, %95 ]
%110 = phi double [ %85, %87 ], [ %103, %95 ]
%111 = icmp ult i64 %92, 3
br i1 %111, label %112, label %114
112: ; preds = %107, %114, %74
%113 = phi double [ %85, %74 ], [ %108, %107 ], [ %139, %114 ]
ret double %113
114: ; preds = %107, %114
%115 = phi i64 [ %140, %114 ], [ %109, %107 ]
%116 = phi double [ %139, %114 ], [ %110, %107 ]
%117 = getelementptr inbounds double, double* %76, i64 %115
%118 = load double, double* %117
%119 = getelementptr inbounds double, double* %77, i64 %115
%120 = load double, double* %119
%121 = tail call double @llvm.fmuladd.f64(double %118, double %120, double %116)
%122 = add nsw i64 %115, 1
%123 = getelementptr inbounds double, double* %76, i64 %122
%124 = load double, double* %123
%125 = getelementptr inbounds double, double* %77, i64 %122
%126 = load double, double* %125
%127 = tail call double @llvm.fmuladd.f64(double %124, double %126, double %121)
%128 = add nsw i64 %115, 2
%129 = getelementptr inbounds double, double* %76, i64 %128
%130 = load double, double* %129
%131 = getelementptr inbounds double, double* %77, i64 %128
%132 = load double, double* %131
%133 = tail call double @llvm.fmuladd.f64(double %130, double %132, double %127)
%134 = add nsw i64 %115, 3
%135 = getelementptr inbounds double, double* %76, i64 %134
%136 = load double, double* %135
%137 = getelementptr inbounds double, double* %77, i64 %134
%138 = load double, double* %137
%139 = tail call double @llvm.fmuladd.f64(double %136, double %138, double %133)
%140 = add nsw i64 %115, 4
%141 = icmp eq i64 %140, %89
br i1 %141, label %112, label %114
}