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

```asm
<main>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x428,%rsp
  movdqa D113(%rip),%xmm2
  movdqa D111(%rip),%xmm7
  lea    0x10(%rsp),%rbp
  movapd D31(%rip),%xmm5
  movapd D51(%rip),%xmm4
  lea    0x410(%rsp),%r13
  movapd D25(%rip),%xmm3
  movdqa D61(%rip),%xmm6
  mov    %fs:0x28,%rax
  mov    %rax,0x418(%rsp)
  xor    %eax,%eax
  mov    %rbp,%rax
  nopl   0x0(%rax,%rax,1)
L205:
  movdqa %xmm2,%xmm0
  movapd %xmm3,%xmm1
  paddd  %xmm7,%xmm2
  add    $0x10,%rax
  cvtdq2pd %xmm0,%xmm8
  mulpd  %xmm5,%xmm8
  pshufd $0xee,%xmm0,%xmm0
  cvtdq2pd %xmm0,%xmm0
  mulpd  %xmm5,%xmm0
  addpd  %xmm4,%xmm8
  addpd  %xmm4,%xmm0
  divpd  %xmm8,%xmm1
  movapd %xmm3,%xmm8
  divpd  %xmm0,%xmm8
  cvttpd2dq %xmm1,%xmm0
  cvttpd2dq %xmm8,%xmm1
  punpcklqdq %xmm1,%xmm0
  psrad  $0x8,%xmm0
  pand   %xmm6,%xmm0
  movaps %xmm0,-0x10(%rax)
  cmp    %rax,%r13
  jne    L205
  lea    D16(%rip),%r12
  nopl   (%rax)
L116:
  xor    %ebx,%ebx
L66:
  mov    0x0(%rbp,%rbx,4),%esi
  mov    %r12,%rdi
  xor    %eax,%eax
  add    $0x1,%rbx
  call   <printf@plt>
  cmp    $0x4,%rbx
  jne    L66
  mov    $0xa,%edi
  add    $0x10,%rbp
  call   <putchar@plt>
  cmp    %rbp,%r13
  jne    L116
  pxor   %xmm9,%xmm9
  xor    %ecx,%ecx
  xor    %edi,%edi
  mov    $0x800000,%r8d
  movsd  D83(%rip),%xmm7
  movsd  D115(%rip),%xmm4
  movapd %xmm9,%xmm0
  movapd %xmm9,%xmm6
  movsd  D51(%rip),%xmm10
  movapd %xmm9,%xmm5
  movsd  D25(%rip),%xmm8
  jmp    L87
  nopw   0x0(%rax,%rax,1)
L27:
  sub    %r9d,%eax
  pxor   %xmm3,%xmm3
  sub    %r9d,%edx
  sub    %r9d,%esi
  pxor   %xmm2,%xmm2
  pxor   %xmm1,%xmm1
  add    $0x1,%edi
  sub    $0xffffff80,%ecx
  cvtsi2sd %eax,%xmm3
  cvtsi2sd %edx,%xmm2
  cvtsi2sd %esi,%xmm1
  mulsd  %xmm4,%xmm3
  mulsd  %xmm4,%xmm2
  mulsd  %xmm4,%xmm1
  mulsd  %xmm3,%xmm3
  mulsd  %xmm2,%xmm2
  mulsd  %xmm1,%xmm1
  addsd  %xmm3,%xmm5
  addsd  %xmm2,%xmm6
  addsd  %xmm1,%xmm0
  cmp    $0x10000,%edi
  je     L129
L87:
  pxor   %xmm1,%xmm1
  movapd %xmm8,%xmm2
  mov    %ecx,%r10d
  mov    %ecx,%r13d
  cvtsi2sd %edi,%xmm1
  sar    $0xc,%r10d
  sar    %r13d
  mov    %ecx,%ebp
  imul   %r10d,%r10d
  mov    %r8d,%eax
  sar    $0x2,%ebp
  mov    %ecx,%ebx
  sub    %r13d,%eax
  sar    $0x3,%ebx
  mov    %ecx,%r11d
  sub    %ebp,%eax
  and    $0x7fff,%r11d
  mulsd  %xmm7,%xmm1
  mov    %r10d,%r12d
  sub    %ebx,%eax
  mov    %r10d,%edx
  sar    %r12d
  sar    $0x2,%edx
  mov    %r11d,%r15d
  mov    %r11d,%r14d
  lea    (%rax,%r12,1),%esi
  sar    $0x5,%r15d
  sar    $0x6,%r11d
  add    %edx,%esi
  mov    %ecx,%edx
  sub    %r15d,%r14d
  sar    $0xf,%edx
  add    %r14d,%r11d
  addsd  %xmm10,%xmm1
  movslq %edx,%rdx
  mov    0x10(%rsp,%rdx,4),%edx
  divsd  %xmm1,%xmm2
  shl    $0x8,%edx
  mov    %edx,%r14d
  and    $0xf00000,%r14d
  lea    -0x700000(%r14),%r15d
  and    $0xffe00000,%r15d
  cvttsd2si %xmm2,%r9d
  je     L157
  cmp    $0x400000,%r14d
  jne    L95
  sar    $0x2,%r11d
L157:
  mov    %r9d,%r14d
  sub    %r11d,%edx
  sar    $0x14,%r14d
  mov    %r14d,%r11d
  and    $0xf,%r11d
  cmp    $0x6,%r11d
  jg     L139
  and    $0xc,%r14d
  mov    %esi,%eax
  jne    L27
L115:
  mov    %r8d,%eax
  sub    %ecx,%eax
  add    %r13d,%eax
  sub    %ebp,%eax
  add    %ebx,%eax
  jmp    L27
  nopl   0x0(%rax)
L139:
  cmp    $0x7,%r11d
  je     L159
  cmp    $0x8,%r11d
  jne    L115
L159:
  sub    %r12d,%eax
  sar    $0x4,%r10d
  add    %r10d,%eax
  jmp    L27
L129:
  mulsd  %xmm7,%xmm0
  ucomisd %xmm0,%xmm9
  ja     L112
  sqrtsd %xmm0,%xmm0
  movapd %xmm0,%xmm2
L48:
  mulsd  %xmm7,%xmm6
  ucomisd %xmm6,%xmm9
  movapd %xmm6,%xmm1
  ja     L10
  sqrtsd %xmm1,%xmm1
L12:
  mulsd  %xmm7,%xmm5
  ucomisd %xmm5,%xmm9
  movapd %xmm5,%xmm0
  ja     L65
  sqrtsd %xmm0,%xmm0
L53:
  lea    D5(%rip),%rdi
  mov    $0x3,%eax
  call   <printf@plt>
  mov    0x418(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L202
  add    $0x428,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L95:
  sar    %r11d
  jmp    L157
L202:
  call   <__stack_chk_fail@plt>
L65:
  movsd  %xmm1,0x8(%rsp)
  movsd  %xmm2,(%rsp)
  call   <sqrt@plt>
  movsd  0x8(%rsp),%xmm1
  movsd  (%rsp),%xmm2
  jmp    L53
L10:
  movapd %xmm6,%xmm0
  movsd  %xmm5,0x8(%rsp)
  movsd  %xmm2,(%rsp)
  call   <sqrt@plt>
  movsd  0x8(%rsp),%xmm5
  movsd  (%rsp),%xmm2
  pxor   %xmm9,%xmm9
  mov    D83(%rip),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm7
  jmp    L12
L112:
  movsd  %xmm6,0x8(%rsp)
  movsd  %xmm5,(%rsp)
  call   <sqrt@plt>
  movsd  0x8(%rsp),%xmm6
  movsd  (%rsp),%xmm5
  pxor   %xmm9,%xmm9
  mov    D83(%rip),%rax
  movapd %xmm0,%xmm2
  movq   %rax,%xmm7
  jmp    L48
```

tool_calls: name: parse_data
arguments: data_label: D113
type: byte[16]
name: parse_data
arguments: data_label: D111
type: byte[16]
name: parse_data
arguments: data_label: D31
type: byte[16]
name: parse_data
arguments: data_label: D51
type: f64[2]
name: parse_data
arguments: data_label: D25
type: f64[2]
name: parse_data
arguments: data_label: D61
type: byte[16]
name: parse_data
arguments: data_label: D16
type: string
name: parse_data
arguments: data_label: D83
type: byte[8]
name: parse_data
arguments: data_label: D115
type: byte[8]
name: parse_data
arguments: data_label: D5
type: string

tool: {"D113": [0x00,0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x02,0x00,0x00,0x00,0x03,0x00,0x00,0x00]}

tool: {"D111": [0x04,0x00,0x00,0x00,0x04,0x00,0x00,0x00,0x04,0x00,0x00,0x00,0x04,0x00,0x00,0x00]}

tool: {"D31": [0x00,0x00,0x00,0x00,0x00,0x00,0x70,0x3f,0x00,0x00,0x00,0x00,0x00,0x00,0x70,0x3f]}

tool: {"D51": [1.0,1.0]}

tool: {"D25": [8388608.0,8388608.0]}

tool: {"D61": [0xfc,0xff,0x00,0x00,0xfc,0xff,0x00,0x00,0xfc,0xff,0x00,0x00,0xfc,0xff,0x00,0x00]}

tool: D16: 16'h%04X, 

tool: {"D83": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3e]}

tool: {"D115": [0x00,0x00,0x00,0x00,0x00,0x00,0x80,0x3e]}

tool: D5: e0=%f e1=%f e2=%f


assistant: ```c
int main() {
  int divtab[256];
  double e0, e1, e2;
  double d0, d1, d2;
  int m1, m2, rf, p0, p1, p2, p3, p4;
  int i0, i1, i2, i3;
  int j0, j1, j2, j3, j4, j5, j6, j7;
  int i, j, k, n;
  e0 = 0;
  e1 = 0;
  e2 = 0;
  for (i = 0; i < 256; i++) {
    rf = 0x800000 / (1.0 + (i / 256.0));
    p0 = (rf + 127) >> 8;
    p0 = rf >> 8;
    p0 &= 0xFFFC;
    divtab[i] = p0;
  }
  for (i = 0; i < 64; i++) {
    for (j = 0; j < 4; j++) {
      printf("16'h%04X, ", divtab[i * 4 + j]);
    }
    printf("\n");
  }
  n = 16;
  for (i = 0; i < (1 << n); i++) {
    m1 = i << (23 - n);
    rf = 0x800000 / (1.0 + (i / ((1 << n) * 1.0)));
    m2 = (m1 >> 12);
    m2 = m2 * m2;
    p0 = 0x800000 - m1 + (m1 >> 1) - (m1 >> 2) + (m1 >> 3);
    p1 = 0x800000 - (m1 >> 1);
    p2 = 0x800000 - (m1 >> 1) - (m1 >> 4) - (m1 >> 6);
    p2 = 0x800000 - (m1 >> 1) - (m1 >> 2) - (m1 >> 3) + (m2 >> 1) + (m2 >> 2);
    i0 = (m1 >> (23 - 8));
    if (i0 < 0)
      i0 = 0;
    if (i0 > 255)
      i0 = 255;
    i1 = i0 - 1;
    if (i1 < 0)
      i1 = 0;
    if (i1 > 255)
      i1 = 255;
    j4 = divtab[i0];
    j5 = divtab[i1];
    j0 = j5 - j4;
    p1 = divtab[i0] << 8;
    j6 = j4 << 8;
    j7 = j5 << 8;
    i2 = (m1 >> (23 - 16)) & 255;
    j2 = divtab[i2];
    j3 = (m1 & 32767);
    j3 = j3 - (j3 >> 5) + (j3 >> 6);
    j1 = j3 >> 1;
    if ((p1 & 0xF00000) == 0x800000)
      j1 = j1 << 1;
    if ((p1 & 0xF00000) == 0x700000)
      j1 = j1 << 1;
    if ((p1 & 0xF00000) == 0x400000)
      j1 = j1 >> 1;
    p1 -= j1;
    switch ((rf >> 20) & 15) {
    case 8:
      p0 = 0x800000 - (m1 >> 1) - (m1 >> 2) - (m1 >> 3) - (m2 >> 1) + (m2 >> 4);
      break;
    case 7:
      p0 = 0x800000 - (m1 >> 1) - (m1 >> 2) - (m1 >> 3) - (m2 >> 1) + (m2 >> 4);
      break;
    case 6:
      p0 = 0x800000 - (m1 >> 1) - (m1 >> 2) - (m1 >> 3) + (m2 >> 1) + (m2 >> 2);
      break;
    case 5:
      p0 = 0x800000 - (m1 >> 1) - (m1 >> 2) - (m1 >> 3) + (m2 >> 1) + (m2 >> 2);
      break;
    case 4:
      p0 = 0x800000 - (m1 >> 1) - (m1 >> 2) - (m1 >> 3) + (m2 >> 1) + (m2 >> 2);
      break;
    }
    d0 = (p0 - rf) / 8388608.0;
    d1 = (p1 - rf) / 8388608.0;
    d2 = (p2 - rf) / 8388608.0;
    e0 += d0 * d0;
    e1 += d1 * d1;
    e2 += d2 * d2;
  }
  printf("e0=%f e1=%f e2=%f\n", sqrt(e0 / (1 << n)), sqrt(e1 / (1 << n)), sqrt(e2 / (1 << n)));
}
```