user: Converte this assembly snippet into its equivalent C code:

```asm
<BestMatrixDistance>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x1000,%rsp
  orq    $0x0,(%rsp)
  sub    $0x998,%rsp
  movsd  D6(%rip),%xmm0
  mov    %rdi,0x18(%rsp)
  mov    %rdi,%rbx
  mov    %rsi,0x48(%rsp)
  mov    %rbx,%rsi
  xor    %ebx,%ebx
  mov    %rcx,0x58(%rsp)
  mov    $0x190,%ecx
  mov    %rdx,0x50(%rsp)
  mov    %r8,0x60(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x1988(%rsp)
  xor    %eax,%eax
  lea    0x80(%rsp),%rax
  mov    %rax,%rdi
  mov    %rax,0x28(%rsp)
  lea    0x78(%rsp),%rax
  rep movsq %ds:(%rsi),%es:(%rdi)
  mov    %rax,0x68(%rsp)
  lea    0xd00(%rsp),%rax
  mov    %rax,0x30(%rsp)
L43:
  mov    0x68(%rsp),%rcx
  mov    0x50(%rsp),%rdx
  lea    0x1(%rbx),%eax
  movsd  %xmm0,0x8(%rsp)
  mov    0x48(%rsp),%rsi
  mov    0x28(%rsp),%rdi
  mov    %eax,0x44(%rsp)
  call   <MatrixLikelihood@plt>
  movsd  0x78(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm0
  comisd %xmm0,%xmm1
  movsd  %xmm1,0x38(%rsp)
  jbe    L59
  mov    0x30(%rsp),%rax
  lea    0x120(%rsp),%rbx
  mov    %rax,0x10(%rsp)
  lea    0xda0(%rsp),%rax
  mov    %rax,0x20(%rsp)
  mov    0x18(%rsp),%rax
  add    $0xa0,%rax
  mov    %rax,0x8(%rsp)
L6:
  mov    0x10(%rsp),%r14
  mov    0x18(%rsp),%rbp
  lea    -0xa0(%rbx),%r12
  mov    %r14,%r13
  nopl   0x0(%rax)
L34:
  movq   $0x0,0x0(%r13)
  mov    %r12,%r15
  mov    %rbp,%r14
  cs nopw 0x0(%rax,%rax,1)
L11:
  movsd  (%r15),%xmm0
  addsd  (%r14),%xmm0
  add    $0x8,%r15
  add    $0xa0,%r14
  call   <exp@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  addsd  0x0(%r13),%xmm0
  movsd  %xmm0,0x0(%r13)
  cmp    %r15,%rbx
  jne    L11
  call   <log@plt>
  add    $0x8,%r13
  add    $0x8,%rbp
  movsd  %xmm0,-0x8(%r13)
  cmp    %rbp,0x8(%rsp)
  jne    L34
  addq   $0xa0,0x10(%rsp)
  add    $0xa0,%rbx
  cmp    0x20(%rsp),%rbx
  jne    L6
  mov    0x28(%rsp),%rdi
  mov    0x30(%rsp),%rsi
  mov    $0x190,%ecx
  rep movsq %ds:(%rsi),%es:(%rdi)
  movsd  0x38(%rsp),%xmm0
  mov    0x44(%rsp),%ebx
  jmp    L43
L59:
  mov    0x60(%rsp),%rax
  movsd  %xmm0,(%rax)
  mov    0x58(%rsp),%rax
  mov    %ebx,(%rax)
  mov    0x1988(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L40
  add    $0x1998,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L40:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: data_label: D6
type_infered: byte[8]

tool: {"D6": [0x00,0x00,0x00,0xe0,0xcf,0x12,0x63,0xc1]}

assistant: ```c
void BestMatrixDistance(double bmx[20][20], double *pri, double *ct, int *ret_dist, double *ret_ll) {
  double smx[20][20];
  double nmx[20][20];
  int c, r, i;
  int dist;
  double ll;
  double best_ll;
  for (r = 0; r < 20; r++)
    for (c = 0; c < 20; c++)
      smx[r][c] = bmx[r][c];
  best_ll = -9999999.0;
  dist = 0;
  while (1) {
    dist++;
    MatrixLikelihood(smx, pri, ct, &ll);
    if (ll > best_ll)
      best_ll = ll;
    else {
      *ret_ll = best_ll;
      *ret_dist = dist - 1;
      return;
    }
    for (r = 0; r < 20; r++)
      for (c = 0; c < 20; c++) {
        nmx[r][c] = 0.0;
        for (i = 0; i < 20; i++)
          nmx[r][c] += exp(smx[r][i] + bmx[i][c]);
        nmx[r][c] = log(nmx[r][c]);
      }
    for (r = 0; r < 20; r++)
      for (c = 0; c < 20; c++)
        smx[r][c] = nmx[r][c];
  }
}
```