user: What would the C source look like for the assembly code provided here:

```asm
<GPS_Math_Miller_EN_To_LatLon>:
  endbr64
  push   %r15
  push   %r14
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x20,%rsp
  mov    %rdi,%rbp
  mov    %rsi,%rbx
  movsd  %xmm5,0x8(%rsp)
  movsd  %xmm6,0x10(%rsp)
  movapd %xmm0,%xmm7
  subsd  %xmm3,%xmm7
  movq   %xmm7,%r15
  subsd  %xmm4,%xmm1
  movq   %xmm1,%r12
  movapd %xmm2,%xmm0
  call   <GPS_Math_Deg_To_Rad@plt>
  movq   %xmm0,%r14
  movsd  0x8(%rsp),%xmm4
  movapd %xmm4,%xmm1
  mulsd  %xmm4,%xmm1
  movsd  0x10(%rsp),%xmm0
  mulsd  %xmm0,%xmm0
  movapd %xmm1,%xmm2
  subsd  %xmm0,%xmm2
  divsd  %xmm1,%xmm2
  movapd %xmm2,%xmm3
  mulsd  %xmm2,%xmm3
  movsd  %xmm3,0x18(%rsp)
  movsd  D179(%rip),%xmm1
  movsd  %xmm2,0x10(%rsp)
  movapd %xmm2,%xmm0
  call   <pow@plt>
  movsd  0x10(%rsp),%xmm7
  movapd %xmm7,%xmm0
  divsd  D195(%rip),%xmm0
  movsd  D149(%rip),%xmm1
  subsd  %xmm0,%xmm1
  movsd  0x18(%rsp),%xmm3
  movapd %xmm3,%xmm0
  mulsd  D41(%rip),%xmm0
  divsd  D13(%rip),%xmm0
  subsd  %xmm0,%xmm1
  mulsd  %xmm3,%xmm7
  movapd %xmm7,%xmm0
  mulsd  D65(%rip),%xmm0
  divsd  D93(%rip),%xmm0
  subsd  %xmm0,%xmm1
  mulsd  0x8(%rsp),%xmm1
  movsd  %xmm1,0x8(%rsp)
  movq   %r14,%xmm5
  comisd D194(%rip),%xmm5
  jbe    L7
  subsd  D206(%rip),%xmm5
  movq   %xmm5,%r14
L7:
  movsd  D87(%rip),%xmm0
  movq   %r12,%xmm6
  mulsd  %xmm6,%xmm0
  divsd  0x8(%rsp),%xmm0
  call   <sinh@plt>
  mov    %eax,%edi
  call   <atan@plt>
  divsd  D87(%rip),%xmm0
  movsd  %xmm0,0x0(%rbp)
  movq   %r15,%xmm0
  divsd  0x8(%rsp),%xmm0
  movq   %r14,%xmm5
  addsd  %xmm5,%xmm0
  movsd  %xmm0,(%rbx)
  movsd  0x0(%rbp),%xmm0
  comisd D91(%rip),%xmm0
  jbe    L17
  mov    D91(%rip),%rax
  mov    %rax,0x0(%rbp)
L4:
  movsd  (%rbx),%xmm0
  comisd D194(%rip),%xmm0
  jbe    L1
  subsd  D206(%rip),%xmm0
  movsd  %xmm0,(%rbx)
L1:
  movsd  (%rbx),%xmm0
  movsd  D16(%rip),%xmm1
  comisd %xmm0,%xmm1
  jbe    L9
  addsd  D206(%rip),%xmm0
  movsd  %xmm0,(%rbx)
L9:
  movsd  (%rbx),%xmm0
  comisd D194(%rip),%xmm0
  jbe    L77
  mov    D194(%rip),%rax
  mov    %rax,(%rbx)
L69:
  movsd  (%rbx),%xmm0
  call   <GPS_Math_Rad_To_Deg@plt>
  movsd  %xmm0,(%rbx)
  movsd  0x0(%rbp),%xmm0
  call   <GPS_Math_Rad_To_Deg@plt>
  movsd  %xmm0,0x0(%rbp)
  add    $0x20,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r14
  pop    %r15
  ret
L17:
  movsd  D177(%rip),%xmm1
  comisd %xmm0,%xmm1
  jbe    L4
  movsd  %xmm1,0x0(%rbp)
  jmp    L4
L77:
  movsd  D16(%rip),%xmm1
  comisd %xmm0,%xmm1
  jbe    L69
  movsd  %xmm1,(%rbx)
  jmp    L69
```

tool_calls: name: retrieve_data
arguments: data_label: D179
type: f64
name: retrieve_data
arguments: data_label: D195
type: f64
name: retrieve_data
arguments: data_label: D149
type: f64
name: retrieve_data
arguments: data_label: D41
type: f64
name: retrieve_data
arguments: data_label: D13
type: f64
name: retrieve_data
arguments: data_label: D65
type: f64
name: retrieve_data
arguments: data_label: D93
type: f64
name: retrieve_data
arguments: data_label: D194
type: f64[2]
name: retrieve_data
arguments: data_label: D206
type: byte[8]
name: retrieve_data
arguments: data_label: D87
type: f64
name: retrieve_data
arguments: data_label: D91
type: byte[16]
name: retrieve_data
arguments: data_label: D16
type: byte[8]
name: retrieve_data
arguments: data_label: D177
type: byte[8]

tool: D179: 0.5

tool: D195: 6

tool: D149: 1

tool: D41: 17

tool: D13: 360

tool: D65: 67

tool: D93: 3024

tool: {"D194": [3.141592653589,6.283185307178]}

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

tool: D87: 0.8

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

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

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

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;
}
```