user: What is the c source code of the assembly code below?

```asm
<int_sim_ad>:
  endbr64
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  movapd %xmm0,%xmm6
  movq   %xmm1,%rbx
  movsd  %xmm2,0x28(%rsp)
  mov    %rdi,%r12
  mov    %rsi,%r13
  movapd %xmm1,%xmm0
  subsd  %xmm6,%xmm0
  mulsd  D84(%rip),%xmm0
  movsd  %xmm0,0x8(%rsp)
  mov    %rsi,%rdi
  movsd  %xmm6,0x10(%rsp)
  movapd %xmm6,%xmm0
  call   *%r12
  movq   %xmm0,%rbp
  mov    %r13,%rdi
  movq   %rbx,%xmm0
  call   *%r12
  movq   %rbp,%xmm1
  addsd  %xmm0,%xmm1
  movsd  %xmm1,(%rsp)
  movsd  0x8(%rsp),%xmm7
  movapd %xmm7,%xmm6
  addsd  %xmm7,%xmm6
  movapd %xmm6,%xmm0
  addsd  0x10(%rsp),%xmm0
  mov    %r13,%rdi
  call   *%r12
  addsd  %xmm0,%xmm0
  addsd  (%rsp),%xmm0
  movq   %xmm0,%rbx
  movsd  0x10(%rsp),%xmm0
  addsd  0x8(%rsp),%xmm0
  mov    %r13,%rdi
  call   *%r12
  movapd %xmm0,%xmm3
  addsd  D97(%rip),%xmm3
  movsd  %xmm3,(%rsp)
  movsd  D5(%rip),%xmm5
  mulsd  0x8(%rsp),%xmm5
  movapd %xmm5,%xmm0
  addsd  0x10(%rsp),%xmm0
  mov    %r13,%rdi
  call   *%r12
  movsd  (%rsp),%xmm1
  addsd  %xmm0,%xmm1
  movq   %rbx,%xmm4
  divsd  D5(%rip),%xmm4
  movsd  %xmm4,0x18(%rsp)
  mulsd  D30(%rip),%xmm1
  movapd %xmm1,%xmm0
  addsd  %xmm1,%xmm0
  addsd  %xmm4,%xmm0
  mulsd  0x8(%rsp),%xmm0
  movsd  %xmm0,0x20(%rsp)
  movsd  0x28(%rsp),%xmm2
  movapd %xmm2,%xmm0
  mulsd  D91(%rip),%xmm0
  andpd  D45(%rip),%xmm0
  comisd %xmm2,%xmm0
  jbe    L15
  mov    $0x4,%ebp
  jmp    L23
L55:
  movsd  %xmm0,0x20(%rsp)
L23:
  add    %ebp,%ebp
  movsd  D0(%rip),%xmm6
  mulsd  0x8(%rsp),%xmm6
  movsd  %xmm6,0x8(%rsp)
  addsd  0x18(%rsp),%xmm1
  movsd  %xmm1,0x18(%rsp)
  cmp    $0x1,%ebp
  jle    L63
  mov    $0x1,%ebx
  pxor   %xmm7,%xmm7
  movsd  %xmm7,(%rsp)
L8:
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebx,%xmm0
  mulsd  0x8(%rsp),%xmm0
  addsd  0x10(%rsp),%xmm0
  mov    %r13,%rdi
  call   *%r12
  addsd  (%rsp),%xmm0
  movsd  %xmm0,(%rsp)
  add    $0x2,%ebx
  cmp    %ebx,%ebp
  jg     L8
L14:
  movsd  D30(%rip),%xmm1
  mulsd  (%rsp),%xmm1
  movapd %xmm1,%xmm0
  addsd  %xmm1,%xmm0
  addsd  0x18(%rsp),%xmm0
  mulsd  0x8(%rsp),%xmm0
  movapd %xmm0,%xmm2
  subsd  0x20(%rsp),%xmm2
  divsd  D50(%rip),%xmm2
  andpd  D45(%rip),%xmm2
  comisd 0x28(%rsp),%xmm2
  ja     L55
L6:
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
L63:
  pxor   %xmm2,%xmm2
  movsd  %xmm2,(%rsp)
  jmp    L14
L15:
  movsd  0x20(%rsp),%xmm0
  jmp    L6
```

tool_calls: name: parse_data
arguments: tag: D84
type_guessed: f32[2]
name: parse_data
arguments: tag: D97
type_guessed: f32[2]
name: parse_data
arguments: tag: D5
type_guessed: f64
name: parse_data
arguments: tag: D30
type_guessed: byte[8]
name: parse_data
arguments: tag: D91
type_guessed: f64
name: parse_data
arguments: tag: D45
type_guessed: qword[2]
name: parse_data
arguments: tag: D0
type_guessed: f64
name: parse_data
arguments: tag: D50
type_guessed: f64

tool: {"D84": [0.0,1.625]}

tool: {"D97": [0.0,0.0]}

tool: D5: 3

tool: {"D30": [0x55,0x55,0x55,0x55,0x55,0x55,0xe5,0x3f]}

tool: D91: 10000000000

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

tool: D0: 0.5

tool: D50: 15

assistant: ```c
double int_sim_ad(double a, double b, double tol, double (*f)(double, void *), void *args) {
  int N1 = 4;
  double S = 0.0;
  double T = 0.0;
  double I = 0.0;
  double I1 = 0.0;
  int i, n;
  int N = N1;
  double h = (b - a) / N;
  S = f(a, args) + f(b, args);
  for (i = 2; i < N; i += 2)
    S += 2 * f(a + i * h, args);
  for (i = 1; i < N; i += 2)
    T += f(a + i * h, args);
  S /= 3.0;
  T *= 2.0 / 3.0;
  I1 = h * (S + 2 * T);
  double err = 1.0e10 * tol;
  while (fabs(err) > tol) {
    N *= 2;
    h *= 0.5;
    I = I1;
    S += T;
    T = 0.0;
    for (i = 1; i < N; i += 2)
      T += f(a + i * h, args);
    T *= 2.0 / 3.0;
    I1 = h * (S + 2 * T);
    err = (I1 - I) / 15.0;
  }
  return I1;
}
```