#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double atan(double);
double atan2(double, double);
double cos(double);
double pow(double, double);
double sin(double);
double sqrt(double);
double tan(double);
void DirectVincenty(double lat_in, double lon_in, double alpha, double s,
                    double *lat_out, double *lon_out, double a, double f) {
  if (a == 0.0 || f == 0.0) {
    a = 6378137.0;
    f = 1.0 / 298.257223563;
  }
  double b = (1.0 - f) * a;
  double sin_a1 = sin(alpha);
  double cos_a1 = cos(alpha);
  double tan_u1 = (1.0 - f) * tan(lat_in);
  double u1 = atan(tan_u1);
  double sin_u1 = sin(u1);
  double cos_u1 = cos(u1);
  double sigma1 = atan2(tan_u1, cos_a1);
  double sin_alpha = cos_u1 * sin_a1;
  double cos2_alpha = (1.0 - sin_alpha) * (1.0 + sin_alpha);
  double u2 = cos2_alpha * (a * a - b * b) / (b * b);
  double sqrt_1pu2 = sqrt(1.0 + u2);
  double k1 = (sqrt_1pu2 - 1.0) / (sqrt_1pu2 + 1.0);
  double k12 = k1 * k1;
  double A = (1.0 + 0.25 * k12) / (1.0 - k1);
  double B = k1 * (1.0 - 3.0 / 8.0 * k12);
  double sigma = s / (b * A);
  double sigmaOld, do_sigma_m, d_sigma;
  do {
    sigmaOld = sigma;
    do_sigma_m = 2.0 * sigma1 + sigma;
    d_sigma =
        B * sin(sigma) *
        (cos(do_sigma_m) +
         0.25 * B *
             (cos(sigma) * (-1.0 + 2.0 * cos(do_sigma_m) * cos(do_sigma_m)) -
              (1.0 / 6.0) * B * cos(do_sigma_m) *
                  (-3.0 + 4.0 * sin(sigma) * sin(sigma)) *
                  (-3.0 + 4.0 * cos(do_sigma_m) * cos(do_sigma_m))));
    sigma = s / (b * A) + d_sigma;
  } while ((sigmaOld - sigma) > 1.0e-6 || (sigmaOld - sigma) < -1.0e-6);
  do_sigma_m = 2.0 * sigma1 + sigma;
  *lat_out = atan2(
      sin_u1 * cos(sigma) + cos_u1 * sin(sigma) * cos_a1,
      (1.0 - f) *
          sqrt(sin_alpha * sin_alpha +
               pow((sin_u1 * sin(sigma) - cos_u1 * cos(sigma) * cos_a1), 2.0)));
  double lambda = atan2(sin(sigma) * sin_a1,
                        cos_u1 * cos(sigma) - sin_u1 * sin(sigma) * cos_a1);
  double C = (f / 16.0) * cos2_alpha * (4.0 + f * (4.0 - 3.0 * cos2_alpha));
  double dLon =
      lambda -
      (1.0 - C) * f * sin_alpha *
          (sigma + C * sin(sigma) *
                       (cos(do_sigma_m) +
                        C * cos(sigma) *
                            (-1.0 + 2.0 * cos(do_sigma_m) * cos(do_sigma_m))));
  *lon_out = lon_in + dLon;
}


x86 assembly (clang, -O0):
<DirectVincenty>:
push   %rbp
mov    %rsp,%rbp
sub    $0x220,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
movsd  %xmm3,-0x20(%rbp)
mov    %rdi,-0x28(%rbp)
mov    %rsi,-0x30(%rbp)
movsd  %xmm4,-0x38(%rbp)
movsd  %xmm5,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11be <DirectVincenty+0x4e>
jp     11be <DirectVincenty+0x4e>
jmp    11d6 <DirectVincenty+0x66>
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11f0 <DirectVincenty+0x80>
jp     11f0 <DirectVincenty+0x80>
movsd  0xe2a(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe15(%rip),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xe20(%rip),%xmm0
subsd  -0x40(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x18(%rbp),%xmm0
call   1040 <sin@plt>
movsd  %xmm0,-0x50(%rbp)
movsd  -0x18(%rbp),%xmm0
call   1070 <cos@plt>
movsd  %xmm0,-0x58(%rbp)
movsd  0xdeb(%rip),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x100(%rbp)
ff
movsd  -0x8(%rbp),%xmm0
call   1060 <tan@plt>
movaps %xmm0,%xmm1
movsd  -0x100(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
call   1050 <atan@plt>
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
call   1040 <sin@plt>
movsd  %xmm0,-0x70(%rbp)
movsd  -0x68(%rbp),%xmm0
call   1070 <cos@plt>
movsd  %xmm0,-0x78(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
call   1030 <atan2@plt>
movsd  %xmm0,-0x80(%rbp)
movsd  -0x78(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  0xd65(%rip),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
movsd  0xd55(%rip),%xmm1
addsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
movsd  -0x38(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm3
movsd  -0x48(%rbp),%xmm2
mulsd  -0x48(%rbp),%xmm2
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x48(%rbp),%xmm1
mulsd  -0x48(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0xce4(%rip),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
call   1090 <sqrt@plt>
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0xa0(%rbp),%xmm0
ff
movsd  0xcbf(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  0xcb3(%rip),%xmm1
addsd  -0xa0(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
mulsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
movsd  0xc7f(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xc6b(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  0xc5f(%rip),%xmm1
subsd  -0xa8(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
movsd  -0xb0(%rbp),%xmm1
ff
movsd  0xc2b(%rip),%xmm2
mulsd  %xmm2,%xmm1
movsd  0xc27(%rip),%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xc0(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
mulsd  -0xb8(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  -0xc8(%rbp),%xmm1
ff
addsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0x160(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x160(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
movsd  -0xd8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movsd  %xmm0,-0x110(%rbp)
ff
movsd  0xb83(%rip),%xmm0
mulsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0x118(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movsd  %xmm0,-0x128(%rbp)
ff
movsd  -0xd8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movaps %xmm0,%xmm1
movsd  0xb6e(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x158(%rbp)
ff
movsd  -0xd8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movaps %xmm0,%xmm1
movsd  -0x158(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  0xb36(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
movsd  0xb1a(%rip),%xmm0
mulsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0x150(%rbp)
ff
movsd  -0xd8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movaps %xmm0,%xmm1
movsd  -0x150(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x140(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  0xac6(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x148(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
call   1040 <sin@plt>
movsd  -0x148(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x140(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm1
movsd  0xa86(%rip),%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
movsd  -0xd8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movaps %xmm0,%xmm1
movsd  0xa66(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x138(%rbp)
ff
movsd  -0xd8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movsd  -0x138(%rbp),%xmm6
ff
movsd  -0x130(%rbp),%xmm4
ff
movsd  -0x128(%rbp),%xmm3
ff
movsd  -0x120(%rbp),%xmm5
ff
movsd  -0x118(%rbp),%xmm1
ff
movsd  -0x110(%rbp),%xmm2
ff
movaps %xmm0,%xmm7
movsd  -0x108(%rbp),%xmm0
ff
mulsd  %xmm7,%xmm6
movsd  0x9fe(%rip),%xmm7
addsd  %xmm7,%xmm6
mulsd  %xmm6,%xmm4
movq   %xmm4,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm4
mulsd  %xmm5,%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  -0x20(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
mulsd  -0xb8(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
addsd  -0xe0(%rbp),%xmm0
ff
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0xd0(%rbp),%xmm0
ff
subsd  -0xc8(%rbp),%xmm0
ff
mov    $0x1,%al
movsd  0x9ab(%rip),%xmm1
ucomisd %xmm1,%xmm0
mov    %al,-0x161(%rbp)
ja     16da <DirectVincenty+0x56a>
movsd  -0xd0(%rbp),%xmm1
ff
subsd  -0xc8(%rbp),%xmm1
ff
movsd  0x98b(%rip),%xmm0
ucomisd %xmm1,%xmm0
seta   %al
mov    %al,-0x161(%rbp)
mov    -0x161(%rbp),%al
test   $0x1,%al
jne    141f <DirectVincenty+0x2af>
movsd  -0x80(%rbp),%xmm0
movsd  -0xc8(%rbp),%xmm1
ff
addsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x210(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movsd  %xmm0,-0x218(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x220(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
call   1040 <sin@plt>
movsd  -0x220(%rbp),%xmm1
ff
movsd  -0x218(%rbp),%xmm2
ff
movaps %xmm0,%xmm3
movsd  -0x210(%rbp),%xmm0
ff
mulsd  %xmm3,%xmm1
mulsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1d8(%rbp)
ff
movsd  0x89b(%rip),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x1e0(%rbp)
ff
movsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x1e8(%rbp)
ff
movsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x1f0(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x1f8(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
call   1040 <sin@plt>
movsd  %xmm0,-0x200(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x208(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movsd  -0x208(%rbp),%xmm1
ff
movsd  -0x200(%rbp),%xmm2
ff
movaps %xmm0,%xmm3
movsd  -0x1f8(%rbp),%xmm0
mulsd  %xmm3,%xmm1
mulsd  -0x58(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0x817(%rip),%xmm1
call   1080 <pow@plt>
movsd  -0x1f0(%rbp),%xmm2
ff
movaps %xmm0,%xmm1
movsd  -0x1e8(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1090 <sqrt@plt>
movsd  -0x1e0(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x1d8(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm1
call   1030 <atan2@plt>
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0xc8(%rbp),%xmm0
ff
call   1040 <sin@plt>
mulsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x1b8(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x1c8(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movsd  %xmm0,-0x1c0(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x1d0(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
call   1040 <sin@plt>
movsd  -0x1d0(%rbp),%xmm2
ff
movsd  -0x1c8(%rbp),%xmm1
ff
movsd  -0x1c0(%rbp),%xmm3
ff
movaps %xmm0,%xmm4
movsd  -0x1b8(%rbp),%xmm0
ff
mulsd  %xmm4,%xmm2
mulsd  -0x58(%rbp),%xmm2
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
call   1030 <atan2@plt>
movsd  %xmm0,-0xe8(%rbp)
ff
movsd  -0x40(%rbp),%xmm0
movsd  0x72b(%rip),%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
movsd  -0x40(%rbp),%xmm1
movsd  -0x90(%rbp),%xmm2
ff
movsd  0x6d2(%rip),%xmm3
mulsd  %xmm3,%xmm2
movsd  0x6ce(%rip),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
movsd  0x6be(%rip),%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  -0xe8(%rbp),%xmm0
ff
movsd  %xmm0,-0x178(%rbp)
ff
movsd  0x67e(%rip),%xmm0
subsd  -0xf0(%rbp),%xmm0
ff
mulsd  -0x40(%rbp),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x170(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
movsd  %xmm0,-0x180(%rbp)
ff
movsd  -0xf0(%rbp),%xmm0
ff
movsd  %xmm0,-0x1b0(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x1b0(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x188(%rbp)
ff
movsd  -0xd8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movsd  %xmm0,-0x190(%rbp)
ff
movsd  -0xf0(%rbp),%xmm0
ff
movsd  %xmm0,-0x1a8(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movaps %xmm0,%xmm1
movsd  -0x1a8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x198(%rbp)
ff
movsd  -0xd8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movaps %xmm0,%xmm1
movsd  0x5ec(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x1a0(%rbp)
ff
movsd  -0xd8(%rbp),%xmm0
ff
call   1070 <cos@plt>
movsd  -0x1a0(%rbp),%xmm6
ff
movsd  -0x198(%rbp),%xmm4
ff
movsd  -0x190(%rbp),%xmm5
ff
movsd  -0x188(%rbp),%xmm2
ff
movsd  -0x180(%rbp),%xmm3
ff
movsd  -0x178(%rbp),%xmm1
ff
movaps %xmm0,%xmm7
movsd  -0x170(%rbp),%xmm0
ff
mulsd  %xmm7,%xmm6
movsd  0x584(%rip),%xmm7
addsd  %xmm7,%xmm6
mulsd  %xmm6,%xmm4
addsd  %xmm5,%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
movsd  -0x10(%rbp),%xmm0
addsd  -0xf8(%rbp),%xmm0
ff
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x220,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<DirectVincenty>:
push   %r14
push   %rbx
sub    $0xb8,%rsp
mov    %rsi,%r14
mov    %rdi,%rbx
movsd  %xmm3,0x48(%rsp)
movsd  %xmm1,0x90(%rsp)
movsd  %xmm0,(%rsp)
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm4
setnp  %al
sete   %cl
and    %al,%cl
ucomisd %xmm0,%xmm5
setnp  %al
sete   %dl
and    %al,%dl
or     %cl,%dl
test   %dl,%dl
je     11ae <DirectVincenty+0x4e>
movsd  0xe52(%rip),%xmm5
je     11b8 <DirectVincenty+0x58>
movsd  0xe50(%rip),%xmm4
movsd  %xmm4,0x10(%rsp)
movsd  0xe4a(%rip),%xmm0
movsd  %xmm5,0x88(%rsp)
subsd  %xmm5,%xmm0
movsd  %xmm0,0x70(%rsp)
mulsd  %xmm0,%xmm4
movsd  %xmm4,0x8(%rsp)
movaps %xmm2,%xmm0
movsd  %xmm2,0x20(%rsp)
call   1040 <sin@plt>
movsd  %xmm0,0x78(%rsp)
movsd  0x20(%rsp),%xmm0
call   1070 <cos@plt>
movsd  %xmm0,0x30(%rsp)
movsd  (%rsp),%xmm0
call   1060 <tan@plt>
mulsd  0x70(%rsp),%xmm0
movsd  %xmm0,(%rsp)
call   1050 <atan@plt>
movsd  %xmm0,0x20(%rsp)
call   1040 <sin@plt>
movsd  %xmm0,0x40(%rsp)
movsd  0x20(%rsp),%xmm0
call   1070 <cos@plt>
movsd  %xmm0,0x28(%rsp)
movsd  (%rsp),%xmm0
movsd  0x30(%rsp),%xmm1
call   1030 <atan2@plt>
movsd  0x10(%rsp),%xmm2
movsd  %xmm0,0x18(%rsp)
movsd  0x78(%rsp),%xmm1
mulsd  0x28(%rsp),%xmm1
movsd  0xd9c(%rip),%xmm3
movapd %xmm3,%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm1,0x80(%rsp)
addsd  %xmm3,%xmm1
mulsd  %xmm0,%xmm1
movsd  0x8(%rsp),%xmm0
mulsd  %xmm0,%xmm0
mulsd  %xmm2,%xmm2
subsd  %xmm0,%xmm2
movsd  %xmm1,0x38(%rsp)
mulsd  %xmm1,%xmm2
divsd  %xmm0,%xmm2
addsd  %xmm3,%xmm2
ucomisd 0xda7(%rip),%xmm2
jb     12c4 <DirectVincenty+0x164>
xorps  %xmm0,%xmm0
sqrtsd %xmm2,%xmm0
jmp    12cd <DirectVincenty+0x16d>
movapd %xmm2,%xmm0
call   1080 <sqrt@plt>
movsd  0xd43(%rip),%xmm1
addsd  %xmm0,%xmm1
movsd  0xd2f(%rip),%xmm3
addsd  %xmm3,%xmm0
divsd  %xmm0,%xmm1
movapd %xmm1,%xmm4
mulsd  %xmm1,%xmm4
movsd  0xd27(%rip),%xmm5
movapd %xmm4,%xmm0
mulsd  %xmm5,%xmm0
addsd  %xmm3,%xmm0
movapd %xmm3,%xmm2
subsd  %xmm1,%xmm2
divsd  %xmm2,%xmm0
mulsd  0xd0f(%rip),%xmm4
addsd  %xmm3,%xmm4
mulsd  %xmm1,%xmm4
movsd  0x8(%rsp),%xmm2
mulsd  %xmm0,%xmm2
movsd  0x48(%rsp),%xmm1
divsd  %xmm2,%xmm1
mulsd  %xmm4,%xmm5
movsd  %xmm5,0xa0(%rsp)
movsd  0xce6(%rip),%xmm0
movsd  %xmm4,0xa8(%rsp)
mulsd  %xmm4,%xmm0
movsd  %xmm0,0x98(%rsp)
movsd  0x18(%rsp),%xmm0
addsd  %xmm0,%xmm0
movsd  %xmm0,0x18(%rsp)
movsd  %xmm1,0x48(%rsp)
cs nopw 0x0(%rax,%rax,1)
movapd %xmm1,%xmm0
movsd  %xmm1,0x10(%rsp)
movsd  0x18(%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
call   1040 <sin@plt>
mulsd  0xa8(%rsp),%xmm0
movsd  %xmm0,0x20(%rsp)
movsd  0x8(%rsp),%xmm0
call   1070 <cos@plt>
movsd  %xmm0,0x68(%rsp)
movsd  0x10(%rsp),%xmm0
call   1070 <cos@plt>
movsd  %xmm0,0x60(%rsp)
movsd  0x8(%rsp),%xmm0
call   1070 <cos@plt>
addsd  %xmm0,%xmm0
movsd  %xmm0,(%rsp)
movsd  0x8(%rsp),%xmm0
call   1070 <cos@plt>
mulsd  (%rsp),%xmm0
addsd  0xc1c(%rip),%xmm0
movsd  %xmm0,(%rsp)
movsd  0x8(%rsp),%xmm0
call   1070 <cos@plt>
movsd  %xmm0,0x58(%rsp)
movsd  0x10(%rsp),%xmm0
call   1040 <sin@plt>
movapd %xmm0,%xmm1
movsd  0xc0f(%rip),%xmm0
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x50(%rsp)
movsd  0x10(%rsp),%xmm0
call   1040 <sin@plt>
movapd %xmm0,%xmm1
mulsd  0x50(%rsp),%xmm1
movsd  0xbf0(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  %xmm1,0x50(%rsp)
movsd  0x8(%rsp),%xmm0
call   1070 <cos@plt>
mulsd  0xbcb(%rip),%xmm0
movsd  %xmm0,0xb0(%rsp)
movsd  0x8(%rsp),%xmm0
call   1070 <cos@plt>
movsd  0x10(%rsp),%xmm2
movsd  (%rsp),%xmm1
mulsd  0xb0(%rsp),%xmm0
addsd  0xba3(%rip),%xmm0
movsd  0x58(%rsp),%xmm3
mulsd  0x98(%rsp),%xmm3
mulsd  0x50(%rsp),%xmm3
mulsd  %xmm0,%xmm3
mulsd  0x60(%rsp),%xmm1
addsd  %xmm3,%xmm1
mulsd  0xa0(%rsp),%xmm1
addsd  0x68(%rsp),%xmm1
mulsd  0x20(%rsp),%xmm1
addsd  0x48(%rsp),%xmm1
subsd  %xmm1,%xmm2
ucomisd 0xb61(%rip),%xmm2
ja     1380 <DirectVincenty+0x220>
movsd  0xb5b(%rip),%xmm0
ucomisd %xmm2,%xmm0
ja     1380 <DirectVincenty+0x220>
movsd  0x18(%rsp),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x18(%rsp)
movapd %xmm1,%xmm0
movsd  %xmm1,(%rsp)
call   1070 <cos@plt>
movsd  %xmm0,0x8(%rsp)
movsd  (%rsp),%xmm0
call   1040 <sin@plt>
mulsd  0x28(%rsp),%xmm0
mulsd  0x30(%rsp),%xmm0
movsd  0x8(%rsp),%xmm1
mulsd  0x40(%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
movsd  (%rsp),%xmm0
call   1040 <sin@plt>
movsd  %xmm0,0x10(%rsp)
movsd  (%rsp),%xmm0
call   1070 <cos@plt>
mulsd  0x28(%rsp),%xmm0
mulsd  0x30(%rsp),%xmm0
movsd  0x10(%rsp),%xmm1
mulsd  0x40(%rsp),%xmm1
subsd  %xmm0,%xmm1
mulsd  %xmm1,%xmm1
movsd  0x80(%rsp),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     15aa <DirectVincenty+0x44a>
sqrtsd %xmm0,%xmm0
jmp    15af <DirectVincenty+0x44f>
call   1080 <sqrt@plt>
movsd  0x70(%rsp),%xmm1
mulsd  %xmm0,%xmm1
movsd  0x8(%rsp),%xmm0
call   1030 <atan2@plt>
movsd  %xmm0,(%rbx)
movsd  (%rsp),%xmm0
call   1040 <sin@plt>
mulsd  0x78(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  (%rsp),%xmm0
call   1070 <cos@plt>
movsd  %xmm0,0x10(%rsp)
movsd  (%rsp),%xmm0
call   1040 <sin@plt>
mulsd  0x40(%rsp),%xmm0
movsd  0x30(%rsp),%xmm2
mulsd  %xmm0,%xmm2
movsd  0x28(%rsp),%xmm1
mulsd  0x10(%rsp),%xmm1
subsd  %xmm2,%xmm1
movsd  0x8(%rsp),%xmm0
call   1030 <atan2@plt>
movsd  %xmm0,0x8(%rsp)
movsd  0xa27(%rip),%xmm0
movsd  0x88(%rsp),%xmm2
mulsd  %xmm2,%xmm0
movsd  0x38(%rsp),%xmm3
mulsd  %xmm3,%xmm0
mulsd  0x9f0(%rip),%xmm3
movsd  0x9e0(%rip),%xmm1
addsd  %xmm1,%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm1,%xmm3
mulsd  %xmm0,%xmm3
movsd  %xmm3,0x38(%rsp)
movsd  0x99a(%rip),%xmm0
subsd  %xmm3,%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  (%rsp),%xmm0
call   1040 <sin@plt>
mulsd  0x38(%rsp),%xmm0
movsd  %xmm0,0x20(%rsp)
movsd  0x18(%rsp),%xmm0
call   1070 <cos@plt>
movsd  %xmm0,0x68(%rsp)
movsd  (%rsp),%xmm0
call   1070 <cos@plt>
mulsd  0x38(%rsp),%xmm0
movsd  %xmm0,0x60(%rsp)
movsd  0x18(%rsp),%xmm0
call   1070 <cos@plt>
addsd  %xmm0,%xmm0
movsd  %xmm0,0x58(%rsp)
movsd  0x18(%rsp),%xmm0
call   1070 <cos@plt>
mulsd  0x58(%rsp),%xmm0
addsd  0x92d(%rip),%xmm0
mulsd  0x60(%rsp),%xmm0
addsd  0x68(%rsp),%xmm0
mulsd  0x20(%rsp),%xmm0
addsd  (%rsp),%xmm0
movsd  0x10(%rsp),%xmm2
mulsd  0x88(%rsp),%xmm2
movsd  0x80(%rsp),%xmm1
mulsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
movsd  0x8(%rsp),%xmm0
subsd  %xmm1,%xmm0
addsd  0x90(%rsp),%xmm0
movsd  %xmm0,(%r14)
add    $0xb8,%rsp
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @DirectVincenty(double %0, double %1, double %2, double %3, double* %4, double* %5, double %6, double %7) {
%9 = alloca double
%10 = alloca double
%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
store double %0, double* %9
store double %1, double* %10
store double %2, double* %11
store double %3, double* %12
store double* %4, double** %13
store double* %5, double** %14
store double %6, double* %15
store double %7, double* %16
%40 = load double, double* %15
%41 = fcmp oeq double %40, 0.000000e+00
br i1 %41, label %45, label %42
42: ; preds = %8
%43 = load double, double* %16
%44 = fcmp oeq double %43, 0.000000e+00
br i1 %44, label %45, label %46
45: ; preds = %42, %8
store double 0x415854A640000000, double* %15
store double 0x3F6B775A84F3E128, double* %16
br label %46
46: ; preds = %45, %42
%47 = load double, double* %16
%48 = fsub double 1.000000e+00, %47
%49 = load double, double* %15
%50 = fmul double %48, %49
store double %50, double* %17
%51 = load double, double* %11
%52 = call double @sin(double %51)
store double %52, double* %18
%53 = load double, double* %11
%54 = call double @cos(double %53)
store double %54, double* %19
%55 = load double, double* %16
%56 = fsub double 1.000000e+00, %55
%57 = load double, double* %9
%58 = call double @tan(double %57)
%59 = fmul double %56, %58
store double %59, double* %20
%60 = load double, double* %20
%61 = call double @atan(double %60)
store double %61, double* %21
%62 = load double, double* %21
%63 = call double @sin(double %62)
store double %63, double* %22
%64 = load double, double* %21
%65 = call double @cos(double %64)
store double %65, double* %23
%66 = load double, double* %20
%67 = load double, double* %19
%68 = call double @atan2(double %66, double %67)
store double %68, double* %24
%69 = load double, double* %23
%70 = load double, double* %18
%71 = fmul double %69, %70
store double %71, double* %25
%72 = load double, double* %25
%73 = fsub double 1.000000e+00, %72
%74 = load double, double* %25
%75 = fadd double 1.000000e+00, %74
%76 = fmul double %73, %75
store double %76, double* %26
%77 = load double, double* %26
%78 = load double, double* %15
%79 = load double, double* %15
%80 = load double, double* %17
%81 = load double, double* %17
%82 = fmul double %80, %81
%83 = fneg double %82
%84 = call double @llvm.fmuladd.f64(double %78, double %79, double %83)
%85 = fmul double %77, %84
%86 = load double, double* %17
%87 = load double, double* %17
%88 = fmul double %86, %87
%89 = fdiv double %85, %88
store double %89, double* %27
%90 = load double, double* %27
%91 = fadd double 1.000000e+00, %90
%92 = call double @sqrt(double %91)
store double %92, double* %28
%93 = load double, double* %28
%94 = fsub double %93, 1.000000e+00
%95 = load double, double* %28
%96 = fadd double %95, 1.000000e+00
%97 = fdiv double %94, %96
store double %97, double* %29
%98 = load double, double* %29
%99 = load double, double* %29
%100 = fmul double %98, %99
store double %100, double* %30
%101 = load double, double* %30
%102 = call double @llvm.fmuladd.f64(double 2.500000e-01, double %101, double 1.000000e+00)
%103 = load double, double* %29
%104 = fsub double 1.000000e+00, %103
%105 = fdiv double %102, %104
store double %105, double* %31
%106 = load double, double* %29
%107 = load double, double* %30
%108 = call double @llvm.fmuladd.f64(double -3.750000e-01, double %107, double 1.000000e+00)
%109 = fmul double %106, %108
store double %109, double* %32
%110 = load double, double* %12
%111 = load double, double* %17
%112 = load double, double* %31
%113 = fmul double %111, %112
%114 = fdiv double %110, %113
store double %114, double* %33
br label %115
115: ; preds = %176, %46
%116 = load double, double* %33
store double %116, double* %34
%117 = load double, double* %24
%118 = load double, double* %33
%119 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %117, double %118)
store double %119, double* %35
%120 = load double, double* %32
%121 = load double, double* %33
%122 = call double @sin(double %121)
%123 = fmul double %120, %122
%124 = load double, double* %35
%125 = call double @cos(double %124)
%126 = load double, double* %32
%127 = fmul double 2.500000e-01, %126
%128 = load double, double* %33
%129 = call double @cos(double %128)
%130 = load double, double* %35
%131 = call double @cos(double %130)
%132 = fmul double 2.000000e+00, %131
%133 = load double, double* %35
%134 = call double @cos(double %133)
%135 = call double @llvm.fmuladd.f64(double %132, double %134, double -1.000000e+00)
%136 = load double, double* %32
%137 = fmul double 0x3FC5555555555555, %136
%138 = load double, double* %35
%139 = call double @cos(double %138)
%140 = fmul double %137, %139
%141 = load double, double* %33
%142 = call double @sin(double %141)
%143 = fmul double 4.000000e+00, %142
%144 = load double, double* %33
%145 = call double @sin(double %144)
%146 = call double @llvm.fmuladd.f64(double %143, double %145, double -3.000000e+00)
%147 = fmul double %140, %146
%148 = load double, double* %35
%149 = call double @cos(double %148)
%150 = fmul double 4.000000e+00, %149
%151 = load double, double* %35
%152 = call double @cos(double %151)
%153 = call double @llvm.fmuladd.f64(double %150, double %152, double -3.000000e+00)
%154 = fmul double %147, %153
%155 = fneg double %154
%156 = call double @llvm.fmuladd.f64(double %129, double %135, double %155)
%157 = call double @llvm.fmuladd.f64(double %127, double %156, double %125)
%158 = fmul double %123, %157
store double %158, double* %36
%159 = load double, double* %12
%160 = load double, double* %17
%161 = load double, double* %31
%162 = fmul double %160, %161
%163 = fdiv double %159, %162
%164 = load double, double* %36
%165 = fadd double %163, %164
store double %165, double* %33
br label %166
166: ; preds = %115
%167 = load double, double* %34
%168 = load double, double* %33
%169 = fsub double %167, %168
%170 = fcmp ogt double %169, 0x3EB0C6F7A0B5ED8D
br i1 %170, label %176, label %171
171: ; preds = %166
%172 = load double, double* %34
%173 = load double, double* %33
%174 = fsub double %172, %173
%175 = fcmp olt double %174, 0xBEB0C6F7A0B5ED8D
br label %176
176: ; preds = %171, %166
%177 = phi i1 [ true, %166 ], [ %175, %171 ]
br i1 %177, label %115, label %178
178: ; preds = %176
%179 = load double, double* %24
%180 = load double, double* %33
%181 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %179, double %180)
store double %181, double* %35
%182 = load double, double* %22
%183 = load double, double* %33
%184 = call double @cos(double %183)
%185 = load double, double* %23
%186 = load double, double* %33
%187 = call double @sin(double %186)
%188 = fmul double %185, %187
%189 = load double, double* %19
%190 = fmul double %188, %189
%191 = call double @llvm.fmuladd.f64(double %182, double %184, double %190)
%192 = load double, double* %16
%193 = fsub double 1.000000e+00, %192
%194 = load double, double* %25
%195 = load double, double* %25
%196 = load double, double* %22
%197 = load double, double* %33
%198 = call double @sin(double %197)
%199 = load double, double* %23
%200 = load double, double* %33
%201 = call double @cos(double %200)
%202 = fmul double %199, %201
%203 = load double, double* %19
%204 = fmul double %202, %203
%205 = fneg double %204
%206 = call double @llvm.fmuladd.f64(double %196, double %198, double %205)
%207 = call double @pow(double %206, double 2.000000e+00)
%208 = call double @llvm.fmuladd.f64(double %194, double %195, double %207)
%209 = call double @sqrt(double %208)
%210 = fmul double %193, %209
%211 = call double @atan2(double %191, double %210)
%212 = load double*, double** %13
store double %211, double* %212
%213 = load double, double* %33
%214 = call double @sin(double %213)
%215 = load double, double* %18
%216 = fmul double %214, %215
%217 = load double, double* %23
%218 = load double, double* %33
%219 = call double @cos(double %218)
%220 = load double, double* %22
%221 = load double, double* %33
%222 = call double @sin(double %221)
%223 = fmul double %220, %222
%224 = load double, double* %19
%225 = fmul double %223, %224
%226 = fneg double %225
%227 = call double @llvm.fmuladd.f64(double %217, double %219, double %226)
%228 = call double @atan2(double %216, double %227)
store double %228, double* %37
%229 = load double, double* %16
%230 = fdiv double %229, 1.600000e+01
%231 = load double, double* %26
%232 = fmul double %230, %231
%233 = load double, double* %16
%234 = load double, double* %26
%235 = call double @llvm.fmuladd.f64(double -3.000000e+00, double %234, double 4.000000e+00)
%236 = call double @llvm.fmuladd.f64(double %233, double %235, double 4.000000e+00)
%237 = fmul double %232, %236
store double %237, double* %38
%238 = load double, double* %37
%239 = load double, double* %38
%240 = fsub double 1.000000e+00, %239
%241 = load double, double* %16
%242 = fmul double %240, %241
%243 = load double, double* %25
%244 = fmul double %242, %243
%245 = load double, double* %33
%246 = load double, double* %38
%247 = load double, double* %33
%248 = call double @sin(double %247)
%249 = fmul double %246, %248
%250 = load double, double* %35
%251 = call double @cos(double %250)
%252 = load double, double* %38
%253 = load double, double* %33
%254 = call double @cos(double %253)
%255 = fmul double %252, %254
%256 = load double, double* %35
%257 = call double @cos(double %256)
%258 = fmul double 2.000000e+00, %257
%259 = load double, double* %35
%260 = call double @cos(double %259)
%261 = call double @llvm.fmuladd.f64(double %258, double %260, double -1.000000e+00)
%262 = call double @llvm.fmuladd.f64(double %255, double %261, double %251)
%263 = call double @llvm.fmuladd.f64(double %249, double %262, double %245)
%264 = fneg double %244
%265 = call double @llvm.fmuladd.f64(double %264, double %263, double %238)
store double %265, double* %39
%266 = load double, double* %10
%267 = load double, double* %39
%268 = fadd double %266, %267
%269 = load double*, double** %14
store double %268, double* %269
ret void
}

LLVM IR (-O2):
define void @DirectVincenty(double %0, double %1, double %2, double %3, double* nocapture writeonly %4, double* nocapture writeonly %5, double %6, double %7) {
%9 = fcmp oeq double %6, 0.000000e+00
%10 = fcmp oeq double %7, 0.000000e+00
%11 = or i1 %9, %10
%12 = select i1 %11, double 0x3F6B775A84F3E128, double %7
%13 = select i1 %11, double 0x415854A640000000, double %6
%14 = fsub double 1.000000e+00, %12
%15 = fmul double %13, %14
%16 = tail call double @sin(double %2)
%17 = tail call double @cos(double %2)
%18 = tail call double @tan(double %0)
%19 = fmul double %14, %18
%20 = tail call double @atan(double %19)
%21 = tail call double @sin(double %20)
%22 = tail call double @cos(double %20)
%23 = tail call double @atan2(double %19, double %17)
%24 = fmul double %16, %22
%25 = fsub double 1.000000e+00, %24
%26 = fadd double %24, 1.000000e+00
%27 = fmul double %25, %26
%28 = fmul double %15, %15
%29 = fneg double %28
%30 = tail call double @llvm.fmuladd.f64(double %13, double %13, double %29)
%31 = fmul double %30, %27
%32 = fdiv double %31, %28
%33 = fadd double %32, 1.000000e+00
%34 = tail call double @sqrt(double %33)
%35 = fadd double %34, -1.000000e+00
%36 = fadd double %34, 1.000000e+00
%37 = fdiv double %35, %36
%38 = fmul double %37, %37
%39 = tail call double @llvm.fmuladd.f64(double %38, double 2.500000e-01, double 1.000000e+00)
%40 = fsub double 1.000000e+00, %37
%41 = fdiv double %39, %40
%42 = tail call double @llvm.fmuladd.f64(double %38, double -3.750000e-01, double 1.000000e+00)
%43 = fmul double %37, %42
%44 = fmul double %15, %41
%45 = fdiv double %3, %44
%46 = fmul double %43, 2.500000e-01
%47 = fmul double %43, 0xBFC5555555555555
br label %48
48: ; preds = %48, %8
%49 = phi double [ %45, %8 ], [ %74, %48 ]
%50 = tail call double @llvm.fmuladd.f64(double %23, double 2.000000e+00, double %49)
%51 = tail call double @sin(double %49)
%52 = fmul double %43, %51
%53 = tail call double @cos(double %50)
%54 = tail call double @cos(double %49)
%55 = tail call double @cos(double %50)
%56 = fmul double %55, 2.000000e+00
%57 = tail call double @cos(double %50)
%58 = tail call double @llvm.fmuladd.f64(double %56, double %57, double -1.000000e+00)
%59 = tail call double @cos(double %50)
%60 = tail call double @sin(double %49)
%61 = fmul double %60, 4.000000e+00
%62 = tail call double @sin(double %49)
%63 = tail call double @llvm.fmuladd.f64(double %61, double %62, double -3.000000e+00)
%64 = tail call double @cos(double %50)
%65 = fmul double %64, 4.000000e+00
%66 = tail call double @cos(double %50)
%67 = tail call double @llvm.fmuladd.f64(double %65, double %66, double -3.000000e+00)
%68 = fmul double %47, %59
%69 = fmul double %68, %63
%70 = fmul double %69, %67
%71 = tail call double @llvm.fmuladd.f64(double %54, double %58, double %70)
%72 = tail call double @llvm.fmuladd.f64(double %46, double %71, double %53)
%73 = fmul double %52, %72
%74 = fadd double %45, %73
%75 = fsub double %49, %74
%76 = fcmp ogt double %75, 0x3EB0C6F7A0B5ED8D
%77 = fcmp olt double %75, 0xBEB0C6F7A0B5ED8D
%78 = or i1 %76, %77
br i1 %78, label %48, label %79
79: ; preds = %48
%80 = tail call double @llvm.fmuladd.f64(double %23, double 2.000000e+00, double %74)
%81 = tail call double @cos(double %74)
%82 = tail call double @sin(double %74)
%83 = fmul double %22, %82
%84 = fmul double %17, %83
%85 = tail call double @llvm.fmuladd.f64(double %21, double %81, double %84)
%86 = tail call double @sin(double %74)
%87 = tail call double @cos(double %74)
%88 = fneg double %22
%89 = fmul double %87, %88
%90 = fmul double %17, %89
%91 = tail call double @llvm.fmuladd.f64(double %21, double %86, double %90)
%92 = fmul double %91, %91
%93 = tail call double @llvm.fmuladd.f64(double %24, double %24, double %92)
%94 = tail call double @sqrt(double %93)
%95 = fmul double %14, %94
%96 = tail call double @atan2(double %85, double %95)
store double %96, double* %4
%97 = tail call double @sin(double %74)
%98 = fmul double %16, %97
%99 = tail call double @cos(double %74)
%100 = tail call double @sin(double %74)
%101 = fneg double %21
%102 = fmul double %100, %101
%103 = fmul double %17, %102
%104 = tail call double @llvm.fmuladd.f64(double %22, double %99, double %103)
%105 = tail call double @atan2(double %98, double %104)
%106 = fmul double %12, 6.250000e-02
%107 = fmul double %106, %27
%108 = tail call double @llvm.fmuladd.f64(double %27, double -3.000000e+00, double 4.000000e+00)
%109 = tail call double @llvm.fmuladd.f64(double %12, double %108, double 4.000000e+00)
%110 = fmul double %107, %109
%111 = fsub double 1.000000e+00, %110
%112 = tail call double @sin(double %74)
%113 = fmul double %110, %112
%114 = tail call double @cos(double %80)
%115 = tail call double @cos(double %74)
%116 = fmul double %110, %115
%117 = tail call double @cos(double %80)
%118 = fmul double %117, 2.000000e+00
%119 = tail call double @cos(double %80)
%120 = tail call double @llvm.fmuladd.f64(double %118, double %119, double -1.000000e+00)
%121 = tail call double @llvm.fmuladd.f64(double %116, double %120, double %114)
%122 = tail call double @llvm.fmuladd.f64(double %113, double %121, double %74)
%123 = fneg double %111
%124 = fmul double %12, %123
%125 = fmul double %24, %124
%126 = tail call double @llvm.fmuladd.f64(double %125, double %122, double %105)
%127 = fadd double %126, %1
store double %127, double* %5
ret void
}