user: Translate the below assembly code into C programming language:

```asm
<spread>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  movsd  %xmm0,-0x28(%rbp)
  mov    %rdi,-0x30(%rbp)
  mov    %esi,-0x34(%rbp)
  movsd  %xmm1,-0x40(%rbp)
  mov    %edx,-0x38(%rbp)
  cmpl   $0xa,-0x38(%rbp)
  jle    L32
  mov    $0x1,%eax
  jmp    L71
L32:
  movsd  -0x40(%rbp),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,-0x14(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x14(%rbp),%xmm0
  ucomisd -0x40(%rbp),%xmm0
  jp     L108
  ucomisd -0x40(%rbp),%xmm0
  jne    L108
  mov    -0x14(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x14(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  addsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,(%rax)
  jmp    L133
L108:
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x38(%rbp),%xmm1
  movsd  D41(%rip),%xmm0
  movapd %xmm1,%xmm2
  mulsd  %xmm0,%xmm2
  movsd  -0x40(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  D119(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  D119(%rip),%xmm1
  comisd %xmm1,%xmm0
  jbe    L12
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x38(%rbp),%xmm1
  movsd  D41(%rip),%xmm0
  movapd %xmm1,%xmm2
  mulsd  %xmm0,%xmm2
  movsd  -0x40(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  D119(%rip),%xmm0
  addsd  %xmm1,%xmm0
  jmp    L33
L12:
  movsd  D119(%rip),%xmm0
L33:
  mov    -0x34(%rbp),%eax
  sub    -0x38(%rbp),%eax
  add    $0x1,%eax
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  comisd %xmm0,%xmm1
  jbe    L143
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x38(%rbp),%xmm1
  movsd  D41(%rip),%xmm0
  movapd %xmm1,%xmm2
  mulsd  %xmm0,%xmm2
  movsd  -0x40(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  D119(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  D119(%rip),%xmm1
  comisd %xmm1,%xmm0
  jbe    L58
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x38(%rbp),%xmm1
  movsd  D41(%rip),%xmm0
  movapd %xmm1,%xmm2
  mulsd  %xmm0,%xmm2
  movsd  -0x40(%rbp),%xmm0
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  movsd  D119(%rip),%xmm0
  addsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%eax
  jmp    L24
L58:
  mov    $0x1,%eax
  jmp    L24
L143:
  mov    -0x34(%rbp),%eax
  sub    -0x38(%rbp),%eax
  add    $0x1,%eax
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  cvttsd2si %xmm0,%eax
L24:
  mov    %eax,-0x10(%rbp)
  mov    -0x10(%rbp),%edx
  mov    -0x38(%rbp),%eax
  add    %edx,%eax
  sub    $0x1,%eax
  mov    %eax,-0xc(%rbp)
  mov    -0x38(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  lea    D78(%rip),%rax
  mov    (%rdx,%rax,1),%eax
  mov    %eax,-0x18(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x10(%rbp),%xmm1
  movsd  -0x40(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0x10(%rbp),%eax
  add    $0x1,%eax
  mov    %eax,-0x1c(%rbp)
  jmp    L120
L16:
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x1c(%rbp),%xmm1
  movsd  -0x40(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movsd  -0x8(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  addl   $0x1,-0x1c(%rbp)
L120:
  mov    -0x1c(%rbp),%eax
  cmp    -0xc(%rbp),%eax
  jle    L16
  mov    -0xc(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  pxor   %xmm3,%xmm3
  cvtsi2sdl -0x18(%rbp),%xmm3
  pxor   %xmm4,%xmm4
  cvtsi2sdl -0xc(%rbp),%xmm4
  movsd  -0x40(%rbp),%xmm1
  subsd  %xmm4,%xmm1
  mulsd  %xmm3,%xmm1
  divsd  %xmm1,%xmm0
  mov    -0xc(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm2,%xmm0
  movsd  %xmm0,(%rax)
  mov    -0xc(%rbp),%eax
  sub    $0x1,%eax
  mov    %eax,-0x1c(%rbp)
  jmp    L89
L19:
  mov    -0x1c(%rbp),%eax
  add    $0x1,%eax
  sub    -0x10(%rbp),%eax
  mov    %eax,%ecx
  mov    -0x18(%rbp),%eax
  cltd
  idiv   %ecx
  mov    %eax,%edx
  mov    -0x1c(%rbp),%eax
  sub    -0xc(%rbp),%eax
  imul   %edx,%eax
  mov    %eax,-0x18(%rbp)
  mov    -0x1c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm2
  movsd  -0x28(%rbp),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  pxor   %xmm3,%xmm3
  cvtsi2sdl -0x18(%rbp),%xmm3
  pxor   %xmm4,%xmm4
  cvtsi2sdl -0x1c(%rbp),%xmm4
  movsd  -0x40(%rbp),%xmm1
  subsd  %xmm4,%xmm1
  mulsd  %xmm3,%xmm1
  divsd  %xmm1,%xmm0
  mov    -0x1c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x8(%rax),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  addsd  %xmm2,%xmm0
  movsd  %xmm0,(%rax)
  subl   $0x1,-0x1c(%rbp)
L89:
  mov    -0x1c(%rbp),%eax
  cmp    -0x10(%rbp),%eax
  jge    L19
L133:
  mov    $0x0,%eax
L71:
  pop    %rbp
  ret
```

tool_calls: name: retrieve_data
arguments: data_tag: D41
datatype: f64
name: retrieve_data
arguments: data_tag: D119
datatype: f64

tool: D41: 0.5

tool: D119: 1

assistant: ```c
int spread(double y, double *yy, int n, double x, int m) {
  int ihi, ilo, ix, j, nden;
  static int nfac[11] = {0, 1, 1, 2, 6, 24, 120, 720, 5040, 40320, 362880};
  double fac;
  if (m > 10)
    return (1);
  ix = (int)x;
  if (x == (double)ix)
    yy[ix - 1] += y;
  else {
    ilo = (((((x - 0.5 * m + 1.0)) > (1) ? ((x - 0.5 * m + 1.0)) : (1))) < (n - m + 1) ? ((((x - 0.5 * m + 1.0)) > (1) ? ((x - 0.5 * m + 1.0)) : (1))) : (n - m + 1));
    ihi = ilo + m - 1;
    nden = nfac[m];
    fac = x - ilo;
    for (j = ilo + 1; j <= ihi; j++)
      fac *= (x - j);
    yy[ihi - 1] += y * fac / (nden * (x - ihi));
    for (j = ihi - 1; j >= ilo; j--) {
      nden = (nden / (j + 1 - ilo)) * (j - ihi);
      yy[j - 1] += y * fac / (nden * (x - j));
    }
  }
  return (0);
}
```