user: 请从给定的汇编指令推导出C源代码：

```asm
<RBR_dist_stats>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  movl   $0x0,(%r8)
  mov    %esi,0x18(%rsp)
  mov    0x80(%rsp),%r12d
  mov    %rdx,0x30(%rsp)
  mov    %rcx,0x38(%rsp)
  movl   $0x42c80000,(%r9)
  test   %esi,%esi
  jle    L121
  cmp    $0x1,%esi
  je     L110
  movslq %esi,%rax
  pxor   %xmm2,%xmm2
  mov    %rdi,%r14
  mov    %r8,%rbx
  mov    %eax,0x24(%rsp)
  mov    $0x1,%r15d
  mov    %r9,%r13
  mov    %rax,0x28(%rsp)
  movl   $0x0,0x20(%rsp)
  movss  %xmm2,0xc(%rsp)
  nopl   0x0(%rax)
L115:
  cmpl   $0x2,0x88(%rsp)
  je     L42
  pxor   %xmm7,%xmm7
  mov    %r12d,0x80(%rsp)
  mov    %r15,%rbp
  mov    %r15,%r12
  cvtsi2ssl 0x88(%rsp),%xmm7
  mov    %r13,%r15
  mov    0x80(%rsp),%r13d
  movss  %xmm7,0x1c(%rsp)
  nopw   0x0(%rax,%rax,1)
L92:
  mov    (%r14,%rbp,8),%rsi
  mov    -0x8(%r14,%r12,8),%rdi
  mov    %r13d,%edx
  movss  %xmm2,0x14(%rsp)
  mov    0x88(%rsp),%ecx
  call   <RBR_ipowdist@plt>
  movaps %xmm0,%xmm1
  movss  D0(%rip),%xmm0
  divss  0x1c(%rsp),%xmm0
  cvttss2si %xmm0,%edi
  pxor   %xmm0,%xmm0
  movss  %xmm1,0x10(%rsp)
  cvtss2sd %xmm1,%xmm0
  call   <pow@plt>
  movss  0x10(%rsp),%xmm1
  comiss (%rbx),%xmm0
  addss  0xc(%rsp),%xmm1
  movss  0x14(%rsp),%xmm2
  addss  %xmm0,%xmm2
  movss  %xmm1,0xc(%rsp)
  jbe    L56
  movss  %xmm0,(%rbx)
L56:
  movss  (%r15),%xmm1
  comiss %xmm0,%xmm1
  jbe    L78
  movss  %xmm0,(%r15)
L78:
  add    $0x1,%rbp
  cmp    %ebp,0x18(%rsp)
  jg     L92
  mov    0x20(%rsp),%edx
  mov    0x24(%rsp),%eax
  mov    %r13d,0x80(%rsp)
  mov    %r15,%r13
  mov    %r12,%r15
  mov    0x80(%rsp),%r12d
  lea    -0x1(%rax,%rdx,1),%eax
  mov    %eax,0x20(%rsp)
L114:
  subl   $0x1,0x24(%rsp)
  add    $0x1,%r15
  cmp    0x28(%rsp),%r15
  jne    L115
L5:
  movaps %xmm2,%xmm1
  pxor   %xmm3,%xmm3
  mulss  %xmm2,%xmm1
  cvtsi2ssl 0x20(%rsp),%xmm3
L60:
  divss  %xmm3,%xmm1
  movss  0xc(%rsp),%xmm0
  mov    0x30(%rsp),%rax
  divss  %xmm3,%xmm2
  subss  %xmm1,%xmm0
  divss  %xmm3,%xmm0
  movss  %xmm2,(%rax)
  cvtss2sd %xmm0,%xmm0
  call   <sqrt@plt>
  mov    0x38(%rsp),%rax
  movss  %xmm0,(%rax)
  add    $0x48,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L42:
  mov    %r12d,0x80(%rsp)
  mov    %r15,%rbp
  mov    %r15,%r12
  mov    %r13,%r15
  mov    0x80(%rsp),%r13d
  nopl   0x0(%rax)
L106:
  mov    (%r14,%rbp,8),%rsi
  mov    -0x8(%r14,%r12,8),%rdi
  mov    $0x2,%ecx
  mov    %r13d,%edx
  movss  %xmm2,0x14(%rsp)
  call   <RBR_ipowdist@plt>
  movaps %xmm0,%xmm1
  cvtss2sd %xmm0,%xmm0
  movss  %xmm1,0x10(%rsp)
  call   <sqrt@plt>
  movss  0x10(%rsp),%xmm1
  comiss (%rbx),%xmm0
  addss  0xc(%rsp),%xmm1
  movss  0x14(%rsp),%xmm2
  addss  %xmm0,%xmm2
  movss  %xmm1,0xc(%rsp)
  jbe    L111
  movss  %xmm0,(%rbx)
L111:
  movss  (%r15),%xmm1
  comiss %xmm0,%xmm1
  jbe    L65
  movss  %xmm0,(%r15)
L65:
  add    $0x1,%rbp
  cmp    %ebp,0x18(%rsp)
  jg     L106
  mov    0x20(%rsp),%ecx
  mov    0x24(%rsp),%eax
  mov    %r13d,0x80(%rsp)
  mov    %r15,%r13
  mov    %r12,%r15
  mov    0x80(%rsp),%r12d
  lea    -0x1(%rax,%rcx,1),%eax
  mov    %eax,0x20(%rsp)
  jmp    L114
L121:
  pxor   %xmm1,%xmm1
  movaps %xmm1,%xmm3
  movss  %xmm1,0xc(%rsp)
  movaps %xmm1,%xmm2
  jmp    L60
L110:
  pxor   %xmm2,%xmm2
  movl   $0x0,0x20(%rsp)
  movss  %xmm2,0xc(%rsp)
  jmp    L5
```

tool_calls: name: retrieve_data
arguments: data_label: D0
type_infered: byte[4]

tool: {"D0": [0x00,0x00,0x80,0x3f]}

assistant: ```c
int RBR_dist_stats(int **coords, int ncoords, float *mean, float *sd, float *longest, float *shortest, int PRECISION, int DIST_POWER) {
  int i, j, k;
  int total;
  float sum, sumpow;
  float powdist, dist;
  sum = 0.0;
  sumpow = 0.0;
  (*longest) = 0.0;
  (*shortest) = 100.0;
  total = 0.0;
  for (j = 0; j < ncoords; ++j) {
    for (k = j + 1; k < ncoords; ++k) {
      total++;
      powdist = RBR_ipowdist(coords[j], coords[k], PRECISION, DIST_POWER);
      if (DIST_POWER == 2)
        dist = sqrt((double)powdist);
      else
        dist = pow((double)powdist, (1 / (float)DIST_POWER));
      sum += dist;
      sumpow += powdist;
      if (dist > (*longest))
        (*longest) = dist;
      if (dist < (*shortest))
        (*shortest) = dist;
    }
  }
  (*mean) = (sum / (float)total);
  (*sd) = sqrt((double)((sumpow - (sum * sum) / (float)total) / (float)total));
  return 0;
}
```