user: Please show me the C implementation that matches the assembly code listed below:

```asm
<frac_pch>:
  endbr64
  push   %r15
  mov    %edx,%eax
  push   %r14
  push   %r13
  push   %r12
  mov    %r9d,%r12d
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rsi,0x28(%rsp)
  mov    %ecx,0x34(%rsp)
  mov    %r8d,0x24(%rsp)
  mov    %fs:0x28,%rdx
  mov    %rdx,0x48(%rsp)
  xor    %edx,%edx
  test   %eax,%eax
  jle    L66
  cvtss2sd %xmm0,%xmm0
  mov    %r8d,%ebx
  mov    %ecx,%edi
  pxor   %xmm1,%xmm1
  addsd  D88(%rip),%xmm0
  cvttsd2si %xmm0,%r8d
  pxor   %xmm0,%xmm0
  mov    %r8d,%edx
  cvtsi2sd %r8d,%xmm0
  lea    (%rax,%r8,1),%ecx
  mulsd  D60(%rip),%xmm0
  cmp    %ebx,%ecx
  cmovg  %ebx,%ecx
  sub    %eax,%edx
  cmp    %edi,%edx
  cmovl  %edi,%edx
  cvtsi2sd %edx,%xmm1
  comisd %xmm1,%xmm0
  ja     L12
L52:
  cmp    %r8d,%r12d
  lea    0x44(%rsp),%rsi
  mov    %rbp,%rdi
  cmovge %r12d,%r8d
  call   <find_pitch@plt>
L66:
  cvtss2sd %xmm0,%xmm0
  addsd  D88(%rip),%xmm0
  mov    0x24(%rsp),%eax
  lea    -0x1(%rax),%edx
  cvttsd2si %xmm0,%r13d
  cmp    %r13d,%eax
  cmovle %edx,%r13d
  cmp    %r12d,%r13d
  cmovge %r13d,%r12d
  lea    (%r12,%r13,1),%ecx
  mov    %r12d,%edx
  mov    %ecx,%r15d
  shr    $0x1f,%r15d
  add    %ecx,%r15d
  sar    %r15d
  neg    %r15d
  lea    0x0(%r13,%r15,1),%ebx
  movslq %r15d,%rcx
  movslq %ebx,%rbx
  lea    0x0(,%rbx,4),%r14
  lea    0x0(%rbp,%rcx,4),%rbx
  lea    0x0(%rbp,%r14,1),%r10
  mov    %rbx,%rdi
  mov    %r10,%rsi
  mov    %r10,0x18(%rsp)
  call   <v_inner@plt>
  mov    %r14,0x10(%rsp)
  mov    %r12d,%edx
  mov    %rbx,%rdi
  lea    0x4(%rbp,%r14,1),%r14
  movss  %xmm0,0x20(%rsp)
  mov    %r14,%rsi
  call   <v_inner@plt>
  mov    0x10(%rsp),%rax
  mov    %r12d,%edx
  mov    %rbx,%rdi
  movss  %xmm0,0xc(%rsp)
  lea    -0x4(%rbp,%rax,1),%rsi
  call   <v_inner@plt>
  comiss 0xc(%rsp),%xmm0
  mov    0x18(%rsp),%r10
  jbe    L30
  sub    $0x1,%r13d
  movss  0x20(%rsp),%xmm7
  movss  %xmm0,0x20(%rsp)
  add    %r13d,%r15d
  movslq %r15d,%r15
  movss  %xmm7,0xc(%rsp)
  lea    0x4(,%r15,4),%rcx
  lea    0x0(%rbp,%rcx,1),%r14
  lea    -0x4(%rbp,%rcx,1),%r10
L30:
  mov    %r10,%rdi
  mov    %r12d,%edx
  mov    %r14,%rsi
  mov    %r10,0x18(%rsp)
  call   <v_inner@plt>
  mov    %r12d,%edx
  mov    %rbx,%rsi
  mov    %rbx,%rdi
  movss  %xmm0,0x10(%rsp)
  call   <v_inner@plt>
  mov    0x18(%rsp),%r10
  mov    %r12d,%edx
  movss  %xmm0,0x30(%rsp)
  mov    %r10,%rsi
  mov    %r10,%rdi
  call   <v_inner@plt>
  mov    %r12d,%edx
  mov    %r14,%rsi
  mov    %r14,%rdi
  movss  %xmm0,0x18(%rsp)
  call   <v_inner@plt>
  movss  0x10(%rsp),%xmm6
  movss  0x18(%rsp),%xmm2
  movss  %xmm0,0x38(%rsp)
  subss  %xmm6,%xmm2
  subss  %xmm6,%xmm0
  mulss  0xc(%rsp),%xmm2
  mulss  0x20(%rsp),%xmm0
  addss  %xmm0,%xmm2
  movaps %xmm2,%xmm0
  movss  %xmm2,0x3c(%rsp)
  call   <fabs@plt>
  comisd D86(%rip),%xmm0
  movss  D30(%rip),%xmm1
  movss  0x38(%rsp),%xmm4
  jbe    L48
  movss  0x20(%rsp),%xmm0
  movss  0xc(%rsp),%xmm1
  mulss  0x10(%rsp),%xmm0
  movss  0x3c(%rsp),%xmm2
  mulss  0x18(%rsp),%xmm1
  subss  %xmm0,%xmm1
  divss  %xmm2,%xmm1
  comiss D67(%rip),%xmm1
  jbe    L17
  movss  D67(%rip),%xmm1
L48:
  pxor   %xmm0,%xmm0
  pxor   %xmm2,%xmm2
  cvtsi2ss %r13d,%xmm0
  cvtsi2ssl 0x24(%rsp),%xmm2
  addss  %xmm0,%xmm1
  minss  %xmm1,%xmm2
  movaps %xmm2,%xmm1
  pxor   %xmm2,%xmm2
  cvtsi2ssl 0x34(%rsp),%xmm2
  maxss  %xmm1,%xmm2
  movaps %xmm2,%xmm3
  movss  %xmm2,0x24(%rsp)
  movss  D10(%rip),%xmm2
  subss  %xmm0,%xmm3
  subss  %xmm3,%xmm2
  movaps %xmm3,%xmm0
  addss  %xmm3,%xmm0
  movaps %xmm2,%xmm5
  movss  %xmm2,0x34(%rsp)
  mulss  %xmm2,%xmm0
  mulss  %xmm2,%xmm5
  mulss  0x18(%rsp),%xmm5
  movss  %xmm3,0x18(%rsp)
  mulss  0x10(%rsp),%xmm0
  addss  %xmm5,%xmm0
  movaps %xmm3,%xmm5
  mulss  %xmm3,%xmm5
  mulss  %xmm5,%xmm4
  addss  %xmm4,%xmm0
  mulss  0x30(%rsp),%xmm0
  call   <sqrt@plt>
  movss  %xmm0,0x10(%rsp)
  call   <fabs@plt>
  movss  0x24(%rsp),%xmm1
  pxor   %xmm2,%xmm2
  comisd D86(%rip),%xmm0
  jbe    L7
  movss  0x34(%rsp),%xmm2
  movss  0x18(%rsp),%xmm3
  mulss  0x20(%rsp),%xmm2
  movss  0x10(%rsp),%xmm4
  mulss  0xc(%rsp),%xmm3
  addss  %xmm3,%xmm2
  divss  %xmm4,%xmm2
L7:
  mov    0x28(%rsp),%rax
  movss  %xmm2,(%rax)
  mov    0x48(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L70
  add    $0x58,%rsp
  movaps %xmm1,%xmm0
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L12:
  cvttsd2si %xmm0,%edx
  jmp    L52
L17:
  movss  D57(%rip),%xmm0
  maxss  %xmm1,%xmm0
  movaps %xmm0,%xmm1
  jmp    L48
L70:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: label: D88
type_guessed: f64
name: extract_symbolic_data
arguments: label: D60
type_guessed: f64
name: extract_symbolic_data
arguments: label: D86
type_guessed: f64
name: extract_symbolic_data
arguments: label: D30
type_guessed: byte[4]
name: extract_symbolic_data
arguments: label: D67
type_guessed: f32
name: extract_symbolic_data
arguments: label: D10
type_guessed: f32
name: extract_symbolic_data
arguments: label: D57
type_guessed: byte[4]

tool: D88: 0.5

tool: D60: 0.75

tool: D86: 0.01

tool: {"D30": [0x00,0x00,0x00,0x3f]}

tool: D67: 2

tool: D10: 1

tool: {"D57": [0x00,0x00,0x80,0xbf]}

assistant: ```c
float frac_pch(float sig_in[], float *pcorr, float fpitch, int range, int pmin, int pmax, int lmin) {
  int length, cbegin, lower, upper, ipitch;
  float c0_0, c0_T, c0_T1, cT_T, cT_T1, cT1_T1, c0_Tm1;
  float frac, frac1;
  float corr, denom;
  if (range > 0) {
    ipitch = fpitch + 0.5;
    lower = ipitch - range;
    upper = ipitch + range;
    if (upper > pmax)
      upper = pmax;
    if (lower < pmin)
      lower = pmin;
    if (lower < 0.75 * ipitch)
      lower = 0.75 * ipitch;
    length = ipitch;
    if (length < lmin)
      length = lmin;
    fpitch = find_pitch(sig_in, &corr, lower, upper, length);
  }
  ipitch = fpitch + 0.5;
  if (ipitch >= pmax)
    ipitch = pmax - 1;
  length = ipitch;
  if (length < lmin)
    length = lmin;
  cbegin = -((length + ipitch) / 2);
  c0_T = v_inner(&sig_in[cbegin], &sig_in[cbegin + ipitch], length);
  c0_T1 = v_inner(&sig_in[cbegin], &sig_in[cbegin + ipitch + 1], length);
  c0_Tm1 = v_inner(&sig_in[cbegin], &sig_in[cbegin + ipitch - 1], length);
  if (c0_Tm1 > c0_T1) {
    c0_T1 = c0_T;
    c0_T = c0_Tm1;
    ipitch--;
  }
  cT_T1 = v_inner(&sig_in[cbegin + ipitch], &sig_in[cbegin + ipitch + 1], length);
  c0_0 = v_inner(&sig_in[cbegin], &sig_in[cbegin], length);
  cT_T = v_inner(&sig_in[cbegin + ipitch], &sig_in[cbegin + ipitch], length);
  cT1_T1 = v_inner(&sig_in[cbegin + ipitch + 1], &sig_in[cbegin + ipitch + 1], length);
  denom = c0_T1 * (cT_T - cT_T1) + c0_T * (cT1_T1 - cT_T1);
  if (fabs(denom) > 0.01)
    frac = (c0_T1 * cT_T - c0_T * cT_T1) / denom;
  else
    frac = 0.5;
  if (frac > 2.0)
    frac = 2.0;
  if (frac < -1.0)
    frac = -1.0;
  fpitch = ipitch + frac;
  if (fpitch > pmax)
    fpitch = pmax;
  if (fpitch < pmin)
    fpitch = pmin;
  frac = fpitch - ipitch;
  frac1 = 1.0 - frac;
  denom = c0_0 * (frac1 * frac1 * cT_T + 2 * frac * frac1 * cT_T1 + frac * frac * cT1_T1);
  denom = sqrt(denom);
  if (fabs(denom) > 0.01)
    *pcorr = (frac1 * c0_T + frac * c0_T1) / denom;
  else
    *pcorr = 0.0;
  return (fpitch);
}
```