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

```asm
<doubleBoxWhiskerCalc>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  mov    %edi,-0x44(%rbp)
  mov    %rsi,-0x50(%rbp)
  mov    %rdx,-0x58(%rbp)
  mov    %rcx,-0x60(%rbp)
  mov    %r8,-0x68(%rbp)
  mov    %r9,-0x70(%rbp)
  cmpl   $0x0,-0x44(%rbp)
  jg     L46
  mov    -0x44(%rbp),%eax
  mov    %eax,%esi
  lea    D24(%rip),%rax
  mov    %rax,%rdi
  call   <errAbort@plt>
L46:
  cmpl   $0x1,-0x44(%rbp)
  jne    L67
  mov    -0x50(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    0x10(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    0x10(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x70(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x70(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x68(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x68(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x60(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x60(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x58(%rbp),%rax
  movsd  %xmm0,(%rax)
  jmp    L50
L67:
  mov    -0x50(%rbp),%rdx
  mov    -0x44(%rbp),%eax
  mov    %rdx,%rsi
  mov    %eax,%edi
  call   <doubleSort@plt>
  mov    -0x50(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x44(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0x44(%rbp),%eax
  sar    %eax
  mov    %eax,-0x34(%rbp)
  mov    -0x44(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L27
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  jmp    L18
L27:
  mov    -0x34(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x34(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  addsd  %xmm0,%xmm1
  movsd  D16(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
L18:
  cmpl   $0x3,-0x44(%rbp)
  jg     L29
  movsd  -0x28(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x10(%rbp),%xmm1
  movsd  D16(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x28(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x8(%rbp),%xmm1
  movsd  D16(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  jmp    L62
L29:
  mov    -0x44(%rbp),%eax
  lea    0x3(%rax),%edx
  test   %eax,%eax
  cmovs  %edx,%eax
  sar    $0x2,%eax
  mov    %eax,-0x30(%rbp)
  mov    -0x44(%rbp),%eax
  sub    $0x1,%eax
  sub    -0x30(%rbp),%eax
  mov    %eax,-0x2c(%rbp)
  mov    -0x2c(%rbp),%ecx
  mov    -0x30(%rbp),%edx
  mov    -0x44(%rbp),%eax
  mov    %ecx,%r8d
  mov    %edx,%ecx
  mov    %eax,%edx
  lea    D30(%rip),%rax
  mov    %rax,%rsi
  mov    $0x4,%edi
  call   <verbose@plt>
  mov    -0x30(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x2c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
L62:
  mov    -0x58(%rbp),%rax
  movsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x60(%rbp),%rax
  movsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x68(%rbp),%rax
  movsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    -0x70(%rbp),%rax
  movsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  mov    0x10(%rbp),%rax
  movsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,(%rax)
L50:
  leave
  ret
```

tool_calls: name: read_data
arguments: label: D24
type: string
name: read_data
arguments: label: D16
type: f64
name: read_data
arguments: label: D30
type: string

tool: D24: doubleBoxWhiskerCalc needs a positive number, not %d for count

tool: D16: 0.5

tool: D30: count %d, q1Ix %d, q3Ix %d


assistant: ```c
void doubleBoxWhiskerCalc(int count, double *array, double *retMin, double *retQ1, double *retMedian, double *retQ3, double *retMax) {
  if (count <= 0)
    errAbort("doubleBoxWhiskerCalc needs a positive number, not %d for count", count);
  if (count == 1) {
    *retMin = *retQ1 = *retMedian = *retQ3 = *retMax = array[0];
    return;
  }
  doubleSort(count, array);
  double min = array[0];
  double max = array[count - 1];
  double median;
  int halfCount = count >> 1;
  if ((count & 1) == 1)
    median = array[halfCount];
  else {
    median = (array[halfCount] + array[halfCount - 1]) * 0.5;
  }
  double q1, q3;
  if (count <= 3) {
    q1 = 0.5 * (median + min);
    q3 = 0.5 * (median + max);
  } else {
    int q1Ix = count / 4;
    int q3Ix = count - 1 - q1Ix;
    verbose(4, "count %d, q1Ix %d, q3Ix %d\n", count, q1Ix, q3Ix);
    q1 = array[q1Ix];
    q3 = array[q3Ix];
  }
  *retMin = min;
  *retQ1 = q1;
  *retMedian = median;
  *retQ3 = q3;
  *retMax = max;
}
```