user: 请从给定的汇编指令推导出C源代码:

```asm
<FFT1>:
  endbr64
  push   %r15
  mov    %rsi,%r15
  push   %r14
  mov    %rdi,%r14
  push   %r13
  push   %r12
  push   %rbp
  mov    %edx,%ebp
  push   %rbx
  sub    $0x88,%rsp
  mov    %edx,0x50(%rsp)
  mov    %ecx,0x6c(%rsp)
  cmp    $0xffffffff,%ecx
  je     L137
L237:
  mov    %ebp,%edi
  call   <Calc_power_of_two@plt>
  mov    %ebp,%edx
  mov    %r15,%rsi
  mov    %r14,%rdi
  mov    %eax,%ecx
  mov    %eax,0x54(%rsp)
  mov    %eax,%ebx
  call   <Make_initial_data@plt>
  pxor   %xmm0,%xmm0
  movsd  D77(%rip),%xmm1
  cvtsi2sd %ebp,%xmm0
  divsd  %xmm0,%xmm1
  movsd  %xmm1,0x48(%rsp)
  test   %ebx,%ebx
  jle    L78
L47:
  mov    0x50(%rsp),%eax
  movl   $0x2,0x30(%rsp)
  movl   $0x0,0x44(%rsp)
  sub    $0x1,%eax
  lea    0x8(,%rax,8),%rdi
  mov    %rax,0x60(%rsp)
  lea    (%r15,%rdi,1),%rax
  mov    %rdi,0x58(%rsp)
  cmp    %rax,%r14
  lea    (%r14,%rdi,1),%rax
  setae  %dl
  cmp    %rax,%r15
  setae  %al
  or     %eax,%edx
  mov    %dl,0x6b(%rsp)
  nop
L173:
  mov    0x50(%rsp),%eax
  mov    0x30(%rsp),%ecx
  pxor   %xmm0,%xmm0
  cltd
  idiv   %ecx
  sar    %ecx
  cvtsi2sd %eax,%xmm0
  mulsd  0x48(%rsp),%xmm0
  mov    %eax,0x40(%rsp)
  movsd  %xmm0,0x20(%rsp)
  test   %eax,%eax
  jle    L77
  movslq 0x30(%rsp),%rax
  movl   $0x0,0x2c(%rsp)
  lea    0x78(%rsp),%r13
  mov    %rax,%rdx
  shl    $0x3,%rax
  mov    %rax,0x38(%rsp)
  lea    D45(%rip),%rax
  mov    %rax,0x10(%rsp)
  lea    D114(%rip),%rax
  mov    %edx,0x28(%rsp)
  mov    %rax,0x18(%rsp)
  mov    %r15,%rax
  mov    %r14,%r15
  mov    %ecx,%r14d
  mov    %rax,%rbx
  nopl   0x0(%rax,%rax,1)
L269:
  mov    0x28(%rsp),%eax
  sub    0x30(%rsp),%eax
  xor    %ebp,%ebp
  mov    %r15,%rdx
  mov    %eax,%r12d
  mov    %eax,0x34(%rsp)
  pxor   %xmm2,%xmm2
  mov    %rbp,%r15
  mov    D107(%rip),%rax
  movapd %xmm2,%xmm5
  mov    %rdx,%rbp
  movq   %rax,%xmm0
  mov    %rbx,%rax
  mov    %r14d,%ebx
  mov    %r12d,%r14d
  mov    %rax,%r12
  jmp    L271
  nopl   0x0(%rax)
L112:
  movapd %xmm5,%xmm0
  lea    0x70(%rsp),%rsi
  mov    %r13,%rdi
  add    $0x1,%r15
  movsd  %xmm5,0x8(%rsp)
  call   <sincos@plt>
  movsd  0x70(%rsp),%xmm0
  movsd  0x78(%rsp),%xmm2
  movsd  0x8(%rsp),%xmm5
L271:
  mov    %r14d,%eax
  lea    (%rbx,%r14,1),%esi
  cmp    %r15d,%ebx
  jg     L227
  mov    %r14d,%esi
  sub    %ebx,%eax
L227:
  movslq %esi,%rsi
  movapd %xmm2,%xmm6
  movapd %xmm0,%xmm4
  cltq
  movsd  (%r12,%rsi,8),%xmm1
  movsd  0x0(%rbp,%rsi,8),%xmm3
  add    $0x1,%r14d
  mov    0x10(%rsp),%rcx
  addsd  0x20(%rsp),%xmm5
  mulsd  %xmm1,%xmm6
  mulsd  %xmm3,%xmm4
  addsd  0x0(%rbp,%rax,8),%xmm4
  mulsd  %xmm3,%xmm2
  mulsd  %xmm1,%xmm0
  addsd  %xmm6,%xmm4
  movsd  %xmm4,(%rcx,%r15,8)
  addsd  (%r12,%rax,8),%xmm0
  mov    0x18(%rsp),%rax
  subsd  %xmm2,%xmm0
  movsd  %xmm0,(%rax,%r15,8)
  cmp    0x28(%rsp),%r14d
  jne    L112
  mov    0x38(%rsp),%rdi
  mov    0x34(%rsp),%edx
  mov    %ebx,%r14d
  mov    %rbp,%r15
  add    %rdi,0x10(%rsp)
  mov    %r12,%rbx
  add    %rdi,0x18(%rsp)
  mov    0x30(%rsp),%edi
  addl   $0x1,0x2c(%rsp)
  mov    0x2c(%rsp),%eax
  lea    (%rdx,%rdi,2),%edx
  mov    %edx,0x28(%rsp)
  cmp    %eax,0x40(%rsp)
  jne    L269
  mov    %rbp,%r14
  mov    %r12,%r15
L77:
  mov    0x50(%rsp),%edx
  test   %edx,%edx
  jle    L44
  cmpb   $0x0,0x6b(%rsp)
  je     L209
  mov    0x58(%rsp),%rbx
  lea    D45(%rip),%r9
  mov    %r14,%rdi
  mov    %r9,%rsi
  mov    %rbx,%rdx
  call   <memcpy@plt>
  lea    D114(%rip),%rsi
  mov    %rbx,%rdx
  mov    %r15,%rdi
  call   <memcpy@plt>
L44:
  addl   $0x1,0x44(%rsp)
  mov    0x44(%rsp),%eax
  shll   0x30(%rsp)
  cmp    0x54(%rsp),%eax
  jl     L173
  cmpl   $0xffffffff,0x6c(%rsp)
  je     L143
L78:
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L209:
  mov    0x60(%rsp),%rcx
  xor    %eax,%eax
  lea    D45(%rip),%r9
  lea    D114(%rip),%rsi
  nopw   0x0(%rax,%rax,1)
L232:
  movsd  (%r9,%rax,8),%xmm0
  mov    %rax,%rdx
  movsd  %xmm0,(%r14,%rax,8)
  movsd  (%rsi,%rax,8),%xmm0
  movsd  %xmm0,(%r15,%rax,8)
  add    $0x1,%rax
  cmp    %rdx,%rcx
  jne    L232
  jmp    L44
L137:
  test   %edx,%edx
  jle    L14
  lea    0xf(%rdi),%rax
  pxor   %xmm1,%xmm1
  sub    %rsi,%rax
  cvtsi2sd %edx,%xmm1
  cmp    $0x1e,%rax
  jbe    L60
  mov    %edx,%eax
  cmp    $0x1,%edx
  je     L60
  shr    %eax
  movapd %xmm1,%xmm3
  movapd D6(%rip),%xmm2
  mov    %eax,%edx
  unpcklpd %xmm3,%xmm3
  xor    %eax,%eax
  shl    $0x4,%rdx
L231:
  movupd (%r14,%rax,1),%xmm0
  divpd  %xmm3,%xmm0
  movups %xmm0,(%r14,%rax,1)
  movupd (%r15,%rax,1),%xmm0
  xorpd  %xmm2,%xmm0
  divpd  %xmm3,%xmm0
  movups %xmm0,(%r15,%rax,1)
  add    $0x10,%rax
  cmp    %rdx,%rax
  jne    L231
  mov    0x50(%rsp),%edx
  mov    %edx,%eax
  and    $0xfffffffe,%eax
  and    $0x1,%edx
  je     L279
  shl    $0x3,%rax
  lea    (%r14,%rax,1),%rdx
  add    %r15,%rax
  movsd  (%rdx),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdx)
  movsd  (%rax),%xmm0
  xorpd  D25(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
L279:
  mov    0x50(%rsp),%ebp
  movsd  %xmm1,0x8(%rsp)
  mov    %ebp,%edi
  call   <Calc_power_of_two@plt>
  mov    %ebp,%edx
  mov    %r15,%rsi
  mov    %r14,%rdi
  mov    %eax,%ecx
  mov    %eax,0x54(%rsp)
  mov    %eax,%ebx
  call   <Make_initial_data@plt>
  movsd  D77(%rip),%xmm0
  movsd  0x8(%rsp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x48(%rsp)
  test   %ebx,%ebx
  jg     L47
L253:
  mov    0x50(%rsp),%edx
  cmp    $0x1,%edx
  je     L252
  shr    %edx
  movapd D6(%rip),%xmm2
  mov    %r15,%rax
  sub    $0x1,%edx
  shl    $0x4,%rdx
  lea    0x10(%r15,%rdx,1),%rdx
L260:
  movupd (%rax),%xmm0
  add    $0x10,%rax
  xorpd  %xmm2,%xmm0
  movups %xmm0,-0x10(%rax)
  cmp    %rdx,%rax
  jne    L260
  mov    0x50(%rsp),%edx
  mov    %edx,%eax
  and    $0xfffffffe,%eax
  and    $0x1,%edx
  je     L78
L239:
  cltq
  lea    (%r15,%rax,8),%rax
  movsd  (%rax),%xmm0
  xorpd  D25(%rip),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L78
L143:
  mov    0x50(%rsp),%eax
  test   %eax,%eax
  jg     L253
  jmp    L78
L60:
  movslq 0x50(%rsp),%rdx
  movq   D25(%rip),%xmm2
  xor    %eax,%eax
L178:
  movsd  (%r14,%rax,8),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%r14,%rax,8)
  movsd  (%r15,%rax,8),%xmm0
  xorpd  %xmm2,%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%r15,%rax,8)
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L178
  jmp    L279
L14:
  mov    0x50(%rsp),%ebp
  jmp    L237
L252:
  xor    %eax,%eax
  jmp    L239
```

tool_calls: name: disassemble_data
arguments: label: D77
type_guessed: byte[8]
name: disassemble_data
arguments: label: D107
type_guessed: f32[2]
name: disassemble_data
arguments: label: D6
type_guessed: f32[4]
name: disassemble_data
arguments: label: D25
type_guessed: f32[4]

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

tool: {"D107": [0.0,1.875]}

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

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

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);
    }
  }
}
```