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

```asm
<pd_add>:
  endbr64
  push   %r13
  mov    %rdi,%r13
  push   %r12
  mov    %rdx,%r12
  push   %rbp
  push   %rbx
  mov    %rsi,%rbx
  sub    $0x48,%rsp
  movsd  0x20(%rsi),%xmm1
  movsd  0x20(%rdi),%xmm0
  mov    %fs:0x28,%rax
  mov    %rax,0x38(%rsp)
  xor    %eax,%eax
  lea    0x30(%rsp),%rbp
  mov    %rbp,%rdi
  call   <dd_two_sum@plt>
  pxor   %xmm5,%xmm5
  movsd  0x18(%rbx),%xmm1
  mov    %rbp,%rdi
  addsd  0x30(%rsp),%xmm5
  movsd  %xmm0,0x10(%rsp)
  movsd  0x18(%r13),%xmm0
  movsd  %xmm5,0x8(%rsp)
  call   <dd_two_sum@plt>
  movsd  0x10(%rsp),%xmm2
  movsd  0x30(%rsp),%xmm1
  mov    %rbp,%rdi
  movsd  %xmm0,(%rsp)
  movapd %xmm2,%xmm0
  call   <dd_two_sum@plt>
  movsd  0x8(%rsp),%xmm5
  movsd  0x10(%rbx),%xmm1
  mov    %rbp,%rdi
  addsd  0x30(%rsp),%xmm5
  movsd  %xmm0,0x18(%rsp)
  movsd  0x10(%r13),%xmm0
  movsd  %xmm5,0x10(%rsp)
  call   <dd_two_sum@plt>
  movsd  (%rsp),%xmm3
  movsd  0x30(%rsp),%xmm1
  mov    %rbp,%rdi
  movsd  %xmm0,0x8(%rsp)
  movapd %xmm3,%xmm0
  call   <dd_two_sum@plt>
  movsd  0x18(%rsp),%xmm4
  movsd  0x30(%rsp),%xmm1
  mov    %rbp,%rdi
  movsd  %xmm0,(%rsp)
  movapd %xmm4,%xmm0
  call   <dd_two_sum@plt>
  movsd  0x10(%rsp),%xmm5
  movsd  0x8(%rbx),%xmm1
  mov    %rbp,%rdi
  addsd  0x30(%rsp),%xmm5
  movsd  %xmm0,0x20(%rsp)
  movsd  0x8(%r13),%xmm0
  movsd  %xmm5,0x18(%rsp)
  call   <dd_two_sum@plt>
  movsd  0x8(%rsp),%xmm2
  movsd  0x30(%rsp),%xmm1
  mov    %rbp,%rdi
  movsd  %xmm0,0x10(%rsp)
  movapd %xmm2,%xmm0
  call   <dd_two_sum@plt>
  movsd  (%rsp),%xmm3
  movsd  0x30(%rsp),%xmm1
  mov    %rbp,%rdi
  movsd  %xmm0,0x8(%rsp)
  movapd %xmm3,%xmm0
  call   <dd_two_sum@plt>
  movsd  0x20(%rsp),%xmm4
  movsd  0x30(%rsp),%xmm1
  mov    %rbp,%rdi
  movsd  %xmm0,(%rsp)
  movapd %xmm4,%xmm0
  call   <dd_two_sum@plt>
  movsd  0x18(%rsp),%xmm5
  movsd  (%rbx),%xmm1
  mov    %rbp,%rdi
  addsd  0x30(%rsp),%xmm5
  movsd  %xmm0,0x28(%rsp)
  movsd  0x0(%r13),%xmm0
  movsd  %xmm5,0x20(%rsp)
  call   <dd_two_sum@plt>
  movsd  0x10(%rsp),%xmm7
  movsd  0x30(%rsp),%xmm1
  mov    %rbp,%rdi
  movsd  %xmm0,0x18(%rsp)
  movapd %xmm7,%xmm0
  call   <dd_two_sum@plt>
  movsd  0x8(%rsp),%xmm2
  movsd  0x30(%rsp),%xmm1
  mov    %rbp,%rdi
  movsd  %xmm0,0x10(%rsp)
  movapd %xmm2,%xmm0
  call   <dd_two_sum@plt>
  movsd  (%rsp),%xmm3
  movsd  0x30(%rsp),%xmm1
  mov    %rbp,%rdi
  movsd  %xmm0,0x8(%rsp)
  movapd %xmm3,%xmm0
  call   <dd_two_sum@plt>
  movsd  0x28(%rsp),%xmm4
  movsd  0x30(%rsp),%xmm1
  mov    %rbp,%rdi
  movsd  %xmm0,(%rsp)
  movapd %xmm4,%xmm0
  call   <dd_two_sum@plt>
  movsd  0x20(%rsp),%xmm5
  movsd  0x10(%rsp),%xmm7
  lea    0x18(%r12),%rcx
  movsd  0x18(%rsp),%xmm6
  addsd  0x30(%rsp),%xmm5
  movapd %xmm0,%xmm4
  lea    0x10(%r12),%rdx
  movsd  (%rsp),%xmm3
  movsd  0x8(%rsp),%xmm2
  lea    0x8(%r12),%rsi
  mov    %r12,%rdi
  lea    0x20(%r12),%r8
  movapd %xmm7,%xmm1
  movapd %xmm6,%xmm0
  call   <pd_fast_renorm@plt>
  mov    0x38(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L0
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
L0:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void pd_add(const double *a, const double *b, double *c) {
  double f0, f1, f2, f3, f4, f5, e;
  f5 = 0.0;
  f4 = dd_two_sum(a[4], b[4], &e);
  f5 = f5 + e;
  f3 = dd_two_sum(a[3], b[3], &e);
  f4 = dd_two_sum(f4, e, &e);
  f5 = f5 + e;
  f2 = dd_two_sum(a[2], b[2], &e);
  f3 = dd_two_sum(f3, e, &e);
  f4 = dd_two_sum(f4, e, &e);
  f5 = f5 + e;
  f1 = dd_two_sum(a[1], b[1], &e);
  f2 = dd_two_sum(f2, e, &e);
  f3 = dd_two_sum(f3, e, &e);
  f4 = dd_two_sum(f4, e, &e);
  f5 = f5 + e;
  f0 = dd_two_sum(a[0], b[0], &e);
  f1 = dd_two_sum(f1, e, &e);
  f2 = dd_two_sum(f2, e, &e);
  f3 = dd_two_sum(f3, e, &e);
  f4 = dd_two_sum(f4, e, &e);
  f5 = f5 + e;
  pd_fast_renorm(f0, f1, f2, f3, f4, f5, &c[0], &c[1], &c[2], &c[3], &c[4]);
}
```