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

```asm
<eval_amdf>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x50,%rsp
  mov    %rdi,-0x28(%rbp)
  mov    %esi,-0x2c(%rbp)
  mov    %rdx,-0x38(%rbp)
  mov    %ecx,-0x30(%rbp)
  mov    %r8d,-0x3c(%rbp)
  mov    %r9,-0x48(%rbp)
  mov    0x10(%rbp),%rax
  movl   $0x0,(%rax)
  mov    0x18(%rbp),%rax
  movl   $0x0,(%rax)
  movl   $0x0,-0x10(%rbp)
  jmp    L59
L31:
  mov    -0x10(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0x3c(%rbp),%eax
  sub    %edx,%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  add    $0x1,%eax
  mov    %eax,-0x8(%rbp)
  mov    -0x8(%rbp),%edx
  mov    -0x2c(%rbp),%eax
  add    %edx,%eax
  sub    $0x1,%eax
  mov    %eax,-0x4(%rbp)
  pxor   %xmm0,%xmm0
  movss  %xmm0,-0x14(%rbp)
  mov    -0x8(%rbp),%eax
  mov    %eax,-0xc(%rbp)
  jmp    L49
L16:
  mov    -0xc(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    -0x10(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x38(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%edx
  mov    -0xc(%rbp),%eax
  add    %edx,%eax
  cltq
  shl    $0x2,%rax
  lea    -0x4(%rax),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  subss  %xmm1,%xmm0
  movd   %xmm0,%eax
  movd   %eax,%xmm0
  call   <fabsf@plt>
  pxor   %xmm0,%xmm0
  cvtsi2ss %rax,%xmm0
  movss  -0x14(%rbp),%xmm1
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x14(%rbp)
  addl   $0x4,-0xc(%rbp)
L49:
  mov    -0xc(%rbp),%eax
  cmp    -0x4(%rbp),%eax
  jle    L16
  mov    -0x10(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  -0x14(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x10(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  mov    0x10(%rbp),%rax
  mov    (%rax),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  comiss %xmm1,%xmm0
  jbe    L34
  mov    0x10(%rbp),%rax
  mov    -0x10(%rbp),%edx
  mov    %edx,(%rax)
L34:
  mov    -0x10(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    0x18(%rbp),%rax
  mov    (%rax),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  comiss %xmm1,%xmm0
  jbe    L43
  mov    0x18(%rbp),%rax
  mov    -0x10(%rbp),%edx
  mov    %edx,(%rax)
L43:
  addl   $0x1,-0x10(%rbp)
L59:
  mov    -0x10(%rbp),%eax
  cmp    -0x30(%rbp),%eax
  jl     L31
  nop
  nop
  leave
  ret
```

assistant: ```c
void eval_amdf(float speech[], int32_t lpita, const int32_t tau[], int32_t ltau, int32_t maxlag, float amdf[], int32_t *minptr, int32_t *maxptr) {
  float sum;
  int i;
  int j;
  int n1;
  int n2;
  *minptr = 0;
  *maxptr = 0;
  for (i = 0; i < ltau; i++) {
    n1 = (maxlag - tau[i]) / 2 + 1;
    n2 = n1 + lpita - 1;
    sum = 0.0f;
    for (j = n1; j <= n2; j += 4)
      sum += fabsf(speech[j - 1] - speech[j + tau[i] - 1]);
    amdf[i] = sum;
    if (amdf[i] < amdf[*minptr])
      *minptr = i;
    if (amdf[i] > amdf[*maxptr])
      *maxptr = i;
  }
}
```