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

```asm
<seventhOrderIIRFilter>:
  endbr64
  push   %r15
  mov    %rsi,%r15
  push   %r14
  mov    %ecx,%r14d
  push   %r13
  mov    %rdx,%r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  sub    $0x18,%rsp
  movsd  (%rsi),%xmm0
  ucomisd D0(%rip),%xmm0
  jp     L4
  jne    L4
L46:
  mov    $0x8,%esi
  mov    $0x6,%edi
  call   <calloc@plt>
  mov    $0x8,%esi
  mov    $0x6,%edi
  mov    %rax,%rbp
  call   <calloc@plt>
  mov    $0x8,%esi
  mov    %r14d,%edi
  mov    %rax,%r12
  call   <calloc@plt>
  pxor   %xmm0,%xmm0
  movups %xmm0,(%r12)
  mov    %rax,%r8
  movups %xmm0,0x10(%r12)
  movups %xmm0,0x20(%r12)
  movups %xmm0,0x0(%rbp)
  movups %xmm0,0x10(%rbp)
  movups %xmm0,0x20(%rbp)
  test   %r14d,%r14d
  jle    L42
  mov    %r14d,%edx
  xor    %esi,%esi
  mov    %rax,%rdi
  shl    $0x3,%rdx
  call   <memset@plt>
  movsd  (%r12),%xmm1
  movsd  0x8(%r12),%xmm6
  mov    %rax,%r8
  movsd  0x10(%r12),%xmm5
  lea    -0x1(%r14),%eax
  movsd  0x18(%r12),%xmm4
  mov    %r8,%rcx
  movsd  0x20(%r12),%xmm3
  lea    0x8(%r13,%rax,8),%rdx
  movsd  0x28(%r12),%xmm2
  sub    %r13,%rcx
L78:
  movsd  (%rbx),%xmm0
  movsd  0x8(%rbx),%xmm7
  lea    (%rcx,%r13,1),%rax
  add    $0x8,%r13
  mulsd  -0x8(%r13),%xmm0
  mulsd  0x0(%rbp),%xmm7
  mulsd  0x10(%r15),%xmm6
  mulsd  0x8(%r15),%xmm1
  mulsd  0x18(%r15),%xmm5
  mulsd  0x20(%r15),%xmm4
  addsd  %xmm7,%xmm0
  movsd  0x10(%rbx),%xmm7
  mulsd  0x8(%rbp),%xmm7
  mulsd  0x28(%r15),%xmm3
  addsd  %xmm6,%xmm1
  mulsd  0x30(%r15),%xmm2
  addsd  %xmm7,%xmm0
  movsd  0x18(%rbx),%xmm7
  mulsd  0x10(%rbp),%xmm7
  addsd  %xmm5,%xmm1
  addsd  %xmm4,%xmm1
  addsd  %xmm7,%xmm0
  movsd  0x20(%rbx),%xmm7
  mulsd  0x18(%rbp),%xmm7
  addsd  %xmm3,%xmm1
  addsd  %xmm2,%xmm1
  addsd  %xmm7,%xmm0
  movsd  0x28(%rbx),%xmm7
  mulsd  0x20(%rbp),%xmm7
  addsd  %xmm7,%xmm0
  movsd  0x30(%rbx),%xmm7
  mulsd  0x28(%rbp),%xmm7
  addsd  %xmm7,%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  movsd  0x20(%rbp),%xmm0
  movsd  %xmm0,0x28(%rbp)
  movsd  0x18(%rbp),%xmm0
  movsd  %xmm0,0x20(%rbp)
  movsd  0x10(%rbp),%xmm0
  movsd  %xmm0,0x18(%rbp)
  movsd  0x8(%rbp),%xmm0
  movsd  %xmm0,0x10(%rbp)
  movsd  0x0(%rbp),%xmm0
  movsd  %xmm0,0x8(%rbp)
  movsd  -0x8(%r13),%xmm0
  movsd  %xmm0,0x0(%rbp)
  movsd  0x20(%r12),%xmm2
  movsd  0x18(%r12),%xmm3
  movsd  0x10(%r12),%xmm4
  movsd  0x8(%r12),%xmm5
  movsd  (%r12),%xmm6
  movsd  %xmm2,0x28(%r12)
  movsd  %xmm3,0x20(%r12)
  movsd  %xmm4,0x18(%r12)
  movsd  %xmm5,0x10(%r12)
  movsd  %xmm6,0x8(%r12)
  movsd  (%rax),%xmm1
  movsd  %xmm1,(%r12)
  cmp    %rdx,%r13
  jne    L78
L42:
  mov    %rbp,%rdi
  mov    %r8,0x8(%rsp)
  call   <free@plt>
  mov    %r12,%rdi
  call   <free@plt>
  mov    0x8(%rsp),%rax
  add    $0x18,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L4:
  divsd  %xmm0,%xmm0
  lea    0x8(%r15),%rdx
  lea    0x18(%rbx),%rax
  cmp    %rax,%rdx
  lea    0x18(%r15),%rcx
  lea    0x8(%rbx),%rdx
  setae  %al
  cmp    %rdx,%rcx
  setbe  %dl
  or     %dl,%al
  movsd  %xmm0,(%r15)
  movsd  (%rbx),%xmm1
  divsd  %xmm0,%xmm1
  movsd  %xmm1,(%rbx)
  je     L44
  lea    0x38(%rbx),%rax
  cmp    %rax,%r15
  setae  %al
  cmp    %rbx,%r15
  setbe  %dl
  or     %dl,%al
  je     L44
  movsd  (%r15),%xmm0
  movupd 0x8(%r15),%xmm1
  unpcklpd %xmm0,%xmm0
  divpd  %xmm0,%xmm1
  movups %xmm1,0x8(%r15)
  movupd 0x8(%rbx),%xmm1
  divpd  %xmm0,%xmm1
  movups %xmm1,0x8(%rbx)
  movupd 0x18(%r15),%xmm1
  divpd  %xmm0,%xmm1
  movups %xmm1,0x18(%r15)
  movupd 0x18(%rbx),%xmm1
  divpd  %xmm0,%xmm1
  movups %xmm1,0x18(%rbx)
  movupd 0x28(%r15),%xmm1
  divpd  %xmm0,%xmm1
  movups %xmm1,0x28(%r15)
  movupd 0x28(%rbx),%xmm1
  divpd  %xmm0,%xmm1
  movups %xmm1,0x28(%rbx)
  jmp    L46
  nopl   0x0(%rax,%rax,1)
L44:
  movsd  (%r15),%xmm1
  movsd  0x8(%r15),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x8(%r15)
  movsd  0x8(%rbx),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x8(%rbx)
  movsd  (%r15),%xmm1
  movsd  0x10(%r15),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x10(%r15)
  movsd  0x10(%rbx),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x10(%rbx)
  movsd  (%r15),%xmm1
  movsd  0x18(%r15),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x18(%r15)
  movsd  0x18(%rbx),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x18(%rbx)
  movsd  (%r15),%xmm1
  movsd  0x20(%r15),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x20(%r15)
  movsd  0x20(%rbx),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x20(%rbx)
  movsd  (%r15),%xmm1
  movsd  0x28(%r15),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x28(%r15)
  movsd  0x28(%rbx),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x28(%rbx)
  movsd  (%r15),%xmm1
  movsd  0x30(%r15),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x30(%r15)
  movsd  0x30(%rbx),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x30(%rbx)
  jmp    L46
```

tool_calls: name: read_data
arguments: data_label: D0
type_guessed: f32[2]

tool: {"D0": [0.0,1.875]}

assistant: ```c
double *seventhOrderIIRFilter(double *B, double *A, double *signal, int signalLength) {
  int i;
  int filterLength = 7;
  double *x, *y, *output;
  double zeros = 0;
  double poles = 0;
  if (A[0] != 1) {
    for (i = 0; i < filterLength; i++) {
      A[i] = A[i] / A[0];
      B[i] = B[i] / A[0];
    }
  }
  x = (double *)calloc((filterLength - 1), sizeof(double));
  y = (double *)calloc((filterLength - 1), sizeof(double));
  output = (double *)calloc((signalLength), sizeof(double));
  for (i = 0; i < filterLength - 1; i++)
    y[i] = 0;
  for (i = 0; i < filterLength - 1; i++)
    x[i] = 0;
  for (i = 0; i < signalLength; i++)
    output[i] = 0;
  for (i = 0; i < signalLength; i++) {
    zeros = B[0] * signal[i] + B[1] * x[0] + B[2] * x[1] + B[3] * x[2] + B[4] * x[3] + B[5] * x[4] + B[6] * x[5];
    poles = A[1] * y[0] + A[2] * y[1] + A[3] * y[2] + A[4] * y[3] + A[5] * y[4] + A[6] * y[5];
    output[i] = zeros - poles;
    x[5] = x[4];
    x[4] = x[3];
    x[3] = x[2];
    x[2] = x[1];
    x[1] = x[0];
    x[0] = signal[i];
    y[5] = y[4];
    y[4] = y[3];
    y[3] = y[2];
    y[2] = y[1];
    y[1] = y[0];
    y[0] = output[i];
  }
  free(x);
  free(y);
  return output;
}
```