user: What is the c source code of the assembly code below:

```asm
<FFT1>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %rdi,-0x58(%rbp)
  mov    %rsi,-0x60(%rbp)
  mov    %edx,-0x64(%rbp)
  mov    %ecx,-0x68(%rbp)
  cmpl   $0xffffffff,-0x68(%rbp)
  jne    L84
  movl   $0x0,-0x50(%rbp)
  jmp    L30
L103:
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x64(%rbp),%xmm1
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movq   D18(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x64(%rbp),%xmm1
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x50(%rbp)
L30:
  mov    -0x50(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L103
L84:
  mov    -0x64(%rbp),%eax
  mov    %eax,%edi
  call   <Calc_power_of_two@plt>
  mov    %eax,-0x38(%rbp)
  mov    -0x38(%rbp),%ecx
  mov    -0x64(%rbp),%edx
  mov    -0x60(%rbp),%rsi
  mov    -0x58(%rbp),%rax
  mov    %rax,%rdi
  call   <Make_initial_data@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x64(%rbp),%xmm1
  movsd  D53(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movl   $0x2,-0x44(%rbp)
  movl   $0x0,-0x50(%rbp)
  jmp    L18
L15:
  mov    -0x64(%rbp),%eax
  cltd
  idivl  -0x44(%rbp)
  mov    %eax,-0x34(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x34(%rbp),%xmm0
  movsd  -0x20(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x44(%rbp),%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  mov    %eax,-0x30(%rbp)
  movl   $0x0,-0x4c(%rbp)
  jmp    L16
L154:
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movl   $0x0,-0x48(%rbp)
  jmp    L52
L51:
  mov    -0x4c(%rbp),%eax
  imul   -0x44(%rbp),%eax
  mov    %eax,%edx
  mov    -0x48(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x2c(%rbp)
  mov    -0x48(%rbp),%eax
  cmp    -0x30(%rbp),%eax
  jge    L108
  mov    -0x2c(%rbp),%eax
  mov    %eax,-0x40(%rbp)
  mov    -0x40(%rbp),%edx
  mov    -0x30(%rbp),%eax
  add    %edx,%eax
  mov    %eax,-0x3c(%rbp)
  jmp    L98
L108:
  mov    -0x2c(%rbp),%eax
  mov    %eax,-0x3c(%rbp)
  mov    -0x2c(%rbp),%eax
  sub    -0x30(%rbp),%eax
  mov    %eax,-0x40(%rbp)
L98:
  mov    -0x3c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x3c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0x40(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm4
  movsd  %xmm4,-0x70(%rbp)
  mov    -0x28(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  mulsd  -0x10(%rbp),%xmm0
  addsd  -0x70(%rbp),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  mov    -0x28(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  movq   %rax,%xmm3
  mulsd  -0x8(%rbp),%xmm3
  movapd %xmm3,%xmm0
  addsd  -0x70(%rbp),%xmm0
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  lea    D57(%rip),%rax
  movsd  %xmm0,(%rdx,%rax,1)
  mov    -0x40(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm5
  movsd  %xmm5,-0x70(%rbp)
  mov    -0x28(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  mulsd  -0x8(%rbp),%xmm0
  addsd  -0x70(%rbp),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  mov    -0x28(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  movq   %rax,%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movsd  -0x70(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  lea    D8(%rip),%rax
  movsd  %xmm0,(%rdx,%rax,1)
  movsd  -0x28(%rbp),%xmm0
  addsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  addl   $0x1,-0x48(%rbp)
L52:
  mov    -0x48(%rbp),%eax
  cmp    -0x44(%rbp),%eax
  jl     L51
  addl   $0x1,-0x4c(%rbp)
L16:
  mov    -0x4c(%rbp),%eax
  cmp    -0x34(%rbp),%eax
  jl     L154
  movl   $0x0,-0x4c(%rbp)
  jmp    L3
L120:
  mov    -0x4c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x58(%rbp),%rax
  add    %rdx,%rax
  mov    -0x4c(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rcx
  lea    D57(%rip),%rdx
  movsd  (%rcx,%rdx,1),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x4c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  mov    -0x4c(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rcx
  lea    D8(%rip),%rdx
  movsd  (%rcx,%rdx,1),%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x4c(%rbp)
L3:
  mov    -0x4c(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L120
  shll   -0x44(%rbp)
  addl   $0x1,-0x50(%rbp)
L18:
  mov    -0x50(%rbp),%eax
  cmp    -0x38(%rbp),%eax
  jl     L15
  cmpl   $0xffffffff,-0x68(%rbp)
  jne    L143
  movl   $0x0,-0x4c(%rbp)
  jmp    L85
L139:
  mov    -0x4c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x4c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x60(%rbp),%rax
  add    %rdx,%rax
  movq   D18(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addl   $0x1,-0x4c(%rbp)
L85:
  mov    -0x4c(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L139
L143:
  nop
  leave
  ret
```

tool_calls: name: parse_data
arguments: data_tag: D18
type_infered: f32[2]
name: parse_data
arguments: data_tag: D53
type_infered: byte[8]

tool: {"D18": [0.0,-0.0]}

tool: {"D53": [0x1e,0x26,0x44,0x54,0xfb,0x21,0x19,0x40]}

assistant: ```c
void FFT1(double *re_part, double *im_part, int number_of_data, int flag) {
  int i, j, k;
  int power;
  double unit_angle, step_angle;
  int dft, half, num_of_dft;
  int num_out, num_in1, num_in2;
  double re_buf, im_buf, angle;
  static double re_part_new[256], im_part_new[256];
  if (flag == -1) {
    for (i = 0; i < number_of_data; i++) {
      *(re_part + i) = *(re_part + i) / number_of_data;
      *(im_part + i) = -*(im_part + i) / number_of_data;
    }
  }
  power = Calc_power_of_two(number_of_data);
  Make_initial_data(re_part, im_part, number_of_data, power);
  unit_angle = 2.0 * 3.141592653589 / number_of_data;
  dft = 2;
  for (i = 0; i < power; i++) {
    num_of_dft = number_of_data / dft;
    step_angle = unit_angle * num_of_dft;
    half = dft / 2;
    for (j = 0; j < num_of_dft; j++) {
      angle = 0.0;
      for (k = 0; k < dft; k++) {
        num_out = j * dft + k;
        if (k < half) {
          num_in1 = num_out;
          num_in2 = num_in1 + half;
        } else {
          num_in2 = num_out;
          num_in1 = num_out - half;
        }
        re_buf = *(re_part + num_in2);
        im_buf = *(im_part + num_in2);
        re_part_new[num_out] = *(re_part + num_in1) + re_buf * cos(angle) + im_buf * sin(angle);
        im_part_new[num_out] = *(im_part + num_in1) + im_buf * cos(angle) - re_buf * sin(angle);
        angle = angle + step_angle;
      }
    }
    for (j = 0; j < number_of_data; j++) {
      *(re_part + j) = re_part_new[j];
      *(im_part + j) = im_part_new[j];
    }
    dft = dft * 2;
  }
  if (flag == -1) {
    for (j = 0; j < number_of_data; j++) {
      *(im_part + j) = -*(im_part + j);
    }
  }
}
```