user: 请将这段汇编代码片段转换成等效的C代码:

```asm
<pow>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x70,%rsp
  movsd  %xmm0,-0x68(%rbp)
  movsd  %xmm1,-0x70(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  lea    -0x68(%rbp),%rax
  mov    %rax,-0x38(%rbp)
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm1
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  jb     L30
  movsd  -0x68(%rbp),%xmm0
  movq   D109(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x68(%rbp)
L30:
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  lea    -0x50(%rbp),%rax
  mov    %rax,-0x30(%rbp)
  mov    -0x38(%rbp),%rax
  mov    0xc(%rax),%eax
  and    $0xffff7ff0,%eax
  mov    %eax,%edx
  mov    -0x30(%rbp),%rax
  mov    %edx,0xc(%rax)
  mov    -0x30(%rbp),%rax
  mov    0xc(%rax),%eax
  shr    $0x4,%eax
  sub    $0x3ff,%eax
  mov    %eax,-0x54(%rbp)
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm0
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm1
  subsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm1,%xmm0
  jp     L66
  pxor   %xmm1,%xmm1
  ucomisd %xmm1,%xmm0
  jne    L66
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x54(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  jmp    L120
L66:
  addl   $0x1,-0x54(%rbp)
  mov    -0x54(%rbp),%eax
  shl    $0x4,%eax
  mov    %eax,%edx
  mov    -0x38(%rbp),%rax
  mov    0xc(%rax),%eax
  mov    %edx,%ecx
  sub    %ecx,%eax
  mov    %eax,%edx
  mov    -0x38(%rbp),%rax
  mov    %edx,0xc(%rax)
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D71(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  mov    -0x38(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  D161(%rip),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x38(%rbp),%rax
  movsd  (%rax),%xmm2
  movsd  D161(%rip),%xmm1
  addsd  %xmm2,%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x28(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x20(%rbp),%xmm1
  movsd  D180(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  D191(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x20(%rbp),%xmm1
  movsd  D145(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x20(%rbp),%xmm1
  movsd  D196(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x20(%rbp),%xmm1
  movsd  D24(%rip),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x28(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  movsd  -0x50(%rbp),%xmm1
  movsd  D177(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x54(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  D20(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
L120:
  movsd  -0x50(%rbp),%xmm0
  mulsd  -0x70(%rbp),%xmm0
  movsd  %xmm0,-0x50(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  lea    -0x48(%rbp),%rax
  mov    %rax,-0x18(%rbp)
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm0
  comisd D214(%rip),%xmm0
  ja     L90
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D74(%rip),%xmm0
  comisd %xmm1,%xmm0
  jbe    L22
L90:
  movsd  D77(%rip),%xmm0
  jmp    L85
L22:
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,-0x54(%rbp)
  mov    -0x18(%rbp),%rax
  mov    0xc(%rax),%edx
  mov    -0x54(%rbp),%eax
  add    $0x3ff,%eax
  shl    $0x4,%eax
  add    %eax,%edx
  mov    -0x18(%rbp),%rax
  mov    %edx,0xc(%rax)
  mov    -0x18(%rbp),%rax
  mov    0xc(%rax),%eax
  and    $0xffff7ff0,%eax
  mov    %eax,%edx
  mov    -0x18(%rbp),%rax
  mov    %edx,0xc(%rax)
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0x54(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x30(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm1,%xmm0
  jp     L115
  pxor   %xmm1,%xmm1
  ucomisd %xmm1,%xmm0
  jne    L47
  mov    -0x18(%rbp),%rax
  movsd  (%rax),%xmm0
  jmp    L85
L47:
  jmp    L115
L31:
  mov    -0x18(%rbp),%rax
  mov    0xc(%rax),%eax
  lea    -0x10(%rax),%edx
  mov    -0x18(%rbp),%rax
  mov    %edx,0xc(%rax)
  mov    -0x18(%rbp),%rax
  mov    0xc(%rax),%eax
  and    $0xffff7ff0,%eax
  mov    %eax,%edx
  mov    -0x18(%rbp),%rax
  mov    %edx,0xc(%rax)
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  D161(%rip),%xmm0
  addsd  %xmm1,%xmm0
  mov    -0x30(%rbp),%rax
  movsd  %xmm0,(%rax)
L115:
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm1
  pxor   %xmm0,%xmm0
  comisd %xmm1,%xmm0
  ja     L31
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm0
  comisd D20(%rip),%xmm0
  jb     L109
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm0
  movsd  D20(%rip),%xmm1
  subsd  %xmm1,%xmm0
  mov    -0x30(%rbp),%rax
  movsd  %xmm0,(%rax)
  movsd  D71(%rip),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  jmp    L39
L109:
  movsd  D161(%rip),%xmm0
  movsd  %xmm0,-0x40(%rbp)
L39:
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm1
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm1
  movsd  -0x20(%rbp),%xmm2
  movsd  D105(%rip),%xmm0
  mulsd  %xmm0,%xmm2
  movsd  D53(%rip),%xmm0
  addsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x20(%rbp),%xmm1
  movsd  D79(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  mov    -0x18(%rbp),%rax
  movsd  (%rax),%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x40(%rbp),%xmm1
  movsd  -0x10(%rbp),%xmm0
  addsd  -0x28(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x10(%rbp),%xmm1
  subsd  -0x28(%rbp),%xmm1
  divsd  %xmm1,%xmm0
  mov    -0x30(%rbp),%rax
  movsd  %xmm0,(%rax)
  mov    -0x30(%rbp),%rax
  movsd  (%rax),%xmm0
L85:
  movq   %xmm0,%rax
  mov    -0x8(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L149
  call   <__stack_chk_fail@plt>
L149:
  movq   %rax,%xmm0
  leave
  ret
```

tool_calls: name: disassemble_data
arguments: data_label: D109
type: f32[2]
name: disassemble_data
arguments: data_label: D71
type: f64
name: disassemble_data
arguments: data_label: D161
type: f64
name: disassemble_data
arguments: data_label: D180
type: f64
name: disassemble_data
arguments: data_label: D191
type: f64
name: disassemble_data
arguments: data_label: D145
type: f64
name: disassemble_data
arguments: data_label: D196
type: f64
name: disassemble_data
arguments: data_label: D24
type: f64
name: disassemble_data
arguments: data_label: D177
type: f64
name: disassemble_data
arguments: data_label: D20
type: f64[2]
name: disassemble_data
arguments: data_label: D214
type: f64[2]
name: disassemble_data
arguments: data_label: D74
type: f32[2]
name: disassemble_data
arguments: data_label: D77
type: f64
name: disassemble_data
arguments: data_label: D105
type: f64
name: disassemble_data
arguments: data_label: D53
type: f64
name: disassemble_data
arguments: data_label: D79
type: f64

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

tool: D71: 1.4142135623730951

tool: D161: 1

tool: D180: 0.2376245609

tool: D191: 0.28525381498

tool: D145: 0.4000059794795

tool: D196: 0.666666635059382

tool: D24: 2.000000000046727

tool: D177: 1.442695040888634

tool: {"D20": [0.5,1023.0]}

tool: {"D214": [1023.0,-1023.0]}

tool: {"D74": [0.0,-4.4990234375]}

tool: D77: 1e+307

tool: D105: 0.0576900723731

tool: D53: 7.2152891511493

tool: D79: 20.8189237930062

assistant: ```c
double pow(double y, double x) {
  typedef union {
    double d;
    unsigned u[4];
  } DBL;
  double z, w, p, p2, q, re;
  unsigned ex;
  int iz;
  DBL *yp, *zp, *wp;
  yp = (DBL *)&y;
  if (yp->d <= 0.0)
    y = -y;
  z = 0.0;
  zp = (DBL *)&z;
  zp->u[3] = yp->u[3] & ~0100017;
  iz = (zp->u[3] >> 4) - 1023;
  if ((yp->d - zp->d) == 0.0)
    z = (double)iz;
  else {
    yp->u[3] -= ++iz << 4;
    yp->d *= 1.4142135623730950488;
    p = (yp->d - 1.0) / (yp->d + 1.0);
    p2 = p * p;
    z = p * (2.000000000046727 + p2 * (0.666666635059382 + p2 * (0.4000059794795 + p2 * (0.28525381498 + p2 * 0.2376245609))));
    z = z * 1.442695040888634 + (double)iz - 0.5;
  }
  z *= x;
  w = 0.0;
  wp = (DBL *)&w;
  if (zp->d > 1023.0 || zp->d < -1023.0)
    return (1E307);
  iz = (int)zp->d;
  wp->u[3] += iz + 1023 << 4;
  wp->u[3] &= ~0100017;
  if ((zp->d -= (double)iz) == 0.0)
    return (wp->d);
  while (zp->d < 0.0) {
    wp->u[3] -= 1 << 4;
    wp->u[3] &= ~0100017;
    zp->d++;
  }
  if (zp->d >= 0.5) {
    zp->d -= 0.5;
    re = 1.41421356237309504880;
  } else
    re = 1.0;
  p2 = zp->d * zp->d;
  p = zp->d * (7.2152891511493 + p2 * 0.0576900723731);
  q = 20.8189237930062 + p2;
  zp->d = re * wp->d * (q + p) / (q - p);
  return (zp->d);
}
```