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

```asm
<bessj>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x60,%rsp
  mov    %edi,-0x54(%rbp)
  movsd  %xmm0,-0x60(%rbp)
  cmpl   $0x1,-0x54(%rbp)
  jg     L111
  lea    D74(%rip),%rax
  mov    %rax,%rdi
  call   <nrerror@plt>
L111:
  movsd  -0x60(%rbp),%xmm0
  movq   D107(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  pxor   %xmm0,%xmm0
  ucomisd -0x10(%rbp),%xmm0
  jp     L72
  pxor   %xmm0,%xmm0
  ucomisd -0x10(%rbp),%xmm0
  jne    L72
  pxor   %xmm0,%xmm0
  jmp    L145
L72:
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x54(%rbp),%xmm1
  movsd  -0x10(%rbp),%xmm0
  comisd %xmm1,%xmm0
  jbe    L29
  movsd  D77(%rip),%xmm0
  divsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  mov    -0x10(%rbp),%rax
  movq   %rax,%xmm0
  call   <bessj0@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x30(%rbp)
  mov    -0x10(%rbp),%rax
  movq   %rax,%xmm0
  call   <bessj1@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x38(%rbp)
  movl   $0x1,-0x44(%rbp)
  jmp    L110
L109:
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x44(%rbp),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  mulsd  -0x38(%rbp),%xmm0
  subsd  -0x30(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  addl   $0x1,-0x44(%rbp)
L110:
  mov    -0x44(%rbp),%eax
  cmp    -0x54(%rbp),%eax
  jl     L109
  movsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  jmp    L78
L29:
  movsd  D77(%rip),%xmm0
  divsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x54(%rbp),%xmm1
  movsd  D43(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movq   %xmm1,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  mov    %eax,%edx
  mov    -0x54(%rbp),%eax
  add    %edx,%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %edx,%eax
  sar    %eax
  add    %eax,%eax
  mov    %eax,-0x3c(%rbp)
  movl   $0x0,-0x40(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  D55(%rip),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0x3c(%rbp),%eax
  mov    %eax,-0x44(%rbp)
  jmp    L96
L77:
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x44(%rbp),%xmm0
  mulsd  -0x8(%rbp),%xmm0
  mulsd  -0x38(%rbp),%xmm0
  subsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x30(%rbp)
  movsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x30(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  movq   D107(%rip),%xmm1
  andpd  %xmm1,%xmm0
  comisd D98(%rip),%xmm0
  jbe    L5
  movsd  -0x38(%rbp),%xmm1
  movsd  D69(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x28(%rbp),%xmm1
  movsd  D69(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x18(%rbp),%xmm1
  movsd  D69(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x20(%rbp),%xmm1
  movsd  D69(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
L5:
  cmpl   $0x0,-0x40(%rbp)
  je     L9
  movsd  -0x20(%rbp),%xmm0
  addsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
L9:
  cmpl   $0x0,-0x40(%rbp)
  sete   %al
  movzbl %al,%eax
  mov    %eax,-0x40(%rbp)
  mov    -0x44(%rbp),%eax
  cmp    -0x54(%rbp),%eax
  jne    L70
  movsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
L70:
  subl   $0x1,-0x44(%rbp)
L96:
  cmpl   $0x0,-0x44(%rbp)
  jg     L77
  movsd  -0x20(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  subsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x18(%rbp),%xmm0
  divsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
L78:
  pxor   %xmm0,%xmm0
  comisd -0x60(%rbp),%xmm0
  jbe    L120
  mov    -0x54(%rbp),%eax
  and    $0x1,%eax
  test   %eax,%eax
  je     L120
  movsd  -0x18(%rbp),%xmm0
  movq   D49(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  jmp    L145
L120:
  movsd  -0x18(%rbp),%xmm0
L145:
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: extract_symbolic_data
arguments: tag: D74
datatype: string
name: extract_symbolic_data
arguments: tag: D107
datatype: qword
name: extract_symbolic_data
arguments: tag: D77
datatype: f64
name: extract_symbolic_data
arguments: tag: D43
datatype: f64
name: extract_symbolic_data
arguments: tag: D55
datatype: f64
name: extract_symbolic_data
arguments: tag: D98
datatype: f64[2]
name: extract_symbolic_data
arguments: tag: D69
datatype: byte[8]
name: extract_symbolic_data
arguments: tag: D49
datatype: f32[2]

tool: D74: Index n less than 2 in bessj

tool: D107: 9223372036854775807

tool: D77: 2

tool: D43: 40

tool: D55: 1

tool: {"D98": [10000000000.0,1e-10]}

tool: {"D69": [0xbb,0xbd,0xd7,0xd9,0xdf,0x7c,0xdb,0x3d]}

tool: {"D49": [0.0,-0.0]}

assistant: ```c
double bessj(int n, double x) {
  double bessj0(double x);
  double bessj1(double x);
  void nrerror(char error_text[]);
  int j, jsum, m;
  double ax, bj, bjm, bjp, sum, tox, ans;
  if (n < 2)
    nrerror("Index n less than 2 in bessj");
  ax = fabs(x);
  if (ax == 0.0)
    return 0.0;
  else if (ax > (double)n) {
    tox = 2.0 / ax;
    bjm = bessj0(ax);
    bj = bessj1(ax);
    for (j = 1; j < n; j++) {
      bjp = j * tox * bj - bjm;
      bjm = bj;
      bj = bjp;
    }
    ans = bj;
  } else {
    tox = 2.0 / ax;
    m = 2 * ((n + (int)sqrt(40.0 * n)) / 2);
    jsum = 0;
    bjp = ans = sum = 0.0;
    bj = 1.0;
    for (j = m; j > 0; j--) {
      bjm = j * tox * bj - bjp;
      bjp = bj;
      bj = bjm;
      if (fabs(bj) > 1.0e10) {
        bj *= 1.0e-10;
        bjp *= 1.0e-10;
        ans *= 1.0e-10;
        sum *= 1.0e-10;
      }
      if (jsum)
        sum += bj;
      jsum = !jsum;
      if (j == n)
        ans = bjp;
    }
    sum = 2.0 * sum - bj;
    ans /= sum;
  }
  return x < 0.0 && (n & 1) ? -ans : ans;
}
```