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

```asm
<estimate_noise_loizou>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  sub    $0xa0,%rsp
  mov    %rdi,-0x58(%rbp)
  mov    %esi,-0x5c(%rbp)
  mov    %rdx,-0x68(%rbp)
  mov    %rcx,-0x70(%rbp)
  mov    %r8,-0x78(%rbp)
  mov    %r9,-0x80(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0x88(%rbp)
  mov    0x18(%rbp),%rax
  mov    %rax,-0x90(%rbp)
  mov    0x20(%rbp),%rax
  mov    %rax,-0x98(%rbp)
  mov    0x28(%rbp),%rax
  mov    %rax,-0xa0(%rbp)
  mov    0x30(%rbp),%rax
  mov    %rax,-0xa8(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x28(%rbp)
  xor    %eax,%eax
  mov    %rsp,%rax
  mov    %rax,%rsi
  pxor   %xmm0,%xmm0
  movss  %xmm0,-0x4c(%rbp)
  mov    -0x5c(%rbp),%eax
  add    $0x1,%eax
  movslq %eax,%rdx
  sub    $0x1,%rdx
  mov    %rdx,-0x48(%rbp)
  movslq %eax,%rdx
  mov    %rdx,-0xc0(%rbp)
  movq   $0x0,-0xb8(%rbp)
  movslq %eax,%rdx
  mov    %rdx,%r14
  mov    $0x0,%r15d
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%edi
  mov    $0x0,%edx
  div    %rdi
  imul   $0x10,%rax,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L54:
  cmp    %rdx,%rsp
  je     L75
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L54
L75:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L98
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L98:
  mov    %rsp,%rax
  add    $0x3,%rax
  shr    $0x2,%rax
  shl    $0x2,%rax
  mov    %rax,-0x40(%rbp)
  mov    -0x5c(%rbp),%eax
  add    $0x1,%eax
  movslq %eax,%rdx
  sub    $0x1,%rdx
  mov    %rdx,-0x38(%rbp)
  movslq %eax,%rdx
  mov    %rdx,%r12
  mov    $0x0,%r13d
  movslq %eax,%rdx
  mov    %rdx,%r10
  mov    $0x0,%r11d
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%edi
  mov    $0x0,%edx
  div    %rdi
  imul   $0x10,%rax,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L16:
  cmp    %rdx,%rsp
  je     L79
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L16
L79:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L115
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L115:
  mov    %rsp,%rax
  add    $0x3,%rax
  shr    $0x2,%rax
  shl    $0x2,%rax
  mov    %rax,-0x30(%rbp)
  movl   $0x0,-0x50(%rbp)
  jmp    L121
L57:
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  movss  D21(%rip),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm2
  movss  D86(%rip),%xmm0
  mulss  %xmm2,%xmm0
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x98(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  comiss %xmm1,%xmm0
  jbe    L114
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x98(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  movss  D88(%rip),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm3
  movss  D97(%rip),%xmm2
  mulss  %xmm2,%xmm3
  movaps %xmm0,%xmm2
  subss  %xmm3,%xmm2
  movss  D19(%rip),%xmm0
  mulss  %xmm2,%xmm0
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  jmp    L94
L114:
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x70(%rbp),%rax
  add    %rax,%rdx
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x90(%rbp),%rax
  add    %rcx,%rax
  movss  (%rdx),%xmm0
  movss  %xmm0,(%rax)
L94:
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x90(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  divss  %xmm1,%xmm0
  movss  %xmm0,-0x4c(%rbp)
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  movss  -0x4c(%rbp),%xmm0
  comiss %xmm1,%xmm0
  jbe    L88
  mov    -0x30(%rbp),%rax
  mov    -0x50(%rbp),%edx
  movslq %edx,%rdx
  movss  D33(%rip),%xmm0
  movss  %xmm0,(%rax,%rdx,4)
  jmp    L93
L88:
  mov    -0x30(%rbp),%rax
  mov    -0x50(%rbp),%edx
  movslq %edx,%rdx
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax,%rdx,4)
L93:
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0xa8(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  movss  D39(%rip),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x30(%rbp),%rax
  mov    -0x50(%rbp),%edx
  movslq %edx,%rdx
  movss  (%rax,%rdx,4),%xmm2
  movss  D97(%rip),%xmm0
  mulss  %xmm2,%xmm0
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0xa0(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm1
  movss  D23(%rip),%xmm0
  mulss  %xmm0,%xmm1
  movss  D72(%rip),%xmm0
  addss  %xmm1,%xmm0
  mov    -0x40(%rbp),%rax
  mov    -0x50(%rbp),%edx
  movslq %edx,%rdx
  movss  %xmm0,(%rax,%rdx,4)
  mov    -0x40(%rbp),%rax
  mov    -0x50(%rbp),%edx
  movslq %edx,%rdx
  movss  (%rax,%rdx,4),%xmm1
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x88(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x40(%rbp),%rax
  mov    -0x50(%rbp),%edx
  movslq %edx,%rdx
  movss  (%rax,%rdx,4),%xmm3
  movss  D33(%rip),%xmm0
  movaps %xmm0,%xmm2
  subss  %xmm3,%xmm2
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x68(%rbp),%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  addl   $0x1,-0x50(%rbp)
L121:
  mov    -0x50(%rbp),%eax
  cmp    -0x5c(%rbp),%eax
  jle    L57
  mov    %rsi,%rsp
  nop
  mov    -0x28(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L17
  call   <__stack_chk_fail@plt>
L17:
  lea    -0x20(%rbp),%rsp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
```

tool_calls: name: read_data
arguments: data_label: D21
type_guessed: byte[4]
name: read_data
arguments: data_label: D86
type_guessed: byte[4]
name: read_data
arguments: data_label: D88
type_guessed: byte[4]
name: read_data
arguments: data_label: D97
type_guessed: byte[4]
name: read_data
arguments: data_label: D19
type_guessed: byte[4]
name: read_data
arguments: data_label: D33
type_guessed: byte[4]
name: read_data
arguments: data_label: D39
type_guessed: byte[4]
name: read_data
arguments: data_label: D23
type_guessed: byte[4]
name: read_data
arguments: data_label: D72
type_guessed: byte[4]

tool: {"D21": [0x33,0x33,0x33,0x3f]}

tool: {"D86": [0x9a,0x99,0x99,0x3e]}

tool: {"D88": [0xee,0x7c,0x7f,0x3f]}

tool: {"D97": [0xcd,0xcc,0x4c,0x3f]}

tool: {"D19": [0x81,0xd6,0x23,0x3c]}

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

tool: {"D39": [0xcd,0xcc,0x4c,0x3e]}

tool: {"D23": [0xd0,0xcc,0x4c,0x3d]}

tool: {"D72": [0x33,0x33,0x73,0x3f]}

assistant: ```c
void estimate_noise_loizou(float *thresh, int fft_size_2, float *p2, float *s_pow_spec, float *prev_s_pow_spec, float *noise_thresholds_p2, float *prev_noise, float *p_min, float *prev_p_min, float *speech_p_p, float *prev_speech_p_p) {
  int k;
  float ratio_ns = 0.f;
  float freq_s[fft_size_2 + 1];
  float speech_p_d[fft_size_2 + 1];
  for (k = 0; k <= fft_size_2; k++) {
    s_pow_spec[k] = 0.7f * prev_s_pow_spec[k] + (1.f - 0.7f) * p2[k];
    if (prev_p_min[k] < s_pow_spec[k]) {
      p_min[k] = 0.998f * prev_p_min[k] + ((1.f - 0.998f) / (1.f - 0.8f)) * (s_pow_spec[k] - 0.8f * prev_s_pow_spec[k]);
    } else {
      p_min[k] = s_pow_spec[k];
    }
    ratio_ns = s_pow_spec[k] / p_min[k];
    if (ratio_ns > thresh[k]) {
      speech_p_d[k] = 1.f;
    } else {
      speech_p_d[k] = 0.f;
    }
    speech_p_p[k] = 0.2f * prev_speech_p_p[k] + (1.f - 0.2f) * speech_p_d[k];
    freq_s[k] = 0.95f + (1.f - 0.95f) * speech_p_p[k];
    noise_thresholds_p2[k] = freq_s[k] * prev_noise[k] + (1.f - freq_s[k]) * p2[k];
  }
}
```