user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<unbiasedest>:
  endbr64
  push   %r15
  xor    %ecx,%ecx
  mov    %rdx,%r15
  push   %r14
  push   %r13
  push   %r12
  mov    %esi,%r12d
  push   %rbp
  push   %rbx
  sub    $0xa8,%rsp
  mov    %rdi,0x28(%rsp)
  lea    0x34(%rsp),%rdx
  mov    %fs:0x28,%rax
  mov    %rax,0x98(%rsp)
  xor    %eax,%eax
  call   <ivmaxmin@plt>
  mov    0x34(%rsp),%eax
  mov    $0x4,%esi
  lea    0x1(%rax),%edi
  call   <calloc@plt>
  mov    %rax,0x20(%rsp)
  test   %rax,%rax
  je     L18
L108:
  mov    0x34(%rsp),%edx
  test   %edx,%edx
  js     L198
  mov    0x20(%rsp),%rcx
  xor    %eax,%eax
  nopl   (%rax)
L162:
  mov    (%r15,%rax,8),%rsi
  add    $0x1,%rax
  add    $0x4,%rcx
  mov    0x4(%rsi),%edx
  add    (%rsi),%edx
  mov    %edx,-0x4(%rcx)
  mov    0x34(%rsp),%edx
  cmp    %eax,%edx
  jge    L162
L198:
  lea    0x1(%rdx),%edi
  pxor   %xmm0,%xmm0
  mov    %r12d,%esi
  call   <initarray_2Ddouble@plt>
  mov    0x34(%rsp),%edx
  mov    %rax,0x38(%rsp)
  test   %edx,%edx
  js     L85
  movsd  D58(%rip),%xmm7
  mov    0x20(%rsp),%rbx
  movq   $0x0,0x10(%rsp)
  movsd  %xmm7,0x18(%rsp)
  nopl   0x0(%rax,%rax,1)
L7:
  mov    0x10(%rsp),%rdi
  movsd  0x18(%rsp),%xmm6
  mov    $0x8,%r13d
  mov    $0x1,%r14d
  mov    (%rax,%rdi,8),%rax
  lea    0x0(,%rdi,8),%rbp
  movsd  %xmm6,(%rax)
  test   %r12d,%r12d
  jg     L192
  jmp    L176
  xchg   %ax,%ax
L257:
  mov    D53(%rip),%rax
  add    $0x1,%r14d
  mov    %rax,(%rcx,%r13,1)
  add    $0x8,%r13
  cmp    %r14d,%r12d
  jl     L176
L192:
  mov    (%r15),%rcx
  mov    %r14d,%esi
  mov    (%rcx),%edi
  call   <ifall@plt>
  mov    (%rbx),%edi
  mov    %r14d,%esi
  movsd  %xmm0,0x8(%rsp)
  call   <ifall@plt>
  movsd  D9(%rip),%xmm5
  mov    0x38(%rsp),%rcx
  comisd %xmm0,%xmm5
  mov    (%rcx,%rbp,1),%rcx
  jae    L257
  movsd  0x8(%rsp),%xmm1
  add    $0x1,%r14d
  divsd  %xmm0,%xmm1
  movsd  %xmm1,(%rcx,%r13,1)
  add    $0x8,%r13
  cmp    %r14d,%r12d
  jge    L192
L176:
  addq   $0x1,0x10(%rsp)
  add    $0x8,%r15
  mov    0x10(%rsp),%rax
  add    $0x4,%rbx
  cmp    %eax,0x34(%rsp)
  jl     L85
  mov    0x38(%rsp),%rax
  jmp    L7
  nopl   0x0(%rax)
L85:
  mov    $0x1,%ebx
  mov    %r12d,%ecx
  shl    %cl,%ebx
  mov    %ebx,%r14d
  test   %ebx,%ebx
  jle    L25
  mov    0x28(%rsp),%rbx
  mov    0x34(%rsp),%esi
  xor    %r13d,%r13d
  pxor   %xmm3,%xmm3
  movsd  D58(%rip),%xmm7
  lea    -0x1(%r12),%eax
  lea    0x4(%rbx,%rax,4),%rbp
  lea    0x40(%rsp),%rbx
  movsd  %xmm7,0x18(%rsp)
L177:
  add    $0x1,%esi
  mov    %rbx,%rdi
  movsd  %xmm3,0x8(%rsp)
  call   <ivzero@plt>
  test   %r12d,%r12d
  mov    0x28(%rsp),%rdx
  mov    %r13d,%ecx
  mov    D36(%rip),%rax
  movsd  0x8(%rsp),%xmm3
  movq   D75(%rip),%xmm4
  movq   %rax,%xmm2
  jle    L209
  nopw   0x0(%rax,%rax,1)
L49:
  mov    %ecx,%eax
  sar    %ecx
  and    $0x1,%eax
  test   %eax,%eax
  je     L186
  mov    (%rdx),%eax
L186:
  cltq
  add    $0x4,%rdx
  addl   $0x1,0x40(%rsp,%rax,4)
  cmp    %rdx,%rbp
  jne    L49
L209:
  mov    0x34(%rsp),%esi
  movsd  0x18(%rsp),%xmm1
  test   %esi,%esi
  js     L154
  mov    0x38(%rsp),%r9
  movsd  0x18(%rsp),%xmm1
  xor    %edx,%edx
  jmp    L182
  xchg   %ax,%ax
L230:
  mov    %ecx,%r8d
  shr    $0x1f,%r8d
  add    %r8d,%eax
  and    $0x1,%eax
  sub    %r8d,%eax
  mov    (%r9,%rdx,8),%r8
  movsd  (%r8,%rcx,8),%xmm0
  comisd %xmm0,%xmm2
  ja     L39
  cmp    $0x1,%eax
  jne    L157
  xorpd  %xmm4,%xmm0
L157:
  mulsd  %xmm0,%xmm1
  add    $0x1,%rdx
  cmp    %edx,%esi
  jl     L154
L182:
  movslq (%rbx,%rdx,4),%rcx
  mov    %rcx,%rax
  test   %rdx,%rdx
  jne    L230
  mov    (%r9),%rax
  movsd  (%rax,%rcx,8),%xmm0
  comisd %xmm0,%xmm2
  jbe    L157
L39:
  mov    0x20(%rsp),%rdi
  call   <free@plt>
  mov    0x34(%rsp),%eax
  lea    0x38(%rsp),%rdi
  lea    0x1(%rax),%esi
  call   <free2D@plt>
  movsd  D53(%rip),%xmm3
L171:
  mov    0x98(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L227
  add    $0xa8,%rsp
  movapd %xmm3,%xmm0
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nop
L154:
  add    $0x1,%r13d
  addsd  %xmm1,%xmm3
  cmp    %r13d,%r14d
  jne    L177
L159:
  movapd %xmm3,%xmm0
  movsd  %xmm3,0x8(%rsp)
  call   <fabs@plt>
  movsd  0x8(%rsp),%xmm3
  cmp    $0x63,%eax
  jle    L71
  movsd  D53(%rip),%xmm3
L71:
  mov    0x20(%rsp),%rdi
  movsd  %xmm3,0x8(%rsp)
  call   <free@plt>
  mov    0x34(%rsp),%eax
  lea    0x38(%rsp),%rdi
  lea    0x1(%rax),%esi
  call   <free2D@plt>
  movsd  0x8(%rsp),%xmm3
  jmp    L171
L18:
  mov    0x34(%rsp),%eax
  lea    D54(%rip),%rdi
  lea    0x1(%rax),%esi
  call   <fatalx@plt>
  jmp    L108
L25:
  pxor   %xmm3,%xmm3
  jmp    L159
L227:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_label: D58
type_infered: f64
name: read_data
arguments: data_label: D53
type_infered: f32[2]
name: read_data
arguments: data_label: D9
type_infered: f64
name: read_data
arguments: data_label: D36
type_infered: f32[2]
name: read_data
arguments: data_label: D75
type_infered: f32[2]
name: read_data
arguments: data_label: D54
type_infered: string

tool: D58: 1

tool: {"D53": [0.0,-6.1103515625]}

tool: D9: 0.1

tool: {"D36": [0.0,-1.875]}

tool: {"D75": [0.0,-0.0]}

tool: D54: Unable to allocate %d unit(s) for item 


assistant: ```c
double unbiasedest(int *ndx, int ndsize, int **counts) {
  double xtop, xbot, yest, y;
  int popind[20];
  int popmax, j, k, n, nmax, a, t, s;
  int *tcounts;
  double **xmomest, yp;
  ivmaxmin(ndx, ndsize, &popmax, NULL);
  if ((tcounts = (int *)calloc((popmax + 1), sizeof(int))) == NULL)
    fatalx("Unable to allocate %d unit(s) for item \n", popmax + 1);
  for (j = 0; j <= popmax; ++j) {
    tcounts[j] = counts[j][0] + counts[j][1];
  }
  xmomest = initarray_2Ddouble(popmax + 1, ndsize, 0.0);
  for (j = 0; j <= popmax; ++j) {
    xmomest[j][0] = 1.0;
    for (k = 1; k <= ndsize; ++k) {
      xtop = ifall(counts[j][0], k);
      xbot = ifall(tcounts[j], k);
      if (xbot <= 0.1)
        xmomest[j][k] = -10000.0;
      else
        xmomest[j][k] = (double)xtop / (double)xbot;
    }
  }
  nmax = (1 << (ndsize)) - 1;
  yest = 0.0;
  for (n = 0; n <= nmax; ++n) {
    t = n;
    ivzero(popind, popmax + 1);
    for (k = 0; k < ndsize; ++k) {
      a = 0;
      s = t & 1;
      t = t >> 1;
      if (s == 1)
        a = ndx[k];
      ++popind[a];
    }
    yp = 1.0;
    for (j = 0; j <= popmax; ++j) {
      t = popind[j];
      s = 0;
      if (j > 0)
        s = t % 2;
      y = xmomest[j][t];
      if (y < -1.0) {
        free(tcounts);
        free2D(&xmomest, popmax + 1);
        return (-10000.0);
      }
      if (s == 1)
        y = -y;
      yp *= y;
    }
    yest += yp;
  }
  if (fabs(yest) >= 100)
    yest = -10000;
  free(tcounts);
  free2D(&xmomest, popmax + 1);
  return (yest);
}
```