user: What would the C source look like for the assembly code provided here:

```asm
<ERI_Interpolate_Linear_Complex>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x60,%rsp
  mov    %rdi,-0x38(%rbp)
  movsd  %xmm0,-0x40(%rbp)
  mov    %rsi,-0x48(%rbp)
  mov    %rdx,-0x50(%rbp)
  mov    %ecx,-0x54(%rbp)
  cmpl   $0x0,-0x54(%rbp)
  setg   %al
  movzbl %al,%eax
  mov    %eax,%edi
  call   <assert@plt>
  movl   $0x0,-0x24(%rbp)
  jmp    L49
L0:
  addl   $0x1,-0x24(%rbp)
L49:
  mov    -0x24(%rbp),%eax
  cmp    -0x54(%rbp),%eax
  jge    L25
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x40(%rbp),%xmm0
  comisd %xmm1,%xmm0
  ja     L0
L25:
  cmpl   $0x0,-0x24(%rbp)
  jne    L5
  mov    -0x50(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x38(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x50(%rbp),%rax
  lea    0x8(%rax),%rdx
  mov    -0x38(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L29
L5:
  mov    -0x24(%rbp),%eax
  cmp    -0x54(%rbp),%eax
  jne    L32
  mov    -0x54(%rbp),%eax
  add    %eax,%eax
  cltq
  shl    $0x3,%rax
  lea    -0x10(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x38(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x54(%rbp),%eax
  add    %eax,%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rax,%rdx
  mov    -0x38(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L29
L32:
  mov    -0x24(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x24(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  mov    -0x24(%rbp),%eax
  add    %eax,%eax
  cltq
  shl    $0x3,%rax
  lea    -0x10(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x24(%rbp),%eax
  add    %eax,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  -0x10(%rbp),%xmm1
  movsd  -0x40(%rbp),%xmm0
  subsd  -0x20(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x18(%rbp),%xmm1
  subsd  -0x20(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  addsd  -0x10(%rbp),%xmm0
  mov    -0x38(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x24(%rbp),%eax
  add    %eax,%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x24(%rbp),%eax
  add    %eax,%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x8(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  -0x10(%rbp),%xmm1
  movsd  -0x40(%rbp),%xmm0
  subsd  -0x20(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x18(%rbp),%xmm1
  subsd  -0x20(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  mov    -0x38(%rbp),%rax
  add    $0x8,%rax
  addsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,(%rax)
L29:
  leave
  ret
```

assistant: ```c
void ERI_Interpolate_Linear_Complex(double out[2], double x, const double *sx, const double *sy, int ns) {
  int i;
  double x0, x1, y0, y1;
  assert(0 < ns);
  for (i = 0; i < ns && sx[i] < x; i++)
    ;
  if (i == 0) {
    out[0] = sy[0];
    out[1] = sy[1];
    return;
  }
  if (i == ns) {
    out[0] = sy[2 * ns - 2];
    out[1] = sy[2 * ns - 1];
    return;
  }
  x0 = sx[i - 1];
  x1 = sx[i];
  y0 = sy[2 * i - 2];
  y1 = sy[2 * i + 0];
  out[0] = y0 + (y1 - y0) * (x - x0) / (x1 - x0);
  y0 = sy[2 * i - 1];
  y1 = sy[2 * i + 1];
  out[1] = y0 + (y1 - y0) * (x - x0) / (x1 - x0);
}
```