user: What is the c source code of the assembly code below:

```asm
<kmeans_it>:
  endbr64
  push   %r15
  movslq %r8d,%r15
  push   %r14
  push   %r13
  mov    %r9,%r13
  push   %r12
  movslq %ecx,%r12
  push   %rbp
  push   %rbx
  sub    $0x88,%rsp
  mov    %rdi,0x60(%rsp)
  lea    0x0(,%r12,4),%rdi
  mov    %edx,0x70(%rsp)
  mov    %ecx,0x2c(%rsp)
  mov    %esi,0x74(%rsp)
  mov    %r15d,0x34(%rsp)
  mov    %r12,0x78(%rsp)
  call   <malloc@plt>
  mov    $0x4,%esi
  mov    %r15,%rdi
  mov    %r15,0x18(%rsp)
  mov    %rax,%r14
  call   <calloc@plt>
  lea    0x0(,%r15,8),%rdi
  mov    %rax,0x48(%rsp)
  call   <malloc@plt>
  mov    %rax,0x50(%rsp)
  test   %r15d,%r15d
  jle    L264
  mov    %rax,%rcx
  mov    %rax,%rbp
  mov    %r15d,%eax
  lea    (%rcx,%rax,8),%rbx
  nopl   0x0(%rax,%rax,1)
L131:
  mov    $0x4,%esi
  mov    %r12,%rdi
  add    $0x8,%rbp
  call   <calloc@plt>
  mov    %rax,-0x8(%rbp)
  cmp    %rbp,%rbx
  jne    L131
  mov    0x74(%rsp),%esi
  test   %esi,%esi
  jle    L123
L98:
  mov    0x74(%rsp),%ebx
  mov    0x70(%rsp),%ecx
  imul   %ecx,%ebx
  test   %ecx,%ecx
  jle    L124
  mov    0x2c(%rsp),%eax
  movslq %ebx,%rbx
  mov    %ecx,0x68(%rsp)
  shl    $0x2,%rbx
  movl   $0x0,0x24(%rsp)
  sub    $0x1,%eax
  movl   $0x0,0x6c(%rsp)
  lea    0x4(,%rax,4),%rax
  mov    %rbx,0x8(%rsp)
  mov    %rax,0x10(%rsp)
  add    %r14,%rax
  mov    %rax,0x58(%rsp)
  nopl   (%rax)
L168:
  mov    0x68(%rsp),%eax
  sub    0x70(%rsp),%eax
  mov    %eax,0x28(%rsp)
  nopl   0x0(%rax)
L241:
  mov    0x2c(%rsp),%r8d
  test   %r8d,%r8d
  jle    L104
  mov    0x60(%rsp),%rax
  movslq 0x28(%rsp),%rsi
  mov    0x58(%rsp),%rdi
  mov    0x8(%rsp),%rbx
  lea    (%rax,%rsi,4),%rdx
  mov    %r14,%rax
  nopl   0x0(%rax)
L105:
  mov    (%rdx),%ecx
  add    $0x4,%rax
  add    %rbx,%rdx
  mov    %ecx,-0x4(%rax)
  cmp    %rax,%rdi
  jne    L105
  mov    0x34(%rsp),%eax
  test   %eax,%eax
  jle    L214
L78:
  mov    0x2c(%rsp),%edx
  mov    0x0(%r13),%rsi
  xor    %eax,%eax
  mov    %r14,%rdi
  mov    $0x1,%r15d
  call   <sqd_distance@plt>
  cmpl   $0x1,0x34(%rsp)
  movl   $0x0,0x24(%rsp)
  mov    %eax,%ebp
  je     L192
  mov    %r13,%rax
  mov    %ebp,%ebx
  mov    %r15,%r13
  mov    0x2c(%rsp),%ebp
  mov    %rax,%r15
  jmp    L125
  nop
L130:
  cmp    %ebx,%eax
  jl     L151
  add    $0x1,%r13
  cmp    %r13,0x18(%rsp)
  je     L108
L125:
  mov    (%r15,%r13,8),%rsi
  xor    %eax,%eax
  mov    %ebp,%edx
  mov    %r14,%rdi
  mov    %r13d,%r12d
  call   <sqd_distance@plt>
  test   %ebx,%ebx
  jns    L130
L151:
  mov    %r12d,0x24(%rsp)
  mov    %eax,%ebx
  add    $0x1,%r13
  cmp    %r13,0x18(%rsp)
  jne    L125
L108:
  mov    %r15,%r13
L192:
  movslq 0x24(%rsp),%rax
  mov    0x48(%rsp),%rcx
  mov    0x2c(%rsp),%r9d
  movslq 0x28(%rsp),%rsi
  lea    (%rcx,%rax,4),%rcx
  mov    (%rcx),%edi
  mov    %rcx,0x38(%rsp)
  mov    %edi,0x30(%rsp)
  test   %r9d,%r9d
  jle    L106
L128:
  mov    0x50(%rsp),%rdx
  lea    0x0(%r13,%rax,8),%r9
  mov    %r13,0x40(%rsp)
  xor    %r15d,%r15d
  mov    %r14,%r13
  mov    %r9,%r14
  mov    (%rdx,%rax,8),%rbp
  mov    0xc8(%rsp),%rax
  mov    %rbp,%rbx
  mov    0x30(%rsp),%ebp
  lea    (%rax,%rsi,4),%r12
  cs nopw 0x0(%rax,%rax,1)
L32:
  mov    (%r14),%rax
  mov    0x0(%r13,%r15,1),%esi
  pxor   %xmm0,%xmm0
  mov    %ebp,%edi
  cvtss2sd (%rbx,%r15,1),%xmm0
  mov    (%rax,%r15,1),%eax
  mov    %eax,(%r12)
  mov    $0x1,%eax
  call   <mean_it@plt>
  pxor   %xmm0,%xmm0
  add    0x8(%rsp),%r12
  cvtsi2ss %eax,%xmm0
  movss  %xmm0,(%rbx,%r15,1)
  add    $0x4,%r15
  cmp    %r15,0x10(%rsp)
  jne    L32
  mov    %r13,%r14
  mov    0x40(%rsp),%r13
L106:
  mov    0x30(%rsp),%ebp
  mov    0x38(%rsp),%rax
  addl   $0x1,0x28(%rsp)
  add    $0x1,%ebp
  mov    %ebp,(%rax)
  mov    0x28(%rsp),%eax
  cmp    %eax,0x68(%rsp)
  jne    L241
  addl   $0x1,0x6c(%rsp)
  mov    0x70(%rsp),%edi
  mov    0x6c(%rsp),%eax
  add    %edi,0x68(%rsp)
  cmp    %eax,0x74(%rsp)
  jg     L168
L124:
  mov    0x34(%rsp),%edi
  test   %edi,%edi
  jle    L266
L123:
  xor    %ebx,%ebx
  xor    %r12d,%r12d
  mov    %r14,0x8(%rsp)
  mov    0x50(%rsp),%rbp
  mov    %rbx,%r14
  mov    0x2c(%rsp),%r15d
  mov    %r12d,%ebx
  mov    0x78(%rsp),%r12
  nopl   0x0(%rax,%rax,1)
L194:
  mov    0x0(%rbp,%r14,8),%rdi
  test   %r15d,%r15d
  jle    L82
  mov    0x0(%r13,%r14,8),%rcx
  xor    %eax,%eax
  nopl   0x0(%rax)
L239:
  cvttss2si (%rdi,%rax,4),%edx
  cmp    %edx,(%rcx,%rax,4)
  je     L43
  mov    %edx,(%rcx,%rax,4)
  mov    $0x1,%ebx
L43:
  add    $0x1,%rax
  cmp    %rax,%r12
  jne    L239
L82:
  call   <free@plt>
  add    $0x1,%r14
  cmp    %r14d,0x34(%rsp)
  jg     L194
  mov    0x8(%rsp),%r14
  mov    %ebx,%r12d
L64:
  mov    0x50(%rsp),%rdi
  call   <free@plt>
  mov    %r14,%rdi
  call   <free@plt>
  add    $0x88,%rsp
  mov    %r12d,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L104:
  mov    0x34(%rsp),%edx
  test   %edx,%edx
  jg     L78
  movslq 0x24(%rsp),%rax
  mov    0x48(%rsp),%rdi
  lea    (%rdi,%rax,4),%rax
  mov    %rax,0x38(%rsp)
  mov    (%rax),%eax
  mov    %eax,0x30(%rsp)
  jmp    L106
  nopw   0x0(%rax,%rax,1)
L214:
  movslq 0x24(%rsp),%rax
  mov    0x48(%rsp),%rcx
  lea    (%rcx,%rax,4),%rcx
  mov    (%rcx),%edi
  mov    %rcx,0x38(%rsp)
  mov    %edi,0x30(%rsp)
  jmp    L128
L264:
  mov    0x74(%rsp),%ecx
  test   %ecx,%ecx
  jg     L98
L266:
  xor    %r12d,%r12d
  jmp    L64
```

assistant: ```c
int kmeans_it(int *image, int width, int height, int depth, int k, int **clusters, int mean_colors, int *ret) {
  int i = 0, j = 0, d = 0, ki = 0;
  int min_dist = -1, temp_dist = -1;
  int cluster = 0;
  int *pixel_it = malloc(depth * sizeof(int));
  int *pix_count = (int *)calloc(k, sizeof(int));
  int changed = 0;
  float **new_clusters = (float **)malloc(k * sizeof(float *));
  for (i = 0; i < k; i++) {
    new_clusters[i] = calloc(depth, sizeof(float));
  }
  for (i = 0; i < width; i++) {
    for (j = 0; j < height; j++) {
      int l = 0;
      min_dist = -1;
      for (l = 0; l < depth; l++) {
        pixel_it[l] = image[((height) * (i) + (j) + (width) * (height) * (l))];
      }
      for (ki = 0; ki < k; ki++) {
        temp_dist = sqd_distance(pixel_it, clusters[ki], depth);
        if (min_dist < 0 || temp_dist < min_dist) {
          min_dist = temp_dist;
          cluster = ki;
        }
      }
      for (l = 0; l < depth; l++) {
        ret[((height) * (i) + (j) + (width) * (height) * (l))] = (mean_colors) ? clusters[cluster][l] : clusters[cluster][l];
        new_clusters[cluster][l] = mean_it(new_clusters[cluster][l], pix_count[cluster], pixel_it[l]);
      }
      pix_count[cluster]++;
    }
  }
  for (i = 0; i < k; i++) {
    for (d = 0; d < depth; d++) {
      if (clusters[i][d] != (int)new_clusters[i][d]) {
        clusters[i][d] = (int)(new_clusters[i][d]);
        changed = 1;
      }
    }
    free(new_clusters[i]);
  }
  free(new_clusters);
  free(pixel_it);
  return changed;
}
```