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

```asm
<rdft>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x30,%rsp
  mov    %edi,-0x14(%rbp)
  mov    %esi,-0x18(%rbp)
  mov    %rdx,-0x20(%rbp)
  mov    %rcx,-0x28(%rbp)
  mov    %r8,-0x30(%rbp)
  mov    -0x28(%rbp),%rax
  mov    (%rax),%eax
  mov    %eax,-0xc(%rbp)
  mov    -0xc(%rbp),%eax
  shl    $0x2,%eax
  cmp    %eax,-0x14(%rbp)
  jle    L22
  mov    -0x14(%rbp),%eax
  sar    $0x2,%eax
  mov    %eax,-0xc(%rbp)
  mov    -0x30(%rbp),%rdx
  mov    -0x28(%rbp),%rcx
  mov    -0xc(%rbp),%eax
  mov    %rcx,%rsi
  mov    %eax,%edi
  call   <makewt@plt>
L22:
  mov    -0x28(%rbp),%rax
  mov    0x4(%rax),%eax
  mov    %eax,-0x8(%rbp)
  mov    -0x8(%rbp),%eax
  shl    $0x2,%eax
  cmp    %eax,-0x14(%rbp)
  jle    L33
  mov    -0x14(%rbp),%eax
  sar    $0x2,%eax
  mov    %eax,-0x8(%rbp)
  mov    -0xc(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rax,%rdx
  mov    -0x28(%rbp),%rcx
  mov    -0x8(%rbp),%eax
  mov    %rcx,%rsi
  mov    %eax,%edi
  call   <makect@plt>
L33:
  cmpl   $0x0,-0x18(%rbp)
  js     L8
  cmpl   $0x4,-0x14(%rbp)
  jle    L61
  mov    -0x30(%rbp),%rdi
  mov    -0xc(%rbp),%ecx
  mov    -0x28(%rbp),%rdx
  mov    -0x20(%rbp),%rsi
  mov    -0x14(%rbp),%eax
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <cftfsub@plt>
  mov    -0xc(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  lea    (%rdx,%rax,1),%rcx
  mov    -0x8(%rbp),%edx
  mov    -0x20(%rbp),%rsi
  mov    -0x14(%rbp),%eax
  mov    %eax,%edi
  call   <rftfsub@plt>
  jmp    L41
L61:
  cmpl   $0x4,-0x14(%rbp)
  jne    L41
  mov    -0x30(%rbp),%rdi
  mov    -0xc(%rbp),%ecx
  mov    -0x28(%rbp),%rdx
  mov    -0x20(%rbp),%rsi
  mov    -0x14(%rbp),%eax
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <cftfsub@plt>
L41:
  mov    -0x20(%rbp),%rax
  movss  (%rax),%xmm0
  mov    -0x20(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm1
  subss  %xmm1,%xmm0
  movss  %xmm0,-0x4(%rbp)
  mov    -0x20(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  addss  %xmm1,%xmm0
  mov    -0x20(%rbp),%rax
  movss  %xmm0,(%rax)
  mov    -0x20(%rbp),%rax
  add    $0x4,%rax
  movss  -0x4(%rbp),%xmm0
  movss  %xmm0,(%rax)
  jmp    L1
L8:
  mov    -0x20(%rbp),%rax
  movss  (%rax),%xmm0
  mov    -0x20(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  movaps %xmm0,%xmm1
  subss  %xmm2,%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x4,%rax
  movss  D0(%rip),%xmm0
  mulss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    -0x20(%rbp),%rax
  movss  (%rax),%xmm0
  mov    -0x20(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm1
  subss  %xmm1,%xmm0
  mov    -0x20(%rbp),%rax
  movss  %xmm0,(%rax)
  cmpl   $0x4,-0x14(%rbp)
  jle    L75
  mov    -0xc(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  lea    (%rdx,%rax,1),%rcx
  mov    -0x8(%rbp),%edx
  mov    -0x20(%rbp),%rsi
  mov    -0x14(%rbp),%eax
  mov    %eax,%edi
  call   <rftbsub@plt>
  mov    -0x30(%rbp),%rdi
  mov    -0xc(%rbp),%ecx
  mov    -0x28(%rbp),%rdx
  mov    -0x20(%rbp),%rsi
  mov    -0x14(%rbp),%eax
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <cftbsub@plt>
  jmp    L1
L75:
  cmpl   $0x4,-0x14(%rbp)
  jne    L1
  mov    -0x30(%rbp),%rdi
  mov    -0xc(%rbp),%ecx
  mov    -0x28(%rbp),%rdx
  mov    -0x20(%rbp),%rsi
  mov    -0x14(%rbp),%eax
  mov    %rdi,%r8
  mov    %eax,%edi
  call   <cftbsub@plt>
L1:
  nop
  leave
  ret
```

tool_calls: name: unpack_data
arguments: tag: D0
type: f32

tool: D0: 0.5

assistant: ```c
void rdft(int n, int isgn, float *a, int *ip, float *w) {
  void makewt(int nw, int *ip, float *w);
  void makect(int nc, int *ip, float *c);
  void cftfsub(int n, float *a, int *ip, int nw, float *w);
  void cftbsub(int n, float *a, int *ip, int nw, float *w);
  void rftfsub(int n, float *a, int nc, float *c);
  void rftbsub(int n, float *a, int nc, float *c);
  int nw, nc;
  float xi;
  nw = ip[0];
  if (n > (nw << 2)) {
    nw = n >> 2;
    makewt(nw, ip, w);
  }
  nc = ip[1];
  if (n > (nc << 2)) {
    nc = n >> 2;
    makect(nc, ip, w + nw);
  }
  if (isgn >= 0) {
    if (n > 4) {
      cftfsub(n, a, ip, nw, w);
      rftfsub(n, a, nc, w + nw);
    } else if (n == 4) {
      cftfsub(n, a, ip, nw, w);
    }
    xi = a[0] - a[1];
    a[0] += a[1];
    a[1] = xi;
  } else {
    a[1] = ((float)0.5) * (a[0] - a[1]);
    a[0] -= a[1];
    if (n > 4) {
      rftbsub(n, a, nc, w + nw);
      cftbsub(n, a, ip, nw, w);
    } else if (n == 4) {
      cftbsub(n, a, ip, nw, w);
    }
  }
}
```