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

```asm
<SubGussFs>:
  endbr64
  push   %r13
  push   %r12
  mov    %rdi,%r12
  push   %rbp
  push   %rbx
  mov    %rsi,%rbx
  sub    $0x48,%rsp
  test   %edx,%edx
  jne    L16
  cvttsd2si %xmm1,%eax
  movslq %eax,%rbp
  sar    $0x1f,%eax
  imul   $0x2aaaaaab,%rbp,%rbp
  shr    $0x20,%rbp
  sub    %eax,%ebp
  add    $0x1,%ebp
  lea    0x0(%rbp,%rbp,2),%eax
  lea    -0x3(%rax,%rax,1),%edx
L6:
  movsd  D27(%rip),%xmm3
  movsd  D131(%rip),%xmm2
  divsd  %xmm3,%xmm0
  movapd %xmm0,%xmm6
  pxor   %xmm0,%xmm0
  cvtsi2sd %edx,%xmm0
  mulsd  %xmm2,%xmm6
  subsd  %xmm0,%xmm1
  cvttsd2si %xmm6,%edi
  movsd  %xmm6,(%rsp)
  divsd  %xmm3,%xmm1
  mulsd  %xmm2,%xmm1
  movsd  %xmm1,0x38(%rsp)
  call   <sin@plt>
  movsd  %xmm0,0x8(%rsp)
  movsd  (%rsp),%xmm0
  call   <cos@plt>
  movsd  0x8(%rsp),%xmm12
  movsd  D248(%rip),%xmm7
  movsd  %xmm0,0x30(%rsp)
  movsd  D149(%rip),%xmm2
  mulsd  %xmm12,%xmm7
  movsd  %xmm12,0x28(%rsp)
  movapd %xmm7,%xmm0
  mulsd  %xmm12,%xmm0
  subsd  %xmm0,%xmm2
  cvttsd2si %xmm2,%r13d
  movsd  %xmm2,0x8(%rsp)
  mov    %r13d,%edi
  call   <sqrt@plt>
  movsd  0x8(%rsp),%xmm2
  movsd  D129(%rip),%xmm3
  mov    %r13d,%edi
  divsd  %xmm2,%xmm3
  divsd  %xmm0,%xmm3
  movsd  %xmm3,0x8(%rsp)
  call   <sqrt@plt>
  movsd  D118(%rip),%xmm5
  movsd  (%rsp),%xmm6
  divsd  %xmm0,%xmm5
  movapd %xmm6,%xmm0
  addsd  %xmm6,%xmm0
  cvttsd2si %xmm0,%edi
  movapd %xmm5,%xmm7
  divsd  0x8(%rsp),%xmm7
  movsd  %xmm5,0x20(%rsp)
  subsd  D149(%rip),%xmm7
  movsd  %xmm7,0x18(%rsp)
  call   <sin@plt>
  movsd  (%rsp),%xmm6
  movsd  D215(%rip),%xmm2
  mulsd  D210(%rip),%xmm0
  mulsd  %xmm6,%xmm2
  subsd  %xmm0,%xmm2
  movsd  D197(%rip),%xmm0
  mulsd  %xmm6,%xmm0
  movsd  %xmm2,0x8(%rsp)
  cvttsd2si %xmm0,%edi
  call   <sin@plt>
  movsd  D172(%rip),%xmm5
  mulsd  (%rsp),%xmm5
  mulsd  D164(%rip),%xmm0
  movsd  0x8(%rsp),%xmm13
  cvttsd2si %xmm5,%edi
  addsd  %xmm0,%xmm13
  movsd  %xmm13,0x10(%rsp)
  call   <sin@plt>
  movsd  %xmm0,0x8(%rsp)
  movsd  (%rsp),%xmm0
  call   <tan@plt>
  movsd  0x38(%rsp),%xmm1
  movsd  0x20(%rsp),%xmm5
  movapd %xmm0,%xmm9
  movsd  0x18(%rsp),%xmm7
  movapd %xmm0,%xmm6
  movsd  0x28(%rsp),%xmm12
  mulsd  %xmm0,%xmm9
  movapd %xmm1,%xmm11
  movsd  D247(%rip),%xmm0
  movsd  0x30(%rsp),%xmm4
  mulsd  %xmm1,%xmm11
  movsd  0x8(%rsp),%xmm14
  movsd  D194(%rip),%xmm8
  mulsd  %xmm7,%xmm0
  movsd  0x10(%rsp),%xmm13
  mulsd  D39(%rip),%xmm14
  movapd %xmm8,%xmm15
  subsd  %xmm9,%xmm15
  movapd %xmm11,%xmm3
  mulsd  %xmm1,%xmm3
  mulsd  %xmm7,%xmm0
  subsd  %xmm14,%xmm13
  mulsd  %xmm5,%xmm11
  mulsd  D97(%rip),%xmm11
  mulsd  D129(%rip),%xmm13
  movapd %xmm3,%xmm2
  mulsd  %xmm1,%xmm2
  addsd  %xmm15,%xmm0
  movsd  D197(%rip),%xmm15
  mulsd  %xmm5,%xmm3
  mulsd  %xmm7,%xmm15
  divsd  D172(%rip),%xmm3
  mulsd  %xmm4,%xmm3
  mulsd  %xmm12,%xmm11
  movapd %xmm2,%xmm10
  mulsd  %xmm5,%xmm10
  divsd  D225(%rip),%xmm10
  mulsd  %xmm12,%xmm10
  addsd  %xmm15,%xmm0
  mulsd  %xmm1,%xmm2
  mulsd  %xmm4,%xmm3
  mulsd  %xmm5,%xmm1
  mulsd  %xmm4,%xmm11
  mulsd  %xmm4,%xmm10
  mulsd  %xmm4,%xmm3
  mulsd  %xmm4,%xmm1
  addsd  %xmm13,%xmm11
  mulsd  %xmm4,%xmm10
  mulsd  %xmm4,%xmm10
  mulsd  %xmm0,%xmm10
  movapd %xmm2,%xmm0
  movsd  D67(%rip),%xmm2
  mulsd  %xmm5,%xmm0
  mulsd  %xmm6,%xmm2
  divsd  D110(%rip),%xmm0
  mulsd  %xmm4,%xmm0
  addsd  %xmm11,%xmm10
  mulsd  %xmm6,%xmm2
  mulsd  %xmm4,%xmm0
  movsd  %xmm10,(%r12)
  subsd  %xmm2,%xmm8
  movapd %xmm9,%xmm2
  mulsd  %xmm6,%xmm2
  mulsd  %xmm4,%xmm0
  mulsd  %xmm6,%xmm2
  mulsd  %xmm4,%xmm0
  addsd  %xmm2,%xmm8
  movsd  D149(%rip),%xmm2
  mulsd  %xmm4,%xmm0
  subsd  %xmm9,%xmm2
  addsd  %xmm7,%xmm2
  mulsd  %xmm8,%xmm0
  mulsd  %xmm2,%xmm3
  addsd  %xmm3,%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  addsd  D3(%rip),%xmm0
  cvtsi2sd %ebp,%xmm1
  mulsd  D95(%rip),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,(%rbx)
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
L16:
  movslq %edx,%rbp
  mov    %edx,%eax
  imul   $0x2aaaaaab,%rbp,%rbp
  sar    $0x1f,%eax
  shr    $0x20,%rbp
  sub    %eax,%ebp
  add    $0x1,%ebp
  jmp    L6
```

tool_calls: name: extract_symbolic_data
arguments: tag: D27
data_type: f64
name: extract_symbolic_data
arguments: tag: D131
data_type: f64
name: extract_symbolic_data
arguments: tag: D248
data_type: byte[8]
name: extract_symbolic_data
arguments: tag: D149
data_type: f64
name: extract_symbolic_data
arguments: tag: D129
data_type: byte[8]
name: extract_symbolic_data
arguments: tag: D118
data_type: f64
name: extract_symbolic_data
arguments: tag: D215
data_type: f64
name: extract_symbolic_data
arguments: tag: D210
data_type: byte[8]
name: extract_symbolic_data
arguments: tag: D197
data_type: f64
name: extract_symbolic_data
arguments: tag: D172
data_type: f64
name: extract_symbolic_data
arguments: tag: D164
data_type: byte[8]
name: extract_symbolic_data
arguments: tag: D247
data_type: byte[8]
name: extract_symbolic_data
arguments: tag: D194
data_type: byte[8]
name: extract_symbolic_data
arguments: tag: D39
data_type: byte[8]
name: extract_symbolic_data
arguments: tag: D97
data_type: byte[8]
name: extract_symbolic_data
arguments: tag: D225
data_type: byte[8]
name: extract_symbolic_data
arguments: tag: D67
data_type: byte[8]
name: extract_symbolic_data
arguments: tag: D110
data_type: byte[8]
name: extract_symbolic_data
arguments: tag: D3
data_type: byte[8]
name: extract_symbolic_data
arguments: tag: D95
data_type: byte[8]

tool: D27: 180

tool: D131: 3.141592653589793

tool: {"D248": [0x40,0x36,0xbf,0x50,0x67,0x6a,0x7b,0x3f]}

tool: D149: 1

tool: {"D129": [0x68,0x51,0x47,0x6b,0x10,0x2b,0x58,0x41]}

tool: D118: 6378245

tool: D215: 1.00505117739

tool: {"D210": [0xbb,0xf8,0x75,0xa6,0x49,0xbc,0x64,0x3f]}

tool: D197: 4

tool: D172: 6

tool: {"D164": [0xe4,0xfe,0x41,0x29,0xfc,0x47,0xc6,0x3e]}

tool: {"D247": [0x00,0x00,0x00,0x00,0x00,0x00,0x22,0x40]}

tool: {"D194": [0x00,0x00,0x00,0x00,0x00,0x00,0x14,0x40]}

tool: {"D39": [0xc5,0xa1,0xa4,0x16,0xf2,0xca,0x2d,0x3e]}

tool: {"D97": [0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0x3f]}

tool: {"D225": [0x00,0x00,0x00,0x00,0x00,0x00,0x38,0x40]}

tool: {"D67": [0x00,0x00,0x00,0x00,0x00,0x00,0x32,0x40]}

tool: {"D110": [0x00,0x00,0x00,0x00,0x00,0x00,0x5e,0x40]}

tool: {"D3": [0x00,0x00,0x00,0x00,0x80,0x84,0x1e,0x41]}

tool: {"D95": [0x00,0x00,0x00,0x00,0x80,0x84,0x2e,0x41]}

assistant: ```c
void SubGussFs(double *X, double *Y, double B, double L, int nCenterLongi) {
  int nzonenum;
  if (nCenterLongi == 0) {
    nzonenum = (int)L / 6 + 1;
    nCenterLongi = nzonenum * 6 - 3;
  } else
    nzonenum = (int)nCenterLongi / 6 + 1;
  double rB = B / 180 * 3.14159265358979323846;
  double rL = (L - nCenterLongi) / 180 * 3.14159265358979323846;
  const double a = 6378245.00;
  const double b = 6356863.50;
  double sqre1 = (a * a - b * b) / (a * a);
  double sinb = sin(rB);
  double cosb = cos(rB);
  double M = a * (1 - sqre1) / (1 - sqre1 * sinb * sinb) / sqrt(1 - sqre1 * sinb * sinb);
  double N = a / sqrt(1 - sqre1 * sinb * sinb);
  double sqrita = N / M - 1;
  double s = a * (1 - sqre1) * (1.00505117739 * rB - 0.00506237764 / 2 * sin(2 * rB) + 0.0000106245 / 4 * sin(4 * rB) - 0.00000002081 / 6 * sin(6 * rB));
  double tanb = tan(rB);
  *X = s + rL * rL * N / 2 * sinb * cosb + rL * rL * rL * rL * N / 24 * sinb * cosb * cosb * cosb * (5 - tanb * tanb + 9 * sqrita * sqrita + 4 * sqrita);
  *Y = rL * N * cosb + rL * rL * rL * N / 6 * cosb * cosb * cosb * (1 - tanb * tanb + sqrita) + rL * rL * rL * rL * rL * N / 120 * cosb * cosb * cosb * cosb * cosb * (5 - 18 * tanb * tanb + tanb * tanb * tanb * tanb);
  *Y = *Y + 500000 + nzonenum * 1.0e+6;
}
```