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

```asm
<GPS_Math_Sinusoid_EN_To_LatLon>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  movq   %xmm0,%r12
  movq   %xmm1,%r14
  mov    %rdi,%r13
  mov    %rsi,%rbx
  movapd %xmm2,%xmm0
  movq   %xmm3,%r15
  movsd  %xmm4,0x20(%rsp)
  movsd  %xmm5,0x10(%rsp)
  movsd  %xmm6,0x8(%rsp)
  call   <GPS_Math_Deg_To_Rad@plt>
  movq   %xmm0,%rbp
  movsd  0x10(%rsp),%xmm4
  movapd %xmm4,%xmm1
  mulsd  %xmm4,%xmm1
  movsd  0x8(%rsp),%xmm0
  mulsd  %xmm0,%xmm0
  movapd %xmm1,%xmm2
  subsd  %xmm0,%xmm2
  divsd  %xmm1,%xmm2
  movapd %xmm2,%xmm8
  mulsd  %xmm2,%xmm8
  movsd  %xmm8,0x8(%rsp)
  movsd  D328(%rip),%xmm0
  movsd  %xmm2,0x18(%rsp)
  subsd  %xmm2,%xmm0
  mov    D357(%rip),%rax
  movq   %rax,%xmm1
  call   <pow@plt>
  movsd  D328(%rip),%xmm6
  subsd  %xmm0,%xmm6
  movapd %xmm6,%xmm1
  addsd  D328(%rip),%xmm0
  divsd  %xmm0,%xmm1
  movapd %xmm1,%xmm6
  mulsd  %xmm1,%xmm6
  movapd %xmm1,%xmm7
  movsd  %xmm6,0x38(%rsp)
  mulsd  %xmm6,%xmm7
  movsd  %xmm1,0x30(%rsp)
  movapd %xmm1,%xmm4
  movsd  %xmm7,0x40(%rsp)
  mulsd  %xmm7,%xmm4
  movsd  %xmm4,0x48(%rsp)
  movq   %r14,%xmm2
  subsd  0x20(%rsp),%xmm2
  movsd  0x18(%rsp),%xmm3
  movapd %xmm3,%xmm1
  mulsd  D279(%rip),%xmm1
  movsd  D328(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movsd  D399(%rip),%xmm9
  movsd  0x8(%rsp),%xmm8
  mulsd  %xmm8,%xmm9
  movapd %xmm9,%xmm1
  mulsd  D237(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movapd %xmm8,%xmm1
  mulsd  %xmm3,%xmm1
  mulsd  D219(%rip),%xmm1
  mulsd  D140(%rip),%xmm1
  subsd  %xmm1,%xmm0
  mulsd  0x10(%rsp),%xmm0
  divsd  %xmm0,%xmm2
  movapd %xmm2,%xmm0
  movsd  %xmm2,0x8(%rsp)
  addsd  %xmm2,%xmm0
  call   <sin@plt>
  movq   %xmm0,%r14
  movsd  0x8(%rsp),%xmm0
  mulsd  D9(%rip),%xmm0
  call   <sin@plt>
  movsd  %xmm0,0x20(%rsp)
  movsd  0x8(%rsp),%xmm0
  mulsd  D285(%rip),%xmm0
  call   <sin@plt>
  movsd  %xmm0,0x28(%rsp)
  movsd  0x8(%rsp),%xmm0
  mulsd  D351(%rip),%xmm0
  call   <sin@plt>
  movapd %xmm0,%xmm2
  movsd  D399(%rip),%xmm1
  mulsd  0x30(%rsp),%xmm1
  movapd %xmm1,%xmm0
  mulsd  D357(%rip),%xmm0
  movsd  0x40(%rsp),%xmm7
  movapd %xmm7,%xmm1
  mulsd  D393(%rip),%xmm1
  movsd  D368(%rip),%xmm4
  mulsd  %xmm4,%xmm1
  subsd  %xmm1,%xmm0
  movq   %r14,%xmm1
  mulsd  %xmm1,%xmm0
  addsd  0x8(%rsp),%xmm0
  movsd  0x38(%rsp),%xmm6
  mulsd  D30(%rip),%xmm6
  movapd %xmm6,%xmm1
  mulsd  D306(%rip),%xmm1
  movsd  0x48(%rsp),%xmm5
  movapd %xmm5,%xmm3
  mulsd  D397(%rip),%xmm3
  mulsd  %xmm4,%xmm3
  subsd  %xmm3,%xmm1
  mulsd  0x20(%rsp),%xmm1
  addsd  %xmm1,%xmm0
  mulsd  D48(%rip),%xmm7
  movapd %xmm7,%xmm1
  divsd  D132(%rip),%xmm1
  mulsd  0x28(%rsp),%xmm1
  addsd  %xmm1,%xmm0
  mulsd  D288(%rip),%xmm5
  movapd %xmm5,%xmm1
  mulsd  D271(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%r13)
  comisd D256(%rip),%xmm0
  ja     L95
  movsd  D80(%rip),%xmm1
  comisd %xmm0,%xmm1
  ja     L90
  movapd %xmm0,%xmm1
  andpd  D121(%rip),%xmm1
  comisd D56(%rip),%xmm1
  jbe    L87
  movsd  D383(%rip),%xmm2
  comisd %xmm1,%xmm2
  ja     L53
L87:
  lea    0x58(%rsp),%rdi
  lea    0x50(%rsp),%rsi
  call   <sincos@plt>
  mov    0x50(%rsp),%r14
  movsd  0x58(%rsp),%xmm0
  movsd  0x18(%rsp),%xmm1
  mulsd  %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  D328(%rip),%xmm0
  subsd  %xmm1,%xmm0
  movsd  D357(%rip),%xmm1
  call   <pow@plt>
  movapd %xmm0,%xmm1
  movq   %r12,%xmm0
  movq   %r15,%xmm6
  subsd  %xmm6,%xmm0
  mulsd  %xmm1,%xmm0
  movq   %r14,%xmm1
  mulsd  0x10(%rsp),%xmm1
  divsd  %xmm1,%xmm0
  movq   %rbp,%xmm4
  addsd  %xmm0,%xmm4
  movq   %xmm4,%rbp
L53:
  mov    %rbp,(%rbx)
  movq   %rbp,%xmm7
  comisd D246(%rip),%xmm7
  jbe    L84
  movq   %rbp,%xmm0
  subsd  D122(%rip),%xmm0
  movsd  %xmm0,(%rbx)
L84:
  movsd  (%rbx),%xmm0
  movsd  D189(%rip),%xmm1
  comisd %xmm0,%xmm1
  jbe    L50
  addsd  D122(%rip),%xmm0
  movsd  %xmm0,(%rbx)
L50:
  movsd  (%rbx),%xmm0
  comisd D246(%rip),%xmm0
  jbe    L22
  mov    D246(%rip),%rax
  mov    %rax,(%rbx)
L58:
  movsd  (%rbx),%xmm0
  call   <GPS_Math_Rad_To_Deg@plt>
  movsd  %xmm0,(%rbx)
  movsd  0x0(%r13),%xmm0
  call   <GPS_Math_Rad_To_Deg@plt>
  movsd  %xmm0,0x0(%r13)
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L22:
  movsd  D189(%rip),%xmm1
  comisd %xmm0,%xmm1
  jbe    L58
  movsd  %xmm1,(%rbx)
  jmp    L58
L95:
  movsd  D256(%rip),%xmm0
L31:
  movsd  %xmm0,0x0(%r13)
  jmp    L53
L90:
  movsd  D80(%rip),%xmm0
  jmp    L31
```

tool_calls: name: unpack_data
arguments: data_label: D328
data_type: f64
name: unpack_data
arguments: data_label: D357
data_type: byte[8]
name: unpack_data
arguments: data_label: D279
data_type: byte[8]
name: unpack_data
arguments: data_label: D399
data_type: f64
name: unpack_data
arguments: data_label: D237
data_type: byte[8]
name: unpack_data
arguments: data_label: D219
data_type: f64
name: unpack_data
arguments: data_label: D140
data_type: byte[8]
name: unpack_data
arguments: data_label: D9
data_type: f64
name: unpack_data
arguments: data_label: D285
data_type: f64
name: unpack_data
arguments: data_label: D351
data_type: f64
name: unpack_data
arguments: data_label: D393
data_type: f64
name: unpack_data
arguments: data_label: D368
data_type: byte[8]
name: unpack_data
arguments: data_label: D30
data_type: f64
name: unpack_data
arguments: data_label: D306
data_type: byte[8]
name: unpack_data
arguments: data_label: D397
data_type: f64
name: unpack_data
arguments: data_label: D48
data_type: f64
name: unpack_data
arguments: data_label: D132
data_type: f64
name: unpack_data
arguments: data_label: D288
data_type: f64
name: unpack_data
arguments: data_label: D271
data_type: byte[8]
name: unpack_data
arguments: data_label: D256
data_type: byte[16]
name: unpack_data
arguments: data_label: D80
data_type: byte[8]
name: unpack_data
arguments: data_label: D121
data_type: qword[2]
name: unpack_data
arguments: data_label: D56
data_type: byte[16]
name: unpack_data
arguments: data_label: D383
data_type: byte[8]
name: unpack_data
arguments: data_label: D246
data_type: f64[2]
name: unpack_data
arguments: data_label: D122
data_type: byte[8]
name: unpack_data
arguments: data_label: D189
data_type: byte[8]

tool: D328: 1

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

tool: {"D279": [0x00,0x00,0x00,0x00,0x00,0x00,0xd0,0x3f]}

tool: D399: 3

tool: {"D237": [0x00,0x00,0x00,0x00,0x00,0x00,0x90,0x3f]}

tool: D219: 5

tool: {"D140": [0x00,0x00,0x00,0x00,0x00,0x00,0x70,0x3f]}

tool: D9: 4

tool: D285: 6

tool: D351: 8

tool: D393: 27

tool: {"D368": [0x00,0x00,0x00,0x00,0x00,0x00,0xa0,0x3f]}

tool: D30: 21

tool: {"D306": [0x00,0x00,0x00,0x00,0x00,0x00,0xb0,0x3f]}

tool: D397: 55

tool: D48: 151

tool: D132: 96

tool: D288: 1097

tool: {"D271": [0x00,0x00,0x00,0x00,0x00,0x00,0x60,0x3f]}

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

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

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

tool: {"D56": [0x42,0xf4,0x94,0x51,0xfb,0x21,0xf9,0x3f,0xfa,0x57,0xf3,0x56,0xfb,0x21,0xf9,0x3f]}

tool: {"D383": [0xfa,0x57,0xf3,0x56,0xfb,0x21,0xf9,0x3f]}

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

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

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

assistant: ```c
void GPS_Math_Sinusoid_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 e2;
  double e4;
  double e6;
  double p2;
  double po2;
  double c0;
  double c1;
  double c2;
  double c3;
  double A1;
  double A0;
  double A2;
  double A3;
  double E1;
  double E2;
  double E3;
  double E4;
  double j;
  double om0;
  double som0;
  double dx;
  double dy;
  double mu;
  double mu2s;
  double mu4s;
  double mu6s;
  double mu8s;
  double phis;
  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;
  j = (double)45. * e6 / (double)1024.;
  c0 = (double)1. - e2 / (double)4. - (double)3. * e4 / (double)64. - (double)5. * e6 / (double)256.;
  c1 = (double)3. * e2 / (double)8. + (double)3. * e4 / (double)32. + j;
  c2 = (double)15. * e4 / (double)256. + j;
  c3 = (double)35. * e6 / (double)3072.;
  om0 = (double)1. - e2;
  som0 = pow(om0, (double).5);
  E1 = ((double)1. - som0) / ((double)1. + som0);
  E2 = E1 * E1;
  E3 = E1 * E2;
  E4 = E1 * E3;
  A0 = (double)3. * E1 / (double)2. - (double)27. * E3 / (double)32.;
  A1 = (double)21. * E2 / (double)16. - (double)55. * E4 / (double)32.;
  A2 = (double)151. * E3 / (double)96.;
  A3 = (double)1097. * E4 / (double)512.;
  dx = E - E0;
  dy = N - N0;
  mu = dy / (c0 * a);
  mu2s = A0 * sin((double)2. * mu);
  mu4s = A1 * sin((double)4. * mu);
  mu6s = A2 * sin((double)6. * mu);
  mu8s = A3 * sin((double)8. * mu);
  *phi = mu + mu2s + mu4s + mu6s + mu8s;
  if (*phi > po2)
    *phi = po2;
  else if (*phi < -po2)
    *phi = -po2;
  if ((((po2 - (double)1.0e-8) < fabs(*phi)) && (fabs(*phi) < (po2 + (double)1.0e-8))))
    *lambda = M0;
  else {
    phis = sin(*phi);
    *lambda = M0 + dx * pow((double)1.0 - e2 * phis * phis, (double).5) / (a * cos(*phi));
  }
  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;
}
```