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

```asm
<makehanning>:
  endbr64
  push   %r15
  mov    %rsi,%r15
  push   %r14
  mov    %rdx,%r14
  push   %r13
  mov    %r9d,%r13d
  push   %r12
  mov    %rdi,%r12
  push   %rbp
  mov    %r8d,%ebp
  push   %rbx
  mov    %ecx,%ebx
  sub    $0x48,%rsp
  mov    0x80(%rsp),%eax
  test   %eax,%eax
  jne    L142
  test   %ecx,%ecx
  jle    L85
  lea    -0x1(%rcx),%eax
  pxor   %xmm7,%xmm7
  movsd  D17(%rip),%xmm4
  movsd  D2(%rip),%xmm3
  cvtsi2sd %eax,%xmm7
  movslq %ecx,%rax
  xor    %ecx,%ecx
  movsd  %xmm4,0x18(%rsp)
  movsd  D21(%rip),%xmm4
  mov    %rax,0x30(%rsp)
  movsd  %xmm3,0x8(%rsp)
  movsd  %xmm7,0x38(%rsp)
  movsd  D0(%rip),%xmm7
  movsd  %xmm4,0x10(%rsp)
  movsd  %xmm7,0x20(%rsp)
  nopl   (%rax)
L189:
  pxor   %xmm0,%xmm0
  mov    %rcx,0x28(%rsp)
  cvtsi2sd %ecx,%xmm0
  mulsd  0x18(%rsp),%xmm0
  divsd  0x38(%rsp),%xmm0
  addsd  0x8(%rsp),%xmm0
  call   <cos@plt>
  pxor   %xmm0,%xmm0
  mov    0x28(%rsp),%rcx
  cvtsi2sd %eax,%xmm0
  addsd  0x20(%rsp),%xmm0
  mulsd  0x10(%rsp),%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%r14,%rcx,4)
  movss  %xmm0,(%r15,%rcx,4)
  movss  %xmm0,(%r12,%rcx,4)
  add    $0x1,%rcx
  cmp    %rcx,0x30(%rsp)
  jne    L189
L10:
  pxor   %xmm3,%xmm3
  cmp    %ebp,%ebx
  jle    L120
  mov    $0x1,%eax
  pxor   %xmm1,%xmm1
  pxor   %xmm3,%xmm3
  xor    %r12d,%r12d
  sub    %ebx,%eax
  cvtsi2sd %eax,%xmm1
  mulsd  0x10(%rsp),%xmm1
  cvtsd2ss %xmm1,%xmm1
  jmp    L63
  nopl   0x0(%rax)
L31:
  add    $0x1,%r12
  addss  D4(%rip),%xmm1
  cmp    %r12d,%ebx
  jle    L120
L63:
  ucomiss %xmm3,%xmm1
  jp     L14
  je     L31
L14:
  pxor   %xmm2,%xmm2
  pxor   %xmm4,%xmm4
  movss  %xmm3,0x20(%rsp)
  cvtss2sd %xmm1,%xmm2
  mulsd  0x8(%rsp),%xmm2
  cvtsi2sd %ebp,%xmm4
  movss  %xmm1,0x18(%rsp)
  movapd %xmm2,%xmm0
  movsd  %xmm2,0x10(%rsp)
  divsd  %xmm4,%xmm0
  call   <sin@plt>
  pxor   %xmm0,%xmm0
  pxor   %xmm4,%xmm4
  movsd  0x10(%rsp),%xmm2
  imul   %ebp,%eax
  cvtss2sd (%r15,%r12,4),%xmm4
  test   %r13d,%r13d
  movss  0x18(%rsp),%xmm1
  movss  0x20(%rsp),%xmm3
  cvtsi2sd %eax,%xmm0
  divsd  %xmm2,%xmm0
  mulsd  %xmm4,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%r15,%r12,4)
  je     L31
  pxor   %xmm4,%xmm4
  movapd %xmm2,%xmm0
  cvtsi2sd %r13d,%xmm4
  divsd  %xmm4,%xmm0
  call   <sin@plt>
  pxor   %xmm0,%xmm0
  movsd  0x10(%rsp),%xmm2
  movss  0x20(%rsp),%xmm3
  imul   %r13d,%eax
  movss  0x18(%rsp),%xmm1
  cvtsi2sd %eax,%xmm0
  divsd  %xmm2,%xmm0
  pxor   %xmm2,%xmm2
  cvtss2sd (%r14,%r12,4),%xmm2
  mulsd  %xmm2,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%r14,%r12,4)
  jmp    L31
  nopl   0x0(%rax,%rax,1)
L120:
  xor    %eax,%eax
  movaps %xmm3,%xmm0
  nopl   (%rax)
L164:
  addss  (%r15,%rax,4),%xmm0
  add    $0x1,%rax
  cmp    %eax,%ebx
  jg     L164
  movss  D52(%rip),%xmm1
  movss  D4(%rip),%xmm4
  xor    %eax,%eax
  divss  %xmm0,%xmm1
  nopl   0x0(%rax)
L66:
  movaps %xmm1,%xmm0
  cmp    %ebp,%ebx
  jle    L62
  movaps %xmm4,%xmm0
  divss  %xmm1,%xmm0
L62:
  movss  (%r15,%rax,4),%xmm2
  mulss  %xmm1,%xmm2
  movss  %xmm2,(%r15,%rax,4)
  mulss  (%r14,%rax,4),%xmm0
  movss  %xmm0,(%r14,%rax,4)
  add    $0x1,%rax
  cmp    %eax,%ebx
  jg     L66
  cmp    %ebp,%ebx
  jg     L85
  test   %r13d,%r13d
  je     L85
  movslq %r13d,%rsi
  mov    %r14,%rcx
  movaps %xmm3,%xmm1
  xor    %eax,%eax
  shl    $0x2,%rsi
  nopw   0x0(%rax,%rax,1)
L9:
  movss  (%rcx),%xmm0
  add    %r13d,%eax
  add    %rsi,%rcx
  mulss  %xmm0,%xmm0
  addss  %xmm0,%xmm1
  cmp    %eax,%ebx
  jg     L9
  movss  D4(%rip),%xmm2
  xor    %eax,%eax
  divss  %xmm1,%xmm2
  nopl   0x0(%rax)
L167:
  movss  (%r14,%rax,4),%xmm0
  mulss  %xmm2,%xmm0
  movss  %xmm0,(%r14,%rax,4)
  add    $0x1,%rax
  cmp    %eax,%ebx
  jg     L167
L85:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nop
L142:
  test   %ecx,%ecx
  jle    L85
  lea    -0x1(%rcx),%eax
  pxor   %xmm3,%xmm3
  movsd  D2(%rip),%xmm7
  movsd  D0(%rip),%xmm4
  cvtsi2sd %eax,%xmm3
  movslq %ecx,%rax
  xor    %ecx,%ecx
  mov    %rax,0x30(%rsp)
  movsd  %xmm7,0x8(%rsp)
  movsd  %xmm4,0x20(%rsp)
  movsd  %xmm3,0x38(%rsp)
  movsd  D17(%rip),%xmm3
  movsd  %xmm3,0x18(%rsp)
  movsd  D21(%rip),%xmm3
  movsd  %xmm3,0x10(%rsp)
  nopw   0x0(%rax,%rax,1)
L100:
  pxor   %xmm0,%xmm0
  mov    %rcx,0x28(%rsp)
  cvtsi2sd %ecx,%xmm0
  mulsd  0x18(%rsp),%xmm0
  divsd  0x38(%rsp),%xmm0
  addsd  0x8(%rsp),%xmm0
  call   <cos@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  addsd  0x20(%rsp),%xmm0
  mulsd  0x10(%rsp),%xmm0
  call   <sqrt@plt>
  mov    0x28(%rsp),%rcx
  movss  %xmm0,(%r14,%rcx,4)
  movss  %xmm0,(%r15,%rcx,4)
  movss  %xmm0,(%r12,%rcx,4)
  add    $0x1,%rcx
  cmp    %rcx,0x30(%rsp)
  jne    L100
  jmp    L10
```

tool_calls: name: read_data
arguments: data_label: D17
type: f64
name: read_data
arguments: data_label: D2
type: f64
name: read_data
arguments: data_label: D21
type: f64
name: read_data
arguments: data_label: D0
type: f64
name: read_data
arguments: data_label: D4
type: byte[4]
name: read_data
arguments: data_label: D52
type: f32

tool: D17: 6.2831853072

tool: D2: 3.14159265358979

tool: D21: 0.5

tool: D0: 1

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

tool: D52: 2

assistant: ```c
void makehanning(float *H, float *A, float *S, int Nw, int N, int I, int odd) {
  int i;
  float sum;
  if (odd) {
    for (i = 0; i < Nw; i++)
      H[i] = A[i] = S[i] = sqrt(0.5 * (1. + cos(3.14159265358979 + 6.2831853072 * i / (Nw - 1))));
  } else {
    for (i = 0; i < Nw; i++)
      H[i] = A[i] = S[i] = 0.5 * (1. + cos(3.14159265358979 + 6.2831853072 * i / (Nw - 1)));
  }
  if (Nw > N) {
    float x;
    x = -(Nw - 1) / 2.;
    for (i = 0; i < Nw; i++, x += 1.)
      if (x != 0.) {
        A[i] *= N * sin(3.14159265358979 * x / N) / (3.14159265358979 * x);
        if (I)
          S[i] *= I * sin(3.14159265358979 * x / I) / (3.14159265358979 * x);
      }
  }
  for (sum = i = 0; i < Nw; i++)
    sum += A[i];
  for (i = 0; i < Nw; i++) {
    float afac = 2. / sum;
    float sfac = Nw > N ? 1. / afac : afac;
    A[i] *= afac;
    S[i] *= sfac;
  }
  if (Nw <= N && I) {
    for (sum = i = 0; i < Nw; i += I)
      sum += S[i] * S[i];
    for (sum = 1. / sum, i = 0; i < Nw; i++)
      S[i] *= sum;
  }
}
```