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

```asm
<print_aft_stats_linear_mba>:
  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    L22
  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    D98(%rip),%rax
  movq   $0x0,0x40(%rsp)
  movl   $0x0,0x48(%rsp)
  movq   $0x0,0x10(%rsp)
  mov    %rax,0x8(%rsp)
  movaps %xmm0,0x30(%rsp)
  jmp    L64
  nopl   0x0(%rax,%rax,1)
L12:
  comisd D98(%rip),%xmm0
  jbe    L98
  addl   $0x1,0x44(%rsp)
L50:
  addsd  0x10(%rsp),%xmm1
  add    $0x10,%rbx
  movsd  %xmm1,0x10(%rsp)
  cmp    %rbx,%r12
  je     L21
L64:
  mov    0x8(%rbx),%ecx
  mov    0x4(%rbx),%edx
  mov    %rbp,%rdi
  mov    0xc(%rbx),%r8d
  mov    (%rbx),%esi
  call   <aft_tetra_qual_mba@plt>
  pxor   %xmm5,%xmm5
  movapd %xmm0,%xmm2
  comisd %xmm0,%xmm5
  minsd  0x8(%rsp),%xmm2
  movapd %xmm0,%xmm1
  movsd  %xmm2,0x8(%rsp)
  jb     L12
  addl   $0x1,0x48(%rsp)
  jmp    L50
  nopl   0x0(%rax)
L21:
  pxor   %xmm7,%xmm7
  movapd %xmm1,%xmm0
  mov    $0x1,%eax
  cvtsi2sd %r13d,%xmm7
  lea    D58(%rip),%rdi
  lea    0x20(%rsp),%rbx
  lea    0x48(%rsp),%r12
  mov    %rbx,%rbp
  lea    D22(%rip),%r13
  divsd  %xmm7,%xmm0
  movsd  %xmm7,0x18(%rsp)
  call   <printf@plt>
  mov    0x48(%rsp),%r14d
  lea    D24(%rip),%rdi
  xor    %eax,%eax
  mov    %r14d,%esi
  call   <printf@plt>
  nopw   0x0(%rax,%rax,1)
L35:
  mov    0x0(%rbp),%esi
  mov    %r13,%rdi
  xor    %eax,%eax
  add    $0x4,%rbp
  call   <printf@plt>
  cmp    %rbp,%r12
  jne    L35
  mov    $0xa,%edi
  lea    D109(%rip),%rbp
  call   <putchar@plt>
  movsd  0x8(%rsp),%xmm0
  mov    $0x1,%eax
  lea    D47(%rip),%rdi
  call   <printf@plt>
  pxor   %xmm0,%xmm0
  mov    $0x1,%eax
  lea    D60(%rip),%rdi
  cvtsi2sd %r14d,%xmm0
  mulsd  D48(%rip),%xmm0
  divsd  0x18(%rsp),%xmm0
  call   <printf@plt>
  xchg   %ax,%ax
L14:
  pxor   %xmm0,%xmm0
  mov    %rbp,%rdi
  mov    $0x1,%eax
  add    $0x4,%rbx
  cvtsi2sdl -0x4(%rbx),%xmm0
  mulsd  D48(%rip),%xmm0
  divsd  0x18(%rsp),%xmm0
  call   <printf@plt>
  cmp    %rbx,%r12
  jne    L14
  mov    $0xa,%edi
  call   <putchar@plt>
L22:
  mov    0x58(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L32
  add    $0x60,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
  cs nopw 0x0(%rax,%rax,1)
L98:
  movsd  D19(%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    L50
L32:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_label: D98
datatype: f64[2]
name: read_data
arguments: data_label: D58
datatype: string
name: read_data
arguments: data_label: D22
datatype: string
name: read_data
arguments: data_label: D24
datatype: string
name: read_data
arguments: data_label: D109
datatype: string
name: read_data
arguments: data_label: D47
datatype: string
name: read_data
arguments: data_label: D60
datatype: string
name: read_data
arguments: data_label: D48
datatype: f64
name: read_data
arguments: data_label: D19
datatype: f32[2]

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

tool: D58: [mba] QS: q_avg: %12.8le  

tool: D22:  %6d 

tool: D24:  ## %6d # 

tool: D109:  %6.2lf 

tool: D47: [mba]     q_min: %12.8le  

tool: D60:  ## %6.2lf # 

tool: D48: 100

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

assistant: ```c
int print_aft_stats_linear_mba(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_mba(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("[mba] QS: q_avg: %12.8le  ", qs);
  printf(" ## %6d # ", stat[10]);
  for (i = 0; i < 10; i++)
    printf(" %6d ", stat[i]);
  printf("\n");
  printf("[mba]     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;
}
```