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

```asm
<GPS_Math_Miller_EN_To_LatLon>:
  endbr64
  push   %r14
  subsd  %xmm3,%xmm0
  subsd  %xmm4,%xmm1
  push   %rbp
  mov    %rsi,%rbp
  push   %rbx
  mov    %rdi,%rbx
  movq   %xmm0,%r14
  movapd %xmm2,%xmm0
  sub    $0x30,%rsp
  movsd  %xmm5,0x8(%rsp)
  movsd  %xmm6,0x18(%rsp)
  movsd  %xmm1,0x28(%rsp)
  call   <GPS_Math_Deg_To_Rad@plt>
  movsd  0x8(%rsp),%xmm4
  movsd  0x18(%rsp),%xmm6
  movsd  D46(%rip),%xmm1
  movsd  %xmm0,0x10(%rsp)
  mulsd  %xmm6,%xmm6
  movapd %xmm4,%xmm7
  mulsd  %xmm4,%xmm7
  movapd %xmm7,%xmm2
  subsd  %xmm6,%xmm2
  divsd  %xmm7,%xmm2
  movapd %xmm2,%xmm5
  movapd %xmm2,%xmm0
  movsd  %xmm2,0x18(%rsp)
  mulsd  %xmm2,%xmm5
  movsd  %xmm5,0x20(%rsp)
  call   <pow@plt>
  movsd  0x20(%rsp),%xmm4
  movsd  D20(%rip),%xmm5
  movsd  0x18(%rsp),%xmm2
  movsd  D108(%rip),%xmm3
  mulsd  %xmm4,%xmm5
  movsd  D10(%rip),%xmm1
  movsd  0x10(%rsp),%xmm7
  movapd %xmm2,%xmm6
  divsd  D166(%rip),%xmm6
  subsd  %xmm6,%xmm3
  mulsd  %xmm4,%xmm2
  comisd %xmm1,%xmm7
  mulsd  D160(%rip),%xmm2
  divsd  D114(%rip),%xmm2
  movapd %xmm5,%xmm0
  divsd  D139(%rip),%xmm0
  subsd  %xmm0,%xmm3
  movsd  0x28(%rsp),%xmm5
  subsd  %xmm2,%xmm3
  mulsd  0x8(%rsp),%xmm3
  jbe    L77
  subsd  D110(%rip),%xmm7
  movsd  %xmm7,0x10(%rsp)
L77:
  mulsd  D128(%rip),%xmm5
  movsd  %xmm3,0x8(%rsp)
  divsd  %xmm3,%xmm5
  movapd %xmm5,%xmm0
  call   <sinh@plt>
  mov    %eax,%edi
  call   <atan@plt>
  movsd  0x8(%rsp),%xmm3
  divsd  D128(%rip),%xmm0
  movsd  %xmm0,(%rbx)
  movq   %r14,%xmm0
  mov    D10(%rip),%rax
  divsd  %xmm3,%xmm0
  movsd  D47(%rip),%xmm3
  movq   %rax,%xmm1
  addsd  0x10(%rsp),%xmm0
  movsd  %xmm0,0x0(%rbp)
  movsd  (%rbx),%xmm2
  comisd %xmm3,%xmm2
  ja     L60
  movsd  D131(%rip),%xmm3
  comisd %xmm2,%xmm3
  ja     L60
L20:
  comisd %xmm1,%xmm0
  jbe    L33
  subsd  D110(%rip),%xmm0
  movsd  %xmm0,0x0(%rbp)
L33:
  movsd  D124(%rip),%xmm2
  comisd %xmm0,%xmm2
  jbe    L75
  addsd  D110(%rip),%xmm0
  movsd  %xmm0,0x0(%rbp)
L75:
  comisd %xmm1,%xmm0
  jbe    L24
  movsd  %xmm1,0x0(%rbp)
  movapd %xmm1,%xmm0
L81:
  call   <GPS_Math_Rad_To_Deg@plt>
  movsd  %xmm0,0x0(%rbp)
  movsd  (%rbx),%xmm0
  call   <GPS_Math_Rad_To_Deg@plt>
  movsd  %xmm0,(%rbx)
  add    $0x30,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r14
  ret
L60:
  movsd  %xmm3,(%rbx)
  movsd  0x0(%rbp),%xmm0
  jmp    L20
  nopl   0x0(%rax,%rax,1)
L24:
  comisd %xmm0,%xmm2
  jbe    L81
  movsd  %xmm2,0x0(%rbp)
  movapd %xmm2,%xmm0
  jmp    L81
```

tool_calls: name: extract_symbolic_data
arguments: label: D46
data_type: f64
name: extract_symbolic_data
arguments: label: D20
data_type: f64
name: extract_symbolic_data
arguments: label: D108
data_type: f64
name: extract_symbolic_data
arguments: label: D10
data_type: f64
name: extract_symbolic_data
arguments: label: D166
data_type: f64
name: extract_symbolic_data
arguments: label: D160
data_type: f64
name: extract_symbolic_data
arguments: label: D114
data_type: f64
name: extract_symbolic_data
arguments: label: D139
data_type: f64
name: extract_symbolic_data
arguments: label: D110
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D128
data_type: f64
name: extract_symbolic_data
arguments: label: D47
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D131
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D124
data_type: byte[8]

tool: D46: 0.5

tool: D20: 17

tool: D108: 1

tool: D10: 3.141592653589

tool: D166: 6

tool: D160: 67

tool: D114: 3024

tool: D139: 360

tool: {"D110": [0x1e,0x26,0x44,0x54,0xfb,0x21,0x19,0x40]}

tool: D128: 0.8

tool: {"D47": [0x1e,0x26,0x44,0x54,0xfb,0x21,0xf9,0x3f]}

tool: {"D131": [0x1e,0x26,0x44,0x54,0xfb,0x21,0xf9,0xbf]}

tool: {"D124": [0x1e,0x26,0x44,0x54,0xfb,0x21,0x09,0xc0]}

assistant: ```c
void GPS_Math_Miller_EN_To_LatLon(double E, double N, double *phi, double *lambda, double M0, double E0, double N0, double a, double b) {
  double a2;
  double b2;
  double R;
  double e;
  double e2;
  double e4;
  double e6;
  double p2;
  double po2;
  double dx;
  double dy;
  dx = E - E0;
  dy = N - N0;
  M0 = GPS_Math_Deg_To_Rad(M0);
  po2 = (double)3.141592653589 / (double)2.0;
  p2 = (double)3.141592653589 * (double)2.0;
  a2 = a * a;
  b2 = b * b;
  e2 = (a2 - b2) / a2;
  e4 = e2 * e2;
  e6 = e4 * e2;
  e = pow(e2, (double)0.5);
  R = a * ((double)1.0 - e2 / (double)6.0 - (double)17.0 * e4 / (double)360.0 - (double)67.0 * e6 / (double)3024.0);
  if (M0 > 3.141592653589)
    M0 -= p2;
  *phi = atan(sinh((double)0.8 * dy / R)) / (double)0.8;
  *lambda = M0 + dx / R;
  if (*phi > po2)
    *phi = po2;
  else if (*phi < -po2)
    *phi = -po2;
  if (*lambda > 3.141592653589)
    *lambda -= p2;
  if (*lambda < -3.141592653589)
    *lambda += p2;
  if (*lambda > 3.141592653589)
    *lambda = 3.141592653589;
  else if (*lambda < -3.141592653589)
    *lambda = -3.141592653589;
  *lambda = GPS_Math_Rad_To_Deg(*lambda);
  *phi = GPS_Math_Rad_To_Deg(*phi);
  return;
}
```