user: 下面的汇编指令对应的C源代码是什么:

```asm
<emalg>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %edi,0x40(%rsp)
  mov    %esi,0x58(%rsp)
  mov    %edx,0x5c(%rsp)
  mov    %rcx,0x18(%rsp)
  mov    %r8,0x20(%rsp)
  mov    %r9,0x28(%rsp)
  mov    0xa0(%rsp),%r14d
  mov    0xb8(%rsp),%r15
  test   %r14d,%r14d
  jne    L206
  mov    %esi,0x44(%rsp)
  mov    %edx,0x30(%rsp)
L94:
  mov    $0x4,%esi
  mov    0x40(%rsp),%edi
  call   <calloc@plt>
  mov    %rax,%r13
  mov    0x44(%rsp),%ebx
  lea    0x0(,%rbx,4),%edi
  call   <malloc@plt>
  mov    %rax,0x50(%rsp)
  test   %ebx,%ebx
  jle    L145
  mov    %r15,%rax
  lea    -0x1(%rbx),%edx
  lea    0x4(%r15,%rdx,4),%rcx
L1:
  movslq (%rax),%rdx
  addl   $0x1,0x0(%r13,%rdx,4)
  add    $0x4,%rax
  cmp    %rcx,%rax
  jne    L1
L145:
  movl   $0xa,0x4c(%rsp)
  movl   $0x0,0x48(%rsp)
  mov    %r15,0xb8(%rsp)
  mov    %r14d,%r15d
  mov    0xb8(%rsp),%r14
  jmp    L19
L206:
  mov    0x5c(%rsp),%eax
  mov    %eax,0x44(%rsp)
  mov    0x58(%rsp),%eax
  mov    %eax,0x30(%rsp)
  jmp    L94
L85:
  mov    0x44(%rsp),%eax
  test   %eax,%eax
  jle    L122
  mov    %eax,%ecx
  mov    $0x0,%eax
  mov    0x50(%rsp),%rsi
L71:
  mov    (%r14,%rax,4),%edx
  mov    %edx,(%rsi,%rax,4)
  add    $0x1,%rax
  cmp    %rax,%rcx
  jne    L71
L122:
  shll   0x4c(%rsp)
  jmp    L149
L180:
  add    $0x1,%rbx
  cmp    %rbx,0x8(%rsp)
  je     L60
L68:
  mov    %ebx,0x10(%rsp)
  cmp    %ebx,%r12d
  je     L180
  sub    $0x8,%rsp
  push   %r14
  push   %rbx
  mov    0x2c(%rsp),%eax
  push   %rax
  mov    0x48(%rsp),%r9
  mov    0xe8(%rsp),%r8
  mov    0x40(%rsp),%rcx
  mov    0xe0(%rsp),%rdx
  mov    0x38(%rsp),%rsi
  mov    0x50(%rsp),%edi
  call   *0xd0(%rsp)
  add    $0x20,%rsp
  movsd  (%rsp),%xmm1
  comisd %xmm0,%xmm1
  jbe    L180
  movslq (%r15,%rbp,4),%rax
  subl   $0x1,0x0(%r13,%rax,4)
  mov    0x10(%rsp),%eax
  mov    %eax,(%r15,%rbp,4)
  addl   $0x1,0x0(%r13,%rbx,4)
  movsd  %xmm0,(%rsp)
  movl   $0x1,0x34(%rsp)
  jmp    L180
L60:
  add    $0x1,%rbp
  cmp    %rbp,0x38(%rsp)
  je     L91
L76:
  mov    %ebp,0x14(%rsp)
  mov    (%r15,%rbp,4),%r12d
  movslq %r12d,%rax
  cmpl   $0x1,0x0(%r13,%rax,4)
  je     L60
  sub    $0x8,%rsp
  push   %r14
  push   %r12
  push   %rbp
  mov    0x48(%rsp),%r9
  mov    0xe8(%rsp),%r8
  mov    0x40(%rsp),%rcx
  mov    0xe0(%rsp),%rdx
  mov    0x38(%rsp),%rsi
  mov    0x50(%rsp),%edi
  call   *0xd0(%rsp)
  movsd  %xmm0,0x20(%rsp)
  add    $0x20,%rsp
  mov    0x40(%rsp),%eax
  test   %eax,%eax
  jle    L60
  mov    %eax,%eax
  mov    %rax,0x8(%rsp)
  mov    $0x0,%ebx
  jmp    L68
L91:
  mov    %r15,0xb8(%rsp)
  mov    %r14d,%r15d
  mov    0xb8(%rsp),%r14
  cmpl   $0x0,0x34(%rsp)
  je     L205
  mov    %r14,%rdx
  mov    0x50(%rsp),%rsi
  mov    0x44(%rsp),%edi
  call   <equal_clusters@plt>
  test   %eax,%eax
  jne    L205
L19:
  mov    0x48(%rsp),%eax
  cltd
  idivl  0x4c(%rsp)
  test   %edx,%edx
  je     L85
L149:
  addl   $0x1,0x48(%rsp)
  sub    $0x8,%rsp
  push   %r15
  push   0xd8(%rsp)
  push   0xd8(%rsp)
  mov    %r14,%r9
  mov    0x40(%rsp),%r8
  mov    0x38(%rsp),%rcx
  mov    0x7c(%rsp),%edx
  mov    0x78(%rsp),%esi
  mov    0x60(%rsp),%edi
  call   *0xc8(%rsp)
  add    $0x20,%rsp
  mov    0x44(%rsp),%eax
  test   %eax,%eax
  jle    L205
  mov    %eax,%eax
  mov    %rax,0x38(%rsp)
  mov    $0x0,%ebp
  movl   $0x0,0x34(%rsp)
  mov    %r14,0xb8(%rsp)
  mov    %r15d,%r14d
  mov    0xb8(%rsp),%r15
  jmp    L76
L205:
  mov    0x50(%rsp),%rdi
  call   <free@plt>
  mov    %r13,%rdi
  call   <free@plt>
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void emalg(int nclusters, int nrows, int ncolumns, double **data, int **mask, double weight[], int transpose, void getclustercenter(int, int, int, double **, int **, int[], double **, int **, int), double metric(int, double **, double **, int **, int **, const double[], int, int, int), int clusterid[], double **cdata, int **cmask) {
  const int nitems = (transpose == 0) ? nrows : ncolumns;
  const int ndata = (transpose == 0) ? ncolumns : nrows;
  int *cn = calloc(nclusters, sizeof(int));
  int *savedids = malloc(nitems * sizeof(int));
  int changed;
  int iteration = 0;
  int period = 10;
  int i, j;
  for (i = 0; i < nitems; i++)
    cn[clusterid[i]]++;
  do {
    if (iteration % period == 0) {
      for (i = 0; i < nitems; i++)
        savedids[i] = clusterid[i];
      period *= 2;
    }
    iteration++;
    getclustercenter(nclusters, nrows, ncolumns, data, mask, clusterid, cdata, cmask, transpose);
    changed = 0;
    for (i = 0; i < nitems; i++) {
      double distance;
      int jnow = clusterid[i];
      if (cn[jnow] == 1)
        continue;
      distance = metric(ndata, data, cdata, mask, cmask, weight, i, jnow, transpose);
      for (j = 0; j < nclusters; j++) {
        double tdistance;
        if (j == jnow)
          continue;
        tdistance = metric(ndata, data, cdata, mask, cmask, weight, i, j, transpose);
        if (tdistance < distance) {
          distance = tdistance;
          cn[clusterid[i]]--;
          clusterid[i] = j;
          cn[j]++;
          changed = 1;
        }
      }
    }
  } while (changed && !equal_clusters(nitems, savedids, clusterid));
  free(savedids);
  free(cn);
  return;
}
```