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

```asm
<rtbis>:
  endbr64
  push   %r14
  movq   %xmm2,%r14
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  sub    $0x20,%rsp
  movsd  %xmm1,0x8(%rsp)
  movsd  %xmm0,0x18(%rsp)
  call   *%rdi
  movsd  %xmm0,0x10(%rsp)
  movsd  0x8(%rsp),%xmm0
  call   *%rbp
  mulsd  0x10(%rsp),%xmm0
  pxor   %xmm5,%xmm5
  movsd  0x18(%rsp),%xmm3
  comisd %xmm5,%xmm0
  jb     L52
  lea    D50(%rip),%rdi
  xor    %eax,%eax
  call   <printf@plt>
  lea    D13(%rip),%rdi
  xor    %eax,%eax
  call   <printf@plt>
  movsd  0x18(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm1
  lea    D1(%rip),%rdi
  mov    $0x2,%eax
  movapd %xmm3,%xmm0
  call   <printf@plt>
  mov    $0x1,%edi
  call   <exit@plt>
  movsd  0x18(%rsp),%xmm3
L52:
  pxor   %xmm6,%xmm6
  comisd 0x10(%rsp),%xmm6
  ja     L60
  movsd  0x8(%rsp),%xmm1
  movapd %xmm3,%xmm2
  subsd  %xmm1,%xmm2
L20:
  mov    $0x28,%ebx
  jmp    L15
  nopl   0x0(%rax,%rax,1)
L79:
  ucomisd %xmm5,%xmm0
  jp     L13
  je     L66
L13:
  sub    $0x1,%ebx
  je     L70
L15:
  mulsd  D23(%rip),%xmm2
  movsd  %xmm1,0x10(%rsp)
  movapd %xmm2,%xmm4
  movsd  %xmm2,0x18(%rsp)
  addsd  %xmm1,%xmm4
  movapd %xmm4,%xmm0
  movsd  %xmm4,0x8(%rsp)
  call   *%rbp
  pxor   %xmm5,%xmm5
  movsd  0x10(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm4
  movapd %xmm0,%xmm3
  movsd  0x18(%rsp),%xmm2
  movq   %r14,%xmm6
  cmpnlesd %xmm5,%xmm3
  andpd  %xmm3,%xmm1
  andnpd %xmm4,%xmm3
  orpd   %xmm3,%xmm1
  movapd %xmm2,%xmm3
  andpd  D64(%rip),%xmm3
  comisd %xmm3,%xmm6
  jbe    L79
L66:
  add    $0x20,%rsp
  movapd %xmm1,%xmm0
  pop    %rbx
  pop    %rbp
  pop    %r14
  ret
  nopl   0x0(%rax,%rax,1)
L60:
  movsd  0x8(%rsp),%xmm2
  movapd %xmm3,%xmm1
  subsd  %xmm3,%xmm2
  jmp    L20
  nopl   0x0(%rax,%rax,1)
L70:
  lea    D50(%rip),%rdi
  xor    %eax,%eax
  call   <printf@plt>
  lea    D41(%rip),%rdi
  call   <puts@plt>
  mov    $0x1,%edi
  call   <exit@plt>
  add    $0x20,%rsp
  pop    %rbx
  movapd %xmm1,%xmm0
  pop    %rbp
  pop    %r14
  ret
```

tool_calls: name: disassemble_data
arguments: tag: D50
type_infered: string
name: disassemble_data
arguments: tag: D13
type_infered: string
name: disassemble_data
arguments: tag: D1
type_infered: string
name: disassemble_data
arguments: tag: D23
type_infered: f64
name: disassemble_data
arguments: tag: D64
type_infered: qword[2]
name: disassemble_data
arguments: tag: D41
type_infered: string

tool: D50: 

Error in function rtbis():

tool: D13: 
root in function must be within

tool: D1: x1 = %f  and x2 = %f


tool: D23: 0.5

tool: {"D64": [0x7fffffffffffffff,0x0000000000000000]}

tool: D41: 
Too many bisections!!!

assistant: ```c
double rtbis(double (*func)(double), double x1, double x2, double xacc) {
  int j;
  double dx, f, fmid, xmid, rtb;
  f = (*func)(x1);
  fmid = (*func)(x2);
  if (f * fmid >= 0.0) {
    printf("\n\nError in function rtbis():");
    printf("\nroot in function must be within");
    printf("x1 = %f  and x2 = %f\n", x1, x2);
    exit(1);
  }
  rtb = f < 0.0 ? (dx = x2 - x1, x1) : (dx = x1 - x2, x2);
  for (j = 0; j < 40; j++) {
    fmid = (*func)(xmid = rtb + (dx *= 0.5));
    if (fmid <= 0.0)
      rtb = xmid;
    if (fabs(dx) < xacc || fmid == 0.0)
      return rtb;
  }
  printf("\n\nError in function rtbis():");
  printf("\nToo many bisections!!!\n");
  exit(1);
}
```