user: 请从给定的汇编指令推导出C源代码:

```asm
<bindingradius>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x60,%rsp
  movsd  %xmm0,-0x38(%rbp)
  movsd  %xmm1,-0x40(%rbp)
  movsd  %xmm2,-0x48(%rbp)
  movsd  %xmm3,-0x50(%rbp)
  mov    %edi,-0x54(%rbp)
  pxor   %xmm0,%xmm0
  comisd -0x38(%rbp),%xmm0
  ja     L158
  pxor   %xmm0,%xmm0
  comisd -0x40(%rbp),%xmm0
  ja     L158
  pxor   %xmm0,%xmm0
  comisd -0x48(%rbp),%xmm0
  jb     L248
L158:
  movsd  D86(%rip),%xmm0
  jmp    L280
L248:
  pxor   %xmm0,%xmm0
  ucomisd -0x38(%rbp),%xmm0
  jp     L22
  pxor   %xmm0,%xmm0
  ucomisd -0x38(%rbp),%xmm0
  jne    L22
  pxor   %xmm0,%xmm0
  jmp    L280
L22:
  pxor   %xmm0,%xmm0
  ucomisd -0x40(%rbp),%xmm0
  jp     L186
  pxor   %xmm0,%xmm0
  ucomisd -0x40(%rbp),%xmm0
  jne    L186
  pxor   %xmm0,%xmm0
  comisd -0x50(%rbp),%xmm0
  jbe    L60
  movsd  -0x48(%rbp),%xmm1
  movsd  D95(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x38(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  jmp    L280
L60:
  cmpl   $0x0,-0x54(%rbp)
  je     L223
  movsd  -0x50(%rbp),%xmm0
  movsd  D5(%rip),%xmm1
  comisd %xmm1,%xmm0
  jbe    L223
  movsd  D5(%rip),%xmm0
  movapd %xmm0,%xmm1
  divsd  -0x50(%rbp),%xmm1
  movsd  D5(%rip),%xmm0
  subsd  %xmm1,%xmm0
  mulsd  -0x38(%rbp),%xmm0
  movsd  -0x48(%rbp),%xmm2
  movsd  D95(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  jmp    L280
L223:
  cmpl   $0x0,-0x54(%rbp)
  je     L42
  movsd  D5(%rip),%xmm0
  comisd -0x50(%rbp),%xmm0
  jb     L42
  movsd  D86(%rip),%xmm0
  jmp    L280
L42:
  movsd  -0x50(%rbp),%xmm0
  pxor   %xmm1,%xmm1
  comisd %xmm1,%xmm0
  jbe    L162
  movsd  -0x48(%rbp),%xmm1
  movsd  D95(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  -0x38(%rbp),%xmm0
  divsd  -0x50(%rbp),%xmm0
  addsd  %xmm0,%xmm1
  movsd  -0x38(%rbp),%xmm0
  divsd  %xmm1,%xmm0
  jmp    L280
L162:
  movsd  D86(%rip),%xmm0
  jmp    L280
L186:
  movsd  -0x48(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  mulsd  -0x40(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x8(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x8(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  jmp    L247
L208:
  movsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x20(%rbp),%xmm0
  addsd  %xmm0,%xmm0
  movsd  %xmm0,-0x20(%rbp)
L247:
  cmpl   $0x0,-0x54(%rbp)
  je     L80
  movsd  -0x50(%rbp),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  jmp    L83
L80:
  movsd  -0x50(%rbp),%xmm0
L83:
  movsd  -0x20(%rbp),%xmm1
  mov    -0x8(%rbp),%rax
  movapd %xmm0,%xmm2
  movq   %rax,%xmm0
  call   <numrxnrate@plt>
  movsd  -0x38(%rbp),%xmm1
  mulsd  -0x40(%rbp),%xmm1
  comisd %xmm0,%xmm1
  ja     L208
  movsd  -0x20(%rbp),%xmm0
  subsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movl   $0x0,-0x24(%rbp)
  jmp    L246
L87:
  movsd  -0x10(%rbp),%xmm1
  movsd  D93(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x18(%rbp),%xmm0
  addsd  -0x10(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  cmpl   $0x0,-0x54(%rbp)
  je     L250
  movsd  -0x50(%rbp),%xmm0
  mulsd  -0x20(%rbp),%xmm0
  jmp    L104
L250:
  movsd  -0x50(%rbp),%xmm0
L104:
  movsd  -0x20(%rbp),%xmm1
  mov    -0x8(%rbp),%rax
  movapd %xmm0,%xmm2
  movq   %rax,%xmm0
  call   <numrxnrate@plt>
  movsd  -0x38(%rbp),%xmm1
  mulsd  -0x40(%rbp),%xmm1
  comisd %xmm0,%xmm1
  jbe    L36
  movsd  -0x20(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
L36:
  addl   $0x1,-0x24(%rbp)
L246:
  cmpl   $0xe,-0x24(%rbp)
  jle    L87
  movsd  -0x10(%rbp),%xmm1
  movsd  D93(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x18(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x20(%rbp),%xmm0
L280:
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: disassemble_data
arguments: data_tag: D86
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D95
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D5
type_guessed: byte[8]
name: disassemble_data
arguments: data_tag: D93
type_guessed: f64

tool: {"D86": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0xbf]}

tool: {"D95": [0x18,0x2d,0x44,0x54,0xfb,0x21,0x29,0x40]}

tool: {"D5": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f]}

tool: D93: 0.5

assistant: ```c
double bindingradius(double rate, double dt, double difc, double b, int rel) {
  double a, lo, dif, step;
  int n;
  if (rate < 0 || dt < 0 || difc <= 0)
    return -1;
  if (rate == 0)
    return 0;
  if (dt == 0) {
    if (b < 0)
      return rate / (4 * 3.14159265358979323846 * difc);
    if (rel && b > 1)
      return rate * (1 - 1 / b) / (4 * 3.14159265358979323846 * difc);
    if (rel && b <= 1)
      return -1;
    if (b > 0)
      return rate / (4 * 3.14159265358979323846 * difc + rate / b);
    return -1;
  }
  step = sqrt(2.0 * difc * dt);
  lo = 0;
  a = step;
  while (numrxnrate(step, a, rel ? b * a : b) < rate * dt) {
    lo = a;
    a *= 2.0;
  }
  dif = a - lo;
  for (n = 0; n < 15; n++) {
    dif *= 0.5;
    a = lo + dif;
    if (numrxnrate(step, a, rel ? b * a : b) < rate * dt)
      lo = a;
  }
  a = lo + 0.5 * dif;
  return a;
}
```