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

```asm
<sum_squares>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x18,%rsp
  mov    %rdi,%r12
  mov    %rsi,%rbp
  mov    %rdx,%rbx
  mov    %rcx,%r13
  mov    %r8,%r14
  mov    %r9,0x8(%rsp)
  mov    0x50(%rsp),%r15
  mov    (%rdx),%eax
  lea    0x0(,%rax,8),%edi
  call   <malloc@plt>
  mov    %rax,%rdi
  cmpl   $0x0,(%rbx)
  jle    L115
  mov    %rax,%rdx
  mov    $0x0,%eax
L84:
  movq   $0x0,(%rdx)
  add    $0x1,%eax
  add    $0x8,%rdx
  cmp    %eax,(%rbx)
  jg     L84
L115:
  mov    $0x0,%esi
  cmpl   $0x0,0x0(%rbp)
  jg     L161
L213:
  cmpl   $0x0,(%rbx)
  jle    L6
  mov    %rdi,%rax
  mov    $0x0,%edx
L236:
  pxor   %xmm1,%xmm1
  cvtsi2sdl 0x0(%rbp),%xmm1
  movsd  (%rax),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  add    $0x1,%edx
  add    $0x8,%rax
  cmp    %edx,(%rbx)
  jg     L236
L6:
  movq   $0x0,(%r15)
  mov    0x0(%rbp),%ecx
  test   %ecx,%ecx
  jle    L37
  mov    %r12,%r8
  mov    $0x0,%esi
  jmp    L128
L206:
  mov    %edx,%eax
  imul   0x0(%rbp),%eax
  add    %esi,%eax
  cltq
  movsd  (%rcx),%xmm0
  addsd  (%r12,%rax,8),%xmm0
  movsd  %xmm0,(%rcx)
  add    $0x1,%edx
  add    $0x8,%rcx
  cmp    %edx,(%rbx)
  jg     L206
L17:
  add    $0x1,%esi
  cmp    %esi,0x0(%rbp)
  jle    L213
L161:
  mov    %rdi,%rcx
  mov    $0x0,%edx
  cmpl   $0x0,(%rbx)
  jg     L206
  jmp    L17
L224:
  mov    %rdi,%rax
  movslq %ecx,%rcx
  shl    $0x3,%rcx
  mov    %edx,%edx
  lea    (%rdi,%rdx,8),%r9
  mov    %r8,%rdx
  pxor   %xmm1,%xmm1
L219:
  movsd  (%rax),%xmm0
  subsd  (%rdx),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  add    $0x8,%rax
  add    %rcx,%rdx
  cmp    %r9,%rax
  jne    L219
L34:
  addsd  (%r15),%xmm1
  movsd  %xmm1,(%r15)
  add    $0x1,%esi
  mov    0x0(%rbp),%ecx
  add    $0x8,%r8
  cmp    %esi,%ecx
  jle    L196
L128:
  mov    (%rbx),%edx
  pxor   %xmm1,%xmm1
  test   %edx,%edx
  jg     L224
  jmp    L34
L196:
  cmpl   $0x0,0x0(%r13)
  jle    L135
L117:
  mov    $0x0,%eax
L148:
  mov    0x58(%rsp),%rsi
  movq   $0x0,(%rsi,%rax,8)
  add    $0x1,%rax
  cmp    %eax,0x0(%r13)
  jg     L148
L135:
  mov    0x0(%rbp),%esi
  test   %esi,%esi
  jle    L209
  mov    $0x0,%r9d
  jmp    L232
L191:
  movslq %esi,%rsi
  shl    $0x3,%rsi
  movslq 0x0(%r13),%r8
  shl    $0x3,%r8
  movslq (%r14,%r9,4),%rax
  mov    0x8(%rsp),%rcx
  lea    (%rcx,%rax,8),%rcx
  mov    %r12,%rdx
  pxor   %xmm1,%xmm1
  mov    $0x0,%eax
L204:
  movsd  (%rdx),%xmm0
  subsd  (%rcx),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  add    $0x1,%eax
  add    %rsi,%rdx
  add    %r8,%rcx
  cmp    %r10d,%eax
  jne    L204
L231:
  movslq (%r14,%r9,4),%rax
  mov    0x58(%rsp),%rsi
  lea    (%rsi,%rax,8),%rax
  addsd  (%rax),%xmm1
  movsd  %xmm1,(%rax)
  mov    0x0(%rbp),%esi
  add    $0x1,%r9
  add    $0x8,%r12
  cmp    %r9d,%esi
  jle    L209
L232:
  mov    (%rbx),%r10d
  pxor   %xmm1,%xmm1
  test   %r10d,%r10d
  jg     L191
  jmp    L231
L37:
  cmpl   $0x0,0x0(%r13)
  jg     L117
L209:
  call   <free@plt>
  add    $0x18,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void sum_squares(double *x, int *nr, int *nc, int *k, int *cluster, double *centers, double *totss, double *withiness) {
  int i, j;
  double *global_average = (double *)malloc(*nc * sizeof(double));
  for (j = 0; j < *nc; ++j) {
    global_average[j] = 0.0;
  }
  for (i = 0; i < *nr; ++i) {
    for (j = 0; j < *nc; ++j) {
      global_average[j] += x[j * (*nr) + i];
    }
  }
  for (j = 0; j < *nc; ++j) {
    global_average[j] = global_average[j] / (*nr);
  }
  *totss = 0.0;
  double tmp_ss, temp;
  for (i = 0; i < (*nr); ++i) {
    tmp_ss = 0.0;
    for (j = 0; j < (*nc); ++j) {
      temp = global_average[j] - x[j * (*nr) + i];
      tmp_ss += temp * temp;
    }
    *totss += tmp_ss;
  }
  int t;
  for (t = 0; t < (*k); ++t) {
    withiness[t] = 0.0;
  }
  for (i = 0; i < (*nr); ++i) {
    tmp_ss = 0;
    for (j = 0; j < (*nc); ++j) {
      temp = x[j * (*nr) + i] - centers[j * (*k) + cluster[i]];
      tmp_ss += temp * temp;
    }
    withiness[cluster[i]] += tmp_ss;
  }
  free(global_average);
  return;
}
```