float **kmeans_clustering(float **feature, int nfeatures, int npoints,
                          int nclusters, float threshold, int *membership) {
  int i, j, n = 0;
  int loop = 0, temp;
  int *new_centers_len;
  float delta;
  float **clusters;
  float **new_centers;
  int *initial;
  int initial_points;
  int c = 0;
  if (nclusters > npoints)
    nclusters = npoints;
  clusters = (float **)malloc(nclusters * sizeof(float *));
  clusters[0] = (float *)malloc(nclusters * nfeatures * sizeof(float));
  for (i = 1; i < nclusters; i++)
    clusters[i] = clusters[i - 1] + nfeatures;
  initial = (int *)malloc(npoints * sizeof(int));
  for (i = 0; i < npoints; i++) {
    initial[i] = i;
  }
  initial_points = npoints;
  for (i = 0; i < nclusters && initial_points >= 0; i++) {
    for (j = 0; j < nfeatures; j++)
      clusters[i][j] = feature[initial[n]][j];
    temp = initial[n];
    initial[n] = initial[initial_points - 1];
    initial[initial_points - 1] = temp;
    initial_points--;
    n++;
  }
  for (i = 0; i < npoints; i++)
    membership[i] = -1;
  new_centers_len = (int *)calloc(nclusters, sizeof(int));
  new_centers = (float **)malloc(nclusters * sizeof(float *));
  new_centers[0] = (float *)calloc(nclusters * nfeatures, sizeof(float));
  for (i = 1; i < nclusters; i++)
    new_centers[i] = new_centers[i - 1] + nfeatures;
  do {
    delta = 0.0;
    delta =
        (float)kmeansCuda(feature, nfeatures, npoints, nclusters, membership,
                          clusters, new_centers_len, new_centers);
    for (i = 0; i < nclusters; i++) {
      for (j = 0; j < nfeatures; j++) {
        if (new_centers_len[i] > 0)
          clusters[i][j] = new_centers[i][j] / new_centers_len[i];
        new_centers[i][j] = 0.0;
      }
      new_centers_len[i] = 0;
    }
    c++;
  } while ((delta > threshold) && (loop++ < 500));
  printf("iterated %d times\n", c);
  free(new_centers[0]);
  free(new_centers);
  free(new_centers_len);
  return clusters;
}


x86 assembly (clang, -O0):
<kmeans_clustering>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
movss  %xmm0,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x0,-0x68(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x10(%rbp),%eax
jle    1198 <kmeans_clustering+0x48>
mov    -0x10(%rbp),%eax
mov    %eax,-0x14(%rbp)
movslq -0x14(%rbp),%rdi
shl    $0x3,%rdi
call   1070 <malloc@plt>
mov    %rax,-0x50(%rbp)
mov    -0x14(%rbp),%eax
imul   -0xc(%rbp),%eax
movslq %eax,%rdi
shl    $0x2,%rdi
call   1070 <malloc@plt>
mov    %rax,%rcx
mov    -0x50(%rbp),%rax
mov    %rcx,(%rax)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    120f <kmeans_clustering+0xbf>
mov    -0x50(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rdx
movslq -0xc(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdx
mov    -0x50(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11cd <kmeans_clustering+0x7d>
movslq -0x10(%rbp),%rdi
shl    $0x2,%rdi
call   1070 <malloc@plt>
mov    %rax,-0x60(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    124f <kmeans_clustering+0xff>
mov    -0x24(%rbp),%edx
mov    -0x60(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1227 <kmeans_clustering+0xd7>
mov    -0x10(%rbp),%eax
mov    %eax,-0x64(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%ecx
xor    %eax,%eax
cmp    -0x14(%rbp),%ecx
mov    %al,-0x69(%rbp)
jge    1277 <kmeans_clustering+0x127>
cmpl   $0x0,-0x64(%rbp)
setge  %al
mov    %al,-0x69(%rbp)
mov    -0x69(%rbp),%al
test   $0x1,%al
jne    1287 <kmeans_clustering+0x137>
jmp    1336 <kmeans_clustering+0x1e6>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12da <kmeans_clustering+0x18a>
mov    -0x8(%rbp),%rax
mov    -0x60(%rbp),%rcx
movslq -0x2c(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    128e <kmeans_clustering+0x13e>
mov    -0x60(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x34(%rbp)
mov    -0x60(%rbp),%rax
mov    -0x64(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x60(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x34(%rbp),%edx
mov    -0x60(%rbp),%rax
mov    -0x64(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x64(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x64(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    125c <kmeans_clustering+0x10c>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1366 <kmeans_clustering+0x216>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movl   $0xffffffff,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    133d <kmeans_clustering+0x1ed>
movslq -0x14(%rbp),%rdi
mov    $0x4,%esi
call   1060 <calloc@plt>
mov    %rax,-0x40(%rbp)
movslq -0x14(%rbp),%rdi
shl    $0x3,%rdi
call   1070 <malloc@plt>
mov    %rax,-0x58(%rbp)
mov    -0x14(%rbp),%eax
imul   -0xc(%rbp),%eax
movslq %eax,%rdi
mov    $0x4,%esi
call   1060 <calloc@plt>
mov    %rax,%rcx
mov    -0x58(%rbp),%rax
mov    %rcx,(%rax)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    13f0 <kmeans_clustering+0x2a0>
mov    -0x58(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rdx
movslq -0xc(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdx
mov    -0x58(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    13ae <kmeans_clustering+0x25e>
jmp    13f5 <kmeans_clustering+0x2a5>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x10(%rbp),%edx
mov    -0x14(%rbp),%ecx
mov    -0x20(%rbp),%r8
mov    -0x50(%rbp),%r9
mov    -0x40(%rbp),%r10
mov    -0x58(%rbp),%rax
mov    %r10,(%rsp)
mov    %rax,0x8(%rsp)
mov    $0x0,%al
call   1050 <kmeansCuda@plt>
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x44(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    14e9 <kmeans_clustering+0x399>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    14cc <kmeans_clustering+0x37c>
mov    -0x40(%rbp),%rax
movslq -0x24(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
jle    14a6 <kmeans_clustering+0x356>
mov    -0x58(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x24(%rbp),%rcx
cvtsi2ssl (%rax,%rcx,4),%xmm1
divss  %xmm1,%xmm0
mov    -0x50(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x58(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    144d <kmeans_clustering+0x2fd>
mov    -0x40(%rbp),%rax
movslq -0x24(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    143a <kmeans_clustering+0x2ea>
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
movss  -0x44(%rbp),%xmm0
xor    %eax,%eax
ucomiss -0x18(%rbp),%xmm0
mov    %al,-0x6a(%rbp)
jbe    151c <kmeans_clustering+0x3cc>
mov    -0x30(%rbp),%eax
mov    %eax,%ecx
add    $0x1,%ecx
mov    %ecx,-0x30(%rbp)
cmp    $0x1f4,%eax
setl   %al
mov    %al,-0x6a(%rbp)
mov    -0x6a(%rbp),%al
test   $0x1,%al
jne    13f5 <kmeans_clustering+0x2a5>
mov    -0x68(%rbp),%esi
lea    0xacf(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
mov    -0x58(%rbp),%rax
mov    (%rax),%rdi
call   1030 <free@plt>
mov    -0x58(%rbp),%rdi
call   1030 <free@plt>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
mov    -0x50(%rbp),%rax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<kmeans_clustering>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %r8,0x28(%rsp)
movd   %xmm0,0x24(%rsp)
mov    %ecx,%ebp
mov    %esi,%r14d
mov    %rdi,0x30(%rsp)
cmp    %edx,%ecx
mov    %edx,0xc(%rsp)
cmovg  %edx,%ebp
movslq %ebp,%r15
lea    0x0(,%r15,8),%rdi
mov    %rdi,0x50(%rsp)
call   1080 <malloc@plt>
mov    %rax,%rbx
mov    %ebp,0x4(%rsp)
mov    %r14d,0x8(%rsp)
imul   %r14d,%ebp
movslq %ebp,%rax
mov    %rax,0x48(%rsp)
lea    0x0(,%rax,4),%rdi
call   1080 <malloc@plt>
mov    %rax,(%rbx)
cmp    $0x2,%r15d
jl     12a0 <kmeans_clustering+0x140>
movslq 0x8(%rsp),%r9
mov    0x4(%rsp),%ecx
lea    -0x1(%rcx),%rdi
add    $0xfffffffffffffffe,%rcx
mov    %edi,%r10d
and    $0x7,%r10d
cmp    $0x7,%rcx
jae    1204 <kmeans_clustering+0xa4>
mov    $0x1,%esi
test   %r10,%r10
jne    127b <kmeans_clustering+0x11b>
jmp    12a0 <kmeans_clustering+0x140>
and    $0xfffffffffffffff8,%rdi
mov    %r9,%r8
shl    $0x5,%r8
lea    0x0(,%r9,4),%rcx
xor    %esi,%esi
nopl   0x0(%rax)
mov    %rax,%rdx
add    %rcx,%rax
mov    %rax,0x8(%rbx,%rsi,8)
add    %rcx,%rax
mov    %rax,0x10(%rbx,%rsi,8)
add    %rcx,%rax
mov    %rax,0x18(%rbx,%rsi,8)
add    %rcx,%rax
mov    %rax,0x20(%rbx,%rsi,8)
add    %rcx,%rax
mov    %rax,0x28(%rbx,%rsi,8)
add    %rcx,%rax
mov    %rax,0x30(%rbx,%rsi,8)
add    %rcx,%rax
mov    %rax,0x38(%rbx,%rsi,8)
add    %rcx,%rax
mov    %rax,0x40(%rbx,%rsi,8)
add    $0x8,%rsi
cmp    %rsi,%rdi
jne    1220 <kmeans_clustering+0xc0>
add    %r8,%rdx
add    $0x1,%rsi
mov    %rdx,%rax
test   %r10,%r10
je     12a0 <kmeans_clustering+0x140>
lea    (%rax,%r9,4),%rax
shl    $0x2,%r9
lea    (%rbx,%rsi,8),%rcx
xor    %edx,%edx
nopl   0x0(%rax)
mov    %rax,(%rcx,%rdx,8)
add    $0x1,%rdx
add    %r9,%rax
cmp    %rdx,%r10
jne    1290 <kmeans_clustering+0x130>
mov    0xc(%rsp),%ebp
movslq %ebp,%r14
lea    0x0(,%r14,4),%rdi
call   1080 <malloc@plt>
mov    %ebp,%edx
test   %ebp,%ebp
mov    0x8(%rsp),%ebp
jle    143c <kmeans_clustering+0x2dc>
cmpl   $0x8,0xc(%rsp)
jae    12d0 <kmeans_clustering+0x170>
xor    %ecx,%ecx
jmp    1430 <kmeans_clustering+0x2d0>
mov    %edx,%ecx
and    $0xfffffff8,%ecx
lea    -0x8(%rcx),%rsi
mov    %rsi,%rdi
shr    $0x3,%rdi
add    $0x1,%rdi
mov    %edi,%r8d
and    $0x3,%r8d
cmp    $0x18,%rsi
jae    1300 <kmeans_clustering+0x1a0>
movdqa 0xd07(%rip),%xmm0
xor    %esi,%esi
jmp    13d8 <kmeans_clustering+0x278>
and    $0xfffffffffffffffc,%rdi
movdqa 0xcf4(%rip),%xmm0
xor    %esi,%esi
movdqa 0xcf9(%rip),%xmm8
movdqa 0xd00(%rip),%xmm9
movdqa 0xd07(%rip),%xmm10
movdqa 0xd0f(%rip),%xmm4
movdqa 0xd17(%rip),%xmm5
movdqa 0xd1f(%rip),%xmm6
movdqa 0xd27(%rip),%xmm7
movdqa 0xd2f(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movdqa %xmm0,%xmm2
paddd  %xmm8,%xmm2
movdqu %xmm0,(%rax,%rsi,4)
movdqu %xmm2,0x10(%rax,%rsi,4)
movdqa %xmm0,%xmm2
paddd  %xmm9,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm10,%xmm3
movdqu %xmm2,0x20(%rax,%rsi,4)
movdqu %xmm3,0x30(%rax,%rsi,4)
movdqa %xmm0,%xmm2
paddd  %xmm4,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm5,%xmm3
movdqu %xmm2,0x40(%rax,%rsi,4)
movdqu %xmm3,0x50(%rax,%rsi,4)
movdqa %xmm0,%xmm2
paddd  %xmm6,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm7,%xmm3
movdqu %xmm2,0x60(%rax,%rsi,4)
movdqu %xmm3,0x70(%rax,%rsi,4)
add    $0x20,%rsi
paddd  %xmm1,%xmm0
add    $0xfffffffffffffffc,%rdi
jne    1360 <kmeans_clustering+0x200>
test   %r8,%r8
je     1420 <kmeans_clustering+0x2c0>
lea    (%rax,%rsi,4),%rsi
add    $0x10,%rsi
shl    $0x5,%r8
xor    %edi,%edi
movdqa 0xc1d(%rip),%xmm1
movdqa 0xc25(%rip),%xmm2
nopl   0x0(%rax,%rax,1)
movdqa %xmm0,%xmm3
paddd  %xmm1,%xmm3
movdqu %xmm0,-0x10(%rsi,%rdi,1)
movdqu %xmm3,(%rsi,%rdi,1)
paddd  %xmm2,%xmm0
add    $0x20,%rdi
cmp    %rdi,%r8
jne    1400 <kmeans_clustering+0x2a0>
cmp    %rdx,%rcx
je     143c <kmeans_clustering+0x2dc>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %ecx,(%rax,%rcx,4)
add    $0x1,%rcx
cmp    %rcx,%rdx
jne    1430 <kmeans_clustering+0x2d0>
mov    %rdx,0x38(%rsp)
mov    %ebp,%ecx
mov    %rcx,0x10(%rsp)
cmpl   $0x0,0x4(%rsp)
jle    16b6 <kmeans_clustering+0x556>
cmpl   $0x0,0xc(%rsp)
js     16b6 <kmeans_clustering+0x556>
mov    0x10(%rsp),%r10
mov    %r10d,%ecx
and    $0xfffffff8,%ecx
mov    %rcx,0x18(%rsp)
add    $0xfffffffffffffff8,%rcx
mov    %rcx,0x60(%rsp)
shr    $0x3,%rcx
add    $0x1,%rcx
mov    %ecx,%r13d
and    $0x3,%r13d
and    $0xfffffffffffffffc,%rcx
mov    %rcx,0x40(%rsp)
mov    %r10d,%edx
and    $0x3,%edx
mov    %r13,0x58(%rsp)
shl    $0x5,%r13
xor    %r11d,%r11d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    (%rax,%r11,4),%r12d
cmpl   $0x0,0x8(%rsp)
jle    168f <kmeans_clustering+0x52f>
movslq %r12d,%rcx
mov    0x30(%rsp),%rsi
mov    (%rsi,%rcx,8),%r8
mov    (%rbx,%r11,8),%r9
cmpl   $0x8,0x8(%rsp)
jae    14f0 <kmeans_clustering+0x390>
xor    %esi,%esi
mov    0x10(%rsp),%rbp
jmp    1600 <kmeans_clustering+0x4a0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x10(%rsp),%rbp
lea    (%r8,%rbp,4),%rcx
cmp    %rcx,%r9
jae    150e <kmeans_clustering+0x3ae>
lea    (%r9,%rbp,4),%rcx
cmp    %rcx,%r8
jae    150e <kmeans_clustering+0x3ae>
xor    %esi,%esi
jmp    1600 <kmeans_clustering+0x4a0>
cmpq   $0x18,0x60(%rsp)
jae    151d <kmeans_clustering+0x3bd>
xor    %esi,%esi
jmp    159c <kmeans_clustering+0x43c>
mov    0x40(%rsp),%rdi
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movups (%r8,%rsi,4),%xmm0
movups 0x10(%r8,%rsi,4),%xmm1
movups %xmm0,(%r9,%rsi,4)
movups %xmm1,0x10(%r9,%rsi,4)
movups 0x20(%r8,%rsi,4),%xmm0
movups 0x30(%r8,%rsi,4),%xmm1
movups %xmm0,0x20(%r9,%rsi,4)
movups %xmm1,0x30(%r9,%rsi,4)
movups 0x40(%r8,%rsi,4),%xmm0
movups 0x50(%r8,%rsi,4),%xmm1
movups %xmm0,0x40(%r9,%rsi,4)
movups %xmm1,0x50(%r9,%rsi,4)
movdqu 0x60(%r8,%rsi,4),%xmm0
movdqu 0x70(%r8,%rsi,4),%xmm1
movdqu %xmm0,0x60(%r9,%rsi,4)
movdqu %xmm1,0x70(%r9,%rsi,4)
add    $0x20,%rsi
add    $0xfffffffffffffffc,%rdi
jne    1530 <kmeans_clustering+0x3d0>
cmpq   $0x0,0x58(%rsp)
je     15df <kmeans_clustering+0x47f>
lea    (%r9,%rsi,4),%rdi
add    $0x10,%rdi
lea    (%r8,%rsi,4),%rsi
add    $0x10,%rsi
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
movdqu -0x10(%rsi,%rcx,1),%xmm0
movdqu (%rsi,%rcx,1),%xmm1
movdqu %xmm0,-0x10(%rdi,%rcx,1)
movdqu %xmm1,(%rdi,%rcx,1)
add    $0x20,%rcx
cmp    %rcx,%r13
jne    15c0 <kmeans_clustering+0x460>
mov    0x18(%rsp),%rcx
mov    %rcx,%rsi
mov    0x10(%rsp),%rbp
cmp    %rbp,%rcx
je     168f <kmeans_clustering+0x52f>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r15,%r10
mov    %rsi,%rdi
not    %rdi
add    %rbp,%rdi
test   %rdx,%rdx
je     1636 <kmeans_clustering+0x4d6>
mov    %rdx,%r15
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movd   (%r8,%rsi,4),%xmm0
movd   %xmm0,(%r9,%rsi,4)
add    $0x1,%rsi
add    $0xffffffffffffffff,%r15
jne    1620 <kmeans_clustering+0x4c0>
cmp    $0x3,%rdi
mov    0x10(%rsp),%rcx
mov    %r10,%r15
jb     168f <kmeans_clustering+0x52f>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%r8,%rsi,4),%xmm0
movss  %xmm0,(%r9,%rsi,4)
movss  0x4(%r8,%rsi,4),%xmm0
movss  %xmm0,0x4(%r9,%rsi,4)
movss  0x8(%r8,%rsi,4),%xmm0
movss  %xmm0,0x8(%r9,%rsi,4)
movd   0xc(%r8,%rsi,4),%xmm0
movd   %xmm0,0xc(%r9,%rsi,4)
add    $0x4,%rsi
cmp    %rsi,%rcx
jne    1650 <kmeans_clustering+0x4f0>
mov    -0x4(%rax,%r14,4),%ecx
mov    %ecx,(%rax,%r11,4)
mov    %r12d,-0x4(%rax,%r14,4)
add    $0x1,%r11
cmp    %r15,%r11
jge    16b6 <kmeans_clustering+0x556>
lea    -0x1(%r14),%rcx
test   %r14,%r14
mov    %rcx,%r14
jg     14b0 <kmeans_clustering+0x350>
cmpl   $0x0,0xc(%rsp)
jle    16d5 <kmeans_clustering+0x575>
mov    0x38(%rsp),%rdx
shl    $0x2,%rdx
mov    0x28(%rsp),%rdi
mov    $0xff,%esi
call   1060 <memset@plt>
mov    $0x4,%esi
mov    %r15,%rdi
call   1070 <calloc@plt>
mov    %rax,%r14
mov    0x50(%rsp),%rdi
call   1080 <malloc@plt>
mov    %rax,%r15
mov    $0x4,%esi
mov    0x48(%rsp),%rdi
call   1070 <calloc@plt>
mov    %rax,(%r15)
cmpl   $0x2,0x4(%rsp)
mov    0x8(%rsp),%r11d
jl     17e0 <kmeans_clustering+0x680>
movslq %r11d,%r9
mov    0x4(%rsp),%ecx
lea    -0x1(%rcx),%rdi
add    $0xfffffffffffffffe,%rcx
mov    %edi,%r10d
and    $0x7,%r10d
cmp    $0x7,%rcx
jae    1743 <kmeans_clustering+0x5e3>
mov    $0x1,%esi
test   %r10,%r10
jne    17bb <kmeans_clustering+0x65b>
jmp    17e0 <kmeans_clustering+0x680>
and    $0xfffffffffffffff8,%rdi
mov    %r9,%r8
shl    $0x5,%r8
lea    0x0(,%r9,4),%rcx
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
mov    %rax,%rdx
add    %rcx,%rax
mov    %rax,0x8(%r15,%rsi,8)
add    %rcx,%rax
mov    %rax,0x10(%r15,%rsi,8)
add    %rcx,%rax
mov    %rax,0x18(%r15,%rsi,8)
add    %rcx,%rax
mov    %rax,0x20(%r15,%rsi,8)
add    %rcx,%rax
mov    %rax,0x28(%r15,%rsi,8)
add    %rcx,%rax
mov    %rax,0x30(%r15,%rsi,8)
add    %rcx,%rax
mov    %rax,0x38(%r15,%rsi,8)
add    %rcx,%rax
mov    %rax,0x40(%r15,%rsi,8)
add    $0x8,%rsi
cmp    %rsi,%rdi
jne    1760 <kmeans_clustering+0x600>
add    %r8,%rdx
add    $0x1,%rsi
mov    %rdx,%rax
test   %r10,%r10
je     17e0 <kmeans_clustering+0x680>
lea    (%rax,%r9,4),%rax
shl    $0x2,%r9
lea    (%r15,%rsi,8),%rcx
xor    %edx,%edx
nopl   0x0(%rax)
mov    %rax,(%rcx,%rdx,8)
add    $0x1,%rdx
add    %r9,%rax
cmp    %rdx,%r10
jne    17d0 <kmeans_clustering+0x670>
mov    0x4(%rsp),%r12d
mov    0x10(%rsp),%r13
and    $0xfffffffe,%r13d
xor    %eax,%eax
mov    %rax,0x18(%rsp)
mov    0x30(%rsp),%rdi
mov    %r11d,%esi
mov    0xc(%rsp),%edx
mov    0x4(%rsp),%ebp
mov    %ebp,%ecx
mov    0x28(%rsp),%r8
mov    %rbx,%r9
xor    %eax,%eax
push   %r15
push   %r14
call   1050 <kmeansCuda@plt>
mov    0x20(%rsp),%r8
mov    0x18(%rsp),%r11d
add    $0x10,%rsp
test   %ebp,%ebp
jle    18f0 <kmeans_clustering+0x790>
xor    %ecx,%ecx
jmp    185c <kmeans_clustering+0x6fc>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movl   $0x0,(%rdx,%rsi,4)
movl   $0x0,(%r14,%rcx,4)
add    $0x1,%rcx
cmp    %r12,%rcx
je     18f0 <kmeans_clustering+0x790>
test   %r11d,%r11d
jle    1847 <kmeans_clustering+0x6e7>
mov    (%r14,%rcx,4),%edi
xorps  %xmm0,%xmm0
cvtsi2ss %edi,%xmm0
mov    (%r15,%rcx,8),%rdx
xor    %esi,%esi
cmp    $0x1,%r11d
jne    18b1 <kmeans_clustering+0x751>
test   $0x1,%r8b
je     1847 <kmeans_clustering+0x6e7>
test   %edi,%edi
jle    1840 <kmeans_clustering+0x6e0>
movss  (%rdx,%rsi,4),%xmm1
divss  %xmm0,%xmm1
mov    (%rbx,%rcx,8),%rdi
movss  %xmm1,(%rdi,%rsi,4)
jmp    1840 <kmeans_clustering+0x6e0>
cs nopw 0x0(%rax,%rax,1)
movl   $0x0,0x4(%rdx,%rsi,4)
add    $0x2,%rsi
cmp    %rsi,%r13
je     1878 <kmeans_clustering+0x718>
test   %edi,%edi
jle    18c8 <kmeans_clustering+0x768>
movss  (%rdx,%rsi,4),%xmm1
divss  %xmm0,%xmm1
mov    (%rbx,%rcx,8),%rbp
movss  %xmm1,0x0(%rbp,%rsi,4)
movl   $0x0,(%rdx,%rsi,4)
test   %edi,%edi
jle    18a0 <kmeans_clustering+0x740>
movss  0x4(%rdx,%rsi,4),%xmm1
divss  %xmm0,%xmm1
mov    (%rbx,%rcx,8),%rbp
movss  %xmm1,0x4(%rbp,%rsi,4)
jmp    18a0 <kmeans_clustering+0x740>
nopl   0x0(%rax)
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
mov    0x18(%rsp),%rax
lea    0x1(%rax),%esi
ucomiss 0x24(%rsp),%xmm0
jbe    1913 <kmeans_clustering+0x7b3>
cmp    $0x1f4,%eax
mov    %esi,%eax
jb     17f0 <kmeans_clustering+0x690>
lea    0x776(%rip),%rdi
xor    %eax,%eax
call   1040 <printf@plt>
mov    (%r15),%rdi
call   1030 <free@plt>
mov    %r15,%rdi
call   1030 <free@plt>
mov    %r14,%rdi
call   1030 <free@plt>
mov    %rbx,%rax
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [19 x i8] c"iterated %d times\0A\00"
define float** @kmeans_clustering(float** %0, i32 %1, i32 %2, i32 %3, float %4, i32* %5) {
%7 = alloca float**
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca float
%12 = alloca i32*
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32*
%19 = alloca float
%20 = alloca float**
%21 = alloca float**
%22 = alloca i32*
%23 = alloca i32
%24 = alloca i32
store float** %0, float*** %7
store i32 %1, i32* %8
store i32 %2, i32* %9
store i32 %3, i32* %10
store float %4, float* %11
store i32* %5, i32** %12
store i32 0, i32* %15
store i32 0, i32* %16
store i32 0, i32* %24
%25 = load i32, i32* %10
%26 = load i32, i32* %9
%27 = icmp sgt i32 %25, %26
br i1 %27, label %28, label %30
28: ; preds = %6
%29 = load i32, i32* %9
store i32 %29, i32* %10
br label %30
30: ; preds = %28, %6
%31 = load i32, i32* %10
%32 = sext i32 %31 to i64
%33 = mul i64 %32, 8
%34 = call i8* @malloc(i64 %33)
%35 = bitcast i8* %34 to float**
store float** %35, float*** %20
%36 = load i32, i32* %10
%37 = load i32, i32* %8
%38 = mul nsw i32 %36, %37
%39 = sext i32 %38 to i64
%40 = mul i64 %39, 4
%41 = call i8* @malloc(i64 %40)
%42 = bitcast i8* %41 to float*
%43 = load float**, float*** %20
%44 = getelementptr inbounds float*, float** %43, i64 0
store float* %42, float** %44
store i32 1, i32* %13
br label %45
45: ; preds = %63, %30
%46 = load i32, i32* %13
%47 = load i32, i32* %10
%48 = icmp slt i32 %46, %47
br i1 %48, label %49, label %66
49: ; preds = %45
%50 = load float**, float*** %20
%51 = load i32, i32* %13
%52 = sub nsw i32 %51, 1
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds float*, float** %50, i64 %53
%55 = load float*, float** %54
%56 = load i32, i32* %8
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds float, float* %55, i64 %57
%59 = load float**, float*** %20
%60 = load i32, i32* %13
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds float*, float** %59, i64 %61
store float* %58, float** %62
br label %63
63: ; preds = %49
%64 = load i32, i32* %13
%65 = add nsw i32 %64, 1
store i32 %65, i32* %13
br label %45
66: ; preds = %45
%67 = load i32, i32* %9
%68 = sext i32 %67 to i64
%69 = mul i64 %68, 4
%70 = call i8* @malloc(i64 %69)
%71 = bitcast i8* %70 to i32*
store i32* %71, i32** %22
store i32 0, i32* %13
br label %72
72: ; preds = %82, %66
%73 = load i32, i32* %13
%74 = load i32, i32* %9
%75 = icmp slt i32 %73, %74
br i1 %75, label %76, label %85
76: ; preds = %72
%77 = load i32, i32* %13
%78 = load i32*, i32** %22
%79 = load i32, i32* %13
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds i32, i32* %78, i64 %80
store i32 %77, i32* %81
br label %82
82: ; preds = %76
%83 = load i32, i32* %13
%84 = add nsw i32 %83, 1
store i32 %84, i32* %13
br label %72
85: ; preds = %72
%86 = load i32, i32* %9
store i32 %86, i32* %23
store i32 0, i32* %13
br label %87
87: ; preds = %152, %85
%88 = load i32, i32* %13
%89 = load i32, i32* %10
%90 = icmp slt i32 %88, %89
br i1 %90, label %91, label %94
91: ; preds = %87
%92 = load i32, i32* %23
%93 = icmp sge i32 %92, 0
br label %94
94: ; preds = %91, %87
%95 = phi i1 [ false, %87 ], [ %93, %91 ]
br i1 %95, label %96, label %155
96: ; preds = %94
store i32 0, i32* %14
br label %97
97: ; preds = %123, %96
%98 = load i32, i32* %14
%99 = load i32, i32* %8
%100 = icmp slt i32 %98, %99
br i1 %100, label %101, label %126
101: ; preds = %97
%102 = load float**, float*** %7
%103 = load i32*, i32** %22
%104 = load i32, i32* %15
%105 = sext i32 %104 to i64
%106 = getelementptr inbounds i32, i32* %103, i64 %105
%107 = load i32, i32* %106
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds float*, float** %102, i64 %108
%110 = load float*, float** %109
%111 = load i32, i32* %14
%112 = sext i32 %111 to i64
%113 = getelementptr inbounds float, float* %110, i64 %112
%114 = load float, float* %113
%115 = load float**, float*** %20
%116 = load i32, i32* %13
%117 = sext i32 %116 to i64
%118 = getelementptr inbounds float*, float** %115, i64 %117
%119 = load float*, float** %118
%120 = load i32, i32* %14
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds float, float* %119, i64 %121
store float %114, float* %122
br label %123
123: ; preds = %101
%124 = load i32, i32* %14
%125 = add nsw i32 %124, 1
store i32 %125, i32* %14
br label %97
126: ; preds = %97
%127 = load i32*, i32** %22
%128 = load i32, i32* %15
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds i32, i32* %127, i64 %129
%131 = load i32, i32* %130
store i32 %131, i32* %17
%132 = load i32*, i32** %22
%133 = load i32, i32* %23
%134 = sub nsw i32 %133, 1
%135 = sext i32 %134 to i64
%136 = getelementptr inbounds i32, i32* %132, i64 %135
%137 = load i32, i32* %136
%138 = load i32*, i32** %22
%139 = load i32, i32* %15
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds i32, i32* %138, i64 %140
store i32 %137, i32* %141
%142 = load i32, i32* %17
%143 = load i32*, i32** %22
%144 = load i32, i32* %23
%145 = sub nsw i32 %144, 1
%146 = sext i32 %145 to i64
%147 = getelementptr inbounds i32, i32* %143, i64 %146
store i32 %142, i32* %147
%148 = load i32, i32* %23
%149 = add nsw i32 %148, -1
store i32 %149, i32* %23
%150 = load i32, i32* %15
%151 = add nsw i32 %150, 1
store i32 %151, i32* %15
br label %152
152: ; preds = %126
%153 = load i32, i32* %13
%154 = add nsw i32 %153, 1
store i32 %154, i32* %13
br label %87
155: ; preds = %94
store i32 0, i32* %13
br label %156
156: ; preds = %165, %155
%157 = load i32, i32* %13
%158 = load i32, i32* %9
%159 = icmp slt i32 %157, %158
br i1 %159, label %160, label %168
160: ; preds = %156
%161 = load i32*, i32** %12
%162 = load i32, i32* %13
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds i32, i32* %161, i64 %163
store i32 -1, i32* %164
br label %165
165: ; preds = %160
%166 = load i32, i32* %13
%167 = add nsw i32 %166, 1
store i32 %167, i32* %13
br label %156
168: ; preds = %156
%169 = load i32, i32* %10
%170 = sext i32 %169 to i64
%171 = call i8* @calloc(i64 %170, i64 4)
%172 = bitcast i8* %171 to i32*
store i32* %172, i32** %18
%173 = load i32, i32* %10
%174 = sext i32 %173 to i64
%175 = mul i64 %174, 8
%176 = call i8* @malloc(i64 %175)
%177 = bitcast i8* %176 to float**
store float** %177, float*** %21
%178 = load i32, i32* %10
%179 = load i32, i32* %8
%180 = mul nsw i32 %178, %179
%181 = sext i32 %180 to i64
%182 = call i8* @calloc(i64 %181, i64 4)
%183 = bitcast i8* %182 to float*
%184 = load float**, float*** %21
%185 = getelementptr inbounds float*, float** %184, i64 0
store float* %183, float** %185
store i32 1, i32* %13
br label %186
186: ; preds = %204, %168
%187 = load i32, i32* %13
%188 = load i32, i32* %10
%189 = icmp slt i32 %187, %188
br i1 %189, label %190, label %207
190: ; preds = %186
%191 = load float**, float*** %21
%192 = load i32, i32* %13
%193 = sub nsw i32 %192, 1
%194 = sext i32 %193 to i64
%195 = getelementptr inbounds float*, float** %191, i64 %194
%196 = load float*, float** %195
%197 = load i32, i32* %8
%198 = sext i32 %197 to i64
%199 = getelementptr inbounds float, float* %196, i64 %198
%200 = load float**, float*** %21
%201 = load i32, i32* %13
%202 = sext i32 %201 to i64
%203 = getelementptr inbounds float*, float** %200, i64 %202
store float* %199, float** %203
br label %204
204: ; preds = %190
%205 = load i32, i32* %13
%206 = add nsw i32 %205, 1
store i32 %206, i32* %13
br label %186
207: ; preds = %186
br label %208
208: ; preds = %291, %207
store float 0.000000e+00, float* %19
%209 = load float**, float*** %7
%210 = load i32, i32* %8
%211 = load i32, i32* %9
%212 = load i32, i32* %10
%213 = load i32*, i32** %12
%214 = load float**, float*** %20
%215 = load i32*, i32** %18
%216 = load float**, float*** %21
%217 = call i32 (float**, i32, i32, i32, i32*, float**, i32*, float**, ...) bitcast (i32 (...)* @kmeansCuda to i32 (float**, i32, i32, i32, i32*, float**, i32*, float**, ...)*)(float** %209, i32 %210, i32 %211, i32 %212, i32* %213, float** %214, i32* %215, float** %216)
%218 = sitofp i32 %217 to float
store float %218, float* %19
store i32 0, i32* %13
br label %219
219: ; preds = %277, %208
%220 = load i32, i32* %13
%221 = load i32, i32* %10
%222 = icmp slt i32 %220, %221
br i1 %222, label %223, label %280
223: ; preds = %219
store i32 0, i32* %14
br label %224
224: ; preds = %269, %223
%225 = load i32, i32* %14
%226 = load i32, i32* %8
%227 = icmp slt i32 %225, %226
br i1 %227, label %228, label %272
228: ; preds = %224
%229 = load i32*, i32** %18
%230 = load i32, i32* %13
%231 = sext i32 %230 to i64
%232 = getelementptr inbounds i32, i32* %229, i64 %231
%233 = load i32, i32* %232
%234 = icmp sgt i32 %233, 0
br i1 %234, label %235, label %260
235: ; preds = %228
%236 = load float**, float*** %21
%237 = load i32, i32* %13
%238 = sext i32 %237 to i64
%239 = getelementptr inbounds float*, float** %236, i64 %238
%240 = load float*, float** %239
%241 = load i32, i32* %14
%242 = sext i32 %241 to i64
%243 = getelementptr inbounds float, float* %240, i64 %242
%244 = load float, float* %243
%245 = load i32*, i32** %18
%246 = load i32, i32* %13
%247 = sext i32 %246 to i64
%248 = getelementptr inbounds i32, i32* %245, i64 %247
%249 = load i32, i32* %248
%250 = sitofp i32 %249 to float
%251 = fdiv float %244, %250
%252 = load float**, float*** %20
%253 = load i32, i32* %13
%254 = sext i32 %253 to i64
%255 = getelementptr inbounds float*, float** %252, i64 %254
%256 = load float*, float** %255
%257 = load i32, i32* %14
%258 = sext i32 %257 to i64
%259 = getelementptr inbounds float, float* %256, i64 %258
store float %251, float* %259
br label %260
260: ; preds = %235, %228
%261 = load float**, float*** %21
%262 = load i32, i32* %13
%263 = sext i32 %262 to i64
%264 = getelementptr inbounds float*, float** %261, i64 %263
%265 = load float*, float** %264
%266 = load i32, i32* %14
%267 = sext i32 %266 to i64
%268 = getelementptr inbounds float, float* %265, i64 %267
store float 0.000000e+00, float* %268
br label %269
269: ; preds = %260
%270 = load i32, i32* %14
%271 = add nsw i32 %270, 1
store i32 %271, i32* %14
br label %224
272: ; preds = %224
%273 = load i32*, i32** %18
%274 = load i32, i32* %13
%275 = sext i32 %274 to i64
%276 = getelementptr inbounds i32, i32* %273, i64 %275
store i32 0, i32* %276
br label %277
277: ; preds = %272
%278 = load i32, i32* %13
%279 = add nsw i32 %278, 1
store i32 %279, i32* %13
br label %219
280: ; preds = %219
%281 = load i32, i32* %24
%282 = add nsw i32 %281, 1
store i32 %282, i32* %24
br label %283
283: ; preds = %280
%284 = load float, float* %19
%285 = load float, float* %11
%286 = fcmp ogt float %284, %285
br i1 %286, label %287, label %291
287: ; preds = %283
%288 = load i32, i32* %16
%289 = add nsw i32 %288, 1
store i32 %289, i32* %16
%290 = icmp slt i32 %288, 500
br label %291
291: ; preds = %287, %283
%292 = phi i1 [ false, %283 ], [ %290, %287 ]
br i1 %292, label %208, label %293
293: ; preds = %291
%294 = load i32, i32* %24
%295 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([19 x i8], [19 x i8]* @.str, i64 0, i64 0), i32 %294)
%296 = load float**, float*** %21
%297 = getelementptr inbounds float*, float** %296, i64 0
%298 = load float*, float** %297
%299 = bitcast float* %298 to i8*
%300 = load float**, float*** %21
%301 = bitcast float** %300 to i8*
%302 = load i32*, i32** %18
%303 = bitcast i32* %302 to i8*
%304 = load float**, float*** %20
ret float** %304
}

LLVM IR (-O2):
@.str = [19 x i8] c"iterated %d times\0A\00"
define float** @kmeans_clustering(float** %0, i32 %1, i32 %2, i32 %3, float %4, i32* %5) {
%7 = bitcast i32* %5 to i8*
%8 = icmp sgt i32 %3, %2
%9 = select i1 %8, i32 %2, i32 %3
%10 = sext i32 %9 to i64
%11 = shl nsw i64 %10, 3
%12 = tail call i8* @malloc(i64 %11)
%13 = bitcast i8* %12 to float**
%14 = mul nsw i32 %9, %1
%15 = sext i32 %14 to i64
%16 = shl nsw i64 %15, 2
%17 = tail call i8* @malloc(i64 %16)
%18 = bitcast i8* %12 to i8**
store i8* %17, i8** %18
%19 = icmp sgt i32 %9, 1
br i1 %19, label %20, label %73
20: ; preds = %6
%21 = bitcast i8* %17 to float*
%22 = sext i32 %1 to i64
%23 = zext i32 %9 to i64
%24 = add nsw i64 %23, -1
%25 = add nsw i64 %23, -2
%26 = and i64 %24, 7
%27 = icmp ult i64 %25, 7
br i1 %27, label %60, label %28
28: ; preds = %20
%29 = and i64 %24, -8
br label %30
30: ; preds = %30, %28
%31 = phi float* [ %21, %28 ], [ %55, %30 ]
%32 = phi i64 [ 1, %28 ], [ %57, %30 ]
%33 = phi i64 [ 0, %28 ], [ %58, %30 ]
%34 = getelementptr inbounds float, float* %31, i64 %22
%35 = getelementptr inbounds float*, float** %13, i64 %32
store float* %34, float** %35
%36 = add nuw nsw i64 %32, 1
%37 = getelementptr inbounds float, float* %34, i64 %22
%38 = getelementptr inbounds float*, float** %13, i64 %36
store float* %37, float** %38
%39 = add nuw nsw i64 %32, 2
%40 = getelementptr inbounds float, float* %37, i64 %22
%41 = getelementptr inbounds float*, float** %13, i64 %39
store float* %40, float** %41
%42 = add nuw nsw i64 %32, 3
%43 = getelementptr inbounds float, float* %40, i64 %22
%44 = getelementptr inbounds float*, float** %13, i64 %42
store float* %43, float** %44
%45 = add nuw nsw i64 %32, 4
%46 = getelementptr inbounds float, float* %43, i64 %22
%47 = getelementptr inbounds float*, float** %13, i64 %45
store float* %46, float** %47
%48 = add nuw nsw i64 %32, 5
%49 = getelementptr inbounds float, float* %46, i64 %22
%50 = getelementptr inbounds float*, float** %13, i64 %48
store float* %49, float** %50
%51 = add nuw nsw i64 %32, 6
%52 = getelementptr inbounds float, float* %49, i64 %22
%53 = getelementptr inbounds float*, float** %13, i64 %51
store float* %52, float** %53
%54 = add nuw nsw i64 %32, 7
%55 = getelementptr inbounds float, float* %52, i64 %22
%56 = getelementptr inbounds float*, float** %13, i64 %54
store float* %55, float** %56
%57 = add nuw nsw i64 %32, 8
%58 = add i64 %33, 8
%59 = icmp eq i64 %58, %29
br i1 %59, label %60, label %30
60: ; preds = %30, %20
%61 = phi float* [ %21, %20 ], [ %55, %30 ]
%62 = phi i64 [ 1, %20 ], [ %57, %30 ]
%63 = icmp eq i64 %26, 0
br i1 %63, label %73, label %64
64: ; preds = %60, %64
%65 = phi float* [ %68, %64 ], [ %61, %60 ]
%66 = phi i64 [ %70, %64 ], [ %62, %60 ]
%67 = phi i64 [ %71, %64 ], [ 0, %60 ]
%68 = getelementptr inbounds float, float* %65, i64 %22
%69 = getelementptr inbounds float*, float** %13, i64 %66
store float* %68, float** %69
%70 = add nuw nsw i64 %66, 1
%71 = add i64 %67, 1
%72 = icmp eq i64 %71, %26
br i1 %72, label %73, label %64
73: ; preds = %60, %64, %6
%74 = sext i32 %2 to i64
%75 = shl nsw i64 %74, 2
%76 = tail call i8* @malloc(i64 %75)
%77 = bitcast i8* %76 to i32*
%78 = icmp sgt i32 %2, 0
br i1 %78, label %79, label %146
79: ; preds = %73
%80 = zext i32 %2 to i64
%81 = icmp ult i32 %2, 8
br i1 %81, label %144, label %82
82: ; preds = %79
%83 = and i64 %80, 4294967288
%84 = add nsw i64 %83, -8
%85 = lshr exact i64 %84, 3
%86 = add nuw nsw i64 %85, 1
%87 = and i64 %86, 3
%88 = icmp ult i64 %84, 24
br i1 %88, label %125, label %89
89: ; preds = %82
%90 = and i64 %86, 4611686018427387900
br label %91
91: ; preds = %91, %89
%92 = phi i64 [ 0, %89 ], [ %121, %91 ]
%93 = phi <4 x i32> [ <i32 0, i32 1, i32 2, i32 3>, %89 ], [ %122, %91 ]
%94 = phi i64 [ 0, %89 ], [ %123, %91 ]
%95 = add <4 x i32> %93, <i32 4, i32 4, i32 4, i32 4>
%96 = getelementptr inbounds i32, i32* %77, i64 %92
%97 = bitcast i32* %96 to <4 x i32>*
store <4 x i32> %93, <4 x i32>* %97
%98 = getelementptr inbounds i32, i32* %96, i64 4
%99 = bitcast i32* %98 to <4 x i32>*
store <4 x i32> %95, <4 x i32>* %99
%100 = or i64 %92, 8
%101 = add <4 x i32> %93, <i32 8, i32 8, i32 8, i32 8>
%102 = add <4 x i32> %93, <i32 12, i32 12, i32 12, i32 12>
%103 = getelementptr inbounds i32, i32* %77, i64 %100
%104 = bitcast i32* %103 to <4 x i32>*
store <4 x i32> %101, <4 x i32>* %104
%105 = getelementptr inbounds i32, i32* %103, i64 4
%106 = bitcast i32* %105 to <4 x i32>*
store <4 x i32> %102, <4 x i32>* %106
%107 = or i64 %92, 16
%108 = add <4 x i32> %93, <i32 16, i32 16, i32 16, i32 16>
%109 = add <4 x i32> %93, <i32 20, i32 20, i32 20, i32 20>
%110 = getelementptr inbounds i32, i32* %77, i64 %107
%111 = bitcast i32* %110 to <4 x i32>*
store <4 x i32> %108, <4 x i32>* %111
%112 = getelementptr inbounds i32, i32* %110, i64 4
%113 = bitcast i32* %112 to <4 x i32>*
store <4 x i32> %109, <4 x i32>* %113
%114 = or i64 %92, 24
%115 = add <4 x i32> %93, <i32 24, i32 24, i32 24, i32 24>
%116 = add <4 x i32> %93, <i32 28, i32 28, i32 28, i32 28>
%117 = getelementptr inbounds i32, i32* %77, i64 %114
%118 = bitcast i32* %117 to <4 x i32>*
store <4 x i32> %115, <4 x i32>* %118
%119 = getelementptr inbounds i32, i32* %117, i64 4
%120 = bitcast i32* %119 to <4 x i32>*
store <4 x i32> %116, <4 x i32>* %120
%121 = add nuw i64 %92, 32
%122 = add <4 x i32> %93, <i32 32, i32 32, i32 32, i32 32>
%123 = add i64 %94, 4
%124 = icmp eq i64 %123, %90
br i1 %124, label %125, label %91
125: ; preds = %91, %82
%126 = phi i64 [ 0, %82 ], [ %121, %91 ]
%127 = phi <4 x i32> [ <i32 0, i32 1, i32 2, i32 3>, %82 ], [ %122, %91 ]
%128 = icmp eq i64 %87, 0
br i1 %128, label %142, label %129
129: ; preds = %125, %129
%130 = phi i64 [ %138, %129 ], [ %126, %125 ]
%131 = phi <4 x i32> [ %139, %129 ], [ %127, %125 ]
%132 = phi i64 [ %140, %129 ], [ 0, %125 ]
%133 = add <4 x i32> %131, <i32 4, i32 4, i32 4, i32 4>
%134 = getelementptr inbounds i32, i32* %77, i64 %130
%135 = bitcast i32* %134 to <4 x i32>*
store <4 x i32> %131, <4 x i32>* %135
%136 = getelementptr inbounds i32, i32* %134, i64 4
%137 = bitcast i32* %136 to <4 x i32>*
store <4 x i32> %133, <4 x i32>* %137
%138 = add nuw i64 %130, 8
%139 = add <4 x i32> %131, <i32 8, i32 8, i32 8, i32 8>
%140 = add i64 %132, 1
%141 = icmp eq i64 %140, %87
br i1 %141, label %142, label %129
142: ; preds = %129, %125
%143 = icmp eq i64 %83, %80
br i1 %143, label %146, label %144
144: ; preds = %79, %142
%145 = phi i64 [ 0, %79 ], [ %83, %142 ]
br label %166
146: ; preds = %166, %142, %73
%147 = icmp sgt i32 %9, 0
%148 = icmp sgt i32 %2, -1
%149 = and i1 %147, %148
br i1 %149, label %150, label %274
150: ; preds = %146
%151 = icmp sgt i32 %1, 0
%152 = zext i32 %1 to i64
%153 = and i64 %152, 4294967288
%154 = add nsw i64 %153, -8
%155 = lshr exact i64 %154, 3
%156 = add nuw nsw i64 %155, 1
%157 = icmp ult i32 %1, 8
%158 = and i64 %152, 4294967288
%159 = and i64 %156, 3
%160 = icmp ult i64 %154, 24
%161 = and i64 %156, 4611686018427387900
%162 = icmp eq i64 %159, 0
%163 = icmp eq i64 %158, %152
%164 = and i64 %152, 3
%165 = icmp eq i64 %164, 0
br label %172
166: ; preds = %144, %166
%167 = phi i64 [ %170, %166 ], [ %145, %144 ]
%168 = getelementptr inbounds i32, i32* %77, i64 %167
%169 = trunc i64 %167 to i32
store i32 %169, i32* %168
%170 = add nuw nsw i64 %167, 1
%171 = icmp eq i64 %170, %80
br i1 %171, label %146, label %166
172: ; preds = %150, %297
%173 = phi i64 [ 0, %150 ], [ %302, %297 ]
%174 = phi i64 [ %74, %150 ], [ %299, %297 ]
%175 = getelementptr inbounds i32, i32* %77, i64 %173
%176 = load i32, i32* %175
br i1 %151, label %177, label %297
177: ; preds = %172
%178 = sext i32 %176 to i64
%179 = getelementptr inbounds float*, float** %0, i64 %178
%180 = load float*, float** %179
%181 = getelementptr inbounds float*, float** %13, i64 %173
%182 = load float*, float** %181
br i1 %157, label %258, label %183
183: ; preds = %177
%184 = getelementptr float, float* %182, i64 %152
%185 = getelementptr float, float* %180, i64 %152
%186 = icmp ult float* %182, %185
%187 = icmp ult float* %180, %184
%188 = and i1 %186, %187
br i1 %188, label %258, label %189
189: ; preds = %183
br i1 %160, label %239, label %190
190: ; preds = %189, %190
%191 = phi i64 [ %236, %190 ], [ 0, %189 ]
%192 = phi i64 [ %237, %190 ], [ 0, %189 ]
%193 = getelementptr inbounds float, float* %180, i64 %191
%194 = bitcast float* %193 to <4 x float>*
%195 = load <4 x float>, <4 x float>* %194
%196 = getelementptr inbounds float, float* %193, i64 4
%197 = bitcast float* %196 to <4 x float>*
%198 = load <4 x float>, <4 x float>* %197
%199 = getelementptr inbounds float, float* %182, i64 %191
%200 = bitcast float* %199 to <4 x float>*
store <4 x float> %195, <4 x float>* %200
%201 = getelementptr inbounds float, float* %199, i64 4
%202 = bitcast float* %201 to <4 x float>*
store <4 x float> %198, <4 x float>* %202
%203 = or i64 %191, 8
%204 = getelementptr inbounds float, float* %180, i64 %203
%205 = bitcast float* %204 to <4 x float>*
%206 = load <4 x float>, <4 x float>* %205
%207 = getelementptr inbounds float, float* %204, i64 4
%208 = bitcast float* %207 to <4 x float>*
%209 = load <4 x float>, <4 x float>* %208
%210 = getelementptr inbounds float, float* %182, i64 %203
%211 = bitcast float* %210 to <4 x float>*
store <4 x float> %206, <4 x float>* %211
%212 = getelementptr inbounds float, float* %210, i64 4
%213 = bitcast float* %212 to <4 x float>*
store <4 x float> %209, <4 x float>* %213
%214 = or i64 %191, 16
%215 = getelementptr inbounds float, float* %180, i64 %214
%216 = bitcast float* %215 to <4 x float>*
%217 = load <4 x float>, <4 x float>* %216
%218 = getelementptr inbounds float, float* %215, i64 4
%219 = bitcast float* %218 to <4 x float>*
%220 = load <4 x float>, <4 x float>* %219
%221 = getelementptr inbounds float, float* %182, i64 %214
%222 = bitcast float* %221 to <4 x float>*
store <4 x float> %217, <4 x float>* %222
%223 = getelementptr inbounds float, float* %221, i64 4
%224 = bitcast float* %223 to <4 x float>*
store <4 x float> %220, <4 x float>* %224
%225 = or i64 %191, 24
%226 = getelementptr inbounds float, float* %180, i64 %225
%227 = bitcast float* %226 to <4 x float>*
%228 = load <4 x float>, <4 x float>* %227
%229 = getelementptr inbounds float, float* %226, i64 4
%230 = bitcast float* %229 to <4 x float>*
%231 = load <4 x float>, <4 x float>* %230
%232 = getelementptr inbounds float, float* %182, i64 %225
%233 = bitcast float* %232 to <4 x float>*
store <4 x float> %228, <4 x float>* %233
%234 = getelementptr inbounds float, float* %232, i64 4
%235 = bitcast float* %234 to <4 x float>*
store <4 x float> %231, <4 x float>* %235
%236 = add nuw i64 %191, 32
%237 = add i64 %192, 4
%238 = icmp eq i64 %237, %161
br i1 %238, label %239, label %190
239: ; preds = %190, %189
%240 = phi i64 [ 0, %189 ], [ %236, %190 ]
br i1 %162, label %257, label %241
241: ; preds = %239, %241
%242 = phi i64 [ %254, %241 ], [ %240, %239 ]
%243 = phi i64 [ %255, %241 ], [ 0, %239 ]
%244 = getelementptr inbounds float, float* %180, i64 %242
%245 = bitcast float* %244 to <4 x float>*
%246 = load <4 x float>, <4 x float>* %245
%247 = getelementptr inbounds float, float* %244, i64 4
%248 = bitcast float* %247 to <4 x float>*
%249 = load <4 x float>, <4 x float>* %248
%250 = getelementptr inbounds float, float* %182, i64 %242
%251 = bitcast float* %250 to <4 x float>*
store <4 x float> %246, <4 x float>* %251
%252 = getelementptr inbounds float, float* %250, i64 4
%253 = bitcast float* %252 to <4 x float>*
store <4 x float> %249, <4 x float>* %253
%254 = add nuw i64 %242, 8
%255 = add i64 %243, 1
%256 = icmp eq i64 %255, %159
br i1 %256, label %257, label %241
257: ; preds = %241, %239
br i1 %163, label %297, label %258
258: ; preds = %183, %177, %257
%259 = phi i64 [ 0, %183 ], [ 0, %177 ], [ %158, %257 ]
%260 = xor i64 %259, -1
%261 = add nsw i64 %260, %152
br i1 %165, label %271, label %262
262: ; preds = %258, %262
%263 = phi i64 [ %268, %262 ], [ %259, %258 ]
%264 = phi i64 [ %269, %262 ], [ 0, %258 ]
%265 = getelementptr inbounds float, float* %180, i64 %263
%266 = load float, float* %265
%267 = getelementptr inbounds float, float* %182, i64 %263
store float %266, float* %267
%268 = add nuw nsw i64 %263, 1
%269 = add i64 %264, 1
%270 = icmp eq i64 %269, %164
br i1 %270, label %271, label %262
271: ; preds = %262, %258
%272 = phi i64 [ %259, %258 ], [ %268, %262 ]
%273 = icmp ult i64 %261, 3
br i1 %273, label %297, label %278
274: ; preds = %297, %146
br i1 %78, label %275, label %306
275: ; preds = %274
%276 = zext i32 %2 to i64
%277 = shl nuw nsw i64 %276, 2
br label %306
278: ; preds = %271, %278
%279 = phi i64 [ %295, %278 ], [ %272, %271 ]
%280 = getelementptr inbounds float, float* %180, i64 %279
%281 = load float, float* %280
%282 = getelementptr inbounds float, float* %182, i64 %279
store float %281, float* %282
%283 = add nuw nsw i64 %279, 1
%284 = getelementptr inbounds float, float* %180, i64 %283
%285 = load float, float* %284
%286 = getelementptr inbounds float, float* %182, i64 %283
store float %285, float* %286
%287 = add nuw nsw i64 %279, 2
%288 = getelementptr inbounds float, float* %180, i64 %287
%289 = load float, float* %288
%290 = getelementptr inbounds float, float* %182, i64 %287
store float %289, float* %290
%291 = add nuw nsw i64 %279, 3
%292 = getelementptr inbounds float, float* %180, i64 %291
%293 = load float, float* %292
%294 = getelementptr inbounds float, float* %182, i64 %291
store float %293, float* %294
%295 = add nuw nsw i64 %279, 4
%296 = icmp eq i64 %295, %152
br i1 %296, label %297, label %278
297: ; preds = %271, %278, %257, %172
%298 = getelementptr inbounds i32, i32* %77, i64 %173
%299 = add nsw i64 %174, -1
%300 = getelementptr inbounds i32, i32* %77, i64 %299
%301 = load i32, i32* %300
store i32 %301, i32* %298
store i32 %176, i32* %300
%302 = add nuw nsw i64 %173, 1
%303 = icmp slt i64 %302, %10
%304 = icmp sgt i64 %174, 0
%305 = select i1 %303, i1 %304, i1 false
br i1 %305, label %172, label %274
306: ; preds = %275, %274
%307 = tail call i8* @calloc(i64 %10, i64 4)
%308 = bitcast i8* %307 to i32*
%309 = tail call i8* @malloc(i64 %11)
%310 = bitcast i8* %309 to float**
%311 = tail call i8* @calloc(i64 %15, i64 4)
%312 = bitcast i8* %309 to i8**
store i8* %311, i8** %312
br i1 %19, label %313, label %336
313: ; preds = %306
%314 = bitcast i8* %311 to float*
%315 = sext i32 %1 to i64
%316 = zext i32 %9 to i64
%317 = add nsw i64 %316, -1
%318 = add nsw i64 %316, -2
%319 = and i64 %317, 7
%320 = icmp ult i64 %318, 7
br i1 %320, label %323, label %321
321: ; preds = %313
%322 = and i64 %317, -8
br label %344
323: ; preds = %344, %313
%324 = phi float* [ %314, %313 ], [ %369, %344 ]
%325 = phi i64 [ 1, %313 ], [ %371, %344 ]
%326 = icmp eq i64 %319, 0
br i1 %326, label %336, label %327
327: ; preds = %323, %327
%328 = phi float* [ %331, %327 ], [ %324, %323 ]
%329 = phi i64 [ %333, %327 ], [ %325, %323 ]
%330 = phi i64 [ %334, %327 ], [ 0, %323 ]
%331 = getelementptr inbounds float, float* %328, i64 %315
%332 = getelementptr inbounds float*, float** %310, i64 %329
store float* %331, float** %332
%333 = add nuw nsw i64 %329, 1
%334 = add i64 %330, 1
%335 = icmp eq i64 %334, %319
br i1 %335, label %336, label %327
336: ; preds = %323, %327, %306
%337 = icmp sgt i32 %1, 0
%338 = zext i32 %9 to i64
%339 = zext i32 %1 to i64
%340 = and i64 %339, 1
%341 = icmp eq i32 %1, 1
%342 = and i64 %339, 4294967294
%343 = icmp eq i64 %340, 0
br label %374
344: ; preds = %344, %321
%345 = phi float* [ %314, %321 ], [ %369, %344 ]
%346 = phi i64 [ 1, %321 ], [ %371, %344 ]
%347 = phi i64 [ 0, %321 ], [ %372, %344 ]
%348 = getelementptr inbounds float, float* %345, i64 %315
%349 = getelementptr inbounds float*, float** %310, i64 %346
store float* %348, float** %349
%350 = add nuw nsw i64 %346, 1
%351 = getelementptr inbounds float, float* %348, i64 %315
%352 = getelementptr inbounds float*, float** %310, i64 %350
store float* %351, float** %352
%353 = add nuw nsw i64 %346, 2
%354 = getelementptr inbounds float, float* %351, i64 %315
%355 = getelementptr inbounds float*, float** %310, i64 %353
store float* %354, float** %355
%356 = add nuw nsw i64 %346, 3
%357 = getelementptr inbounds float, float* %354, i64 %315
%358 = getelementptr inbounds float*, float** %310, i64 %356
store float* %357, float** %358
%359 = add nuw nsw i64 %346, 4
%360 = getelementptr inbounds float, float* %357, i64 %315
%361 = getelementptr inbounds float*, float** %310, i64 %359
store float* %360, float** %361
%362 = add nuw nsw i64 %346, 5
%363 = getelementptr inbounds float, float* %360, i64 %315
%364 = getelementptr inbounds float*, float** %310, i64 %362
store float* %363, float** %364
%365 = add nuw nsw i64 %346, 6
%366 = getelementptr inbounds float, float* %363, i64 %315
%367 = getelementptr inbounds float*, float** %310, i64 %365
store float* %366, float** %367
%368 = add nuw nsw i64 %346, 7
%369 = getelementptr inbounds float, float* %366, i64 %315
%370 = getelementptr inbounds float*, float** %310, i64 %368
store float* %369, float** %370
%371 = add nuw nsw i64 %346, 8
%372 = add i64 %347, 8
%373 = icmp eq i64 %372, %322
br i1 %373, label %323, label %344
374: ; preds = %336, %426
%375 = phi i32 [ %427, %426 ], [ 0, %336 ]
%376 = tail call i32 (float**, i32, i32, i32, i32*, float**, i32*, float**, ...) bitcast (i32 (...)* @kmeansCuda to i32 (float**, i32, i32, i32, i32*, float**, i32*, float**, ...)*)(float** %0, i32 %1, i32 %2, i32 %9, i32* %5, float** %13, i32* %308, float** %310)
%377 = sitofp i32 %376 to float
br i1 %147, label %378, label %426
378: ; preds = %374, %422
%379 = phi i64 [ %424, %422 ], [ 0, %374 ]
br i1 %337, label %380, label %422
380: ; preds = %378
%381 = getelementptr inbounds i32, i32* %308, i64 %379
%382 = load i32, i32* %381
%383 = icmp sgt i32 %382, 0
%384 = getelementptr inbounds float*, float** %310, i64 %379
%385 = sitofp i32 %382 to float
%386 = getelementptr inbounds float*, float** %13, i64 %379
%387 = load float*, float** %384
br i1 %341, label %411, label %388
388: ; preds = %380, %406
%389 = phi i64 [ %408, %406 ], [ 0, %380 ]
%390 = phi i64 [ %409, %406 ], [ 0, %380 ]
br i1 %383, label %391, label %397
391: ; preds = %388
%392 = getelementptr inbounds float, float* %387, i64 %389
%393 = load float, float* %392
%394 = fdiv float %393, %385
%395 = load float*, float** %386
%396 = getelementptr inbounds float, float* %395, i64 %389
store float %394, float* %396
br label %397
397: ; preds = %391, %388
%398 = getelementptr inbounds float, float* %387, i64 %389
store float 0.000000e+00, float* %398
%399 = or i64 %389, 1
br i1 %383, label %400, label %406
400: ; preds = %397
%401 = getelementptr inbounds float, float* %387, i64 %399
%402 = load float, float* %401
%403 = fdiv float %402, %385
%404 = load float*, float** %386
%405 = getelementptr inbounds float, float* %404, i64 %399
store float %403, float* %405
br label %406
406: ; preds = %400, %397
%407 = getelementptr inbounds float, float* %387, i64 %399
store float 0.000000e+00, float* %407
%408 = add nuw nsw i64 %389, 2
%409 = add i64 %390, 2
%410 = icmp eq i64 %409, %342
br i1 %410, label %411, label %388
411: ; preds = %406, %380
%412 = phi i64 [ 0, %380 ], [ %408, %406 ]
br i1 %343, label %422, label %413
413: ; preds = %411
br i1 %383, label %414, label %420
414: ; preds = %413
%415 = getelementptr inbounds float, float* %387, i64 %412
%416 = load float, float* %415
%417 = fdiv float %416, %385
%418 = load float*, float** %386
%419 = getelementptr inbounds float, float* %418, i64 %412
store float %417, float* %419
br label %420
420: ; preds = %414, %413
%421 = getelementptr inbounds float, float* %387, i64 %412
store float 0.000000e+00, float* %421
br label %422
422: ; preds = %420, %411, %378
%423 = getelementptr inbounds i32, i32* %308, i64 %379
store i32 0, i32* %423
%424 = add nuw nsw i64 %379, 1
%425 = icmp eq i64 %424, %338
br i1 %425, label %426, label %378
426: ; preds = %422, %374
%427 = add nuw nsw i32 %375, 1
%428 = fcmp ogt float %377, %4
%429 = icmp ult i32 %375, 500
%430 = select i1 %428, i1 %429, i1 false
br i1 %430, label %374, label %431
431: ; preds = %426
%432 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([19 x i8], [19 x i8]* @.str, i64 0, i64 0), i32 %427)
%433 = load i8*, i8** %312
tail call void @free(i8* %433)
tail call void @free(i8* %309)
tail call void @free(i8* %307)
ret float** %13
}