user: Please show me the C implementation that matches the assembly code listed below:

```asm
<GPS_Math_Polycon_EN_To_LatLon>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x88,%rsp
  movq   %xmm0,%r13
  movq   %xmm1,%rbx
  mov    %rdi,%r15
  mov    %rsi,%r14
  movapd %xmm2,%xmm0
  movq   %xmm3,%rbp
  movsd  %xmm4,0x18(%rsp)
  movsd  %xmm5,0x20(%rsp)
  movsd  %xmm6,0x30(%rsp)
  movsd  %xmm7,0x8(%rsp)
  call   <GPS_Math_Deg_To_Rad@plt>
  movsd  %xmm0,(%rsp)
  movq   %rbp,%xmm0
  call   <GPS_Math_Deg_To_Rad@plt>
  movsd  %xmm0,0x60(%rsp)
  comisd D378(%rip),%xmm0
  jbe    L57
  movapd %xmm0,%xmm4
  subsd  D412(%rip),%xmm4
  movsd  %xmm4,0x60(%rsp)
L57:
  movsd  0x8(%rsp),%xmm1
  mulsd  %xmm1,%xmm1
  movsd  0x30(%rsp),%xmm4
  movapd %xmm4,%xmm0
  mulsd  %xmm4,%xmm0
  movapd %xmm0,%xmm2
  subsd  %xmm1,%xmm2
  divsd  %xmm0,%xmm2
  movapd %xmm2,%xmm13
  movsd  %xmm2,0x10(%rsp)
  mulsd  %xmm2,%xmm2
  movapd %xmm13,%xmm3
  mulsd  %xmm2,%xmm3
  movapd %xmm3,%xmm0
  mulsd  D218(%rip),%xmm0
  mulsd  D331(%rip),%xmm0
  movsd  D203(%rip),%xmm1
  movapd %xmm2,%xmm7
  mulsd  %xmm1,%xmm7
  movapd %xmm7,%xmm4
  movapd %xmm13,%xmm5
  mulsd  D220(%rip),%xmm5
  movsd  D130(%rip),%xmm6
  subsd  %xmm5,%xmm6
  movapd %xmm7,%xmm5
  mulsd  D399(%rip),%xmm5
  subsd  %xmm5,%xmm6
  movapd %xmm3,%xmm5
  mulsd  D58(%rip),%xmm5
  movsd  D358(%rip),%xmm7
  mulsd  %xmm7,%xmm5
  subsd  %xmm5,%xmm6
  movsd  %xmm6,0x38(%rsp)
  mulsd  %xmm13,%xmm1
  mulsd  D305(%rip),%xmm1
  mulsd  D21(%rip),%xmm4
  addsd  %xmm4,%xmm1
  addsd  %xmm0,%xmm1
  movsd  %xmm1,0x8(%rsp)
  mulsd  D282(%rip),%xmm2
  mulsd  %xmm7,%xmm2
  addsd  %xmm0,%xmm2
  movsd  %xmm2,0x40(%rsp)
  mulsd  D319(%rip),%xmm3
  divsd  D442(%rip),%xmm3
  movsd  %xmm3,0x48(%rsp)
  movsd  (%rsp),%xmm4
  movapd %xmm4,%xmm0
  addsd  %xmm4,%xmm0
  call   <sin@plt>
  movq   %xmm0,%r12
  movsd  (%rsp),%xmm0
  mulsd  D126(%rip),%xmm0
  call   <sin@plt>
  movq   %xmm0,%rbp
  movsd  (%rsp),%xmm0
  mulsd  D389(%rip),%xmm0
  call   <sin@plt>
  movapd %xmm0,%xmm1
  movsd  (%rsp),%xmm0
  mulsd  0x38(%rsp),%xmm0
  movq   %r12,%xmm6
  mulsd  0x8(%rsp),%xmm6
  subsd  %xmm6,%xmm0
  movq   %rbp,%xmm6
  mulsd  0x40(%rsp),%xmm6
  addsd  %xmm6,%xmm0
  mulsd  0x48(%rsp),%xmm1
  subsd  %xmm1,%xmm0
  movsd  0x30(%rsp),%xmm4
  mulsd  %xmm4,%xmm0
  movq   %rbx,%xmm1
  subsd  0x20(%rsp),%xmm1
  movq   %r13,%xmm2
  subsd  0x18(%rsp),%xmm2
  movapd %xmm2,%xmm7
  divsd  %xmm4,%xmm7
  movsd  %xmm7,0x68(%rsp)
  movapd %xmm0,%xmm2
  xorpd  D128(%rip),%xmm2
  subsd  D130(%rip),%xmm2
  comisd %xmm2,%xmm1
  jbe    L107
  movsd  D130(%rip),%xmm2
  subsd  %xmm0,%xmm2
  comisd %xmm1,%xmm2
  jbe    L107
  movq   $0x0,(%r15)
  movsd  0x60(%rsp),%xmm3
  addsd  %xmm7,%xmm3
  movsd  %xmm3,0x60(%rsp)
  jmp    L79
L107:
  addsd  %xmm1,%xmm0
  divsd  0x30(%rsp),%xmm0
  movapd %xmm0,%xmm7
  movsd  %xmm0,0x50(%rsp)
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  0x68(%rsp),%xmm6
  movapd %xmm6,%xmm0
  mulsd  %xmm6,%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x58(%rsp)
  movsd  %xmm7,(%rsp)
L35:
  movsd  (%rsp),%xmm0
  call   <sin@plt>
  movapd %xmm0,%xmm1
  movsd  0x10(%rsp),%xmm4
  mulsd  %xmm0,%xmm4
  movapd %xmm4,%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D130(%rip),%xmm2
  subsd  %xmm0,%xmm2
  mov    D423(%rip),%rax
  movq   %rax,%xmm1
  movapd %xmm2,%xmm0
  call   <pow@plt>
  movq   %xmm0,%rbx
  movsd  (%rsp),%xmm0
  call   <tan@plt>
  movq   %rbx,%xmm6
  mulsd  %xmm0,%xmm6
  movsd  %xmm6,0x18(%rsp)
  lea    0x78(%rsp),%rbp
  lea    0x70(%rsp),%r12
  movsd  (%rsp),%xmm7
  movapd %xmm7,%xmm0
  addsd  %xmm7,%xmm0
  mov    %r12,%rsi
  mov    %rbp,%rdi
  call   <sincos@plt>
  movsd  0x70(%rsp),%xmm2
  movsd  %xmm2,0x20(%rsp)
  movsd  0x8(%rsp),%xmm8
  mulsd  0x78(%rsp),%xmm8
  movsd  %xmm8,0x28(%rsp)
  movsd  D126(%rip),%xmm0
  mulsd  (%rsp),%xmm0
  mov    %r12,%rsi
  mov    %rbp,%rdi
  call   <sincos@plt>
  mov    0x70(%rsp),%r13
  mov    0x78(%rsp),%rbx
  movsd  D389(%rip),%xmm0
  mulsd  (%rsp),%xmm0
  mov    %r12,%rsi
  mov    %rbp,%rdi
  call   <sincos@plt>
  movsd  0x8(%rsp),%xmm5
  movapd %xmm5,%xmm0
  addsd  %xmm5,%xmm0
  mulsd  0x20(%rsp),%xmm0
  movsd  0x38(%rsp),%xmm5
  movapd %xmm5,%xmm1
  subsd  %xmm0,%xmm1
  movsd  D126(%rip),%xmm2
  movsd  0x40(%rsp),%xmm10
  mulsd  %xmm10,%xmm2
  movq   %r13,%xmm11
  mulsd  %xmm11,%xmm2
  addsd  %xmm1,%xmm2
  movsd  D389(%rip),%xmm0
  movsd  0x48(%rsp),%xmm9
  mulsd  %xmm9,%xmm0
  mulsd  0x70(%rsp),%xmm0
  subsd  %xmm0,%xmm2
  movsd  (%rsp),%xmm7
  mulsd  %xmm7,%xmm5
  movapd %xmm5,%xmm1
  movsd  0x28(%rsp),%xmm8
  subsd  %xmm8,%xmm1
  movq   %rbx,%xmm0
  mulsd  %xmm10,%xmm0
  addsd  %xmm0,%xmm1
  movapd %xmm9,%xmm0
  mulsd  0x78(%rsp),%xmm0
  subsd  %xmm0,%xmm1
  movsd  0x30(%rsp),%xmm5
  mulsd  %xmm5,%xmm1
  divsd  %xmm5,%xmm1
  movsd  0x50(%rsp),%xmm5
  movapd %xmm5,%xmm3
  mulsd  %xmm1,%xmm3
  movapd %xmm1,%xmm0
  mulsd  %xmm1,%xmm0
  addsd  0x58(%rsp),%xmm0
  subsd  %xmm1,%xmm5
  movsd  0x18(%rsp),%xmm6
  movapd %xmm6,%xmm1
  mulsd  %xmm3,%xmm1
  addsd  %xmm5,%xmm1
  movsd  D423(%rip),%xmm4
  mulsd  %xmm0,%xmm4
  mulsd  %xmm6,%xmm4
  subsd  %xmm4,%xmm1
  addsd  %xmm3,%xmm3
  subsd  %xmm3,%xmm0
  movsd  0x10(%rsp),%xmm3
  mulsd  %xmm8,%xmm3
  mulsd  %xmm3,%xmm0
  mulsd  D220(%rip),%xmm0
  mulsd  %xmm6,%xmm0
  mulsd  %xmm2,%xmm6
  movapd %xmm6,%xmm3
  movsd  D414(%rip),%xmm4
  divsd  %xmm8,%xmm4
  subsd  %xmm4,%xmm3
  mulsd  %xmm5,%xmm3
  addsd  %xmm3,%xmm0
  subsd  %xmm2,%xmm0
  divsd  %xmm0,%xmm1
  subsd  %xmm1,%xmm7
  movsd  %xmm7,(%rsp)
  andpd  D439(%rip),%xmm1
  comisd D359(%rip),%xmm1
  ja     L35
  movsd  %xmm7,(%r15)
  comisd D68(%rip),%xmm7
  ja     L42
  movsd  D236(%rip),%xmm0
  comisd %xmm7,%xmm0
  ja     L97
  movapd %xmm7,%xmm0
  andpd  D439(%rip),%xmm0
  comisd D396(%rip),%xmm0
  jbe    L66
  movsd  D40(%rip),%xmm1
  comisd %xmm0,%xmm1
  ja     L79
L66:
  movsd  0x68(%rsp),%xmm0
  mulsd  0x18(%rsp),%xmm0
  call   <asin@plt>
  movq   %xmm0,%rbx
  movsd  (%rsp),%xmm0
  call   <sin@plt>
  movapd %xmm0,%xmm1
  movq   %rbx,%xmm0
  divsd  %xmm1,%xmm0
  addsd  0x60(%rsp),%xmm0
  movsd  %xmm0,0x60(%rsp)
L79:
  movsd  0x60(%rsp),%xmm4
  movsd  %xmm4,(%r14)
  comisd D378(%rip),%xmm4
  jbe    L46
  movapd %xmm4,%xmm0
  subsd  D412(%rip),%xmm0
  movsd  %xmm0,(%r14)
L46:
  movsd  (%r14),%xmm0
  movsd  D368(%rip),%xmm1
  comisd %xmm0,%xmm1
  jbe    L21
  addsd  D412(%rip),%xmm0
  movsd  %xmm0,(%r14)
L21:
  movsd  (%r14),%xmm0
  comisd D378(%rip),%xmm0
  jbe    L128
  mov    D378(%rip),%rax
  mov    %rax,(%r14)
L94:
  movsd  (%r14),%xmm0
  call   <GPS_Math_Rad_To_Deg@plt>
  movsd  %xmm0,(%r14)
  movsd  (%r15),%xmm0
  call   <GPS_Math_Rad_To_Deg@plt>
  movsd  %xmm0,(%r15)
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L128:
  movsd  D368(%rip),%xmm1
  comisd %xmm0,%xmm1
  jbe    L94
  movsd  %xmm1,(%r14)
  jmp    L94
L42:
  movsd  D68(%rip),%xmm0
L150:
  movsd  %xmm0,(%r15)
  jmp    L79
L97:
  movsd  D236(%rip),%xmm0
  jmp    L150
```

tool_calls: name: retrieve_data
arguments: data_tag: D378
type_guessed: f64[2]
name: retrieve_data
arguments: data_tag: D412
type_guessed: byte[8]
name: retrieve_data
arguments: data_tag: D218
type_guessed: f64
name: retrieve_data
arguments: data_tag: D331
type_guessed: f32[2]
name: retrieve_data
arguments: data_tag: D203
type_guessed: f64
name: retrieve_data
arguments: data_tag: D220
type_guessed: f32[2]
name: retrieve_data
arguments: data_tag: D130
type_guessed: f64
name: retrieve_data
arguments: data_tag: D399
type_guessed: f32[2]
name: retrieve_data
arguments: data_tag: D58
type_guessed: f64
name: retrieve_data
arguments: data_tag: D358
type_guessed: f32[2]
name: retrieve_data
arguments: data_tag: D305
type_guessed: f32[2]
name: retrieve_data
arguments: data_tag: D21
type_guessed: f32[2]
name: retrieve_data
arguments: data_tag: D282
type_guessed: f64
name: retrieve_data
arguments: data_tag: D319
type_guessed: f64
name: retrieve_data
arguments: data_tag: D442
type_guessed: f64
name: retrieve_data
arguments: data_tag: D126
type_guessed: f64
name: retrieve_data
arguments: data_tag: D389
type_guessed: f64
name: retrieve_data
arguments: data_tag: D128
type_guessed: f32[4]
name: retrieve_data
arguments: data_tag: D423
type_guessed: f32[2]
name: retrieve_data
arguments: data_tag: D414
type_guessed: f64
name: retrieve_data
arguments: data_tag: D439
type_guessed: qword[2]
name: retrieve_data
arguments: data_tag: D359
type_guessed: byte[16]
name: retrieve_data
arguments: data_tag: D68
type_guessed: byte[16]
name: retrieve_data
arguments: data_tag: D236
type_guessed: byte[8]
name: retrieve_data
arguments: data_tag: D396
type_guessed: byte[16]
name: retrieve_data
arguments: data_tag: D40
type_guessed: byte[8]
name: retrieve_data
arguments: data_tag: D368
type_guessed: byte[8]

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

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

tool: D218: 45

tool: {"D331": [0.0,0.8125]}

tool: D203: 3

tool: {"D220": [0.0,1.625]}

tool: D130: 1

tool: {"D399": [0.0,1.125]}

tool: D58: 5

tool: {"D358": [0.0,0.9375]}

tool: {"D305": [0.0,1.5]}

tool: {"D21": [0.0,1.25]}

tool: D282: 15

tool: D319: 35

tool: D442: 3072

tool: D126: 4

tool: D389: 6

tool: {"D128": [0.0,-0.0,0.0,0.0]}

tool: {"D423": [0.0,1.75]}

tool: D414: 2

tool: {"D439": [0x7fffffffffffffff,0x0000000000000000]}

tool: {"D359": [0x11,0xea,0x2d,0x81,0x99,0x97,0x71,0x3d,0xac,0x61,0xe9,0xd7,0xf0,0x21,0xf9,0x3f]}

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

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

tool: {"D396": [0xac,0x61,0xe9,0xd7,0xf0,0x21,0xf9,0x3f,0x90,0xea,0x9e,0xd0,0x05,0x22,0xf9,0x3f]}

tool: {"D40": [0x90,0xea,0x9e,0xd0,0x05,0x22,0xf9,0x3f]}

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

assistant: ```c
void GPS_Math_Polycon_EN_To_LatLon(double E, double N, double *phi, double *lambda, double phi0, double M0, double E0, double N0, double a, double b) {
  double p2;
  double po2;
  double a2;
  double b2;
  double e2;
  double e4;
  double e6;
  double AM0;
  double c0;
  double c1;
  double c2;
  double c3;
  double j;
  double te4;
  double phi0s2;
  double phi0s4;
  double phi0s6;
  double dx;
  double dy;
  double dxoa;
  double AA;
  double BB;
  double CC = (double)0.;
  double PHIn;
  double PHId;
  double PHIs;
  double PHI;
  double PHIs2;
  double PHIs4;
  double PHIs6;
  double Mn;
  double Mnp;
  double Ma;
  double AAMa;
  double mpb;
  double AAmin;
  double tol;
  double lat;
  phi0 = GPS_Math_Deg_To_Rad(phi0);
  M0 = GPS_Math_Deg_To_Rad(M0);
  p2 = (double)3.141592653589 * (double)2.0;
  po2 = (double)3.141592653589 / (double)2.0;
  if (M0 > 3.141592653589)
    M0 -= p2;
  a2 = a * a;
  b2 = b * b;
  e2 = (a2 - b2) / a2;
  e4 = e2 * e2;
  e6 = e2 * e4;
  j = (double)45.0 * e6 / (double)1024.0;
  te4 = (double)3.0 * e4;
  c0 = (double)1.0 - e2 / (double)4.0 - te4 / (double)64.0 - (double)5.0 * e6 / (double)256.0;
  c1 = (double)3.0 * e2 / (double)8.0 + te4 / (double)32.0 + j;
  c2 = (double)15.0 * e4 / (double)256.0 + j;
  c3 = (double)35.0 * e6 / (double)3072.0;
  lat = c0 * phi0;
  phi0s2 = c1 * sin((double)2.0 * phi0);
  phi0s4 = c2 * sin((double)4.0 * phi0);
  phi0s6 = c3 * sin((double)6.0 * phi0);
  AM0 = a * (lat - phi0s2 + phi0s4 - phi0s6);
  tol = (double)1.0e-12;
  dx = E - E0;
  dy = N - N0;
  dxoa = dx / a;
  if ((((-AM0 - (double)1.) < dy) && (dy < (-AM0 + (double)1.)))) {
    *phi = (double)0.;
    *lambda = dxoa + M0;
  } else {
    AA = (AM0 + dy) / a;
    BB = dxoa * dxoa + (AA * AA);
    PHIn = AA;
    PHId = (double)1.;
    while (fabs(PHId) > tol) {
      PHIs = sin(PHIn);
      CC = pow((double)1. - e2 * PHIs * PHIs, (double).5) * tan(PHIn);
      PHI = PHIn * c0;
      PHIs2 = c1 * sin((double)2.0 * PHIn);
      PHIs4 = c2 * sin((double)4.0 * PHIn);
      PHIs6 = c3 * sin((double)6.0 * PHIn);
      Mn = a * (PHI - PHIs2 + PHIs4 - PHIs6);
      Mnp = c0 - (double)2. * c1 * cos((double)2. * PHIn) + (double)4. * c2 * cos((double)4. * PHIn) - (double)6. * c3 * cos((double)6. * PHIn);
      Ma = Mn / a;
      AAMa = AA * Ma;
      mpb = Ma * Ma + BB;
      AAmin = AA - Ma;
      PHId = (AAMa * CC + AAmin - (double).5 * mpb * CC) / (e2 * PHIs2 * (mpb - (double)2. * AAMa) / (double)4. * CC + AAmin * (CC * Mnp - (double)2. / PHIs2) - Mnp);
      PHIn -= PHId;
    }
    *phi = PHIn;
    if (*phi > po2)
      *phi = po2;
    else if (*phi < -po2)
      *phi = -po2;
    if ((((po2 - (double).00001) < fabs(*phi)) && (fabs(*phi) < (po2 + (double).00001))))
      *lambda = M0;
    else
      *lambda = (asin(dxoa * CC)) / sin(*phi) + M0;
  }
  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;
}
```