user: What is the c source code of the assembly code below?

```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  D51(%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  D168(%rip),%xmm5
  movsd  0x18(%rsp),%xmm2
  movsd  D131(%rip),%xmm3
  mulsd  %xmm4,%xmm5
  movsd  D145(%rip),%xmm1
  movsd  0x10(%rsp),%xmm7
  movapd %xmm2,%xmm6
  divsd  D32(%rip),%xmm6
  subsd  %xmm6,%xmm3
  mulsd  %xmm4,%xmm2
  comisd %xmm1,%xmm7
  mulsd  D124(%rip),%xmm2
  divsd  D3(%rip),%xmm2
  movapd %xmm5,%xmm0
  divsd  D67(%rip),%xmm0
  subsd  %xmm0,%xmm3
  movsd  0x28(%rsp),%xmm5
  subsd  %xmm2,%xmm3
  mulsd  0x8(%rsp),%xmm3
  jbe    L44
  subsd  D21(%rip),%xmm7
  movsd  %xmm7,0x10(%rsp)
L44:
  mulsd  D20(%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  D20(%rip),%xmm0
  movsd  %xmm0,(%rbx)
  movq   %r14,%xmm0
  mov    D145(%rip),%rax
  divsd  %xmm3,%xmm0
  movsd  D106(%rip),%xmm3
  movq   %rax,%xmm1
  addsd  0x10(%rsp),%xmm0
  movsd  %xmm0,0x0(%rbp)
  movsd  (%rbx),%xmm2
  comisd %xmm3,%xmm2
  ja     L18
  movsd  D165(%rip),%xmm3
  comisd %xmm2,%xmm3
  ja     L18
L49:
  comisd %xmm1,%xmm0
  jbe    L29
  subsd  D21(%rip),%xmm0
  movsd  %xmm0,0x0(%rbp)
L29:
  movsd  D126(%rip),%xmm2
  comisd %xmm0,%xmm2
  jbe    L37
  addsd  D21(%rip),%xmm0
  movsd  %xmm0,0x0(%rbp)
L37:
  comisd %xmm1,%xmm0
  jbe    L45
  movsd  %xmm1,0x0(%rbp)
  movapd %xmm1,%xmm0
L8:
  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
L18:
  movsd  %xmm3,(%rbx)
  movsd  0x0(%rbp),%xmm0
  jmp    L49
  nopl   0x0(%rax,%rax,1)
L45:
  comisd %xmm0,%xmm2
  jbe    L8
  movsd  %xmm2,0x0(%rbp)
  movapd %xmm2,%xmm0
  jmp    L8
```

tool_calls: name: read_data
arguments: label: D51
type_guessed: f64
name: read_data
arguments: label: D168
type_guessed: f64
name: read_data
arguments: label: D131
type_guessed: f64
name: read_data
arguments: label: D145
type_guessed: f64
name: read_data
arguments: label: D32
type_guessed: f64
name: read_data
arguments: label: D124
type_guessed: f64
name: read_data
arguments: label: D3
type_guessed: f64
name: read_data
arguments: label: D67
type_guessed: f64
name: read_data
arguments: label: D21
type_guessed: byte[8]
name: read_data
arguments: label: D20
type_guessed: f64
name: read_data
arguments: label: D106
type_guessed: byte[8]
name: read_data
arguments: label: D165
type_guessed: byte[8]
name: read_data
arguments: label: D126
type_guessed: byte[8]

tool: D51: 0.5

tool: D168: 17

tool: D131: 1

tool: D145: 3.141592653589

tool: D32: 6

tool: D124: 67

tool: D3: 3024

tool: D67: 360

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

tool: D20: 0.8

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

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

tool: {"D126": [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;
}
```