user: Converte this assembly snippet into its equivalent C code:

```asm
<gtpt>:
  endbr64
  push   %r12
  movapd %xmm2,%xmm8
  movapd %xmm0,%xmm10
  mov    %rdi,%r12
  push   %rbp
  movapd %xmm1,%xmm9
  mov    %rsi,%rbp
  movapd %xmm7,%xmm2
  push   %rbx
  mov    %rdx,%rbx
  sub    $0x70,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x68(%rsp)
  xor    %eax,%eax
  ucomisd D10(%rip),%xmm6
  jp     L62
  je     L46
L62:
  ucomisd D19(%rip),%xmm6
  jp     L41
  jne    L41
  movapd %xmm10,%xmm6
  movapd %xmm9,%xmm7
  lea    0x48(%rsp),%rdx
  movsd  %xmm10,0x28(%rsp)
  subsd  %xmm3,%xmm6
  movapd %xmm8,%xmm3
  subsd  %xmm4,%xmm7
  movsd  %xmm9,0x18(%rsp)
  subsd  %xmm5,%xmm3
  lea    0x60(%rsp),%rsi
  movsd  %xmm8,0x10(%rsp)
  lea    0x58(%rsp),%rdi
  movapd %xmm7,%xmm0
  movsd  %xmm6,0x20(%rsp)
  movapd %xmm3,%xmm1
  movsd  %xmm3,0x8(%rsp)
  movsd  %xmm7,(%rsp)
  call   <gt_agany@plt>
  movsd  0x48(%rsp),%xmm0
  lea    0x50(%rsp),%rdi
  call   <gt_fullrot@plt>
  movsd  0x8(%rsp),%xmm3
  movsd  (%rsp),%xmm7
  lea    0x40(%rsp),%rsi
  movsd  0x50(%rsp),%xmm0
  lea    0x38(%rsp),%rdi
  movapd %xmm3,%xmm2
  movapd %xmm7,%xmm1
  call   <rt@plt>
  movsd  0x20(%rsp),%xmm6
  movsd  0x38(%rsp),%xmm0
  movsd  (%rsp),%xmm7
  movsd  0x28(%rsp),%xmm10
  subsd  %xmm6,%xmm6
  movsd  0x18(%rsp),%xmm9
  movsd  0x8(%rsp),%xmm3
  subsd  %xmm7,%xmm0
  movsd  0x10(%rsp),%xmm8
  addsd  %xmm10,%xmm6
  addsd  %xmm9,%xmm0
  movsd  %xmm6,(%r12)
  movsd  %xmm0,0x0(%rbp)
  movsd  0x40(%rsp),%xmm0
  subsd  %xmm3,%xmm0
  addsd  %xmm8,%xmm0
  movsd  %xmm0,(%rbx)
L63:
  mov    0x68(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L67
  add    $0x70,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  ret
  nopw   0x0(%rax,%rax,1)
L41:
  mov    0x68(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L67
  add    $0x70,%rsp
  lea    D23(%rip),%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  jmp    10b0 <printf@plt>
  nopl   0x0(%rax,%rax,1)
L46:
  movapd %xmm0,%xmm7
  movapd %xmm1,%xmm6
  lea    0x48(%rsp),%rdx
  movsd  %xmm0,0x28(%rsp)
  subsd  %xmm3,%xmm7
  subsd  %xmm4,%xmm6
  movapd %xmm8,%xmm3
  movsd  %xmm1,0x20(%rsp)
  subsd  %xmm5,%xmm3
  lea    0x60(%rsp),%rsi
  movsd  %xmm8,0x18(%rsp)
  lea    0x58(%rsp),%rdi
  movapd %xmm6,%xmm1
  movapd %xmm7,%xmm0
  movsd  %xmm6,0x8(%rsp)
  movsd  %xmm3,0x10(%rsp)
  movsd  %xmm7,(%rsp)
  call   <gt_agany@plt>
  movsd  0x48(%rsp),%xmm0
  lea    0x50(%rsp),%rdi
  call   <gt_fullrot@plt>
  movsd  0x8(%rsp),%xmm6
  movsd  (%rsp),%xmm7
  lea    0x38(%rsp),%rsi
  movsd  0x50(%rsp),%xmm0
  lea    0x40(%rsp),%rdi
  movapd %xmm6,%xmm2
  movapd %xmm7,%xmm1
  call   <rt@plt>
  movsd  0x40(%rsp),%xmm0
  movsd  (%rsp),%xmm7
  movsd  0x28(%rsp),%xmm10
  movsd  0x8(%rsp),%xmm6
  subsd  %xmm7,%xmm0
  movsd  0x20(%rsp),%xmm9
  movsd  0x10(%rsp),%xmm3
  movsd  0x18(%rsp),%xmm8
  subsd  %xmm3,%xmm3
  addsd  %xmm10,%xmm0
  movsd  %xmm0,(%r12)
  movsd  0x38(%rsp),%xmm0
  subsd  %xmm6,%xmm0
  addsd  %xmm9,%xmm0
  movsd  %xmm0,0x0(%rbp)
  movapd %xmm3,%xmm0
  addsd  %xmm8,%xmm0
  movsd  %xmm0,(%rbx)
  jmp    L63
L67:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: label: D10
data_type: byte[8]
name: parse_data
arguments: label: D19
data_type: byte[8]
name: parse_data
arguments: label: D23
data_type: string

tool: {"D10": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

tool: {"D19": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f]}

tool: D23: rotate Error!


assistant: ```c
void gtpt(double x1, double y1, double z1, double x2, double y2, double z2, double *xr, double *yr, double *zr, double vh, double fnr) {
  double brrt_x, brrt_y, brrt_z;
  double agxy, agz, bsl, drl;
  double rot_x, rot_y, rot_z;
  double afrt_x, afrt_y, afrt_z;
  double ag_a, frt;
  double ag_vh, bsl_vh;
  if (vh == 0) {
    brrt_x = x1 - x2;
    brrt_y = y1 - y2;
    brrt_z = z1 - z2;
    gt_agany(brrt_x, brrt_y, &ag_vh, &bsl_vh, fnr, &ag_a);
    gt_fullrot(ag_a, &frt);
    rt(frt, brrt_x, brrt_y, &afrt_x, &afrt_y);
    afrt_z = brrt_z;
    *xr = x1 + (afrt_x - brrt_x);
    *yr = y1 + (afrt_y - brrt_y);
    *zr = z1 + (afrt_z - brrt_z);
  } else if (vh == 1) {
    brrt_x = x1 - x2;
    brrt_y = y1 - y2;
    brrt_z = z1 - z2;
    gt_agany(brrt_y, brrt_z, &ag_vh, &bsl_vh, fnr, &ag_a);
    gt_fullrot(ag_a, &frt);
    rt(frt, brrt_y, brrt_z, &afrt_y, &afrt_z);
    afrt_x = brrt_x;
    *xr = x1 + (afrt_x - brrt_x);
    *yr = y1 + (afrt_y - brrt_y);
    *zr = z1 + (afrt_z - brrt_z);
  } else {
    printf("rotate Error!\n");
  }
}
```