user: What would the C source look like for the assembly code provided here:

```asm
<estimate_noise_loizou>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x28,%rsp
  mov    %rdi,%rax
  mov    %esi,%ebx
  mov    %r9,%r15
  mov    0x10(%rbp),%rdi
  mov    0x18(%rbp),%rsi
  mov    0x20(%rbp),%r9
  mov    0x28(%rbp),%r10
  mov    0x30(%rbp),%r11
  mov    %fs:0x28,%r12
  mov    %r12,-0x38(%rbp)
  xor    %r12d,%r12d
  lea    0x1(%rbx),%r12d
  movslq %r12d,%r12
  shl    $0x2,%r12
  lea    0xf(%r12),%r13
  mov    %r13,%r14
  and    $0xfffffffffffffff0,%r14
  mov    %r14,-0x48(%rbp)
  and    $0xfffffffffffff000,%r13
  mov    %rsp,%r14
  sub    %r13,%r14
L73:
  cmp    %r14,%rsp
  je     L10
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L73
L10:
  mov    -0x48(%rbp),%r14
  mov    %r14,%r13
  and    $0xfff,%r13d
  sub    %r13,%rsp
  test   %r13,%r13
  je     L59
  orq    $0x0,-0x8(%rsp,%r13,1)
L59:
  mov    %rsp,%r13
  add    $0xf,%r12
  mov    %r12,%r14
  and    $0xfffffffffffffff0,%r14
  mov    %r14,-0x48(%rbp)
  and    $0xfffffffffffff000,%r12
  mov    %rsp,%r14
  sub    %r12,%r14
  mov    %r14,%r12
  mov    -0x48(%rbp),%r14
L77:
  cmp    %r12,%rsp
  je     L65
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L77
L65:
  mov    %r14,-0x48(%rbp)
  and    $0xfff,%r14d
  sub    %r14,%rsp
  test   %r14,%r14
  je     L88
  orq    $0x0,-0x8(%rsp,%r14,1)
L88:
  mov    %rsp,%r12
  test   %ebx,%ebx
  js     L119
  lea    0x1(%rbx),%r14d
  shl    $0x2,%r14
  mov    $0x0,%ebx
  movss  D1(%rip),%xmm8
  movss  D86(%rip),%xmm7
  movss  D89(%rip),%xmm11
  movss  D10(%rip),%xmm10
  movss  D97(%rip),%xmm9
  movss  D75(%rip),%xmm6
  pxor   %xmm5,%xmm5
  movss  D64(%rip),%xmm4
  movaps %xmm11,%xmm3
  movss  D36(%rip),%xmm2
  jmp    L84
L36:
  movss  %xmm0,(%rsi,%rbx,1)
L49:
  movss  (%rcx,%rbx,1),%xmm12
  divss  (%rsi,%rbx,1),%xmm12
  movaps %xmm6,%xmm1
  movss  (%rax,%rbx,1),%xmm0
  cmpltss %xmm12,%xmm0
  andps  %xmm0,%xmm1
  andnps %xmm5,%xmm0
  orps   %xmm0,%xmm1
  movss  %xmm1,(%r12,%rbx,1)
  movaps %xmm4,%xmm0
  mulss  (%r11,%rbx,1),%xmm0
  mulss  %xmm3,%xmm1
  addss  %xmm1,%xmm0
  movss  %xmm0,(%r10,%rbx,1)
  mulss  %xmm2,%xmm0
  addss  D47(%rip),%xmm0
  movss  %xmm0,0x0(%r13,%rbx,1)
  movss  D75(%rip),%xmm1
  subss  %xmm0,%xmm1
  mulss  (%rdx,%rbx,1),%xmm1
  mulss  (%rdi,%rbx,1),%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,(%r15,%rbx,1)
  add    $0x4,%rbx
  cmp    %r14,%rbx
  je     L119
L84:
  movaps %xmm8,%xmm0
  mulss  (%r8,%rbx,1),%xmm0
  movaps %xmm7,%xmm1
  mulss  (%rdx,%rbx,1),%xmm1
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rcx,%rbx,1)
  movss  (%r9,%rbx,1),%xmm1
  comiss %xmm1,%xmm0
  jbe    L36
  movaps %xmm11,%xmm12
  mulss  (%r8,%rbx,1),%xmm12
  subss  %xmm12,%xmm0
  mulss  %xmm10,%xmm0
  mulss  %xmm9,%xmm1
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rsi,%rbx,1)
  jmp    L49
L119:
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L101
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
L101:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: data_label: D1
datatype: byte[4]
name: unpack_data
arguments: data_label: D86
datatype: byte[4]
name: unpack_data
arguments: data_label: D89
datatype: byte[4]
name: unpack_data
arguments: data_label: D10
datatype: byte[4]
name: unpack_data
arguments: data_label: D97
datatype: byte[4]
name: unpack_data
arguments: data_label: D75
datatype: byte[4]
name: unpack_data
arguments: data_label: D64
datatype: byte[4]
name: unpack_data
arguments: data_label: D36
datatype: byte[4]
name: unpack_data
arguments: data_label: D47
datatype: byte[4]

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

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

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

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

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

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

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

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

tool: {"D47": [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];
  }
}
```