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

```asm
<DFT64>:
  endbr64
  lea    0x100(%rdi),%rdx
  mov    %rdi,%rax
  xchg   %ax,%ax
L53:
  movss  (%rax),%xmm0
  movss  0x4(%rax),%xmm3
  add    $0x10,%rax
  movss  -0x8(%rax),%xmm1
  movss  -0x4(%rax),%xmm4
  movaps %xmm0,%xmm2
  subss  %xmm3,%xmm0
  addss  %xmm3,%xmm2
  movaps %xmm1,%xmm3
  subss  %xmm4,%xmm1
  addss  %xmm4,%xmm3
  movaps %xmm2,%xmm4
  subss  %xmm3,%xmm2
  addss  %xmm3,%xmm4
  movss  %xmm2,-0x8(%rax)
  movaps %xmm0,%xmm2
  subss  %xmm1,%xmm0
  addss  %xmm1,%xmm2
  movss  %xmm4,-0x10(%rax)
  movss  %xmm0,-0x4(%rax)
  movss  %xmm2,-0xc(%rax)
  cmp    %rax,%rdx
  jne    L53
  movss  D4(%rip),%xmm4
  mov    %rdi,%rax
  nopw   0x0(%rax,%rax,1)
L24:
  movss  0x10(%rax),%xmm1
  movss  (%rax),%xmm0
  add    $0x20,%rax
  movss  -0x4(%rax),%xmm3
  movaps %xmm0,%xmm2
  subss  %xmm1,%xmm0
  addss  %xmm1,%xmm2
  movss  -0x8(%rax),%xmm1
  movss  %xmm0,-0x10(%rax)
  movss  -0x18(%rax),%xmm0
  movss  %xmm2,-0x20(%rax)
  movaps %xmm0,%xmm2
  subss  %xmm1,%xmm0
  addss  %xmm1,%xmm2
  movss  -0x1c(%rax),%xmm1
  movss  %xmm0,-0x8(%rax)
  movss  -0xc(%rax),%xmm0
  movss  %xmm2,-0x18(%rax)
  movaps %xmm0,%xmm2
  subss  %xmm3,%xmm0
  addss  %xmm3,%xmm2
  movaps %xmm1,%xmm3
  mulss  %xmm4,%xmm0
  mulss  %xmm4,%xmm2
  subss  %xmm2,%xmm1
  addss  %xmm2,%xmm3
  movaps %xmm0,%xmm2
  movss  %xmm1,-0xc(%rax)
  movss  -0x14(%rax),%xmm1
  movss  %xmm3,-0x1c(%rax)
  addss  %xmm1,%xmm2
  subss  %xmm0,%xmm1
  movss  %xmm2,-0x14(%rax)
  movss  %xmm1,-0x4(%rax)
  cmp    %rax,%rdx
  jne    L24
  mov    %rdi,%r9
  xor    %r10d,%r10d
  lea    0x50(%rsi),%r11
  lea    0x80(%rsi),%r8
L62:
  movss  0x20(%r9),%xmm1
  mov    %r11,%rcx
  mov    %r9,%rdx
  mov    %r9,%rax
  movss  (%r9),%xmm0
  movaps %xmm0,%xmm2
  subss  %xmm1,%xmm0
  addss  %xmm1,%xmm2
  movss  0x30(%r9),%xmm1
  movss  %xmm0,0x20(%r9)
  movss  0x10(%r9),%xmm0
  movss  %xmm2,(%r9)
  movaps %xmm0,%xmm2
  subss  %xmm1,%xmm0
  addss  %xmm1,%xmm2
  movss  %xmm0,0x30(%r9)
  movss  %xmm2,0x10(%r9)
L1:
  movss  0x24(%rax),%xmm0
  movss  (%rcx),%xmm5
  add    $0x10,%rcx
  add    $0x4,%rax
  movss  0x3c(%rdx),%xmm4
  movss  -0x50(%rcx),%xmm6
  sub    $0x4,%rdx
  movaps %xmm0,%xmm3
  movss  (%rax),%xmm2
  movss  0x20(%rdx),%xmm1
  mulss  %xmm5,%xmm3
  movaps %xmm4,%xmm7
  mulss  %xmm6,%xmm7
  mulss  %xmm5,%xmm4
  mulss  %xmm6,%xmm0
  addss  %xmm7,%xmm3
  subss  %xmm4,%xmm0
  movaps %xmm2,%xmm4
  subss  %xmm3,%xmm2
  addss  %xmm3,%xmm4
  movss  %xmm2,0x20(%rax)
  movaps %xmm1,%xmm2
  subss  %xmm0,%xmm1
  addss  %xmm0,%xmm2
  movss  %xmm4,(%rax)
  movss  %xmm1,0x40(%rdx)
  movss  %xmm2,0x20(%rdx)
  cmp    %rcx,%r8
  jne    L1
  add    $0x10,%r10d
  add    $0x40,%r9
  cmp    $0x40,%r10d
  jne    L62
  mov    %rdi,%r9
  xor    %r10d,%r10d
  lea    0x48(%rsi),%r11
L31:
  movss  0x40(%r9),%xmm1
  mov    %r11,%rcx
  mov    %r9,%rdx
  mov    %r9,%rax
  movss  (%r9),%xmm0
  movaps %xmm0,%xmm2
  subss  %xmm1,%xmm0
  addss  %xmm1,%xmm2
  movss  0x60(%r9),%xmm1
  movss  %xmm0,0x40(%r9)
  movss  0x20(%r9),%xmm0
  movss  %xmm2,(%r9)
  movaps %xmm0,%xmm2
  subss  %xmm1,%xmm0
  addss  %xmm1,%xmm2
  movss  %xmm0,0x60(%r9)
  movss  %xmm2,0x20(%r9)
  cs nopw 0x0(%rax,%rax,1)
L34:
  movss  0x44(%rax),%xmm0
  movss  (%rcx),%xmm5
  add    $0x8,%rcx
  add    $0x4,%rax
  movss  0x7c(%rdx),%xmm4
  movss  -0x48(%rcx),%xmm6
  sub    $0x4,%rdx
  movaps %xmm0,%xmm3
  movss  (%rax),%xmm2
  movss  0x40(%rdx),%xmm1
  mulss  %xmm5,%xmm3
  movaps %xmm4,%xmm7
  mulss  %xmm6,%xmm7
  mulss  %xmm5,%xmm4
  mulss  %xmm6,%xmm0
  addss  %xmm7,%xmm3
  subss  %xmm4,%xmm0
  movaps %xmm2,%xmm4
  subss  %xmm3,%xmm2
  addss  %xmm3,%xmm4
  movss  %xmm2,0x40(%rax)
  movaps %xmm1,%xmm2
  subss  %xmm0,%xmm1
  addss  %xmm0,%xmm2
  movss  %xmm4,(%rax)
  movss  %xmm1,0x80(%rdx)
  movss  %xmm2,0x40(%rdx)
  cmp    %rcx,%r8
  jne    L34
  add    $0x20,%r10d
  sub    $0xffffffffffffff80,%r9
  cmp    $0x40,%r10d
  jne    L31
  movss  0x80(%rdi),%xmm1
  movss  (%rdi),%xmm0
  lea    0x7c(%rdi),%rax
  lea    0x84(%rdi),%rdx
  add    $0x44,%rsi
  add    $0x40,%rdi
  movaps %xmm0,%xmm2
  subss  %xmm1,%xmm0
  addss  %xmm1,%xmm2
  movss  0x80(%rdi),%xmm1
  movss  %xmm0,0x40(%rdi)
  movss  (%rdi),%xmm0
  movss  %xmm2,-0x40(%rdi)
  movaps %xmm0,%xmm2
  subss  %xmm1,%xmm0
  addss  %xmm1,%xmm2
  movss  %xmm0,0x80(%rdi)
  movss  %xmm2,(%rdi)
L4:
  movss  (%rdx),%xmm0
  movss  (%rsi),%xmm5
  sub    $0x4,%rax
  add    $0x4,%rsi
  movss  -0x44(%rsi),%xmm6
  movss  -0x80(%rdx),%xmm2
  add    $0x4,%rdx
  movss  0x84(%rax),%xmm4
  movaps %xmm0,%xmm3
  movss  0x4(%rax),%xmm1
  mulss  %xmm5,%xmm3
  mulss  %xmm6,%xmm0
  movaps %xmm4,%xmm7
  mulss  %xmm6,%xmm7
  mulss  %xmm5,%xmm4
  addss  %xmm7,%xmm3
  subss  %xmm4,%xmm0
  movaps %xmm2,%xmm4
  subss  %xmm3,%xmm2
  addss  %xmm3,%xmm4
  movss  %xmm2,-0x4(%rdx)
  movaps %xmm1,%xmm2
  subss  %xmm0,%xmm1
  addss  %xmm0,%xmm2
  movss  %xmm4,-0x84(%rdx)
  movss  %xmm1,0x84(%rax)
  movss  %xmm2,0x4(%rax)
  cmp    %rax,%rdi
  jne    L4
  ret
```

tool_calls: name: disassemble_data
arguments: tag: D4
type_guessed: byte[4]

tool: {"D4": [0xf3,0x04,0x35,0x3f]}

assistant: ```c
void DFT64(float *A, const float *sinTab) {
  int i, j, theta;
  float alpha, beta, beta1, beta2, y0, y1, y2, y3;
  for (i = 0; i < 64; i += 4) {
    alpha = A[i];
    beta = A[i + 1];
    beta1 = A[i + 2];
    beta2 = A[i + 3];
    y0 = alpha + beta;
    y1 = alpha - beta;
    y2 = beta1 + beta2;
    y3 = beta1 - beta2;
    A[i] = y0 + y2;
    A[i + 2] = y0 - y2;
    A[i + 1] = y1 + y3;
    A[i + 3] = y1 - y3;
  }
  for (i = 0; i < 64; i += 8) {
    alpha = A[i];
    beta = A[i + 4];
    A[i] = alpha + beta;
    A[i + 4] = alpha - beta;
    alpha = A[i + 2];
    beta = A[i + 6];
    A[i + 2] = alpha + beta;
    A[i + 6] = alpha - beta;
    alpha = A[i + 1];
    beta1 = 0.70710678118654752440084436210485f * (A[i + 5] + A[i + 7]);
    beta2 = 0.70710678118654752440084436210485f * (A[i + 5] - A[i + 7]);
    A[i + 1] = alpha + beta1;
    A[i + 5] = alpha - beta1;
    alpha = A[i + 3];
    A[i + 3] = alpha + beta2;
    A[i + 7] = alpha - beta2;
  }
  for (i = 0; i < 64; i += 16) {
    theta = 4;
    alpha = A[i];
    beta = A[i + 8];
    A[i] = alpha + beta;
    A[i + 8] = alpha - beta;
    alpha = A[i + 4];
    beta = A[i + 12];
    A[i + 4] = alpha + beta;
    A[i + 12] = alpha - beta;
    for (j = 1; j < 4; j++) {
      y0 = sinTab[theta];
      y1 = sinTab[theta + 16];
      alpha = A[i + j];
      beta = A[i - j + 8];
      beta1 = A[i + j + 8] * y1 + A[i - j + 16] * y0;
      beta2 = A[i + j + 8] * y0 - A[i - j + 16] * y1;
      theta += 4;
      A[i + j] = alpha + beta1;
      A[i + j + 8] = alpha - beta1;
      A[i - j + 8] = beta + beta2;
      A[i - j + 16] = beta - beta2;
    }
  }
  for (i = 0; i < 64; i += 32) {
    theta = 2;
    alpha = A[i];
    beta = A[i + 16];
    A[i] = alpha + beta;
    A[i + 16] = alpha - beta;
    alpha = A[i + 8];
    beta = A[i + 24];
    A[i + 8] = alpha + beta;
    A[i + 24] = alpha - beta;
    for (j = 1; j < 8; j++) {
      y0 = sinTab[theta];
      y1 = sinTab[theta + 16];
      alpha = A[i + j];
      beta = A[i - j + 16];
      beta1 = A[i + j + 16] * y1 + A[i - j + 32] * y0;
      beta2 = A[i + j + 16] * y0 - A[i - j + 32] * y1;
      theta += 2;
      A[i + j] = alpha + beta1;
      A[i + j + 16] = alpha - beta1;
      A[i - j + 16] = beta + beta2;
      A[i - j + 32] = beta - beta2;
    }
  }
  theta = 1;
  alpha = A[0];
  beta = A[32];
  A[0] = alpha + beta;
  A[32] = alpha - beta;
  alpha = A[16];
  beta = A[48];
  A[16] = alpha + beta;
  A[48] = alpha - beta;
  for (j = 1; j < 16; j++) {
    y0 = sinTab[theta];
    y1 = sinTab[theta + 16];
    alpha = A[j];
    beta = A[-j + 32];
    beta1 = A[j + 32] * y1 + A[-j + 64] * y0;
    beta2 = A[j + 32] * y0 - A[-j + 64] * y1;
    theta += 1;
    A[j] = alpha + beta1;
    A[j + 32] = alpha - beta1;
    A[-j + 32] = beta + beta2;
    A[-j + 64] = beta - beta2;
  }
}
```