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

```asm
<print_aft_stats_linear>:
  endbr64
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x60,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x58(%rsp)
  xor    %eax,%eax
  test   %esi,%esi
  jle    L27
  lea    -0x1(%rsi),%eax
  pxor   %xmm0,%xmm0
  mov    %rdi,%rbp
  mov    %esi,%r13d
  shl    $0x4,%rax
  movaps %xmm0,0x20(%rsp)
  mov    %rdx,%rbx
  lea    0x10(%rdx,%rax,1),%r12
  mov    D82(%rip),%rax
  movq   $0x0,0x40(%rsp)
  movl   $0x0,0x48(%rsp)
  movq   $0x0,0x10(%rsp)
  mov    %rax,0x8(%rsp)
  movaps %xmm0,0x30(%rsp)
  jmp    L56
  nopl   0x0(%rax,%rax,1)
L50:
  comisd D82(%rip),%xmm0
  jbe    L21
  addl   $0x1,0x44(%rsp)
L19:
  addsd  0x10(%rsp),%xmm1
  add    $0x10,%rbx
  movsd  %xmm1,0x10(%rsp)
  cmp    %r12,%rbx
  je     L20
L56:
  mov    0x8(%rbx),%ecx
  mov    0x4(%rbx),%edx
  mov    %rbp,%rdi
  mov    0xc(%rbx),%r8d
  mov    (%rbx),%esi
  call   <aft_tetra_qual@plt>
  pxor   %xmm5,%xmm5
  movapd %xmm0,%xmm2
  comisd %xmm0,%xmm5
  minsd  0x8(%rsp),%xmm2
  movapd %xmm0,%xmm1
  movsd  %xmm2,0x8(%rsp)
  jb     L50
  addl   $0x1,0x48(%rsp)
  jmp    L19
  nopl   0x0(%rax)
L20:
  pxor   %xmm7,%xmm7
  movapd %xmm1,%xmm0
  mov    $0x1,%eax
  cvtsi2sd %r13d,%xmm7
  lea    D74(%rip),%rdi
  lea    0x20(%rsp),%rbx
  lea    0x48(%rsp),%r12
  mov    %rbx,%rbp
  lea    D71(%rip),%r13
  divsd  %xmm7,%xmm0
  movsd  %xmm7,0x18(%rsp)
  call   <printf@plt>
  mov    0x48(%rsp),%r14d
  lea    D42(%rip),%rdi
  xor    %eax,%eax
  mov    %r14d,%esi
  call   <printf@plt>
  nopw   0x0(%rax,%rax,1)
L63:
  mov    0x0(%rbp),%esi
  mov    %r13,%rdi
  xor    %eax,%eax
  add    $0x4,%rbp
  call   <printf@plt>
  cmp    %rbp,%r12
  jne    L63
  mov    $0xa,%edi
  lea    D96(%rip),%rbp
  call   <putchar@plt>
  movsd  0x8(%rsp),%xmm0
  mov    $0x1,%eax
  lea    D14(%rip),%rdi
  call   <printf@plt>
  pxor   %xmm0,%xmm0
  mov    $0x1,%eax
  lea    D34(%rip),%rdi
  cvtsi2sd %r14d,%xmm0
  mulsd  D62(%rip),%xmm0
  divsd  0x18(%rsp),%xmm0
  call   <printf@plt>
  xchg   %ax,%ax
L77:
  pxor   %xmm0,%xmm0
  mov    %rbp,%rdi
  mov    $0x1,%eax
  add    $0x4,%rbx
  cvtsi2sdl -0x4(%rbx),%xmm0
  mulsd  D62(%rip),%xmm0
  divsd  0x18(%rsp),%xmm0
  call   <printf@plt>
  cmp    %rbx,%r12
  jne    L77
  mov    $0xa,%edi
  call   <putchar@plt>
L27:
  mov    0x58(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L30
  add    $0x60,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
  cs nopw 0x0(%rax,%rax,1)
L21:
  movsd  D1(%rip),%xmm0
  movsd  %xmm1,0x18(%rsp)
  mulsd  %xmm1,%xmm0
  call   <ceil@plt>
  movsd  0x18(%rsp),%xmm1
  sub    $0x1,%eax
  cltq
  addl   $0x1,0x20(%rsp,%rax,4)
  jmp    L19
L30:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: tag: D82
data_type: f64[2]
name: parse_data
arguments: tag: D74
data_type: string
name: parse_data
arguments: tag: D71
data_type: string
name: parse_data
arguments: tag: D42
data_type: string
name: parse_data
arguments: tag: D96
data_type: string
name: parse_data
arguments: tag: D14
data_type: string
name: parse_data
arguments: tag: D34
data_type: string
name: parse_data
arguments: tag: D62
data_type: f64
name: parse_data
arguments: tag: D1
data_type: f32[2]

tool: {"D82": [1.0,10.0]}

tool: D74: QS: q_avg: %12.8le  

tool: D71:  %6d 

tool: D42:  ## %6d # 

tool: D96:  %6.2lf 

tool: D14:     q_min: %12.8le  

tool: D34:  ## %6.2lf # 

tool: D62: 100

tool: {"D1": [0.0,2.5625]}

assistant: ```c
int print_aft_stats_linear(double *vertex, int nt, int *tetra) {
  double q, worst, qs;
  int stat[11], i;
  if (nt < 1)
    return 0;
  for (i = 0; i < 11; i++)
    stat[i] = 0;
  worst = 1.0, qs = 0.0;
  for (i = 0; i < nt; i++) {
    q = aft_tetra_qual(vertex, tetra[4 * i + 0], tetra[4 * i + 1], tetra[4 * i + 2], tetra[4 * i + 3]);
    if (q < worst)
      worst = q;
    if (q <= 0.0)
      stat[10]++;
    else if (q > 1.0)
      stat[9]++;
    else
      stat[(int)ceil(q * 10) - 1]++;
    qs += q;
  }
  qs /= nt;
  printf("QS: q_avg: %12.8le  ", qs);
  printf(" ## %6d # ", stat[10]);
  for (i = 0; i < 10; i++)
    printf(" %6d ", stat[i]);
  printf("\n");
  printf("    q_min: %12.8le  ", worst);
  printf(" ## %6.2lf # ", 100.0 * stat[10] / nt);
  for (i = 0; i < 10; i++)
    printf(" %6.2lf ", 100.0 * stat[i] / nt);
  printf("\n");
  return 0;
}
```