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

```asm
<step1>:
  endbr64
  push   %r15
  mov    %r9,%rax
  movss  D1(%rip),%xmm0
  mov    %rdi,%r10
  push   %r14
  mov    %rsi,%r11
  mov    %rcx,%rdi
  push   %r13
  push   %r12
  mov    %rdx,%r12
  mov    %r8,%rdx
  push   %rbp
  mov    %r12,%r8
  xor    %ebp,%ebp
  push   %rbx
  mov    0x60(%rsp),%rbx
  mov    0x58(%rsp),%r9
  nopw   0x0(%rax,%rax,1)
L36:
  lea    0x0(,%rbp,4),%rsi
  xor    %ecx,%ecx
  nopw   0x0(%rax,%rax,1)
L14:
  mov    (%r10,%rsi,1),%r12d
  mov    %r12d,(%rdx,%rcx,1)
  mov    (%r11,%rsi,1),%r12d
  add    $0x100,%rsi
  mov    %r12d,(%rax,%rcx,1)
  add    $0x4,%rcx
  cmp    $0x20,%rcx
  jne    L14
  mov    (%rdx),%esi
  mov    0x10(%rdx),%r12d
  pxor   %xmm1,%xmm1
  mov    (%rax),%ecx
  add    %esi,%r12d
  mov    %r12d,(%rdx)
  mov    0x10(%rax),%r12d
  add    %ecx,%r12d
  mov    %r12d,(%rax)
  sub    0x10(%rdx),%esi
  mov    %esi,0x10(%rdx)
  sub    0x10(%rax),%ecx
  mov    %ecx,0x10(%rax)
  mov    0x4(%rdx),%esi
  mov    0x14(%rdx),%r12d
  mov    0x4(%rax),%ecx
  add    %esi,%r12d
  mov    %r12d,0x4(%rdx)
  mov    0x14(%rax),%r12d
  add    %ecx,%r12d
  mov    %r12d,0x4(%rax)
  sub    0x14(%rdx),%esi
  mov    %esi,0x14(%rdx)
  sub    0x14(%rax),%ecx
  mov    %ecx,0x14(%rax)
  mov    0x8(%rdx),%esi
  mov    0x18(%rdx),%r12d
  mov    0x8(%rax),%ecx
  add    %esi,%r12d
  mov    %r12d,0x8(%rdx)
  mov    0x18(%rax),%r12d
  add    %ecx,%r12d
  mov    %r12d,0x8(%rax)
  sub    0x18(%rdx),%esi
  mov    %esi,0x18(%rdx)
  sub    0x18(%rax),%ecx
  mov    %ecx,0x18(%rax)
  mov    0xc(%rax),%ecx
  mov    0xc(%rdx),%esi
  mov    0x1c(%rdx),%r12d
  add    %esi,%r12d
  mov    %r12d,0xc(%rdx)
  mov    0x1c(%rax),%r12d
  add    %ecx,%r12d
  mov    %r12d,0xc(%rax)
  sub    0x1c(%rdx),%esi
  mov    %esi,0x1c(%rdx)
  sub    0x1c(%rax),%ecx
  mov    %ecx,0x1c(%rax)
  mov    0x14(%rax),%ecx
  mov    0x14(%rdx),%esi
  add    %esi,%ecx
  cvtsi2ss %ecx,%xmm1
  mulss  %xmm0,%xmm1
  cvttss2si %xmm1,%ecx
  pxor   %xmm1,%xmm1
  mov    %ecx,0x14(%rdx)
  mov    0x14(%rax),%ecx
  sub    %esi,%ecx
  mov    0x18(%rax),%esi
  cvtsi2ss %ecx,%xmm1
  mulss  %xmm0,%xmm1
  cvttss2si %xmm1,%ecx
  pxor   %xmm1,%xmm1
  mov    %ecx,0x14(%rax)
  mov    0x18(%rdx),%ecx
  mov    %esi,0x18(%rdx)
  mov    0x1c(%rax),%esi
  neg    %ecx
  mov    %ecx,0x18(%rax)
  mov    0x1c(%rdx),%ecx
  sub    %ecx,%esi
  neg    %ecx
  cvtsi2ss %esi,%xmm1
  mulss  %xmm0,%xmm1
  cvttss2si %xmm1,%esi
  pxor   %xmm1,%xmm1
  mov    %esi,0x1c(%rdx)
  sub    0x1c(%rax),%ecx
  cvtsi2ss %ecx,%xmm1
  mulss  %xmm0,%xmm1
  cvttss2si %xmm1,%ecx
  mov    %ecx,0x1c(%rax)
  mov    (%rdx),%esi
  mov    0x8(%rdx),%r12d
  mov    (%rax),%ecx
  add    %esi,%r12d
  mov    %r12d,(%rdx)
  mov    0x8(%rax),%r12d
  add    %ecx,%r12d
  mov    %r12d,(%rax)
  sub    0x8(%rdx),%esi
  mov    %esi,0x8(%rdx)
  sub    0x8(%rax),%ecx
  mov    %ecx,0x8(%rax)
  mov    0x4(%rdx),%esi
  mov    0xc(%rdx),%r12d
  mov    0x4(%rax),%ecx
  add    %esi,%r12d
  mov    %r12d,0x4(%rdx)
  mov    0xc(%rax),%r12d
  add    %ecx,%r12d
  mov    %r12d,0x4(%rax)
  sub    0xc(%rdx),%esi
  mov    %esi,0xc(%rdx)
  sub    0xc(%rax),%ecx
  mov    %ecx,0xc(%rax)
  mov    0xc(%rdx),%esi
  mov    %ecx,0xc(%rdx)
  mov    (%rax),%ecx
  neg    %esi
  mov    %esi,0xc(%rax)
  mov    (%rdx),%esi
  mov    0x4(%rdx),%r12d
  add    %esi,%r12d
  mov    %r12d,(%rdx)
  mov    0x4(%rax),%r12d
  add    %ecx,%r12d
  mov    %r12d,(%rax)
  sub    0x4(%rdx),%esi
  mov    %esi,0x4(%rdx)
  sub    0x4(%rax),%ecx
  mov    %ecx,0x4(%rax)
  mov    0x8(%rdx),%esi
  mov    0xc(%rdx),%r12d
  mov    0x8(%rax),%ecx
  add    %esi,%r12d
  mov    %r12d,0x8(%rdx)
  mov    0xc(%rax),%r12d
  add    %ecx,%r12d
  mov    %r12d,0x8(%rax)
  sub    0xc(%rdx),%esi
  mov    %esi,0xc(%rdx)
  sub    0xc(%rax),%ecx
  mov    %ecx,0xc(%rax)
  mov    0x10(%rdx),%esi
  mov    0x18(%rdx),%r12d
  mov    0x10(%rax),%ecx
  add    %esi,%r12d
  mov    %r12d,0x10(%rdx)
  mov    0x18(%rax),%r12d
  add    %ecx,%r12d
  mov    %r12d,0x10(%rax)
  sub    0x18(%rdx),%esi
  mov    %esi,0x18(%rdx)
  sub    0x18(%rax),%ecx
  mov    %ecx,0x18(%rax)
  mov    0x14(%rdx),%esi
  mov    0x1c(%rdx),%r12d
  mov    0x14(%rax),%ecx
  add    %esi,%r12d
  mov    %r12d,0x14(%rdx)
  mov    0x1c(%rax),%r12d
  add    %ecx,%r12d
  mov    %r12d,0x14(%rax)
  sub    0x1c(%rdx),%esi
  mov    %esi,0x1c(%rdx)
  sub    0x1c(%rax),%ecx
  mov    %ecx,0x1c(%rax)
  mov    0x1c(%rdx),%esi
  mov    %ecx,0x1c(%rdx)
  mov    0x10(%rax),%ecx
  neg    %esi
  mov    %esi,0x1c(%rax)
  mov    0x10(%rdx),%esi
  mov    0x14(%rdx),%r12d
  add    %esi,%r12d
  mov    %r12d,0x10(%rdx)
  mov    0x14(%rax),%r12d
  add    %ecx,%r12d
  mov    %r12d,0x10(%rax)
  sub    0x14(%rdx),%esi
  mov    %esi,0x14(%rdx)
  sub    0x14(%rax),%ecx
  mov    %ecx,0x14(%rax)
  mov    0x18(%rdx),%esi
  mov    0x1c(%rdx),%r12d
  mov    0x18(%rax),%ecx
  add    %esi,%r12d
  mov    %r12d,0x18(%rdx)
  mov    0x1c(%rax),%r12d
  add    %ecx,%r12d
  mov    %r12d,0x18(%rax)
  sub    0x1c(%rdx),%esi
  mov    %esi,0x1c(%rdx)
  sub    0x1c(%rax),%ecx
  mov    %ecx,0x1c(%rax)
  mov    $0x4,%ecx
  nopl   0x0(%rax)
L5:
  cvttss2si -0x4(%rbx,%rcx,1),%esi
  mov    (%rdx,%rcx,1),%r14d
  cvttss2si -0x4(%r9,%rcx,1),%r12d
  mov    (%rax,%rcx,1),%r15d
  mov    %esi,%r13d
  imul   %r12d,%r15d
  imul   %r14d,%r13d
  imul   %r14d,%r12d
  sub    %r15d,%r13d
  mov    %r13d,(%rdx,%rcx,1)
  imul   (%rax,%rcx,1),%esi
  add    %r12d,%esi
  mov    %esi,(%rax,%rcx,1)
  add    $0x4,%rcx
  cmp    $0x20,%rcx
  jne    L5
  xor    %ecx,%ecx
  nopl   0x0(%rax)
L0:
  mov    (%rdx,%rcx,1),%esi
  mov    %esi,(%r8,%rcx,1)
  mov    (%rax,%rcx,1),%esi
  mov    %esi,(%rdi,%rcx,1)
  add    $0x4,%rcx
  cmp    $0x20,%rcx
  jne    L0
  add    $0x1,%rbp
  add    $0x20,%r8
  add    $0x20,%rdi
  add    $0x1c,%rbx
  add    $0x1c,%r9
  cmp    $0x40,%rbp
  jne    L36
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: extract_symbolic_data
arguments: label: D1
data_type: byte[4]

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

assistant: ```c
void step1(int work_x[], int work_y[], int DATA_x[], int DATA_y[], int data_x[], int data_y[], int smem[], int reversed[], float sin_64[448], float cos_64[448], float sin_512[448], float cos_512[448]) {
  int tid, hi, lo, i, j, stride;
  stride = 64;
  for (tid = 0; tid < 64; tid++) {
    for (i = 0; i < 8; i++) {
      data_x[i] = work_x[i * stride + tid];
      data_y[i] = work_y[i * stride + tid];
    }
    {
      int exp_1_8_x, exp_1_4_x, exp_3_8_x;
      int exp_1_8_y, exp_1_4_y, exp_3_8_y;
      int tmp_1, tmp_2;
      exp_1_8_x = 1;
      exp_1_8_y = -1;
      exp_1_4_x = 0;
      exp_1_4_y = -1;
      exp_3_8_x = -1;
      exp_3_8_y = -1;
      {
        int c0_x = *&data_x[0];
        int c0_y = *&data_y[0];
        *&data_x[0] = (c0_x + *&data_x[4]);
        *&data_y[0] = (c0_y + *&data_y[4]);
        *&data_x[4] = (c0_x - *&data_x[4]);
        *&data_y[4] = (c0_y - *&data_y[4]);
      };
      {
        int c0_x = *&data_x[1];
        int c0_y = *&data_y[1];
        *&data_x[1] = (c0_x + *&data_x[5]);
        *&data_y[1] = (c0_y + *&data_y[5]);
        *&data_x[5] = (c0_x - *&data_x[5]);
        *&data_y[5] = (c0_y - *&data_y[5]);
      };
      {
        int c0_x = *&data_x[2];
        int c0_y = *&data_y[2];
        *&data_x[2] = (c0_x + *&data_x[6]);
        *&data_y[2] = (c0_y + *&data_y[6]);
        *&data_x[6] = (c0_x - *&data_x[6]);
        *&data_y[6] = (c0_y - *&data_y[6]);
      };
      {
        int c0_x = *&data_x[3];
        int c0_y = *&data_y[3];
        *&data_x[3] = (c0_x + *&data_x[7]);
        *&data_y[3] = (c0_y + *&data_y[7]);
        *&data_x[7] = (c0_x - *&data_x[7]);
        *&data_y[7] = (c0_y - *&data_y[7]);
      };
      tmp_1 = data_x[5];
      data_x[5] = (0.70710678118654752440f * (data_x[5] * exp_1_8_x - data_y[5] * exp_1_8_y));
      data_y[5] = (0.70710678118654752440f * (tmp_1 * exp_1_8_y + data_y[5] * exp_1_8_x));
      tmp_1 = data_x[6];
      data_x[6] = (data_x[6] * exp_1_4_x - data_y[6] * exp_1_4_y);
      data_y[6] = (tmp_1 * exp_1_4_y + data_y[6] * exp_1_4_x);
      tmp_1 = data_x[7];
      data_x[7] = (0.70710678118654752440f * (data_x[7] * exp_3_8_x - data_y[7] * exp_3_8_y));
      data_y[7] = (0.70710678118654752440f * (tmp_1 * exp_3_8_y + data_y[7] * exp_3_8_x));
      {
        int exp_1_44_x;
        int exp_1_44_y;
        int tmp;
        exp_1_44_x = 0.0;
        exp_1_44_y = -1.0;
        {
          int c0_x = *&data_x[0];
          int c0_y = *&data_y[0];
          *&data_x[0] = (c0_x + *&data_x[2]);
          *&data_y[0] = (c0_y + *&data_y[2]);
          *&data_x[2] = (c0_x - *&data_x[2]);
          *&data_y[2] = (c0_y - *&data_y[2]);
        };
        {
          int c0_x = *&data_x[1];
          int c0_y = *&data_y[1];
          *&data_x[1] = (c0_x + *&data_x[3]);
          *&data_y[1] = (c0_y + *&data_y[3]);
          *&data_x[3] = (c0_x - *&data_x[3]);
          *&data_y[3] = (c0_y - *&data_y[3]);
        };
        tmp = *&data_x[3];
        *&data_x[3] = *&data_x[3] * exp_1_44_x - *&data_y[3] * exp_1_44_y;
        *&data_y[3] = tmp * exp_1_44_y - *&data_y[3] * exp_1_44_x;
        {
          int c0_x = *&data_x[0];
          int c0_y = *&data_y[0];
          *&data_x[0] = (c0_x + *&data_x[1]);
          *&data_y[0] = (c0_y + *&data_y[1]);
          *&data_x[1] = (c0_x - *&data_x[1]);
          *&data_y[1] = (c0_y - *&data_y[1]);
        };
        {
          int c0_x = *&data_x[2];
          int c0_y = *&data_y[2];
          *&data_x[2] = (c0_x + *&data_x[3]);
          *&data_y[2] = (c0_y + *&data_y[3]);
          *&data_x[3] = (c0_x - *&data_x[3]);
          *&data_y[3] = (c0_y - *&data_y[3]);
        };
      };
      {
        int exp_1_44_x;
        int exp_1_44_y;
        int tmp;
        exp_1_44_x = 0.0;
        exp_1_44_y = -1.0;
        {
          int c0_x = *&data_x[4];
          int c0_y = *&data_y[4];
          *&data_x[4] = (c0_x + *&data_x[6]);
          *&data_y[4] = (c0_y + *&data_y[6]);
          *&data_x[6] = (c0_x - *&data_x[6]);
          *&data_y[6] = (c0_y - *&data_y[6]);
        };
        {
          int c0_x = *&data_x[5];
          int c0_y = *&data_y[5];
          *&data_x[5] = (c0_x + *&data_x[7]);
          *&data_y[5] = (c0_y + *&data_y[7]);
          *&data_x[7] = (c0_x - *&data_x[7]);
          *&data_y[7] = (c0_y - *&data_y[7]);
        };
        tmp = *&data_x[7];
        *&data_x[7] = *&data_x[7] * exp_1_44_x - *&data_y[7] * exp_1_44_y;
        *&data_y[7] = tmp * exp_1_44_y - *&data_y[7] * exp_1_44_x;
        {
          int c0_x = *&data_x[4];
          int c0_y = *&data_y[4];
          *&data_x[4] = (c0_x + *&data_x[5]);
          *&data_y[4] = (c0_y + *&data_y[5]);
          *&data_x[5] = (c0_x - *&data_x[5]);
          *&data_y[5] = (c0_y - *&data_y[5]);
        };
        {
          int c0_x = *&data_x[6];
          int c0_y = *&data_y[6];
          *&data_x[6] = (c0_x + *&data_x[7]);
          *&data_y[6] = (c0_y + *&data_y[7]);
          *&data_x[7] = (c0_x - *&data_x[7]);
          *&data_y[7] = (c0_y - *&data_y[7]);
        };
      };
    };
    for (j = 1; j < 8; j++) {
      int A_x, A_y, tmp;
      A_x = cos_512[tid * 7 + j - 1];
      A_y = sin_512[tid * 7 + j - 1];
      tmp = data_x[j];
      data_x[j] = (data_x[j] * A_x - data_y[j] * A_y);
      data_y[j] = (tmp * A_y + data_y[j] * A_x);
    }
    for (i = 0; i < 8; i++) {
      DATA_x[tid * 8 + i] = data_x[i];
      DATA_y[tid * 8 + i] = data_y[i];
    }
  }
}
```