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

```asm
<assign_class_by_pam>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %edi,%r12d
  mov    %esi,%r13d
  mov    %rdx,0x40(%rsp)
  mov    %rcx,%rbx
  mov    %r8,%rbp
  mov    $0x0,%eax
  call   <allocate_double_1D@plt>
  mov    %eax,0x4c(%rsp)
  movslq %eax,%rsi
  mov    %rsi,(%rsp)
  mov    %r12d,%edi
  mov    $0x0,%eax
  call   <allocate_double_1D@plt>
  mov    %eax,0x50(%rsp)
  movslq %eax,%rdi
  mov    %rdi,0x8(%rsp)
  mov    %r12d,%edi
  mov    $0x0,%eax
  call   <allocate_double_1D@plt>
  mov    %eax,0x54(%rsp)
  movslq %eax,%r11
  mov    %r11,0x10(%rsp)
  mov    $0x1,%edi
  mov    $0x0,%eax
  call   <allocate_double_1D@plt>
  movslq %eax,%r11
  mov    %r11,0x18(%rsp)
  mov    $0x2,%edi
  mov    $0x0,%eax
  call   <allocate_double_1D@plt>
  movslq %eax,%r15
  mov    %r12d,%edi
  mov    $0x0,%eax
  call   <allocate_int_1D@plt>
  mov    %eax,0x58(%rsp)
  movslq %eax,%r10
  mov    %r10,0x20(%rsp)
  mov    %r12d,%edi
  mov    $0x0,%eax
  call   <allocate_int_1D@plt>
  mov    %eax,0x5c(%rsp)
  movslq %eax,%r10
  mov    %r10,0x28(%rsp)
  mov    %r12d,%edi
  mov    $0x0,%eax
  call   <allocate_int_1D@plt>
  mov    %eax,%r14d
  movslq %eax,%r9
  mov    %r9,0x30(%rsp)
  mov    $0x1,%edi
  mov    $0x0,%eax
  call   <allocate_int_1D@plt>
  movslq %eax,%r10
  mov    %r10,0x38(%rsp)
  mov    0x40(%rsp),%rdx
  mov    (%rdx),%rax
  lea    -0x8(%rax),%rdx
  mov    %rdx,0x40(%rsp)
  test   %r12d,%r12d
  jle    L97
  movslq %r14d,%r14
  movslq 0x54(%rsp),%r11
  mov    %r12d,%eax
  lea    (%r14,%rax,4),%r10
  mov    %r11,%rdx
  mov    %r14,%rax
  mov    %rbp,%r9
  sub    %r14,%r9
  movslq 0x5c(%rsp),%r8
  sub    %r14,%r8
  movslq 0x58(%rsp),%rdi
  sub    %r14,%rdi
  movslq 0x50(%rsp),%rsi
  sub    %r11,%rsi
  movslq 0x4c(%rsp),%rcx
  sub    %r11,%rcx
L77:
  movl   $0x0,(%r9,%rax,1)
  movl   $0x0,(%rax)
  movl   $0x0,(%r8,%rax,1)
  movl   $0x0,(%rdi,%rax,1)
  movq   $0x0,(%rdx)
  movq   $0x0,(%rsi,%rdx,1)
  movq   $0x0,(%rcx,%rdx,1)
  add    $0x4,%rax
  add    $0x8,%rdx
  cmp    %r10,%rax
  jne    L77
  test   %r13d,%r13d
  jle    L67
L50:
  mov    %rbx,%rax
  lea    -0x1(%r13),%edx
  lea    0x4(%rbx,%rdx,4),%rdx
L38:
  movl   $0x0,(%rax)
  add    $0x4,%rax
  cmp    %rdx,%rax
  jne    L38
  mov    0x18(%rsp),%rax
  movq   $0x0,(%rax)
  movq   $0x0,0x8(%r15)
  movq   $0x0,(%r15)
  mov    0x38(%rsp),%rcx
  movl   $0x1,(%rcx)
  push   %rcx
  push   %rbp
  push   %rbx
  push   %r15
  push   %rax
  push   0x38(%rsp)
  push   0x38(%rsp)
  push   0x38(%rsp)
  mov    0x70(%rsp),%r9
  mov    0x68(%rsp),%r8
  mov    0x60(%rsp),%rcx
  mov    0x80(%rsp),%rdx
  mov    %r13d,%esi
  mov    %r12d,%edi
  mov    $0x0,%eax
  call   <phyclust_pam@plt>
  add    $0x40,%rsp
  test   %r12d,%r12d
  jle    L40
L45:
  mov    $0x0,%eax
L78:
  subl   $0x1,0x0(%rbp,%rax,4)
  add    $0x1,%rax
  cmp    %eax,%r12d
  jg     L78
  test   %r13d,%r13d
  jle    L83
L40:
  mov    $0x0,%eax
L99:
  subl   $0x1,(%rbx,%rax,4)
  add    $0x1,%rax
  cmp    %eax,%r13d
  jg     L99
L83:
  mov    (%rsp),%rdi
  call   <free@plt>
  mov    0x8(%rsp),%rdi
  call   <free@plt>
  mov    0x10(%rsp),%rdi
  call   <free@plt>
  mov    0x18(%rsp),%rdi
  call   <free@plt>
  mov    %r15,%rdi
  call   <free@plt>
  mov    0x20(%rsp),%rdi
  call   <free@plt>
  mov    0x28(%rsp),%rdi
  call   <free@plt>
  mov    0x30(%rsp),%rdi
  call   <free@plt>
  mov    0x38(%rsp),%rdi
  call   <free@plt>
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L67:
  mov    0x18(%rsp),%rax
  movq   $0x0,(%rax)
  movq   $0x0,0x8(%r15)
  movq   $0x0,(%r15)
  mov    0x38(%rsp),%rcx
  movl   $0x1,(%rcx)
  push   %rcx
  push   %rbp
  push   %rbx
  push   %r15
  push   %rax
  push   0x38(%rsp)
  push   0x38(%rsp)
  push   0x38(%rsp)
  mov    0x70(%rsp),%r9
  mov    0x68(%rsp),%r8
  mov    0x60(%rsp),%rcx
  mov    0x80(%rsp),%rdx
  mov    %r13d,%esi
  mov    %r12d,%edi
  mov    $0x0,%eax
  call   <phyclust_pam@plt>
  add    $0x40,%rsp
  jmp    L45
L97:
  test   %r13d,%r13d
  jg     L50
  mov    0x18(%rsp),%rax
  movq   $0x0,(%rax)
  movq   $0x0,0x8(%r15)
  movq   $0x0,(%r15)
  mov    0x38(%rsp),%rcx
  movl   $0x1,(%rcx)
  push   %rcx
  push   %rbp
  push   %rbx
  push   %r15
  push   %rax
  push   0x38(%rsp)
  push   0x38(%rsp)
  push   0x38(%rsp)
  mov    0x70(%rsp),%r9
  mov    0x68(%rsp),%r8
  mov    0x60(%rsp),%rcx
  mov    0x80(%rsp),%rdx
  mov    %r13d,%esi
  mov    %r12d,%edi
  mov    $0x0,%eax
  call   <phyclust_pam@plt>
  add    $0x40,%rsp
  jmp    L83
```

assistant: ```c
void assign_class_by_pam(int N_X, int K, double **EDM_LT_pam, int *center_id, int *class_id) {
  int i;
  double *dys, *radus, *damer, *avsyl, *ttsyl, *obj;
  int *nsend, *nrepr, *nelem, *nisol;
  radus = allocate_double_1D(N_X);
  damer = allocate_double_1D(N_X);
  avsyl = allocate_double_1D(N_X);
  ttsyl = allocate_double_1D(1);
  obj = allocate_double_1D(2);
  nsend = allocate_int_1D(N_X);
  nrepr = allocate_int_1D(N_X);
  nelem = allocate_int_1D(N_X);
  nisol = allocate_int_1D(1);
  dys = EDM_LT_pam[0] - 1;
  for (i = 0; i < N_X; i++) {
    nsend[i] = nrepr[i] = nelem[i] = class_id[i] = 0;
    radus[i] = damer[i] = avsyl[i] = 0;
  }
  for (i = 0; i < K; i++) {
    center_id[i] = 0;
  }
  ttsyl[0] = 0;
  obj[0] = obj[1] = 0;
  nisol[0] = 1;
  phyclust_pam(N_X, K, dys, nsend, nrepr, nelem, radus, damer, avsyl, ttsyl, obj, center_id, class_id, nisol);
  for (i = 0; i < N_X; i++) {
    class_id[i]--;
  }
  for (i = 0; i < K; i++) {
    center_id[i]--;
  }
  free(radus);
  free(damer);
  free(avsyl);
  free(ttsyl);
  free(obj);
  free(nsend);
  free(nrepr);
  free(nelem);
  free(nisol);
}
```