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

```asm
<istanpt>:
  endbr64
  push   %r12
  movapd %xmm0,%xmm8
  movapd %xmm2,%xmm9
  movapd %xmm4,%xmm2
  push   %rbp
  movapd %xmm6,%xmm4
  sub    $0x88,%rsp
  andpd  D7(%rip),%xmm0
  movsd  %xmm3,0x18(%rsp)
  movapd %xmm5,%xmm3
  movapd %xmm7,%xmm5
  comisd %xmm9,%xmm0
  ja     L85
  pxor   %xmm7,%xmm7
  movsd  %xmm7,0x8(%rsp)
L78:
  comisd %xmm1,%xmm7
  movapd %xmm8,%xmm0
  movsd  %xmm5,0x48(%rsp)
  divsd  %xmm9,%xmm0
  movsd  %xmm4,0x40(%rsp)
  movsd  %xmm3,0x38(%rsp)
  movsd  %xmm2,0x30(%rsp)
  movsd  %xmm9,0x28(%rsp)
  movsd  %xmm1,0x20(%rsp)
  movsd  %xmm8,0x10(%rsp)
  ja     L27
  call   <acos@plt>
  movsd  0x48(%rsp),%xmm5
  movsd  0x40(%rsp),%xmm4
  movsd  0x38(%rsp),%xmm3
  movsd  0x30(%rsp),%xmm2
  movapd %xmm0,%xmm10
  movsd  0x28(%rsp),%xmm9
  movsd  0x20(%rsp),%xmm1
  movsd  0x10(%rsp),%xmm8
L24:
  mulsd  %xmm8,%xmm8
  movsd  0x8(%rsp),%xmm7
  mulsd  %xmm1,%xmm1
  addsd  %xmm1,%xmm8
  ucomisd %xmm8,%xmm7
  ja     L52
  sqrtsd %xmm8,%xmm8
L34:
  movsd  D33(%rip),%xmm0
  movsd  D26(%rip),%xmm6
  comisd %xmm8,%xmm0
  jbe    L59
L19:
  movapd %xmm10,%xmm0
  lea    0x78(%rsp),%rbp
  lea    0x70(%rsp),%r12
  movsd  %xmm5,0x68(%rsp)
  addsd  %xmm6,%xmm0
  mov    %r12,%rsi
  mov    %rbp,%rdi
  movsd  %xmm4,0x60(%rsp)
  movsd  %xmm3,0x58(%rsp)
  movsd  %xmm2,0x50(%rsp)
  movsd  %xmm9,0x30(%rsp)
  movsd  %xmm10,0x28(%rsp)
  movsd  %xmm6,0x20(%rsp)
  call   <sincos@plt>
  movsd  0x20(%rsp),%xmm6
  mov    %r12,%rsi
  movsd  0x28(%rsp),%xmm10
  movsd  0x70(%rsp),%xmm7
  movsd  0x78(%rsp),%xmm1
  mov    %rbp,%rdi
  subsd  %xmm6,%xmm10
  movsd  %xmm7,0x10(%rsp)
  movsd  %xmm1,0x48(%rsp)
  movapd %xmm10,%xmm0
  call   <sincos@plt>
  movsd  0x18(%rsp),%xmm1
  movsd  0x78(%rsp),%xmm10
  movsd  0x30(%rsp),%xmm9
  movsd  0x70(%rsp),%xmm8
  mulsd  %xmm1,%xmm10
  movapd %xmm1,%xmm11
  movsd  0x10(%rsp),%xmm0
  movsd  0x48(%rsp),%xmm1
  mulsd  %xmm9,%xmm8
  movsd  0x68(%rsp),%xmm5
  movsd  0x60(%rsp),%xmm4
  mulsd  %xmm11,%xmm1
  movsd  0x58(%rsp),%xmm3
  movsd  0x50(%rsp),%xmm2
  mulsd  %xmm9,%xmm0
  movsd  %xmm5,0x30(%rsp)
  movsd  0xa8(%rsp),%xmm7
  movsd  0xa0(%rsp),%xmm6
  movsd  %xmm4,0x28(%rsp)
  movsd  %xmm10,0x38(%rsp)
  movsd  %xmm8,0x40(%rsp)
  movsd  %xmm3,0x20(%rsp)
  movsd  %xmm2,0x18(%rsp)
  call   <ellipse2tr@plt>
  movsd  0x38(%rsp),%xmm10
  movsd  0x40(%rsp),%xmm8
  movsd  0xa8(%rsp),%xmm7
  movsd  0x30(%rsp),%xmm5
  movsd  %xmm0,0x10(%rsp)
  movsd  0x28(%rsp),%xmm4
  movapd %xmm10,%xmm1
  movapd %xmm8,%xmm0
  movsd  0xa0(%rsp),%xmm6
  movsd  0x20(%rsp),%xmm3
  movsd  0x18(%rsp),%xmm2
  call   <ellipse2tr@plt>
  mulsd  0x10(%rsp),%xmm0
  xor    %eax,%eax
  comisd 0x8(%rsp),%xmm0
  seta   %al
  add    $0x88,%rsp
  add    $0x71,%eax
  pop    %rbp
  pop    %r12
  ret
  nopl   (%rax)
L85:
  pxor   %xmm6,%xmm6
  comisd %xmm8,%xmm6
  movsd  %xmm6,0x8(%rsp)
  movapd %xmm9,%xmm8
  movapd %xmm6,%xmm7
  jbe    L78
  xorpd  D6(%rip),%xmm8
  jmp    L78
  nopl   0x0(%rax)
L27:
  call   <acos@plt>
  movsd  0x10(%rsp),%xmm8
  movsd  D12(%rip),%xmm10
  movsd  0x20(%rsp),%xmm1
  movsd  0x28(%rsp),%xmm9
  subsd  %xmm0,%xmm10
  movsd  0x30(%rsp),%xmm2
  movsd  0x38(%rsp),%xmm3
  movsd  0x40(%rsp),%xmm4
  movsd  0x48(%rsp),%xmm5
  jmp    L24
  nopl   0x0(%rax)
L59:
  divsd  %xmm8,%xmm6
  movsd  %xmm5,0x40(%rsp)
  movsd  %xmm4,0x38(%rsp)
  movsd  %xmm3,0x30(%rsp)
  movsd  %xmm2,0x28(%rsp)
  movsd  %xmm9,0x20(%rsp)
  movsd  %xmm10,0x10(%rsp)
  movapd %xmm6,%xmm0
  call   <asin@plt>
  movsd  0x40(%rsp),%xmm5
  movsd  0x38(%rsp),%xmm4
  movsd  0x30(%rsp),%xmm3
  movsd  0x28(%rsp),%xmm2
  movapd %xmm0,%xmm6
  movsd  0x20(%rsp),%xmm9
  movsd  0x10(%rsp),%xmm10
  jmp    L19
L52:
  movapd %xmm8,%xmm0
  movsd  %xmm5,0x40(%rsp)
  movsd  %xmm4,0x38(%rsp)
  movsd  %xmm3,0x30(%rsp)
  movsd  %xmm2,0x28(%rsp)
  movsd  %xmm9,0x20(%rsp)
  movsd  %xmm10,0x10(%rsp)
  call   <sqrt@plt>
  movsd  0x40(%rsp),%xmm5
  movsd  0x38(%rsp),%xmm4
  movsd  0x30(%rsp),%xmm3
  movsd  0x28(%rsp),%xmm2
  movapd %xmm0,%xmm8
  movsd  0x20(%rsp),%xmm9
  movsd  0x10(%rsp),%xmm10
  jmp    L34
```

tool_calls: name: unpack_data
arguments: label: D7
datatype: qword[2]
name: unpack_data
arguments: label: D33
datatype: byte[8]
name: unpack_data
arguments: label: D26
datatype: byte[8]
name: unpack_data
arguments: label: D6
datatype: f32[4]
name: unpack_data
arguments: label: D12
datatype: byte[8]

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

tool: {"D33": [0x2c,0x43,0x1c,0xeb,0xe2,0x36,0x1a,0x3f]}

tool: {"D26": [0x8d,0xed,0xb5,0xa0,0xf7,0xc6,0xc0,0x3e]}

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

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

assistant: ```c
int istanpt(double x, double y, double A1, double B1, double AA, double BB, double CC, double DD, double EE, double FF) {
  double x1, y1, x2, y2, theta, test1, test2, branch, eps_radian;
  if (fabs(x) > A1)
    x = (x < 0) ? -A1 : A1;
  if (y < 0.0)
    theta = (2.0 * (2.0 * asin(1.0))) - acos(x / A1);
  else
    theta = acos(x / A1);
  branch = sqrt(x * x + y * y);
  if (branch < 100.0 * +1.0E-06)
    eps_radian = 2.0 * +1.0E-06;
  else
    eps_radian = asin(2.0 * +1.0E-06 / branch);
  x1 = A1 * cos(theta + eps_radian);
  y1 = B1 * sin(theta + eps_radian);
  x2 = A1 * cos(theta - eps_radian);
  y2 = B1 * sin(theta - eps_radian);
  test1 = ellipse2tr(x1, y1, AA, BB, CC, DD, EE, FF);
  test2 = ellipse2tr(x2, y2, AA, BB, CC, DD, EE, FF);
  if ((test1 * test2) > 0.0)
    return 114;
  else
    return 113;
}
```