void clusterCenter(double *x, double *c, int *assign, int N, int K, int D,
                   int *count) {
  for (int k = 0; k < K; ++k) {
    for (int d = 0; d < D; ++d) {
      c[k * D + d] = 0;
    }
    count[k] = 0;
  }
  for (int n = 0; n < N; ++n) {
    for (int d = 0; d < D; ++d) {
      c[assign[n] * D + d] += x[n * D + d];
    }
    count[assign[n]]++;
  }
  for (int k = 0; k < K; k++) {
    for (int d = 0; d < D; ++d) {
      c[k * D + d] /= count[k];
    }
  }
}


x86 assembly (clang, -O0):
<clusterCenter>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1189 <clusterCenter+0x89>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    116c <clusterCenter+0x6c>
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x24(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1139 <clusterCenter+0x39>
mov    0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1126 <clusterCenter+0x26>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    121d <clusterCenter+0x11d>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11f6 <clusterCenter+0xf6>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x24(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x30(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
imul   -0x24(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11a3 <clusterCenter+0xa3>
mov    0x10(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x30(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    %edx,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1190 <clusterCenter+0x90>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1290 <clusterCenter+0x190>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    127d <clusterCenter+0x17d>
mov    0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x24(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1237 <clusterCenter+0x137>
jmp    1282 <clusterCenter+0x182>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1224 <clusterCenter+0x124>
pop    %rbp
ret

x86 assembly (clang, -O2):
<clusterCenter>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %r9d,%ebx
mov    %rdx,0x40(%rsp)
mov    %rsi,0x8(%rsp)
mov    %rdi,%r15
mov    %r9d,%r13d
mov    %r8d,%r14d
mov    %r8d,0x1c(%rsp)
test   %r8d,%r8d
mov    %r14,0x38(%rsp)
mov    %rdi,0x20(%rsp)
jle    11ae <clusterCenter+0x9e>
lea    0x0(,%r13,8),%rax
mov    %rax,0x10(%rsp)
xor    %ebp,%ebp
cmpl   $0x1,0x1c(%rsp)
mov    %ecx,0x28(%rsp)
jne    1575 <clusterCenter+0x465>
xor    %r12d,%r12d
testb  $0x1,0x38(%rsp)
mov    0x20(%rsp),%r15
mov    0x28(%rsp),%ecx
je     11ae <clusterCenter+0x9e>
mov    %ecx,%ebp
test   %ebx,%ebx
jle    119c <clusterCenter+0x8c>
mov    %r12d,%eax
imul   %ebx,%eax
cltq
mov    0x8(%rsp),%rcx
lea    (%rcx,%rax,8),%rdi
xor    %esi,%esi
mov    0x10(%rsp),%rdx
call   1030 <memset@plt>
mov    0xb0(%rsp),%rax
movl   $0x0,(%rax,%r12,4)
mov    %ebp,%ecx
movslq %ebx,%r12
test   %ecx,%ecx
mov    0x8(%rsp),%rbp
mov    %r12,0x30(%rsp)
jle    142c <clusterCenter+0x31c>
mov    %ecx,%ecx
lea    0x0(,%r13,8),%rax
add    %rbp,%rax
mov    %rax,0x70(%rsp)
mov    %r13d,%eax
and    $0xfffffffc,%eax
mov    %rax,0x68(%rsp)
add    $0xfffffffffffffffc,%rax
mov    %rax,0x60(%rsp)
mov    %rax,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    %rdx,%rax
mov    %rdx,0x58(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,0x50(%rsp)
mov    %r13,%rax
neg    %rax
mov    %rax,0x10(%rsp)
lea    0x30(%r15),%rdi
lea    0x0(,%r12,8),%rsi
lea    0x30(%rbp),%rax
mov    %rax,0x48(%rsp)
lea    0x8(%rbp),%rax
mov    %rax,0x28(%rsp)
lea    0x8(%r15),%r8
xor    %r10d,%r10d
mov    0x20(%rsp),%r14
jmp    126b <clusterCenter+0x15b>
nopl   0x0(%rax)
mov    0x40(%rsp),%rax
mov    (%rax,%r10,4),%r9d
movslq %r9d,%rax
mov    0xb0(%rsp),%rdx
addl   $0x1,(%rdx,%rax,4)
add    $0x1,%r10
add    %rsi,%rdi
add    %rsi,%r8
cmp    %rcx,%r10
je     142c <clusterCenter+0x31c>
test   %ebx,%ebx
jle    1240 <clusterCenter+0x130>
mov    %r10,%r15
imul   %r12,%r15
mov    0x40(%rsp),%rax
mov    (%rax,%r10,4),%r9d
mov    %r9d,%eax
imul   %ebx,%eax
movslq %eax,%r11
cmp    $0x4,%ebx
jb     12bf <clusterCenter+0x1af>
lea    (%r15,%r13,1),%rax
lea    (%r14,%rax,8),%rax
lea    0x0(,%r11,8),%rdx
add    %rbp,%rdx
cmp    %rax,%rdx
jae    1341 <clusterCenter+0x231>
lea    (%r14,%r15,8),%rax
mov    0x70(%rsp),%rdx
lea    (%rdx,%r11,8),%rdx
cmp    %rdx,%rax
jae    1341 <clusterCenter+0x231>
xor    %edx,%edx
mov    %rdx,%rax
test   $0x1,%r13b
je     12ef <clusterCenter+0x1df>
add    %rdx,%r15
mov    0x20(%rsp),%rax
movsd  (%rax,%r15,8),%xmm0
lea    (%rdx,%r11,1),%rax
addsd  0x0(%rbp,%rax,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
mov    %rdx,%rax
or     $0x1,%rax
not    %rdx
cmp    0x10(%rsp),%rdx
je     1249 <clusterCenter+0x139>
mov    0x28(%rsp),%rdx
lea    (%rdx,%r11,8),%rdx
cs nopw 0x0(%rax,%rax,1)
movsd  -0x8(%r8,%rax,8),%xmm0
addsd  -0x8(%rdx,%rax,8),%xmm0
movsd  %xmm0,-0x8(%rdx,%rax,8)
movsd  (%r8,%rax,8),%xmm0
addsd  (%rdx,%rax,8),%xmm0
movsd  %xmm0,(%rdx,%rax,8)
add    $0x2,%rax
cmp    %rax,%r13
jne    1310 <clusterCenter+0x200>
jmp    1249 <clusterCenter+0x139>
cmpq   $0x0,0x60(%rsp)
je     13c3 <clusterCenter+0x2b3>
mov    0x48(%rsp),%rax
lea    (%rax,%r11,8),%rbp
mov    0x50(%rsp),%rdx
xor    %eax,%eax
nopl   0x0(%rax)
movupd -0x30(%rdi,%rax,8),%xmm0
movupd -0x20(%rdi,%rax,8),%xmm1
movupd -0x30(%rbp,%rax,8),%xmm2
addpd  %xmm0,%xmm2
movupd -0x20(%rbp,%rax,8),%xmm0
addpd  %xmm1,%xmm0
movupd -0x10(%rbp,%rax,8),%xmm1
movupd 0x0(%rbp,%rax,8),%xmm3
movupd %xmm2,-0x30(%rbp,%rax,8)
movupd %xmm0,-0x20(%rbp,%rax,8)
movupd -0x10(%rdi,%rax,8),%xmm0
addpd  %xmm1,%xmm0
movupd (%rdi,%rax,8),%xmm1
addpd  %xmm3,%xmm1
movupd %xmm0,-0x10(%rbp,%rax,8)
movupd %xmm1,0x0(%rbp,%rax,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%rdx
jne    1360 <clusterCenter+0x250>
jmp    13c5 <clusterCenter+0x2b5>
xor    %eax,%eax
testb  $0x1,0x58(%rsp)
mov    0x30(%rsp),%r12
je     1411 <clusterCenter+0x301>
lea    (%rax,%r15,1),%rdx
mov    %rcx,%rbp
mov    0x20(%rsp),%rcx
movupd (%rcx,%rdx,8),%xmm0
movupd 0x10(%rcx,%rdx,8),%xmm1
mov    %rbp,%rcx
add    %r11,%rax
mov    0x8(%rsp),%rdx
movupd (%rdx,%rax,8),%xmm2
addpd  %xmm0,%xmm2
movupd 0x10(%rdx,%rax,8),%xmm0
addpd  %xmm1,%xmm0
movupd %xmm2,(%rdx,%rax,8)
movupd %xmm0,0x10(%rdx,%rax,8)
mov    0x68(%rsp),%rax
mov    %rax,%rdx
cmp    %r13,%rax
mov    0x8(%rsp),%rbp
je     1249 <clusterCenter+0x139>
jmp    12c1 <clusterCenter+0x1b1>
cmpl   $0x0,0x1c(%rsp)
mov    0x38(%rsp),%r14
jle    1566 <clusterCenter+0x456>
mov    %r13d,%r11d
and    $0xfffffffe,%r11d
lea    -0x2(%r11),%r9
mov    %r9,%r10
shr    %r10
add    $0x1,%r10
mov    %r10,%r8
and    $0xfffffffffffffffe,%r8
mov    0x8(%rsp),%rdx
lea    0x10(%rdx),%rbp
mov    0x30(%rsp),%rax
lea    0x0(,%rax,8),%rdi
xor    %eax,%eax
jmp    1493 <clusterCenter+0x383>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%rax
add    %rdi,%rbp
add    %rdi,%rdx
cmp    %r14,%rax
je     1566 <clusterCenter+0x456>
test   %ebx,%ebx
jle    1480 <clusterCenter+0x370>
mov    0xb0(%rsp),%rcx
xorps  %xmm0,%xmm0
cvtsi2sdl (%rcx,%rax,4),%xmm0
cmp    $0x1,%ebx
jne    14c0 <clusterCenter+0x3b0>
xor    %ecx,%ecx
jmp    1540 <clusterCenter+0x430>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
test   %r9,%r9
je     155c <clusterCenter+0x44c>
mov    %r8,%rcx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
movupd -0x10(%rbp,%rsi,8),%xmm2
movupd 0x0(%rbp,%rsi,8),%xmm3
divpd  %xmm1,%xmm2
movupd %xmm2,-0x10(%rbp,%rsi,8)
divpd  %xmm1,%xmm3
movupd %xmm3,0x0(%rbp,%rsi,8)
add    $0x4,%rsi
add    $0xfffffffffffffffe,%rcx
jne    14e0 <clusterCenter+0x3d0>
test   $0x1,%r10b
je     152f <clusterCenter+0x41f>
mov    %rax,%rcx
imul   0x30(%rsp),%rcx
add    %rcx,%rsi
mov    0x8(%rsp),%rcx
movupd (%rcx,%rsi,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rcx,%rsi,8)
mov    %r11,%rcx
cmp    %r13,%r11
je     1480 <clusterCenter+0x370>
nopl   0x0(%rax,%rax,1)
movsd  (%rdx,%rcx,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%r13
jne    1540 <clusterCenter+0x430>
jmp    1480 <clusterCenter+0x370>
xor    %esi,%esi
test   $0x1,%r10b
jne    1510 <clusterCenter+0x400>
jmp    152f <clusterCenter+0x41f>
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
and    $0xfffffffe,%r14d
lea    (%rbx,%rbx,1),%r15d
xor    %r12d,%r12d
jmp    15b1 <clusterCenter+0x4a1>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0xb0(%rsp),%rax
movl   $0x0,0x4(%rax,%r12,4)
add    $0x2,%r12
add    %r15d,%ebp
cmp    %r12,%r14
je     1169 <clusterCenter+0x59>
test   %ebx,%ebx
jle    15cd <clusterCenter+0x4bd>
movslq %ebp,%rax
mov    0x8(%rsp),%rcx
lea    (%rcx,%rax,8),%rdi
xor    %esi,%esi
mov    0x10(%rsp),%rdx
call   1030 <memset@plt>
mov    0xb0(%rsp),%rax
movl   $0x0,(%rax,%r12,4)
test   %ebx,%ebx
jle    1590 <clusterCenter+0x480>
lea    (%rbx,%rbp,1),%eax
cltq
mov    0x8(%rsp),%rcx
lea    (%rcx,%rax,8),%rdi
xor    %esi,%esi
mov    0x10(%rsp),%rdx
call   1030 <memset@plt>
jmp    1590 <clusterCenter+0x480>

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

LLVM IR (-O2):
define void @clusterCenter(double* nocapture %0, double* nocapture %1, i32* nocapture %2, i32 %3, i32 %4, i32 %5, i32* nocapture %6) {
%8 = icmp sgt i32 %4, 0
br i1 %8, label %9, label %39
9: ; preds = %7
%10 = icmp sgt i32 %5, 0
%11 = zext i32 %5 to i64
%12 = shl nuw nsw i64 %11, 3
%13 = zext i32 %4 to i64
%14 = and i64 %13, 1
%15 = icmp eq i32 %4, 1
br i1 %15, label %27, label %16
16: ; preds = %9
%17 = and i64 %13, 4294967294
br label %18
18: ; preds = %70, %16
%19 = phi i64 [ 0, %16 ], [ %72, %70 ]
%20 = phi i64 [ 0, %16 ], [ %73, %70 ]
br i1 %10, label %21, label %61
21: ; preds = %18
%22 = trunc i64 %19 to i32
%23 = mul i32 %22, %5
%24 = sext i32 %23 to i64
%25 = getelementptr double, double* %1, i64 %24
%26 = bitcast double* %25 to i8*
br label %61
27: ; preds = %70, %9
%28 = phi i64 [ 0, %9 ], [ %72, %70 ]
%29 = icmp eq i64 %14, 0
br i1 %29, label %39, label %30
30: ; preds = %27
br i1 %10, label %31, label %37
31: ; preds = %30
%32 = trunc i64 %28 to i32
%33 = mul i32 %32, %5
%34 = sext i32 %33 to i64
%35 = getelementptr double, double* %1, i64 %34
%36 = bitcast double* %35 to i8*
br label %37
37: ; preds = %31, %30
%38 = getelementptr inbounds i32, i32* %6, i64 %28
store i32 0, i32* %38
br label %39
39: ; preds = %37, %27, %7
%40 = icmp sgt i32 %3, 0
br i1 %40, label %41, label %177
41: ; preds = %39
%42 = icmp sgt i32 %5, 0
%43 = sext i32 %5 to i64
%44 = zext i32 %3 to i64
%45 = zext i32 %5 to i64
%46 = getelementptr double, double* %1, i64 %45
%47 = and i64 %45, 4294967292
%48 = add nsw i64 %47, -4
%49 = lshr exact i64 %48, 2
%50 = add nuw nsw i64 %49, 1
%51 = icmp ult i32 %5, 4
%52 = and i64 %45, 4294967292
%53 = and i64 %50, 1
%54 = icmp eq i64 %48, 0
%55 = and i64 %50, 9223372036854775806
%56 = icmp eq i64 %53, 0
%57 = icmp eq i64 %52, %45
%58 = and i64 %45, 1
%59 = icmp eq i64 %58, 0
%60 = sub nsw i64 0, %45
br label %75
61: ; preds = %21, %18
%62 = getelementptr inbounds i32, i32* %6, i64 %19
store i32 0, i32* %62
%63 = or i64 %19, 1
br i1 %10, label %64, label %70
64: ; preds = %61
%65 = trunc i64 %63 to i32
%66 = mul i32 %65, %5
%67 = sext i32 %66 to i64
%68 = getelementptr double, double* %1, i64 %67
%69 = bitcast double* %68 to i8*
br label %70
70: ; preds = %64, %61
%71 = getelementptr inbounds i32, i32* %6, i64 %63
store i32 0, i32* %71
%72 = add nuw nsw i64 %19, 2
%73 = add i64 %20, 2
%74 = icmp eq i64 %73, %17
br i1 %74, label %27, label %18
75: ; preds = %41, %194
%76 = phi i64 [ 0, %41 ], [ %200, %194 ]
%77 = mul i64 %76, %43
%78 = getelementptr double, double* %0, i64 %77
%79 = add i64 %77, %45
%80 = getelementptr double, double* %0, i64 %79
br i1 %42, label %84, label %81
81: ; preds = %75
%82 = getelementptr inbounds i32, i32* %2, i64 %76
%83 = load i32, i32* %82
br label %194
84: ; preds = %75
%85 = mul nsw i64 %76, %43
%86 = getelementptr inbounds i32, i32* %2, i64 %76
%87 = load i32, i32* %86
%88 = mul nsw i32 %87, %5
%89 = sext i32 %88 to i64
br i1 %51, label %162, label %90
90: ; preds = %84
%91 = getelementptr double, double* %1, i64 %89
%92 = getelementptr double, double* %46, i64 %89
%93 = icmp ult double* %91, %80
%94 = icmp ult double* %78, %92
%95 = and i1 %93, %94
br i1 %95, label %162, label %96
96: ; preds = %90
br i1 %54, label %140, label %97
97: ; preds = %96, %97
%98 = phi i64 [ %137, %97 ], [ 0, %96 ]
%99 = phi i64 [ %138, %97 ], [ 0, %96 ]
%100 = add nsw i64 %98, %85
%101 = getelementptr inbounds double, double* %0, i64 %100
%102 = bitcast double* %101 to <2 x double>*
%103 = load <2 x double>, <2 x double>* %102
%104 = getelementptr inbounds double, double* %101, i64 2
%105 = bitcast double* %104 to <2 x double>*
%106 = load <2 x double>, <2 x double>* %105
%107 = add nsw i64 %98, %89
%108 = getelementptr inbounds double, double* %1, i64 %107
%109 = bitcast double* %108 to <2 x double>*
%110 = load <2 x double>, <2 x double>* %109
%111 = getelementptr inbounds double, double* %108, i64 2
%112 = bitcast double* %111 to <2 x double>*
%113 = load <2 x double>, <2 x double>* %112
%114 = fadd <2 x double> %103, %110
%115 = fadd <2 x double> %106, %113
%116 = bitcast double* %108 to <2 x double>*
store <2 x double> %114, <2 x double>* %116
%117 = bitcast double* %111 to <2 x double>*
store <2 x double> %115, <2 x double>* %117
%118 = or i64 %98, 4
%119 = add nsw i64 %118, %85
%120 = getelementptr inbounds double, double* %0, i64 %119
%121 = bitcast double* %120 to <2 x double>*
%122 = load <2 x double>, <2 x double>* %121
%123 = getelementptr inbounds double, double* %120, i64 2
%124 = bitcast double* %123 to <2 x double>*
%125 = load <2 x double>, <2 x double>* %124
%126 = add nsw i64 %118, %89
%127 = getelementptr inbounds double, double* %1, i64 %126
%128 = bitcast double* %127 to <2 x double>*
%129 = load <2 x double>, <2 x double>* %128
%130 = getelementptr inbounds double, double* %127, i64 2
%131 = bitcast double* %130 to <2 x double>*
%132 = load <2 x double>, <2 x double>* %131
%133 = fadd <2 x double> %122, %129
%134 = fadd <2 x double> %125, %132
%135 = bitcast double* %127 to <2 x double>*
store <2 x double> %133, <2 x double>* %135
%136 = bitcast double* %130 to <2 x double>*
store <2 x double> %134, <2 x double>* %136
%137 = add nuw i64 %98, 8
%138 = add i64 %99, 2
%139 = icmp eq i64 %138, %55
br i1 %139, label %140, label %97
140: ; preds = %97, %96
%141 = phi i64 [ 0, %96 ], [ %137, %97 ]
br i1 %56, label %161, label %142
142: ; preds = %140
%143 = add nsw i64 %141, %85
%144 = getelementptr inbounds double, double* %0, i64 %143
%145 = bitcast double* %144 to <2 x double>*
%146 = load <2 x double>, <2 x double>* %145
%147 = getelementptr inbounds double, double* %144, i64 2
%148 = bitcast double* %147 to <2 x double>*
%149 = load <2 x double>, <2 x double>* %148
%150 = add nsw i64 %141, %89
%151 = getelementptr inbounds double, double* %1, i64 %150
%152 = bitcast double* %151 to <2 x double>*
%153 = load <2 x double>, <2 x double>* %152
%154 = getelementptr inbounds double, double* %151, i64 2
%155 = bitcast double* %154 to <2 x double>*
%156 = load <2 x double>, <2 x double>* %155
%157 = fadd <2 x double> %146, %153
%158 = fadd <2 x double> %149, %156
%159 = bitcast double* %151 to <2 x double>*
store <2 x double> %157, <2 x double>* %159
%160 = bitcast double* %154 to <2 x double>*
store <2 x double> %158, <2 x double>* %160
br label %161
161: ; preds = %140, %142
br i1 %57, label %194, label %162
162: ; preds = %90, %84, %161
%163 = phi i64 [ 0, %90 ], [ 0, %84 ], [ %52, %161 ]
%164 = xor i64 %163, -1
br i1 %59, label %174, label %165
165: ; preds = %162
%166 = add nsw i64 %163, %85
%167 = getelementptr inbounds double, double* %0, i64 %166
%168 = load double, double* %167
%169 = add nsw i64 %163, %89
%170 = getelementptr inbounds double, double* %1, i64 %169
%171 = load double, double* %170
%172 = fadd double %168, %171
store double %172, double* %170
%173 = or i64 %163, 1
br label %174
174: ; preds = %165, %162
%175 = phi i64 [ %163, %162 ], [ %173, %165 ]
%176 = icmp eq i64 %164, %60
br i1 %176, label %194, label %202
177: ; preds = %194, %39
br i1 %8, label %178, label %262
178: ; preds = %177
%179 = icmp sgt i32 %5, 0
%180 = sext i32 %5 to i64
%181 = zext i32 %4 to i64
%182 = zext i32 %5 to i64
%183 = and i64 %182, 4294967294
%184 = add nsw i64 %183, -2
%185 = lshr exact i64 %184, 1
%186 = add nuw i64 %185, 1
%187 = icmp eq i32 %5, 1
%188 = and i64 %182, 4294967294
%189 = and i64 %186, 1
%190 = icmp eq i64 %184, 0
%191 = and i64 %186, -2
%192 = icmp eq i64 %189, 0
%193 = icmp eq i64 %188, %182
br label %221
194: ; preds = %174, %202, %161, %81
%195 = phi i32 [ %83, %81 ], [ %87, %161 ], [ %87, %202 ], [ %87, %174 ]
%196 = sext i32 %195 to i64
%197 = getelementptr inbounds i32, i32* %6, i64 %196
%198 = load i32, i32* %197
%199 = add nsw i32 %198, 1
store i32 %199, i32* %197
%200 = add nuw nsw i64 %76, 1
%201 = icmp eq i64 %200, %44
br i1 %201, label %177, label %75
202: ; preds = %174, %202
%203 = phi i64 [ %219, %202 ], [ %175, %174 ]
%204 = add nsw i64 %203, %85
%205 = getelementptr inbounds double, double* %0, i64 %204
%206 = load double, double* %205
%207 = add nsw i64 %203, %89
%208 = getelementptr inbounds double, double* %1, i64 %207
%209 = load double, double* %208
%210 = fadd double %206, %209
store double %210, double* %208
%211 = add nuw nsw i64 %203, 1
%212 = add nsw i64 %211, %85
%213 = getelementptr inbounds double, double* %0, i64 %212
%214 = load double, double* %213
%215 = add nsw i64 %211, %89
%216 = getelementptr inbounds double, double* %1, i64 %215
%217 = load double, double* %216
%218 = fadd double %214, %217
store double %218, double* %216
%219 = add nuw nsw i64 %203, 2
%220 = icmp eq i64 %219, %45
br i1 %220, label %194, label %202
221: ; preds = %178, %263
%222 = phi i64 [ 0, %178 ], [ %264, %263 ]
br i1 %179, label %223, label %263
223: ; preds = %221
%224 = getelementptr inbounds i32, i32* %6, i64 %222
%225 = load i32, i32* %224
%226 = sitofp i32 %225 to double
%227 = mul nsw i64 %222, %180
br i1 %187, label %260, label %228
228: ; preds = %223
%229 = insertelement <2 x double> poison, double %226, i64 0
%230 = shufflevector <2 x double> %229, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %190, label %250, label %231
231: ; preds = %228, %231
%232 = phi i64 [ %247, %231 ], [ 0, %228 ]
%233 = phi i64 [ %248, %231 ], [ 0, %228 ]
%234 = add nsw i64 %232, %227
%235 = getelementptr inbounds double, double* %1, i64 %234
%236 = bitcast double* %235 to <2 x double>*
%237 = load <2 x double>, <2 x double>* %236
%238 = fdiv <2 x double> %237, %230
%239 = bitcast double* %235 to <2 x double>*
store <2 x double> %238, <2 x double>* %239
%240 = or i64 %232, 2
%241 = add nsw i64 %240, %227
%242 = getelementptr inbounds double, double* %1, i64 %241
%243 = bitcast double* %242 to <2 x double>*
%244 = load <2 x double>, <2 x double>* %243
%245 = fdiv <2 x double> %244, %230
%246 = bitcast double* %242 to <2 x double>*
store <2 x double> %245, <2 x double>* %246
%247 = add nuw i64 %232, 4
%248 = add nuw i64 %233, 2
%249 = icmp eq i64 %248, %191
br i1 %249, label %250, label %231
250: ; preds = %231, %228
%251 = phi i64 [ 0, %228 ], [ %247, %231 ]
br i1 %192, label %259, label %252
252: ; preds = %250
%253 = add nsw i64 %251, %227
%254 = getelementptr inbounds double, double* %1, i64 %253
%255 = bitcast double* %254 to <2 x double>*
%256 = load <2 x double>, <2 x double>* %255
%257 = fdiv <2 x double> %256, %230
%258 = bitcast double* %254 to <2 x double>*
store <2 x double> %257, <2 x double>* %258
br label %259
259: ; preds = %250, %252
br i1 %193, label %263, label %260
260: ; preds = %223, %259
%261 = phi i64 [ 0, %223 ], [ %188, %259 ]
br label %266
262: ; preds = %263, %177
ret void
263: ; preds = %266, %259, %221
%264 = add nuw nsw i64 %222, 1
%265 = icmp eq i64 %264, %181
br i1 %265, label %262, label %221
266: ; preds = %260, %266
%267 = phi i64 [ %272, %266 ], [ %261, %260 ]
%268 = add nsw i64 %267, %227
%269 = getelementptr inbounds double, double* %1, i64 %268
%270 = load double, double* %269
%271 = fdiv double %270, %226
store double %271, double* %269
%272 = add nuw nsw i64 %267, 1
%273 = icmp eq i64 %272, %182
br i1 %273, label %263, label %266
}