void GPS_Math_TCylEA_LatLon_To_EN(double phi, double lambda, double *E,
                                  double *N, double phi0, double M0, double E0,
                                  double N0, double a, double b) {
  double p2;
  double po2;
  double a2;
  double b2;
  double e;
  double e2;
  double e4;
  double e6;
  double AM0;
  double qp;
  double om;
  double oo;
  double c0;
  double c1;
  double c2;
  double c3;
  double b0;
  double b1;
  double B2;
  double b3;
  double A0;
  double A1;
  double A2;
  double sf;
  double x;
  double som;
  double phis;
  double j;
  double te4;
  double lat;
  double phi0s2;
  double phi0s4;
  double phi0s6;
  double E1;
  double E2;
  double E3;
  double E4;
  double dlam;
  double qq;
  double qqo;
  double bt;
  double btc;
  double PHI;
  double PHIs2;
  double PHIs4;
  double PHIs6;
  double bts2;
  double bts4;
  double bts6;
  double PHIc;
  double PHIcs;
  double Mc;
  sf = (double)1.0;
  lambda = GPS_Math_Deg_To_Rad(lambda);
  phi = GPS_Math_Deg_To_Rad(phi);
  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;
  e = pow(e2, (double).5);
  om = (double)1. - e2;
  som = pow(om, (double).5);
  oo = (double)1. / ((double)2. * e);
  phis = sin(po2);
  x = e * phis;
  qp = om * (phis / ((double)1. - e2 * phis * phis) -
             oo * log(((double)1. - x) / ((double)1. + x)));
  A0 = e2 / (double)3. + (double)31. * e4 / (double)180. +
       (double)517. * e6 / (double)5040.;
  A1 = (double)23. * e4 / (double)360. + (double)251. * e6 / (double)3780.;
  A2 = (double)761. * e6 / (double)45360.;
  E1 = ((double)1.0 - som) / ((double)1.0 + som);
  E2 = E1 * E1;
  E3 = E2 * E1;
  E4 = E3 * E1;
  b0 = (double)3. * E1 / (double)2. - (double)27. * E3 / (double)32.;
  b1 = (double)21. * E2 / (double)16. - (double)55. * E4 / (double)32.;
  B2 = (double)151. * E3 / (double)96.;
  b3 = (double)1097. * E4 / (double)512.;
  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);
  dlam = lambda - M0;
  if (dlam > 3.141592653589)
    dlam -= p2;
  if (dlam < -3.141592653589)
    dlam += p2;
  phis = sin(phi);
  if (phi == po2) {
    qq = qp;
    qqo = (double)1.;
  } else {
    x = e * phis;
    qq = om * (phis / ((double)1. - e2 * phis * phis) -
               oo * log(((double)1. - x) / ((double)1. + x)));
    qqo = qq / qp;
  }
  if (qqo > (double)1.)
    qqo = (double)1.;
  else if (qqo < (double)-1.)
    qqo = (double)-1.;
  bt = asin(qqo);
  btc = atan(tan(bt) / cos(dlam));
  if ((fabs(btc) - po2) > (double)1.0e-8)
    PHIc = btc;
  else {
    bts2 = A0 * sin((double)2.0 * btc);
    bts4 = A1 * sin((double)4.0 * btc);
    bts6 = A2 * sin((double)6.0 * btc);
    PHIc = btc + bts2 + bts4 + bts6;
  }
  PHIcs = sin(PHIc);
  *E = a * cos(bt) * cos(PHIc) * sin(dlam) /
           (sf * cos(btc) * pow((double)1. - e2 * PHIcs * PHIcs, (double).5)) +
       E0;
  PHI = c0 * PHIc;
  PHIs2 = c1 * sin((double)2.0 * PHIc);
  PHIs4 = c2 * sin((double)4.0 * PHIc);
  PHIs6 = c3 * sin((double)6.0 * PHIc);
  Mc = a * (PHI - PHIs2 + PHIs4 - PHIs6);
  *N = sf * (Mc - AM0) + N0;
  return;
}


x86 assembly (clang, -O0):
<GPS_Math_TCylEA_LatLon_To_EN>:
push   %rbp
mov    %rsp,%rbp
sub    $0x2a0,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
mov    %rdi,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
movsd  %xmm2,-0x28(%rbp)
movsd  %xmm3,-0x30(%rbp)
movsd  %xmm4,-0x38(%rbp)
movsd  %xmm5,-0x40(%rbp)
movsd  %xmm6,-0x48(%rbp)
movsd  %xmm7,-0x50(%rbp)
movsd  0xe55(%rip),%xmm0
movsd  %xmm0,-0x110(%rbp)
ff
movsd  -0x10(%rbp),%xmm0
mov    $0x1,%al
call   1030 <GPS_Math_Deg_To_Rad@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x10(%rbp)
movsd  -0x8(%rbp),%xmm0
mov    $0x1,%al
call   1030 <GPS_Math_Deg_To_Rad@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    $0x1,%al
call   1030 <GPS_Math_Deg_To_Rad@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x30(%rbp),%xmm0
mov    $0x1,%al
call   1030 <GPS_Math_Deg_To_Rad@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xde9(%rip),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  0xdd4(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  0xdba(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1263 <GPS_Math_TCylEA_LatLon_To_EN+0xe3>
movsd  -0x58(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x48(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x50(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x68(%rbp),%xmm0
subsd  -0x70(%rbp),%xmm0
divsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
mulsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  0xe57(%rip),%xmm1
call   10a0 <pow@plt>
movsd  %xmm0,-0x78(%rbp)
movsd  0xd3d(%rip),%xmm0
subsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
movsd  0xe28(%rip),%xmm1
call   10a0 <pow@plt>
movsd  %xmm0,-0x120(%rbp)
ff
movsd  0xd23(%rip),%xmm1
mulsd  -0x78(%rbp),%xmm1
movsd  0xcfe(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
call   1040 <sin@plt>
movsd  %xmm0,-0x128(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
mulsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0x118(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0x210(%rbp)
ff
movsd  -0x128(%rbp),%xmm0
ff
movsd  -0x80(%rbp),%xmm1
mulsd  -0x128(%rbp),%xmm1
ff
movsd  -0x128(%rbp),%xmm2
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm1
movsd  0xc7b(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x218(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
movsd  %xmm0,-0x220(%rbp)
ff
movsd  0xc53(%rip),%xmm0
subsd  -0x118(%rbp),%xmm0
ff
movsd  0xc43(%rip),%xmm1
addsd  -0x118(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
call   1090 <log@plt>
movsd  -0x220(%rbp),%xmm1
ff
movsd  -0x218(%rbp),%xmm2
ff
movaps %xmm0,%xmm3
movsd  -0x210(%rbp),%xmm0
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  0xc2f(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  0xcd3(%rip),%xmm1
mulsd  -0x88(%rbp),%xmm1
ff
movsd  0xcbb(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  0xca3(%rip),%xmm1
mulsd  -0x90(%rbp),%xmm1
ff
movsd  0xc8b(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
movsd  0xc6b(%rip),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
movsd  0xc53(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  0xc3f(%rip),%xmm1
mulsd  -0x90(%rbp),%xmm1
ff
movsd  0xc27(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
movsd  0xc07(%rip),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
movsd  0xbef(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
movsd  0xb23(%rip),%xmm0
subsd  -0x120(%rbp),%xmm0
ff
movsd  0xb13(%rip),%xmm1
addsd  -0x120(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x160(%rbp)
ff
movsd  -0x160(%rbp),%xmm0
ff
mulsd  -0x160(%rbp),%xmm0
ff
movsd  %xmm0,-0x168(%rbp)
ff
movsd  -0x168(%rbp),%xmm0
ff
mulsd  -0x160(%rbp),%xmm0
ff
movsd  %xmm0,-0x170(%rbp)
ff
movsd  -0x170(%rbp),%xmm0
ff
mulsd  -0x160(%rbp),%xmm0
ff
movsd  %xmm0,-0x178(%rbp)
ff
movsd  0xaff(%rip),%xmm0
mulsd  -0x160(%rbp),%xmm0
ff
movsd  0xab7(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  0xb43(%rip),%xmm1
mulsd  -0x170(%rbp),%xmm1
ff
movsd  0xac3(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movsd  0xb13(%rip),%xmm0
mulsd  -0x168(%rbp),%xmm0
ff
movsd  0xafb(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  0xae7(%rip),%xmm1
mulsd  -0x178(%rbp),%xmm1
ff
movsd  0xa7f(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  0xab7(%rip),%xmm0
mulsd  -0x170(%rbp),%xmm0
ff
movsd  0xa9f(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
movsd  0xa83(%rip),%xmm0
mulsd  -0x178(%rbp),%xmm0
ff
movsd  0xa6b(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  0xa4f(%rip),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
movsd  0xa37(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
movsd  0xa0b(%rip),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x138(%rbp)
ff
movsd  -0x80(%rbp),%xmm1
movsd  0x9ae(%rip),%xmm0
divsd  %xmm0,%xmm1
movsd  0x992(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  -0x138(%rbp),%xmm1
ff
movsd  0x9de(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  0x9c6(%rip),%xmm1
mulsd  -0x90(%rbp),%xmm1
ff
movsd  0x98e(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  0x996(%rip),%xmm0
mulsd  -0x80(%rbp),%xmm0
movsd  0x981(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  -0x138(%rbp),%xmm1
ff
movsd  0x965(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
addsd  -0x130(%rbp),%xmm0
ff
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  0x93d(%rip),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
movsd  0x925(%rip),%xmm1
divsd  %xmm1,%xmm0
addsd  -0x130(%rbp),%xmm0
ff
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  0x901(%rip),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
movsd  0x8e9(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
ff
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x140(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0x208(%rbp)
ff
movsd  0x8a8(%rip),%xmm0
mulsd  -0x28(%rbp),%xmm0
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x208(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x148(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
movsd  %xmm0,-0x200(%rbp)
ff
movsd  0x867(%rip),%xmm0
mulsd  -0x28(%rbp),%xmm0
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x200(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x150(%rbp)
ff
movsd  -0xd0(%rbp),%xmm0
ff
movsd  %xmm0,-0x1f8(%rbp)
ff
movsd  0x826(%rip),%xmm0
mulsd  -0x28(%rbp),%xmm0
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x1f8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x158(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
movsd  -0x140(%rbp),%xmm1
ff
subsd  -0x148(%rbp),%xmm1
ff
addsd  -0x150(%rbp),%xmm1
ff
subsd  -0x158(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x10(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x180(%rbp)
ff
movsd  -0x180(%rbp),%xmm0
ff
movsd  0x792(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1891 <GPS_Math_TCylEA_LatLon_To_EN+0x711>
movsd  -0x58(%rbp),%xmm1
movsd  -0x180(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x180(%rbp)
ff
movsd  0x88f(%rip),%xmm0
ucomisd -0x180(%rbp),%xmm0
ff
jbe    18bc <GPS_Math_TCylEA_LatLon_To_EN+0x73c>
movsd  -0x58(%rbp),%xmm0
addsd  -0x180(%rbp),%xmm0
ff
movsd  %xmm0,-0x180(%rbp)
ff
movsd  -0x8(%rbp),%xmm0
call   1040 <sin@plt>
movsd  %xmm0,-0x128(%rbp)
ff
movsd  -0x8(%rbp),%xmm0
ucomisd -0x60(%rbp),%xmm0
jne    1909 <GPS_Math_TCylEA_LatLon_To_EN+0x789>
jp     1909 <GPS_Math_TCylEA_LatLon_To_EN+0x789>
movsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0x188(%rbp)
ff
movsd  0x71c(%rip),%xmm0
movsd  %xmm0,-0x190(%rbp)
ff
jmp    1a15 <GPS_Math_TCylEA_LatLon_To_EN+0x895>
movsd  -0x78(%rbp),%xmm0
mulsd  -0x128(%rbp),%xmm0
ff
movsd  %xmm0,-0x118(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0x228(%rbp)
ff
movsd  -0x128(%rbp),%xmm0
ff
movsd  -0x80(%rbp),%xmm1
mulsd  -0x128(%rbp),%xmm1
ff
movsd  -0x128(%rbp),%xmm2
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm1
movsd  0x6aa(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x230(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
movsd  %xmm0,-0x238(%rbp)
ff
movsd  0x682(%rip),%xmm0
subsd  -0x118(%rbp),%xmm0
ff
movsd  0x672(%rip),%xmm1
addsd  -0x118(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
call   1090 <log@plt>
movsd  -0x238(%rbp),%xmm1
ff
movsd  -0x230(%rbp),%xmm2
ff
movaps %xmm0,%xmm3
movsd  -0x228(%rbp),%xmm0
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x188(%rbp)
ff
movsd  -0x188(%rbp),%xmm0
ff
divsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0x190(%rbp)
ff
movsd  -0x190(%rbp),%xmm0
ff
movsd  0x5f3(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1a44 <GPS_Math_TCylEA_LatLon_To_EN+0x8c4>
movsd  0x5e1(%rip),%xmm0
movsd  %xmm0,-0x190(%rbp)
ff
jmp    1a6f <GPS_Math_TCylEA_LatLon_To_EN+0x8ef>
movsd  0x6e4(%rip),%xmm0
ucomisd -0x190(%rbp),%xmm0
ff
jbe    1a6a <GPS_Math_TCylEA_LatLon_To_EN+0x8ea>
movsd  0x6ce(%rip),%xmm0
movsd  %xmm0,-0x190(%rbp)
ff
jmp    1a6f <GPS_Math_TCylEA_LatLon_To_EN+0x8ef>
movsd  -0x190(%rbp),%xmm0
ff
call   1060 <asin@plt>
movsd  %xmm0,-0x198(%rbp)
ff
movsd  -0x198(%rbp),%xmm0
ff
call   1070 <tan@plt>
movsd  %xmm0,-0x240(%rbp)
ff
movsd  -0x180(%rbp),%xmm0
ff
call   1080 <cos@plt>
movaps %xmm0,%xmm1
movsd  -0x240(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
call   1050 <atan@plt>
movsd  %xmm0,-0x1a0(%rbp)
ff
movsd  -0x1a0(%rbp),%xmm0
ff
movaps 0x66f(%rip),%xmm1
pand   %xmm1,%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  0x656(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1b01 <GPS_Math_TCylEA_LatLon_To_EN+0x981>
movsd  -0x1a0(%rbp),%xmm0
ff
movsd  %xmm0,-0x1e0(%rbp)
ff
jmp    1bdd <GPS_Math_TCylEA_LatLon_To_EN+0xa5d>
movsd  -0xf8(%rbp),%xmm0
ff
movsd  %xmm0,-0x258(%rbp)
ff
movsd  0x517(%rip),%xmm0
mulsd  -0x1a0(%rbp),%xmm0
ff
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x258(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x1c8(%rbp)
ff
movsd  -0x100(%rbp),%xmm0
ff
movsd  %xmm0,-0x250(%rbp)
ff
movsd  0x4d3(%rip),%xmm0
mulsd  -0x1a0(%rbp),%xmm0
ff
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x250(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x1d0(%rbp)
ff
movsd  -0x108(%rbp),%xmm0
ff
movsd  %xmm0,-0x248(%rbp)
ff
movsd  0x48f(%rip),%xmm0
mulsd  -0x1a0(%rbp),%xmm0
ff
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x248(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x1d8(%rbp)
ff
movsd  -0x1a0(%rbp),%xmm0
ff
addsd  -0x1c8(%rbp),%xmm0
ff
addsd  -0x1d0(%rbp),%xmm0
ff
addsd  -0x1d8(%rbp),%xmm0
ff
movsd  %xmm0,-0x1e0(%rbp)
ff
movsd  -0x1e0(%rbp),%xmm0
ff
call   1040 <sin@plt>
movsd  %xmm0,-0x1e8(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x2a0(%rbp)
ff
movsd  -0x198(%rbp),%xmm0
ff
call   1080 <cos@plt>
movaps %xmm0,%xmm1
movsd  -0x2a0(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x298(%rbp)
ff
movsd  -0x1e0(%rbp),%xmm0
ff
call   1080 <cos@plt>
movaps %xmm0,%xmm1
movsd  -0x298(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x290(%rbp)
ff
movsd  -0x180(%rbp),%xmm0
ff
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x290(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x278(%rbp)
ff
movsd  -0x110(%rbp),%xmm0
ff
movsd  %xmm0,-0x288(%rbp)
ff
movsd  -0x1a0(%rbp),%xmm0
ff
call   1080 <cos@plt>
movaps %xmm0,%xmm1
movsd  -0x288(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x280(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
mulsd  -0x1e8(%rbp),%xmm0
ff
movsd  -0x1e8(%rbp),%xmm1
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0x341(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  0x43d(%rip),%xmm1
call   10a0 <pow@plt>
movsd  -0x280(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x278(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
addsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0xb8(%rbp),%xmm0
ff
mulsd  -0x1e0(%rbp),%xmm0
ff
movsd  %xmm0,-0x1a8(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0x270(%rbp)
ff
movsd  0x2f0(%rip),%xmm0
mulsd  -0x1e0(%rbp),%xmm0
ff
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x270(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x1b0(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
movsd  %xmm0,-0x268(%rbp)
ff
movsd  0x2ac(%rip),%xmm0
mulsd  -0x1e0(%rbp),%xmm0
ff
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x268(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x1b8(%rbp)
ff
movsd  -0xd0(%rbp),%xmm0
ff
movsd  %xmm0,-0x260(%rbp)
ff
movsd  0x268(%rip),%xmm0
mulsd  -0x1e0(%rbp),%xmm0
ff
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x260(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x1c0(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
movsd  -0x1a8(%rbp),%xmm1
ff
subsd  -0x1b0(%rbp),%xmm1
ff
addsd  -0x1b8(%rbp),%xmm1
ff
subsd  -0x1c0(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x1f0(%rbp)
ff
movsd  -0x110(%rbp),%xmm0
ff
movsd  -0x1f0(%rbp),%xmm2
ff
subsd  -0x98(%rbp),%xmm2
ff
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x2a0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<GPS_Math_TCylEA_LatLon_To_EN>:
push   %r14
push   %rbx
sub    $0x118,%rsp
movsd  %xmm7,(%rsp)
movsd  %xmm6,0x50(%rsp)
movsd  %xmm5,0xe8(%rsp)
movsd  %xmm4,0xe0(%rsp)
movsd  %xmm3,0x10(%rsp)
movsd  %xmm2,0x8(%rsp)
mov    %rsi,%r14
mov    %rdi,%rbx
movsd  %xmm0,0x70(%rsp)
movaps %xmm1,%xmm0
mov    $0x1,%al
call   1030 <GPS_Math_Deg_To_Rad@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movapd %xmm0,0xb0(%rsp)
movsd  0x70(%rsp),%xmm0
mov    $0x1,%al
call   1030 <GPS_Math_Deg_To_Rad@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x80(%rsp)
movsd  0x8(%rsp),%xmm0
mov    $0x1,%al
call   1030 <GPS_Math_Deg_To_Rad@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x10(%rsp),%xmm0
mov    $0x1,%al
call   1030 <GPS_Math_Deg_To_Rad@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0xdd4(%rip),%xmm0
cmp    $0x4,%eax
jl     1235 <GPS_Math_TCylEA_LatLon_To_EN+0xb5>
addsd  %xmm0,%xmm1
movsd  %xmm1,0x98(%rsp)
movsd  0x50(%rsp),%xmm1
mulsd  %xmm1,%xmm1
movsd  (%rsp),%xmm2
mulsd  %xmm2,%xmm2
movapd %xmm1,%xmm0
subsd  %xmm2,%xmm0
divsd  %xmm1,%xmm0
movapd %xmm0,0x10(%rsp)
movapd %xmm0,%xmm2
mulsd  %xmm0,%xmm2
movapd %xmm2,0x60(%rsp)
movapd %xmm0,%xmm1
mulsd  %xmm2,%xmm1
movapd %xmm1,0x30(%rsp)
movsd  0xd81(%rip),%xmm1
call   10a0 <pow@plt>
movsd  %xmm0,0xa0(%rsp)
movsd  0xd73(%rip),%xmm0
subsd  0x10(%rsp),%xmm0
movapd %xmm0,0x20(%rsp)
movsd  0xd57(%rip),%xmm1
call   10a0 <pow@plt>
movsd  0xa0(%rsp),%xmm1
movapd %xmm1,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rsp)
movsd  0xd3c(%rip),%xmm3
movapd %xmm3,%xmm0
subsd  %xmm1,%xmm0
addsd  %xmm3,%xmm1
divsd  %xmm1,%xmm0
call   1090 <log@plt>
movapd %xmm0,0x100(%rsp)
movsd  0xd1e(%rip),%xmm0
divsd  (%rsp),%xmm0
movapd %xmm0,0xf0(%rsp)
movsd  0xd68(%rip),%xmm0
movapd 0x30(%rsp),%xmm6
mulsd  %xmm6,%xmm0
mulsd  0xd5e(%rip),%xmm0
movsd  %xmm0,0xc8(%rsp)
movsd  0xced(%rip),%xmm3
movapd 0x60(%rsp),%xmm7
movapd %xmm7,%xmm2
mulsd  %xmm3,%xmm2
movsd  0xd3f(%rip),%xmm0
movapd 0x10(%rsp),%xmm5
mulsd  %xmm5,%xmm0
addsd  0xcb5(%rip),%xmm0
movsd  0xd2d(%rip),%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movsd  0xd20(%rip),%xmm4
mulsd  %xmm6,%xmm4
movsd  0xd1c(%rip),%xmm0
movapd %xmm5,%xmm1
mulsd  %xmm3,%xmm1
mulsd  0xd14(%rip),%xmm1
mulsd  0xd14(%rip),%xmm2
mulsd  %xmm0,%xmm4
movsd  %xmm4,0xc0(%rsp)
addsd  %xmm1,%xmm2
movapd %xmm2,0x70(%rsp)
movsd  0xcfd(%rip),%xmm1
mulsd  %xmm7,%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x40(%rsp)
movsd  0xcef(%rip),%xmm0
mulsd  %xmm6,%xmm0
movsd  %xmm0,0x48(%rsp)
movsd  0x8(%rsp),%xmm0
addsd  %xmm0,%xmm0
call   1040 <sin@plt>
movsd  %xmm0,0xd8(%rsp)
movsd  0xcd5(%rip),%xmm0
mulsd  0x8(%rsp),%xmm0
call   1040 <sin@plt>
movsd  %xmm0,0x58(%rsp)
movsd  0xcc4(%rip),%xmm0
mulsd  0x8(%rsp),%xmm0
call   1040 <sin@plt>
movsd  %xmm0,0xd0(%rsp)
movapd 0xb0(%rsp),%xmm1
subsd  0x98(%rsp),%xmm1
movsd  0xbbe(%rip),%xmm2
addsd  %xmm1,%xmm2
movsd  0xc92(%rip),%xmm0
cmpltsd %xmm1,%xmm0
andpd  %xmm0,%xmm2
andnpd %xmm1,%xmm0
orpd   %xmm2,%xmm0
movsd  0xc81(%rip),%xmm1
movapd %xmm0,%xmm2
cmpltsd 0xc7c(%rip),%xmm2
addsd  %xmm0,%xmm1
andpd  %xmm2,%xmm1
andnpd %xmm0,%xmm2
orpd   %xmm1,%xmm2
movapd %xmm2,0xb0(%rsp)
movsd  0x80(%rsp),%xmm0
call   1040 <sin@plt>
movsd  0xa0(%rsp),%xmm3
mulsd  %xmm0,%xmm3
movapd 0xc60(%rip),%xmm2
movapd 0x10(%rsp),%xmm1
xorpd  %xmm1,%xmm2
movapd %xmm2,0xa0(%rsp)
movapd %xmm0,%xmm4
mulsd  %xmm1,%xmm4
mulsd  %xmm0,%xmm4
movsd  0xb39(%rip),%xmm2
movapd %xmm2,%xmm1
subsd  %xmm4,%xmm1
shufpd $0x2,0xc38(%rip),%xmm0
unpcklpd 0x20(%rsp),%xmm1
divpd  %xmm1,%xmm0
movapd %xmm0,0x80(%rsp)
movapd %xmm2,%xmm0
movapd %xmm3,%xmm1
subsd  %xmm3,%xmm0
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
call   1090 <log@plt>
movsd  0xaf4(%rip),%xmm1
movapd 0xf0(%rsp),%xmm2
unpcklpd %xmm2,%xmm2
unpcklpd 0x100(%rsp),%xmm0
mulpd  %xmm2,%xmm0
addpd  0x80(%rsp),%xmm0
movapd 0x20(%rsp),%xmm2
unpcklpd %xmm2,%xmm2
mulpd  %xmm0,%xmm2
movapd %xmm2,%xmm0
unpckhpd %xmm2,%xmm0
divsd  %xmm0,%xmm2
ucomisd %xmm1,%xmm2
ja     157d <GPS_Math_TCylEA_LatLon_To_EN+0x3fd>
movapd %xmm2,%xmm0
movsd  0xaad(%rip),%xmm1
ucomisd %xmm2,%xmm1
movapd %xmm2,%xmm1
jbe    157d <GPS_Math_TCylEA_LatLon_To_EN+0x3fd>
movsd  0xa9b(%rip),%xmm1
movsd  (%rsp),%xmm0
addsd  0xc0(%rsp),%xmm0
movsd  %xmm0,(%rsp)
movapd 0x70(%rsp),%xmm0
movsd  0xc8(%rsp),%xmm2
addsd  %xmm2,%xmm0
movapd %xmm0,0x70(%rsp)
movsd  0x40(%rsp),%xmm0
addsd  %xmm2,%xmm0
movsd  %xmm0,0x40(%rsp)
movsd  0x48(%rsp),%xmm0
divsd  0xb01(%rip),%xmm0
movsd  %xmm0,0x48(%rsp)
movapd %xmm1,%xmm0
call   1060 <asin@plt>
movsd  %xmm0,0x80(%rsp)
call   1070 <tan@plt>
movsd  %xmm0,0x20(%rsp)
movapd 0xb0(%rsp),%xmm0
call   1080 <cos@plt>
movsd  0x20(%rsp),%xmm1
divsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
call   1050 <atan@plt>
movapd %xmm0,%xmm2
movapd 0xb19(%rip),%xmm0
andpd  %xmm2,%xmm0
addsd  0xad5(%rip),%xmm0
ucomisd 0xad5(%rip),%xmm0
movapd %xmm2,0x20(%rsp)
ja     1734 <GPS_Math_TCylEA_LatLon_To_EN+0x5b4>
movsd  0x9e9(%rip),%xmm0
movapd 0x60(%rsp),%xmm1
mulsd  %xmm1,%xmm0
movapd 0x10(%rsp),%xmm2
divsd  0x9c9(%rip),%xmm2
divsd  0x9d1(%rip),%xmm0
addsd  %xmm2,%xmm0
movsd  0x9cd(%rip),%xmm4
mulsd  0x9d5(%rip),%xmm1
movapd 0x30(%rsp),%xmm2
mulsd  %xmm2,%xmm4
movsd  0x9d3(%rip),%xmm3
divsd  0x9b3(%rip),%xmm4
mulsd  %xmm2,%xmm3
addsd  %xmm0,%xmm4
movsd  %xmm4,0x10(%rsp)
divsd  0x9ad(%rip),%xmm1
divsd  0x9b5(%rip),%xmm3
mulsd  0x9b5(%rip),%xmm2
addsd  %xmm1,%xmm3
movsd  %xmm3,0x60(%rsp)
divsd  0x9ab(%rip),%xmm2
movapd %xmm2,0x30(%rsp)
movapd 0x20(%rsp),%xmm0
addsd  %xmm0,%xmm0
call   1040 <sin@plt>
mulsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  0x9e2(%rip),%xmm0
mulsd  0x20(%rsp),%xmm0
call   1040 <sin@plt>
mulsd  0x60(%rsp),%xmm0
movsd  %xmm0,0x60(%rsp)
movsd  0x9cb(%rip),%xmm0
mulsd  0x20(%rsp),%xmm0
call   1040 <sin@plt>
movsd  0x10(%rsp),%xmm2
mulsd  0x30(%rsp),%xmm0
addsd  0x20(%rsp),%xmm2
addsd  0x60(%rsp),%xmm2
addsd  %xmm0,%xmm2
movsd  %xmm2,0x10(%rsp)
movsd  0x8(%rsp),%xmm0
mulsd  (%rsp),%xmm0
movsd  0xd8(%rsp),%xmm1
mulsd  0x70(%rsp),%xmm1
subsd  %xmm1,%xmm0
movsd  0x58(%rsp),%xmm1
mulsd  0x40(%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  0xd0(%rsp),%xmm0
mulsd  0x48(%rsp),%xmm0
subsd  %xmm0,%xmm1
mulsd  0x50(%rsp),%xmm1
movsd  %xmm1,0x58(%rsp)
movapd %xmm2,%xmm0
call   1040 <sin@plt>
movsd  %xmm0,0x8(%rsp)
movsd  0x80(%rsp),%xmm0
call   1080 <cos@plt>
mulsd  0x50(%rsp),%xmm0
movsd  %xmm0,0x30(%rsp)
movsd  0x10(%rsp),%xmm0
call   1080 <cos@plt>
mulsd  0x30(%rsp),%xmm0
movsd  %xmm0,0x30(%rsp)
movapd 0xb0(%rsp),%xmm0
call   1040 <sin@plt>
mulsd  0x30(%rsp),%xmm0
movsd  %xmm0,0x30(%rsp)
movaps 0x20(%rsp),%xmm0
call   1080 <cos@plt>
movsd  %xmm0,0x20(%rsp)
movapd 0xa0(%rsp),%xmm0
movsd  0x8(%rsp),%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm1,%xmm0
addsd  0x800(%rip),%xmm0
movsd  0x7f0(%rip),%xmm1
call   10a0 <pow@plt>
mulsd  0x20(%rsp),%xmm0
movsd  0x30(%rsp),%xmm1
divsd  %xmm0,%xmm1
addsd  0xe0(%rsp),%xmm1
movsd  %xmm1,(%rbx)
movsd  (%rsp),%xmm0
movsd  0x10(%rsp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,(%rsp)
movapd %xmm1,%xmm0
addsd  %xmm1,%xmm0
call   1040 <sin@plt>
mulsd  0x70(%rsp),%xmm0
movsd  (%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movsd  0x859(%rip),%xmm0
mulsd  0x10(%rsp),%xmm0
call   1040 <sin@plt>
mulsd  0x40(%rsp),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x10(%rsp),%xmm0
mulsd  0x838(%rip),%xmm0
call   1040 <sin@plt>
mulsd  0x48(%rsp),%xmm0
movsd  (%rsp),%xmm1
subsd  %xmm0,%xmm1
mulsd  0x50(%rsp),%xmm1
subsd  0x58(%rsp),%xmm1
addsd  0xe8(%rsp),%xmm1
movsd  %xmm1,(%r14)
add    $0x118,%rsp
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @GPS_Math_TCylEA_LatLon_To_EN(double %0, double %1, double* %2, double* %3, double %4, double %5, double %6, double %7, double %8, double %9) {
%11 = alloca double
%12 = alloca double
%13 = alloca double*
%14 = alloca double*
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double
%26 = alloca double
%27 = alloca double
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca double
%32 = alloca double
%33 = alloca double
%34 = alloca double
%35 = alloca double
%36 = alloca double
%37 = alloca double
%38 = alloca double
%39 = alloca double
%40 = alloca double
%41 = alloca double
%42 = alloca double
%43 = alloca double
%44 = alloca double
%45 = alloca double
%46 = alloca double
%47 = alloca double
%48 = alloca double
%49 = alloca double
%50 = alloca double
%51 = alloca double
%52 = alloca double
%53 = alloca double
%54 = alloca double
%55 = alloca double
%56 = alloca double
%57 = alloca double
%58 = alloca double
%59 = alloca double
%60 = alloca double
%61 = alloca double
%62 = alloca double
%63 = alloca double
%64 = alloca double
%65 = alloca double
%66 = alloca double
%67 = alloca double
%68 = alloca double
%69 = alloca double
%70 = alloca double
%71 = alloca double
%72 = alloca double
store double %0, double* %11
store double %1, double* %12
store double* %2, double** %13
store double* %3, double** %14
store double %4, double* %15
store double %5, double* %16
store double %6, double* %17
store double %7, double* %18
store double %8, double* %19
store double %9, double* %20
store double 1.000000e+00, double* %44
%73 = load double, double* %12
%74 = call i32 (double, ...) bitcast (i32 (...)* @GPS_Math_Deg_To_Rad to i32 (double, ...)*)(double %73)
%75 = sitofp i32 %74 to double
store double %75, double* %12
%76 = load double, double* %11
%77 = call i32 (double, ...) bitcast (i32 (...)* @GPS_Math_Deg_To_Rad to i32 (double, ...)*)(double %76)
%78 = sitofp i32 %77 to double
store double %78, double* %11
%79 = load double, double* %15
%80 = call i32 (double, ...) bitcast (i32 (...)* @GPS_Math_Deg_To_Rad to i32 (double, ...)*)(double %79)
%81 = sitofp i32 %80 to double
store double %81, double* %15
%82 = load double, double* %16
%83 = call i32 (double, ...) bitcast (i32 (...)* @GPS_Math_Deg_To_Rad to i32 (double, ...)*)(double %82)
%84 = sitofp i32 %83 to double
store double %84, double* %16
store double 0x401921FB5444261E, double* %21
store double 0x3FF921FB5444261E, double* %22
%85 = load double, double* %16
%86 = fcmp ogt double %85, 0x400921FB5444261E
br i1 %86, label %87, label %91
87: ; preds = %10
%88 = load double, double* %21
%89 = load double, double* %16
%90 = fsub double %89, %88
store double %90, double* %16
br label %91
91: ; preds = %87, %10
%92 = load double, double* %19
%93 = load double, double* %19
%94 = fmul double %92, %93
store double %94, double* %23
%95 = load double, double* %20
%96 = load double, double* %20
%97 = fmul double %95, %96
store double %97, double* %24
%98 = load double, double* %23
%99 = load double, double* %24
%100 = fsub double %98, %99
%101 = load double, double* %23
%102 = fdiv double %100, %101
store double %102, double* %26
%103 = load double, double* %26
%104 = load double, double* %26
%105 = fmul double %103, %104
store double %105, double* %27
%106 = load double, double* %26
%107 = load double, double* %27
%108 = fmul double %106, %107
store double %108, double* %28
%109 = load double, double* %26
%110 = call double @pow(double %109, double 5.000000e-01)
store double %110, double* %25
%111 = load double, double* %26
%112 = fsub double 1.000000e+00, %111
store double %112, double* %31
%113 = load double, double* %31
%114 = call double @pow(double %113, double 5.000000e-01)
store double %114, double* %46
%115 = load double, double* %25
%116 = fmul double 2.000000e+00, %115
%117 = fdiv double 1.000000e+00, %116
store double %117, double* %32
%118 = load double, double* %22
%119 = call double @sin(double %118)
store double %119, double* %47
%120 = load double, double* %25
%121 = load double, double* %47
%122 = fmul double %120, %121
store double %122, double* %45
%123 = load double, double* %31
%124 = load double, double* %47
%125 = load double, double* %26
%126 = load double, double* %47
%127 = fmul double %125, %126
%128 = load double, double* %47
%129 = fneg double %127
%130 = call double @llvm.fmuladd.f64(double %129, double %128, double 1.000000e+00)
%131 = fdiv double %124, %130
%132 = load double, double* %32
%133 = load double, double* %45
%134 = fsub double 1.000000e+00, %133
%135 = load double, double* %45
%136 = fadd double 1.000000e+00, %135
%137 = fdiv double %134, %136
%138 = call double @log(double %137)
%139 = fneg double %132
%140 = call double @llvm.fmuladd.f64(double %139, double %138, double %131)
%141 = fmul double %123, %140
store double %141, double* %30
%142 = load double, double* %26
%143 = fdiv double %142, 3.000000e+00
%144 = load double, double* %27
%145 = fmul double 3.100000e+01, %144
%146 = fdiv double %145, 1.800000e+02
%147 = fadd double %143, %146
%148 = load double, double* %28
%149 = fmul double 5.170000e+02, %148
%150 = fdiv double %149, 5.040000e+03
%151 = fadd double %147, %150
store double %151, double* %41
%152 = load double, double* %27
%153 = fmul double 2.300000e+01, %152
%154 = fdiv double %153, 3.600000e+02
%155 = load double, double* %28
%156 = fmul double 2.510000e+02, %155
%157 = fdiv double %156, 3.780000e+03
%158 = fadd double %154, %157
store double %158, double* %42
%159 = load double, double* %28
%160 = fmul double 7.610000e+02, %159
%161 = fdiv double %160, 4.536000e+04
store double %161, double* %43
%162 = load double, double* %46
%163 = fsub double 1.000000e+00, %162
%164 = load double, double* %46
%165 = fadd double 1.000000e+00, %164
%166 = fdiv double %163, %165
store double %166, double* %54
%167 = load double, double* %54
%168 = load double, double* %54
%169 = fmul double %167, %168
store double %169, double* %55
%170 = load double, double* %55
%171 = load double, double* %54
%172 = fmul double %170, %171
store double %172, double* %56
%173 = load double, double* %56
%174 = load double, double* %54
%175 = fmul double %173, %174
store double %175, double* %57
%176 = load double, double* %54
%177 = fmul double 3.000000e+00, %176
%178 = fdiv double %177, 2.000000e+00
%179 = load double, double* %56
%180 = fmul double 2.700000e+01, %179
%181 = fdiv double %180, 3.200000e+01
%182 = fsub double %178, %181
store double %182, double* %37
%183 = load double, double* %55
%184 = fmul double 2.100000e+01, %183
%185 = fdiv double %184, 1.600000e+01
%186 = load double, double* %57
%187 = fmul double 5.500000e+01, %186
%188 = fdiv double %187, 3.200000e+01
%189 = fsub double %185, %188
store double %189, double* %38
%190 = load double, double* %56
%191 = fmul double 1.510000e+02, %190
%192 = fdiv double %191, 9.600000e+01
store double %192, double* %39
%193 = load double, double* %57
%194 = fmul double 1.097000e+03, %193
%195 = fdiv double %194, 5.120000e+02
store double %195, double* %40
%196 = load double, double* %28
%197 = fmul double 4.500000e+01, %196
%198 = fdiv double %197, 1.024000e+03
store double %198, double* %48
%199 = load double, double* %27
%200 = fmul double 3.000000e+00, %199
store double %200, double* %49
%201 = load double, double* %26
%202 = fdiv double %201, 4.000000e+00
%203 = fsub double 1.000000e+00, %202
%204 = load double, double* %49
%205 = fdiv double %204, 6.400000e+01
%206 = fsub double %203, %205
%207 = load double, double* %28
%208 = fmul double 5.000000e+00, %207
%209 = fdiv double %208, 2.560000e+02
%210 = fsub double %206, %209
store double %210, double* %33
%211 = load double, double* %26
%212 = fmul double 3.000000e+00, %211
%213 = fdiv double %212, 8.000000e+00
%214 = load double, double* %49
%215 = fdiv double %214, 3.200000e+01
%216 = fadd double %213, %215
%217 = load double, double* %48
%218 = fadd double %216, %217
store double %218, double* %34
%219 = load double, double* %27
%220 = fmul double 1.500000e+01, %219
%221 = fdiv double %220, 2.560000e+02
%222 = load double, double* %48
%223 = fadd double %221, %222
store double %223, double* %35
%224 = load double, double* %28
%225 = fmul double 3.500000e+01, %224
%226 = fdiv double %225, 3.072000e+03
store double %226, double* %36
%227 = load double, double* %33
%228 = load double, double* %15
%229 = fmul double %227, %228
store double %229, double* %50
%230 = load double, double* %34
%231 = load double, double* %15
%232 = fmul double 2.000000e+00, %231
%233 = call double @sin(double %232)
%234 = fmul double %230, %233
store double %234, double* %51
%235 = load double, double* %35
%236 = load double, double* %15
%237 = fmul double 4.000000e+00, %236
%238 = call double @sin(double %237)
%239 = fmul double %235, %238
store double %239, double* %52
%240 = load double, double* %36
%241 = load double, double* %15
%242 = fmul double 6.000000e+00, %241
%243 = call double @sin(double %242)
%244 = fmul double %240, %243
store double %244, double* %53
%245 = load double, double* %19
%246 = load double, double* %50
%247 = load double, double* %51
%248 = fsub double %246, %247
%249 = load double, double* %52
%250 = fadd double %248, %249
%251 = load double, double* %53
%252 = fsub double %250, %251
%253 = fmul double %245, %252
store double %253, double* %29
%254 = load double, double* %12
%255 = load double, double* %16
%256 = fsub double %254, %255
store double %256, double* %58
%257 = load double, double* %58
%258 = fcmp ogt double %257, 0x400921FB5444261E
br i1 %258, label %259, label %263
259: ; preds = %91
%260 = load double, double* %21
%261 = load double, double* %58
%262 = fsub double %261, %260
store double %262, double* %58
br label %263
263: ; preds = %259, %91
%264 = load double, double* %58
%265 = fcmp olt double %264, 0xC00921FB5444261E
br i1 %265, label %266, label %270
266: ; preds = %263
%267 = load double, double* %21
%268 = load double, double* %58
%269 = fadd double %268, %267
store double %269, double* %58
br label %270
270: ; preds = %266, %263
%271 = load double, double* %11
%272 = call double @sin(double %271)
store double %272, double* %47
%273 = load double, double* %11
%274 = load double, double* %22
%275 = fcmp oeq double %273, %274
br i1 %275, label %276, label %278
276: ; preds = %270
%277 = load double, double* %30
store double %277, double* %59
store double 1.000000e+00, double* %60
br label %304
278: ; preds = %270
%279 = load double, double* %25
%280 = load double, double* %47
%281 = fmul double %279, %280
store double %281, double* %45
%282 = load double, double* %31
%283 = load double, double* %47
%284 = load double, double* %26
%285 = load double, double* %47
%286 = fmul double %284, %285
%287 = load double, double* %47
%288 = fneg double %286
%289 = call double @llvm.fmuladd.f64(double %288, double %287, double 1.000000e+00)
%290 = fdiv double %283, %289
%291 = load double, double* %32
%292 = load double, double* %45
%293 = fsub double 1.000000e+00, %292
%294 = load double, double* %45
%295 = fadd double 1.000000e+00, %294
%296 = fdiv double %293, %295
%297 = call double @log(double %296)
%298 = fneg double %291
%299 = call double @llvm.fmuladd.f64(double %298, double %297, double %290)
%300 = fmul double %282, %299
store double %300, double* %59
%301 = load double, double* %59
%302 = load double, double* %30
%303 = fdiv double %301, %302
store double %303, double* %60
br label %304
304: ; preds = %278, %276
%305 = load double, double* %60
%306 = fcmp ogt double %305, 1.000000e+00
br i1 %306, label %307, label %308
307: ; preds = %304
store double 1.000000e+00, double* %60
br label %313
308: ; preds = %304
%309 = load double, double* %60
%310 = fcmp olt double %309, -1.000000e+00
br i1 %310, label %311, label %312
311: ; preds = %308
store double -1.000000e+00, double* %60
br label %312
312: ; preds = %311, %308
br label %313
313: ; preds = %312, %307
%314 = load double, double* %60
%315 = call double @asin(double %314)
store double %315, double* %61
%316 = load double, double* %61
%317 = call double @tan(double %316)
%318 = load double, double* %58
%319 = call double @cos(double %318)
%320 = fdiv double %317, %319
%321 = call double @atan(double %320)
store double %321, double* %62
%322 = load double, double* %62
%323 = call double @llvm.fabs.f64(double %322)
%324 = load double, double* %22
%325 = fsub double %323, %324
%326 = fcmp ogt double %325, 1.000000e-08
br i1 %326, label %327, label %329
327: ; preds = %313
%328 = load double, double* %62
store double %328, double* %70
br label %352
329: ; preds = %313
%330 = load double, double* %41
%331 = load double, double* %62
%332 = fmul double 2.000000e+00, %331
%333 = call double @sin(double %332)
%334 = fmul double %330, %333
store double %334, double* %67
%335 = load double, double* %42
%336 = load double, double* %62
%337 = fmul double 4.000000e+00, %336
%338 = call double @sin(double %337)
%339 = fmul double %335, %338
store double %339, double* %68
%340 = load double, double* %43
%341 = load double, double* %62
%342 = fmul double 6.000000e+00, %341
%343 = call double @sin(double %342)
%344 = fmul double %340, %343
store double %344, double* %69
%345 = load double, double* %62
%346 = load double, double* %67
%347 = fadd double %345, %346
%348 = load double, double* %68
%349 = fadd double %347, %348
%350 = load double, double* %69
%351 = fadd double %349, %350
store double %351, double* %70
br label %352
352: ; preds = %329, %327
%353 = load double, double* %70
%354 = call double @sin(double %353)
store double %354, double* %71
%355 = load double, double* %19
%356 = load double, double* %61
%357 = call double @cos(double %356)
%358 = fmul double %355, %357
%359 = load double, double* %70
%360 = call double @cos(double %359)
%361 = fmul double %358, %360
%362 = load double, double* %58
%363 = call double @sin(double %362)
%364 = fmul double %361, %363
%365 = load double, double* %44
%366 = load double, double* %62
%367 = call double @cos(double %366)
%368 = fmul double %365, %367
%369 = load double, double* %26
%370 = load double, double* %71
%371 = fmul double %369, %370
%372 = load double, double* %71
%373 = fneg double %371
%374 = call double @llvm.fmuladd.f64(double %373, double %372, double 1.000000e+00)
%375 = call double @pow(double %374, double 5.000000e-01)
%376 = fmul double %368, %375
%377 = fdiv double %364, %376
%378 = load double, double* %17
%379 = fadd double %377, %378
%380 = load double*, double** %13
store double %379, double* %380
%381 = load double, double* %33
%382 = load double, double* %70
%383 = fmul double %381, %382
store double %383, double* %63
%384 = load double, double* %34
%385 = load double, double* %70
%386 = fmul double 2.000000e+00, %385
%387 = call double @sin(double %386)
%388 = fmul double %384, %387
store double %388, double* %64
%389 = load double, double* %35
%390 = load double, double* %70
%391 = fmul double 4.000000e+00, %390
%392 = call double @sin(double %391)
%393 = fmul double %389, %392
store double %393, double* %65
%394 = load double, double* %36
%395 = load double, double* %70
%396 = fmul double 6.000000e+00, %395
%397 = call double @sin(double %396)
%398 = fmul double %394, %397
store double %398, double* %66
%399 = load double, double* %19
%400 = load double, double* %63
%401 = load double, double* %64
%402 = fsub double %400, %401
%403 = load double, double* %65
%404 = fadd double %402, %403
%405 = load double, double* %66
%406 = fsub double %404, %405
%407 = fmul double %399, %406
store double %407, double* %72
%408 = load double, double* %44
%409 = load double, double* %72
%410 = load double, double* %29
%411 = fsub double %409, %410
%412 = load double, double* %18
%413 = call double @llvm.fmuladd.f64(double %408, double %411, double %412)
%414 = load double*, double** %14
store double %413, double* %414
ret void
}

LLVM IR (-O2):
define void @GPS_Math_TCylEA_LatLon_To_EN(double %0, double %1, double* nocapture writeonly %2, double* nocapture writeonly %3, double %4, double %5, double %6, double %7, double %8, double %9) {
%11 = tail call i32 (double, ...) bitcast (i32 (...)* @GPS_Math_Deg_To_Rad to i32 (double, ...)*)(double %1)
%12 = sitofp i32 %11 to double
%13 = tail call i32 (double, ...) bitcast (i32 (...)* @GPS_Math_Deg_To_Rad to i32 (double, ...)*)(double %0)
%14 = sitofp i32 %13 to double
%15 = tail call i32 (double, ...) bitcast (i32 (...)* @GPS_Math_Deg_To_Rad to i32 (double, ...)*)(double %4)
%16 = sitofp i32 %15 to double
%17 = tail call i32 (double, ...) bitcast (i32 (...)* @GPS_Math_Deg_To_Rad to i32 (double, ...)*)(double %5)
%18 = sitofp i32 %17 to double
%19 = icmp sgt i32 %17, 3
%20 = fadd double %18, 0xC01921FB5444261E
%21 = select i1 %19, double %20, double %18
%22 = fmul double %8, %8
%23 = fmul double %9, %9
%24 = fsub double %22, %23
%25 = fdiv double %24, %22
%26 = fmul double %25, %25
%27 = fmul double %25, %26
%28 = tail call double @pow(double %25, double 5.000000e-01)
%29 = fsub double 1.000000e+00, %25
%30 = tail call double @pow(double %29, double 5.000000e-01)
%31 = fmul double %28, 2.000000e+00
%32 = fsub double 1.000000e+00, %28
%33 = fadd double %28, 1.000000e+00
%34 = fdiv double %32, %33
%35 = tail call double @log(double %34)
%36 = fdiv double -1.000000e+00, %31
%37 = fdiv double %25, 3.000000e+00
%38 = fmul double %26, 3.100000e+01
%39 = fdiv double %38, 1.800000e+02
%40 = fadd double %37, %39
%41 = fmul double %27, 5.170000e+02
%42 = fdiv double %41, 5.040000e+03
%43 = fadd double %40, %42
%44 = fmul double %26, 2.300000e+01
%45 = fdiv double %44, 3.600000e+02
%46 = fmul double %27, 2.510000e+02
%47 = fdiv double %46, 3.780000e+03
%48 = fadd double %45, %47
%49 = fmul double %27, 7.610000e+02
%50 = fdiv double %49, 4.536000e+04
%51 = fmul double %27, 4.500000e+01
%52 = fmul double %51, 0x3F50000000000000
%53 = fmul double %26, 3.000000e+00
%54 = fmul double %25, 2.500000e-01
%55 = fsub double 1.000000e+00, %54
%56 = fmul double %53, 1.562500e-02
%57 = fsub double %55, %56
%58 = fmul double %27, 5.000000e+00
%59 = fmul double %58, 3.906250e-03
%60 = fsub double %57, %59
%61 = fmul double %25, 3.000000e+00
%62 = fmul double %61, 1.250000e-01
%63 = fmul double %53, 3.125000e-02
%64 = fadd double %62, %63
%65 = fadd double %64, %52
%66 = fmul double %26, 1.500000e+01
%67 = fmul double %66, 3.906250e-03
%68 = fadd double %67, %52
%69 = fmul double %27, 3.500000e+01
%70 = fdiv double %69, 3.072000e+03
%71 = fmul double %16, 2.000000e+00
%72 = tail call double @sin(double %71)
%73 = fmul double %16, 4.000000e+00
%74 = tail call double @sin(double %73)
%75 = fmul double %16, 6.000000e+00
%76 = tail call double @sin(double %75)
%77 = fsub double %12, %21
%78 = fcmp ogt double %77, 0x400921FB5444261E
%79 = fadd double %77, 0xC01921FB5444261E
%80 = select i1 %78, double %79, double %77
%81 = fcmp olt double %80, 0xC00921FB5444261E
%82 = fadd double %80, 0x401921FB5444261E
%83 = select i1 %81, double %82, double %80
%84 = tail call double @sin(double %14)
%85 = fmul double %28, %84
%86 = fneg double %25
%87 = fmul double %84, %86
%88 = tail call double @llvm.fmuladd.f64(double %87, double %84, double 1.000000e+00)
%89 = insertelement <2 x double> <double poison, double 1.000000e+00>, double %84, i64 0
%90 = insertelement <2 x double> poison, double %88, i64 0
%91 = insertelement <2 x double> %90, double %29, i64 1
%92 = fdiv <2 x double> %89, %91
%93 = fsub double 1.000000e+00, %85
%94 = fadd double %85, 1.000000e+00
%95 = fdiv double %93, %94
%96 = tail call double @log(double %95)
%97 = insertelement <2 x double> poison, double %36, i64 0
%98 = shufflevector <2 x double> %97, <2 x double> poison, <2 x i32> zeroinitializer
%99 = insertelement <2 x double> poison, double %96, i64 0
%100 = insertelement <2 x double> %99, double %35, i64 1
%101 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %98, <2 x double> %100, <2 x double> %92)
%102 = insertelement <2 x double> poison, double %29, i64 0
%103 = shufflevector <2 x double> %102, <2 x double> poison, <2 x i32> zeroinitializer
%104 = fmul <2 x double> %103, %101
%105 = extractelement <2 x double> %104, i64 0
%106 = extractelement <2 x double> %104, i64 1
%107 = fdiv double %105, %106
%108 = fcmp ogt double %107, 1.000000e+00
br i1 %108, label %112, label %109
109: ; preds = %10
%110 = fcmp olt double %107, -1.000000e+00
br i1 %110, label %111, label %112
111: ; preds = %109
br label %112
112: ; preds = %10, %109, %111
%113 = phi double [ -1.000000e+00, %111 ], [ %107, %109 ], [ 1.000000e+00, %10 ]
%114 = tail call double @asin(double %113)
%115 = tail call double @tan(double %114)
%116 = tail call double @cos(double %83)
%117 = fdiv double %115, %116
%118 = tail call double @atan(double %117)
%119 = tail call double @llvm.fabs.f64(double %118)
%120 = fadd double %119, 0xBFF921FB5444261E
%121 = fcmp ogt double %120, 1.000000e-08
br i1 %121, label %135, label %122
122: ; preds = %112
%123 = fmul double %118, 2.000000e+00
%124 = tail call double @sin(double %123)
%125 = fmul double %43, %124
%126 = fmul double %118, 4.000000e+00
%127 = tail call double @sin(double %126)
%128 = fmul double %48, %127
%129 = fmul double %118, 6.000000e+00
%130 = tail call double @sin(double %129)
%131 = fmul double %50, %130
%132 = fadd double %118, %125
%133 = fadd double %132, %128
%134 = fadd double %133, %131
br label %135
135: ; preds = %112, %122
%136 = phi double [ %134, %122 ], [ %118, %112 ]
%137 = fmul double %60, %16
%138 = fmul double %65, %72
%139 = fsub double %137, %138
%140 = fmul double %68, %74
%141 = fadd double %139, %140
%142 = fmul double %70, %76
%143 = fsub double %141, %142
%144 = fmul double %143, %8
%145 = tail call double @sin(double %136)
%146 = tail call double @cos(double %114)
%147 = fmul double %146, %8
%148 = tail call double @cos(double %136)
%149 = fmul double %147, %148
%150 = tail call double @sin(double %83)
%151 = fmul double %149, %150
%152 = tail call double @cos(double %118)
%153 = fmul double %145, %86
%154 = tail call double @llvm.fmuladd.f64(double %153, double %145, double 1.000000e+00)
%155 = tail call double @pow(double %154, double 5.000000e-01)
%156 = fmul double %152, %155
%157 = fdiv double %151, %156
%158 = fadd double %157, %6
store double %158, double* %2
%159 = fmul double %60, %136
%160 = fmul double %136, 2.000000e+00
%161 = tail call double @sin(double %160)
%162 = fmul double %65, %161
%163 = fmul double %136, 4.000000e+00
%164 = tail call double @sin(double %163)
%165 = fmul double %68, %164
%166 = fmul double %136, 6.000000e+00
%167 = tail call double @sin(double %166)
%168 = fmul double %70, %167
%169 = fsub double %159, %162
%170 = fadd double %169, %165
%171 = fsub double %170, %168
%172 = fmul double %171, %8
%173 = fsub double %172, %144
%174 = fadd double %173, %7
store double %174, double* %3
ret void
}