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

```asm
<FFT_7>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xc8,%rsp
  mov    %edx,%ebx
  mov    %fs:0x28,%rdx
  mov    %rdx,0xb8(%rsp)
  xor    %edx,%edx
  test   %ebx,%ebx
  jle    L67
  mov    %esi,%eax
  mov    %ecx,%esi
  mov    %edi,%r15d
  imul   %eax,%r15d
  lea    0x20(%rsp),%r13
  lea    -0x1(%rax),%edx
  lea    0x24(%rsp,%rdx,4),%rcx
  mov    $0x0,%r10d
  movsd  D2(%rip),%xmm15
  movsd  D138(%rip),%xmm14
  movsd  D78(%rip),%xmm13
  jmp    L5
L38:
  movsd  0x48(%rsp),%xmm1
  movsd  0x70(%rsp),%xmm0
  movapd %xmm1,%xmm12
  addsd  %xmm0,%xmm12
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rsp)
  movsd  0x88(%rsp),%xmm7
  movsd  0xb0(%rsp),%xmm4
  movapd %xmm7,%xmm1
  subsd  %xmm4,%xmm1
  addsd  %xmm4,%xmm7
  movsd  0x50(%rsp),%xmm4
  movsd  0x68(%rsp),%xmm5
  movapd %xmm4,%xmm11
  addsd  %xmm5,%xmm11
  subsd  %xmm4,%xmm5
  movsd  %xmm5,0x8(%rsp)
  movsd  0x90(%rsp),%xmm6
  movsd  0xa8(%rsp),%xmm4
  movapd %xmm6,%xmm9
  subsd  %xmm4,%xmm9
  addsd  %xmm4,%xmm6
  movsd  0x58(%rsp),%xmm4
  movsd  0x60(%rsp),%xmm5
  movapd %xmm4,%xmm10
  addsd  %xmm5,%xmm10
  subsd  %xmm4,%xmm5
  movsd  %xmm5,0x10(%rsp)
  movsd  0x98(%rsp),%xmm5
  movsd  0xa0(%rsp),%xmm4
  movapd %xmm5,%xmm8
  subsd  %xmm4,%xmm8
  addsd  %xmm4,%xmm5
  movsd  D165(%rip),%xmm4
  mulsd  %xmm1,%xmm4
  movsd  D74(%rip),%xmm0
  movsd  %xmm8,0x18(%rsp)
  mulsd  %xmm8,%xmm0
  subsd  %xmm0,%xmm4
  movsd  D115(%rip),%xmm0
  mulsd  %xmm9,%xmm0
  subsd  %xmm0,%xmm4
  movapd %xmm10,%xmm0
  mulsd  %xmm15,%xmm0
  movapd %xmm0,%xmm8
  movapd %xmm11,%xmm0
  mulsd  %xmm14,%xmm0
  subsd  %xmm0,%xmm8
  movapd %xmm12,%xmm0
  mulsd  %xmm13,%xmm0
  subsd  %xmm0,%xmm8
  addsd  %xmm3,%xmm8
  movslq 0x34(%rsp),%rbp
  movapd %xmm8,%xmm0
  subsd  %xmm4,%xmm0
  movsd  %xmm0,(%r8,%rbp,8)
  movslq 0x28(%rsp),%r12
  addsd  %xmm8,%xmm4
  movsd  %xmm4,(%r8,%r12,8)
  movsd  D165(%rip),%xmm0
  mulsd  (%rsp),%xmm0
  movapd %xmm0,%xmm8
  movsd  D74(%rip),%xmm0
  mulsd  0x10(%rsp),%xmm0
  subsd  %xmm0,%xmm8
  movsd  D115(%rip),%xmm0
  mulsd  0x8(%rsp),%xmm0
  subsd  %xmm0,%xmm8
  movapd %xmm5,%xmm0
  mulsd  %xmm15,%xmm0
  movapd %xmm0,%xmm4
  movapd %xmm6,%xmm0
  mulsd  %xmm14,%xmm0
  subsd  %xmm0,%xmm4
  movapd %xmm7,%xmm0
  mulsd  %xmm13,%xmm0
  subsd  %xmm0,%xmm4
  addsd  %xmm2,%xmm4
  movapd %xmm8,%xmm0
  addsd  %xmm4,%xmm0
  movsd  %xmm0,(%r9,%r12,8)
  subsd  %xmm8,%xmm4
  movsd  %xmm4,(%r9,%rbp,8)
  movsd  D74(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm0,%xmm8
  movsd  D165(%rip),%xmm0
  mulsd  %xmm9,%xmm0
  movapd %xmm0,%xmm4
  addsd  %xmm8,%xmm4
  movsd  D115(%rip),%xmm0
  mulsd  0x18(%rsp),%xmm0
  addsd  %xmm0,%xmm4
  movapd %xmm12,%xmm0
  mulsd  %xmm15,%xmm0
  movapd %xmm0,%xmm8
  movapd %xmm11,%xmm0
  mulsd  %xmm13,%xmm0
  subsd  %xmm0,%xmm8
  movapd %xmm10,%xmm0
  mulsd  %xmm14,%xmm0
  subsd  %xmm0,%xmm8
  addsd  %xmm3,%xmm8
  movslq 0x38(%rsp),%rbp
  movapd %xmm8,%xmm0
  subsd  %xmm4,%xmm0
  movsd  %xmm0,(%r8,%rbp,8)
  movslq 0x24(%rsp),%r12
  addsd  %xmm8,%xmm4
  movsd  %xmm4,(%r8,%r12,8)
  movsd  D74(%rip),%xmm0
  mulsd  (%rsp),%xmm0
  movapd %xmm0,%xmm4
  movsd  D165(%rip),%xmm0
  mulsd  0x8(%rsp),%xmm0
  movapd %xmm0,%xmm8
  addsd  %xmm4,%xmm8
  movsd  D115(%rip),%xmm0
  mulsd  0x10(%rsp),%xmm0
  addsd  %xmm0,%xmm8
  movapd %xmm7,%xmm0
  mulsd  %xmm15,%xmm0
  movapd %xmm0,%xmm4
  movapd %xmm6,%xmm0
  mulsd  %xmm13,%xmm0
  subsd  %xmm0,%xmm4
  movapd %xmm5,%xmm0
  mulsd  %xmm14,%xmm0
  subsd  %xmm0,%xmm4
  addsd  %xmm2,%xmm4
  movapd %xmm8,%xmm0
  addsd  %xmm4,%xmm0
  movsd  %xmm0,(%r9,%r12,8)
  subsd  %xmm8,%xmm4
  movsd  %xmm4,(%r9,%rbp,8)
  mulsd  D115(%rip),%xmm1
  movsd  D165(%rip),%xmm0
  mulsd  0x18(%rsp),%xmm0
  addsd  %xmm0,%xmm1
  mulsd  D74(%rip),%xmm9
  subsd  %xmm9,%xmm1
  movapd %xmm11,%xmm0
  mulsd  %xmm15,%xmm0
  movapd %xmm0,%xmm4
  movapd %xmm10,%xmm0
  mulsd  %xmm13,%xmm0
  subsd  %xmm0,%xmm4
  movapd %xmm12,%xmm0
  mulsd  %xmm14,%xmm0
  subsd  %xmm0,%xmm4
  addsd  %xmm3,%xmm4
  movslq 0x30(%rsp),%rbp
  movapd %xmm4,%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%r8,%rbp,8)
  movslq 0x2c(%rsp),%r12
  addsd  %xmm4,%xmm1
  movsd  %xmm1,(%r8,%r12,8)
  movsd  (%rsp),%xmm0
  mulsd  D115(%rip),%xmm0
  movsd  D165(%rip),%xmm1
  mulsd  0x10(%rsp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  D74(%rip),%xmm1
  mulsd  0x8(%rsp),%xmm1
  subsd  %xmm1,%xmm0
  movapd %xmm6,%xmm1
  mulsd  %xmm15,%xmm1
  movapd %xmm5,%xmm4
  mulsd  %xmm13,%xmm4
  subsd  %xmm4,%xmm1
  movapd %xmm7,%xmm4
  mulsd  %xmm14,%xmm4
  subsd  %xmm4,%xmm1
  addsd  %xmm2,%xmm1
  movapd %xmm0,%xmm4
  addsd  %xmm1,%xmm4
  movsd  %xmm4,(%r9,%r12,8)
  subsd  %xmm0,%xmm1
  movsd  %xmm1,(%r9,%rbp,8)
  addsd  %xmm12,%xmm3
  addsd  %xmm11,%xmm3
  addsd  %xmm10,%xmm3
  movsd  %xmm3,(%r11)
  addsd  %xmm7,%xmm2
  addsd  %xmm6,%xmm2
  addsd  %xmm5,%xmm2
  movsd  %xmm2,(%rdx)
  add    %r15d,%esi
  add    $0x1,%r10d
  cmp    %r10d,%ebx
  je     L67
L5:
  mov    %r13,%rdx
  mov    %esi,%r11d
  test   %eax,%eax
  jle    L15
L68:
  mov    %r11d,(%rdx)
  add    %edi,%r11d
  add    $0x4,%rdx
  cmp    %rcx,%rdx
  jne    L68
L15:
  movslq 0x20(%rsp),%rdx
  shl    $0x3,%rdx
  lea    (%r8,%rdx,1),%r11
  movsd  (%r11),%xmm3
  add    %r9,%rdx
  movsd  (%rdx),%xmm2
  cmp    $0x1,%eax
  jle    L38
  mov    $0x1,%ebp
L39:
  movslq 0x20(%rsp,%rbp,4),%r12
  movsd  (%r8,%r12,8),%xmm0
  movsd  %xmm0,0x40(%rsp,%rbp,8)
  movsd  (%r9,%r12,8),%xmm0
  movsd  %xmm0,0x80(%rsp,%rbp,8)
  add    $0x1,%rbp
  cmp    %ebp,%eax
  jg     L39
  jmp    L38
L67:
  mov    0xb8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L57
  add    $0xc8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L57:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_tag: D2
type_guessed: f64
name: read_data
arguments: data_tag: D138
type_guessed: f64
name: read_data
arguments: data_tag: D78
type_guessed: f64
name: read_data
arguments: data_tag: D165
type_guessed: f64
name: read_data
arguments: data_tag: D74
type_guessed: f64
name: read_data
arguments: data_tag: D115
type_guessed: f64

tool: D2: 0.6234898018587335

tool: D138: 0.9009688679024191

tool: D78: 0.2225209339563144

tool: D165: 0.9749279121818236

tool: D74: 0.7818314824680298

tool: D115: 0.4338837391175581

assistant: ```c
void FFT_7(int sofarRadix, int radix, int remainRadix, int groupOffset, double yRe[], double yIm[]) {
  int yIndx[7];
  int groupNo, blockNo, adr, i;
  double are0, aim0, aRe[7], aIm[7];
  double are16add, are61sub, aim16sub, aim16add, are25add, are52sub, aim25sub, aim25add;
  double are34add, are43sub, aim34sub, aim34add;
  double ax, ay;
  adr = groupOffset;
  for (groupNo = 0; groupNo < remainRadix; groupNo++) {
    for (blockNo = 0; blockNo < radix; blockNo++) {
      yIndx[blockNo] = adr;
      adr = adr + sofarRadix;
    }
    are0 = yRe[yIndx[0]];
    aim0 = yIm[yIndx[0]];
    i = 1;
    for (blockNo = 1; blockNo < radix; blockNo += 1) {
      aRe[blockNo] = yRe[yIndx[blockNo]];
      aIm[blockNo] = yIm[yIndx[blockNo]];
    }
    ax = aRe[1];
    ay = aRe[6];
    are16add = ax + ay;
    are61sub = ay - ax;
    ax = aIm[1];
    ay = aIm[6];
    aim16sub = ax - ay;
    aim16add = ax + ay;
    ax = aRe[2];
    ay = aRe[5];
    are25add = ax + ay;
    are52sub = ay - ax;
    ax = aIm[2];
    ay = aIm[5];
    aim25sub = ax - ay;
    aim25add = ax + ay;
    ax = aRe[3];
    ay = aRe[4];
    are34add = ax + ay;
    are43sub = ay - ax;
    ax = aIm[3];
    ay = aIm[4];
    aim34sub = ax - ay;
    aim34add = ax + ay;
    ax = 0.974927912181823607018131682993931217232785801 * aim16sub - 0.781831482468029808708444526674057750232334519 * aim34sub - 0.433883739117558120475768332848358754609990728 * aim25sub;
    ay = 0.623489801858733530525004884004239810632274731 * are34add - 0.900968867902419126236102319507445051165919162 * are25add - 0.222520933956314404288902564496794759466355569 * are16add + are0;
    yRe[yIndx[5]] = ay - ax;
    yRe[yIndx[2]] = ay + ax;
    ax = 0.974927912181823607018131682993931217232785801 * are61sub - 0.781831482468029808708444526674057750232334519 * are43sub - 0.433883739117558120475768332848358754609990728 * are52sub;
    ay = 0.623489801858733530525004884004239810632274731 * aim34add - 0.900968867902419126236102319507445051165919162 * aim25add - 0.222520933956314404288902564496794759466355569 * aim16add + aim0;
    yIm[yIndx[2]] = ax + ay;
    yIm[yIndx[5]] = ay - ax;
    ax = 0.781831482468029808708444526674057750232334519 * aim16sub + 0.974927912181823607018131682993931217232785801 * aim25sub + 0.433883739117558120475768332848358754609990728 * aim34sub;
    ay = 0.623489801858733530525004884004239810632274731 * are16add - 0.222520933956314404288902564496794759466355569 * are25add - 0.900968867902419126236102319507445051165919162 * are34add + are0;
    yRe[yIndx[6]] = ay - ax;
    yRe[yIndx[1]] = ay + ax;
    ax = 0.781831482468029808708444526674057750232334519 * are61sub + 0.974927912181823607018131682993931217232785801 * are52sub + 0.433883739117558120475768332848358754609990728 * are43sub;
    ay = 0.623489801858733530525004884004239810632274731 * aim16add - 0.222520933956314404288902564496794759466355569 * aim25add - 0.900968867902419126236102319507445051165919162 * aim34add + aim0;
    yIm[yIndx[1]] = ax + ay;
    yIm[yIndx[6]] = ay - ax;
    ax = 0.433883739117558120475768332848358754609990728 * aim16sub + 0.974927912181823607018131682993931217232785801 * aim34sub - 0.781831482468029808708444526674057750232334519 * aim25sub;
    ay = 0.623489801858733530525004884004239810632274731 * are25add - 0.222520933956314404288902564496794759466355569 * are34add - 0.900968867902419126236102319507445051165919162 * are16add + are0;
    yRe[yIndx[4]] = ay - ax;
    yRe[yIndx[3]] = ay + ax;
    ax = 0.433883739117558120475768332848358754609990728 * are61sub + 0.974927912181823607018131682993931217232785801 * are43sub - 0.781831482468029808708444526674057750232334519 * are52sub;
    ay = 0.623489801858733530525004884004239810632274731 * aim25add - 0.222520933956314404288902564496794759466355569 * aim34add - 0.900968867902419126236102319507445051165919162 * aim16add + aim0;
    yIm[yIndx[3]] = ax + ay;
    yIm[yIndx[4]] = ay - ax;
    yRe[yIndx[0]] = are0 + are16add + are25add + are34add;
    yIm[yIndx[0]] = aim0 + aim16add + aim25add + aim34add;
    groupOffset = groupOffset + sofarRadix * radix;
    adr = groupOffset;
  }
}
```