user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<GeogStep>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  add    $0xffffffffffffff80,%rsp
  movsd  %xmm0,-0x48(%rbp)
  movsd  %xmm1,-0x50(%rbp)
  movsd  %xmm2,-0x58(%rbp)
  movsd  %xmm3,-0x60(%rbp)
  mov    %rdi,-0x68(%rbp)
  mov    %rsi,-0x70(%rbp)
  mov    -0x60(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x38(%rbp)
  mov    -0x58(%rbp),%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x30(%rbp)
  mov    -0x58(%rbp),%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm0,-0x28(%rbp)
  movsd  -0x50(%rbp),%xmm0
  addsd  -0x60(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  -0x28(%rbp),%xmm2
  movsd  D58(%rip),%xmm1
  addsd  %xmm2,%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x78(%rbp)
  movsd  -0x50(%rbp),%xmm0
  movapd %xmm0,%xmm4
  subsd  -0x60(%rbp),%xmm4
  movq   %xmm4,%rax
  movq   %rax,%xmm0
  call   <sin@plt>
  movsd  D58(%rip),%xmm1
  subsd  -0x28(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  -0x78(%rbp),%xmm1
  addsd  %xmm0,%xmm1
  movsd  D15(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movsd  -0x20(%rbp),%xmm0
  movsd  D58(%rip),%xmm1
  comisd %xmm1,%xmm0
  jbe    L21
  mov    D71(%rip),%rax
  jmp    L1
L21:
  movsd  D17(%rip),%xmm0
  comisd -0x20(%rbp),%xmm0
  jbe    L39
  mov    D89(%rip),%rax
  jmp    L1
L39:
  mov    -0x20(%rbp),%rax
  movq   %rax,%xmm0
  call   <asin@plt>
  movq   %xmm0,%rax
L1:
  mov    -0x70(%rbp),%rdx
  mov    %rax,(%rdx)
  movsd  -0x38(%rbp),%xmm0
  mulsd  -0x30(%rbp),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x50(%rbp),%xmm0
  addsd  -0x60(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  -0x28(%rbp),%xmm2
  movsd  D58(%rip),%xmm0
  addsd  %xmm2,%xmm0
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,-0x78(%rbp)
  movsd  -0x50(%rbp),%xmm0
  subsd  -0x60(%rbp),%xmm0
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  call   <cos@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %eax,%xmm1
  movsd  D58(%rip),%xmm0
  subsd  -0x28(%rbp),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x78(%rbp),%xmm1
  addsd  %xmm0,%xmm1
  movsd  D15(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm0
  mov    -0x18(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <atan2@plt>
  movq   %xmm0,%rax
  mov    %rax,-0x8(%rbp)
  movsd  -0x48(%rbp),%xmm0
  addsd  -0x8(%rbp),%xmm0
  movq   %xmm0,%rax
  pxor   %xmm1,%xmm1
  movq   %rax,%xmm0
  call   <GeogLonR@plt>
  movq   %xmm0,%rax
  mov    -0x68(%rbp),%rdx
  mov    %rax,(%rdx)
  nop
  leave
  ret
```

tool_calls: name: disassemble_data
arguments: label: D58
data_type: f64
name: disassemble_data
arguments: label: D15
data_type: f64
name: disassemble_data
arguments: label: D71
data_type: f64
name: disassemble_data
arguments: label: D17
data_type: f32[2]
name: disassemble_data
arguments: label: D89
data_type: byte[8]

tool: D58: 1

tool: D15: 0.5

tool: D71: 1.5707963267948966

tool: {"D17": [0.0,-1.875]}

tool: {"D89": [0x18,0x2d,0x44,0x54,0xfb,0x21,0xf9,0xbf]}

assistant: ```c
void GeogStep(const double o1, const double a1, const double d, const double s, double *o2, double *a2) {
  double sin_s, sin_d, cos_d, dlon, a, x, y;
  sin_s = sin(s);
  sin_d = sin(d);
  cos_d = cos(d);
  a = 0.5 * (sin(a1 + s) * (1.0 + cos_d) + sin(a1 - s) * (1.0 - cos_d));
  *a2 = (a > 1.0) ? 1.570796326794896619231 : (a < -1.0) ? -1.570796326794896619231 : asin(a);
  y = sin_s * sin_d;
  x = 0.5 * (cos(a1 + s) * (1 + cos_d) + cos(a1 - s) * (1 - cos_d));
  dlon = atan2(y, x);
  *o2 = GeogLonR(o1 + dlon, 0.0);
}
```