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

```asm
<algdiv>:
  endbr64
  movapd %xmm0,%xmm7
  sub    $0x38,%rsp
  movsd  D17(%rip),%xmm6
  divsd  %xmm1,%xmm0
  comisd %xmm7,%xmm1
  jb     L27
  movapd %xmm0,%xmm2
  movapd %xmm0,%xmm12
  movapd %xmm6,%xmm4
  addsd  %xmm6,%xmm2
  divsd  %xmm2,%xmm12
  divsd  %xmm2,%xmm4
  movapd %xmm7,%xmm2
  subsd  D62(%rip),%xmm2
  addsd  %xmm1,%xmm2
  movsd  %xmm2,0x8(%rsp)
L10:
  movapd %xmm4,%xmm2
  movapd %xmm4,%xmm9
  movapd %xmm6,%xmm8
  movsd  %xmm7,0x28(%rsp)
  mulsd  %xmm4,%xmm2
  movsd  %xmm6,0x20(%rsp)
  movsd  %xmm1,0x18(%rsp)
  addsd  %xmm2,%xmm9
  movapd %xmm2,%xmm10
  movapd %xmm2,%xmm3
  movapd %xmm2,%xmm11
  movapd %xmm2,%xmm5
  addsd  %xmm6,%xmm9
  mulsd  %xmm9,%xmm10
  divsd  %xmm1,%xmm8
  mulsd  D61(%rip),%xmm9
  addsd  %xmm4,%xmm10
  addsd  %xmm6,%xmm10
  mulsd  %xmm10,%xmm3
  mulsd  D65(%rip),%xmm10
  addsd  %xmm4,%xmm3
  addsd  %xmm6,%xmm3
  mulsd  %xmm3,%xmm11
  mulsd  D56(%rip),%xmm3
  addsd  %xmm4,%xmm11
  mulsd  %xmm8,%xmm8
  addsd  %xmm6,%xmm11
  mulsd  %xmm11,%xmm5
  mulsd  D5(%rip),%xmm11
  addsd  %xmm5,%xmm4
  addsd  %xmm6,%xmm4
  mulsd  D53(%rip),%xmm4
  mulsd  %xmm8,%xmm4
  addsd  %xmm11,%xmm4
  mulsd  %xmm8,%xmm4
  addsd  %xmm4,%xmm3
  mulsd  %xmm8,%xmm3
  addsd  %xmm10,%xmm3
  mulsd  %xmm8,%xmm3
  movapd %xmm3,%xmm2
  movapd %xmm12,%xmm3
  divsd  %xmm1,%xmm3
  addsd  %xmm9,%xmm2
  mulsd  %xmm8,%xmm2
  addsd  D42(%rip),%xmm2
  mulsd  %xmm2,%xmm3
  movsd  %xmm3,0x10(%rsp)
  call   <alnrel@plt>
  movsd  0x8(%rsp),%xmm2
  movsd  0x18(%rsp),%xmm1
  mulsd  %xmm0,%xmm2
  movapd %xmm1,%xmm0
  movsd  %xmm2,0x8(%rsp)
  call   <log@plt>
  movsd  0x20(%rsp),%xmm6
  movsd  0x28(%rsp),%xmm7
  movsd  0x8(%rsp),%xmm2
  movsd  0x10(%rsp),%xmm3
  subsd  %xmm6,%xmm0
  mulsd  %xmm7,%xmm0
  comisd %xmm2,%xmm0
  jae    L3
  subsd  %xmm0,%xmm3
  add    $0x38,%rsp
  subsd  %xmm2,%xmm3
  movapd %xmm3,%xmm0
  ret
  nopl   (%rax)
L27:
  movapd %xmm1,%xmm4
  movapd %xmm6,%xmm12
  divsd  %xmm7,%xmm4
  movapd %xmm4,%xmm2
  addsd  %xmm6,%xmm2
  divsd  %xmm2,%xmm12
  divsd  %xmm2,%xmm4
  movapd %xmm1,%xmm2
  subsd  D62(%rip),%xmm2
  addsd  %xmm7,%xmm2
  movsd  %xmm2,0x8(%rsp)
  jmp    L10
  nopl   0x0(%rax)
L3:
  subsd  %xmm2,%xmm3
  add    $0x38,%rsp
  subsd  %xmm0,%xmm3
  movapd %xmm3,%xmm0
  ret
```

tool_calls: name: parse_data
arguments: data_label: D17
data_type: f64
name: parse_data
arguments: data_label: D62
data_type: f64
name: parse_data
arguments: data_label: D61
data_type: byte[8]
name: parse_data
arguments: data_label: D65
data_type: byte[8]
name: parse_data
arguments: data_label: D56
data_type: byte[8]
name: parse_data
arguments: data_label: D5
data_type: byte[8]
name: parse_data
arguments: data_label: D53
data_type: byte[8]
name: parse_data
arguments: data_label: D42
data_type: byte[8]

tool: D17: 1

tool: D62: 0.5

tool: {"D61": [0x12,0x84,0xbb,0x16,0x6c,0xc1,0x66,0xbf]}

tool: {"D65": [0x60,0x89,0x48,0xd4,0x9f,0x01,0x4a,0x3f]}

tool: {"D56": [0xdf,0x30,0xf4,0x8f,0xec,0x80,0x43,0xbf]}

tool: {"D5": [0x6b,0x4a,0xb9,0x4f,0xd9,0x6f,0x4b,0x3f]}

tool: {"D53": [0x38,0x18,0xca,0xcb,0x25,0x16,0x5b,0xbf]}

tool: {"D42": [0x53,0x55,0x55,0x55,0x55,0x55,0xb5,0x3f]}

assistant: ```c
double algdiv(double a, double b) {
  static double c0 = .0833333333333333;
  static double c1 = -.00277777777760991;
  static double c2 = 7.9365066682539e-4;
  static double c3 = -5.9520293135187e-4;
  static double c4 = 8.37308034031215e-4;
  static double c5 = -.00165322962780713;
  double ret_val, r__1;
  double c__, d__, h__, t, u, v, w, x, s3, s5, x2, s7, s9, s11;
  if (a <= b) {
    goto L10;
  }
  h__ = b / a;
  c__ = 1.0 / (h__ + 1.0);
  x = h__ / (h__ + 1.0);
  d__ = a + (b - 0.5);
  goto L20;
L10:
  h__ = a / b;
  c__ = h__ / (h__ + 1.0);
  x = 1.0 / (h__ + 1.0);
  d__ = b + (a - 0.5);
L20:
  x2 = x * x;
  s3 = x + x2 + 1.0;
  s5 = x + x2 * s3 + 1.0;
  s7 = x + x2 * s5 + 1.0;
  s9 = x + x2 * s7 + 1.0;
  s11 = x + x2 * s9 + 1.0;
  r__1 = 1.0 / b;
  t = r__1 * r__1;
  w = ((((c5 * s11 * t + c4 * s9) * t + c3 * s7) * t + c2 * s5) * t + c1 * s3) * t + c0;
  w *= c__ / b;
  r__1 = a / b;
  u = d__ * alnrel(r__1);
  v = a * (log(b) - 1.0);
  if (u <= v) {
    goto L30;
  }
  ret_val = w - v - u;
  return ret_val;
L30:
  ret_val = w - u - v;
  return ret_val;
}
```