#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 const cos(double const);
int fabs(double const);
double fmod(double const, double);
double const sin(double const);
double sqrt(double);
void EllipticToRectangular(const double a, const double n, const double elem[6],
                           const double dt, double xyz[]) {
  const double L = fmod(elem[1], 2.0 * 3.14159265358979323846);
  double Le = L - elem[2] * sin(L) + elem[3] * cos(L);
  for (;;) {
    const double cLe = cos(Le);
    const double sLe = sin(Le);
    const double dLe = (L - Le + elem[2] * sLe - elem[3] * cLe) /
                       (1.0 - elem[2] * cLe - elem[3] * sLe);
    Le += dLe;
    if (fabs(dLe) <= 1e-14)
      break;
  }
  {
    const double cLe = cos(Le);
    const double sLe = sin(Le);
    const double dlf = -elem[2] * sLe + elem[3] * cLe;
    const double phi = sqrt(1.0 - elem[2] * elem[2] - elem[3] * elem[3]);
    const double psi = 1.0 / (1.0 + phi);
    const double x1 = a * (cLe - elem[2] - psi * dlf * elem[3]);
    const double y1 = a * (sLe - elem[3] + psi * dlf * elem[2]);
    const double elem_4q = elem[4] * elem[4];
    const double elem_5q = elem[5] * elem[5];
    const double dwho = 2.0 * sqrt(1.0 - elem_4q - elem_5q);
    const double rtp = 1.0 - elem_5q - elem_5q;
    const double rtq = 1.0 - elem_4q - elem_4q;
    const double rdg = 2.0 * elem[5] * elem[4];
    xyz[0] = x1 * rtp + y1 * rdg;
    xyz[1] = x1 * rdg + y1 * rtq;
    xyz[2] = (-x1 * elem[5] + y1 * elem[4]) * dwho;
    const double rsam1 = -elem[2] * cLe - elem[3] * sLe;
    const double h = a * n / (1.0 + rsam1);
    const double vx1 = h * (-sLe - psi * rsam1 * elem[3]);
    const double vy1 = h * (cLe + psi * rsam1 * elem[2]);
    xyz[3] = vx1 * rtp + vy1 * rdg;
    xyz[4] = vx1 * rdg + vy1 * rtq;
    xyz[5] = (-vx1 * elem[5] + vy1 * elem[4]) * dwho;
  }
}


x86 assembly (clang, -O0):
<EllipticToRectangular>:
push   %rbp
mov    %rsp,%rbp
sub    $0x100,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
mov    %rdi,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
mov    %rsi,-0x28(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  0xe7d(%rip),%xmm1
call   1030 <fmod@plt>
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
call   1050 <sin@plt>
movsd  -0xf8(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xf0(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
call   1060 <cos@plt>
movsd  -0xe8(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xe0(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1060 <cos@plt>
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1050 <sin@plt>
movsd  %xmm0,-0x48(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm1
movsd  -0x48(%rbp),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm0
movsd  -0x40(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm2
movsd  -0x40(%rbp),%xmm1
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm1,%xmm2
movsd  0xd4a(%rip),%xmm1
addsd  %xmm1,%xmm2
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm1
movsd  -0x48(%rbp),%xmm3
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x50(%rbp),%xmm0
call   1040 <fabs@plt>
cvtsi2sd %eax,%xmm1
movsd  0xce3(%rip),%xmm0
ucomisd %xmm1,%xmm0
jb     1334 <EllipticToRectangular+0x1e4>
jmp    1339 <EllipticToRectangular+0x1e9>
jmp    122a <EllipticToRectangular+0xda>
movsd  -0x38(%rbp),%xmm0
call   1060 <cos@plt>
movsd  %xmm0,-0x58(%rbp)
movsd  -0x38(%rbp),%xmm0
call   1050 <sin@plt>
movsd  %xmm0,-0x60(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x60(%rbp),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm1
mulsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm1
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm0,%xmm1
movsd  0xc44(%rip),%xmm0
addsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1070 <sqrt@plt>
movsd  %xmm0,-0x70(%rbp)
movsd  0xbfd(%rip),%xmm1
addsd  -0x70(%rbp),%xmm1
movsd  0xbf0(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
mov    -0x18(%rbp),%rax
subsd  0x10(%rax),%xmm2
movsd  -0x78(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm3
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm2
mov    -0x18(%rbp),%rax
subsd  0x18(%rax),%xmm2
movsd  -0x78(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x18(%rbp),%rax
movsd  0x20(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulsd  0x20(%rax),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x18(%rbp),%rax
movsd  0x28(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulsd  0x28(%rax),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0xb23(%rip),%xmm0
subsd  -0x90(%rbp),%xmm0
ff
subsd  -0x98(%rbp),%xmm0
ff
call   1070 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xb0b(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  0xaef(%rip),%xmm0
subsd  -0x98(%rbp),%xmm0
ff
subsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  0xacf(%rip),%xmm0
subsd  -0x90(%rbp),%xmm0
ff
subsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x18(%rbp),%rax
movsd  0xab3(%rip),%xmm0
mulsd  0x28(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulsd  0x20(%rax),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  -0xa8(%rbp),%xmm2
ff
movsd  -0x88(%rbp),%xmm1
ff
mulsd  -0xb8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x80(%rbp),%xmm0
movsd  -0xb8(%rbp),%xmm2
ff
movsd  -0x88(%rbp),%xmm1
ff
mulsd  -0xb0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0x80(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x18(%rbp),%rax
movsd  0x28(%rax),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mov    -0x18(%rbp),%rax
mulsd  0x20(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
mov    -0x28(%rbp),%rax
movsd  %xmm0,0x10(%rax)
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x58(%rbp),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm1
mulsd  -0x60(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x8(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  0x97f(%rip),%xmm1
addsd  -0xc0(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
movsd  -0x60(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
movsd  -0x78(%rbp),%xmm1
mulsd  -0xc0(%rbp),%xmm1
ff
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm3
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  -0xc8(%rbp),%xmm0
ff
movsd  -0x58(%rbp),%xmm2
movsd  -0x78(%rbp),%xmm1
mulsd  -0xc0(%rbp),%xmm1
ff
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movsd  -0xd0(%rbp),%xmm0
ff
movsd  -0xa8(%rbp),%xmm2
ff
movsd  -0xd8(%rbp),%xmm1
ff
mulsd  -0xb8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,0x18(%rax)
movsd  -0xd0(%rbp),%xmm0
ff
movsd  -0xb8(%rbp),%xmm2
ff
movsd  -0xd8(%rbp),%xmm1
ff
mulsd  -0xb0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,0x20(%rax)
movsd  -0xd0(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x18(%rbp),%rax
movsd  0x28(%rax),%xmm2
movsd  -0xd8(%rbp),%xmm1
ff
mov    -0x18(%rbp),%rax
mulsd  0x20(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
mov    -0x28(%rbp),%rax
movsd  %xmm0,0x28(%rax)
add    $0x100,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<EllipticToRectangular>:
push   %r14
push   %rbx
sub    $0x68,%rsp
mov    %rsi,%r14
mov    %rdi,%rbx
movsd  %xmm1,0x48(%rsp)
movsd  %xmm0,0x20(%rsp)
movsd  0x8(%rdi),%xmm0
movsd  0xe8a(%rip),%xmm1
call   1030 <fmod@plt>
movsd  %xmm0,0x8(%rsp)
movsd  0x10(%rbx),%xmm1
movsd  %xmm1,(%rsp)
call   1050 <sin@plt>
mulsd  (%rsp),%xmm0
movsd  0x8(%rsp),%xmm1
movapd %xmm1,%xmm2
subsd  %xmm0,%xmm2
movsd  %xmm2,(%rsp)
movsd  0x18(%rbx),%xmm0
movsd  %xmm0,0x10(%rsp)
movapd %xmm1,%xmm0
call   1060 <cos@plt>
mulsd  0x10(%rsp),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
nopl   0x0(%rax)
movsd  (%rsp),%xmm0
call   1060 <cos@plt>
movaps %xmm0,0x10(%rsp)
movsd  (%rsp),%xmm0
call   1050 <sin@plt>
movsd  0x8(%rsp),%xmm2
movsd  (%rsp),%xmm5
subsd  %xmm5,%xmm2
shufpd $0x2,0xe1f(%rip),%xmm2
movsd  0x10(%rbx),%xmm3
movsd  0x18(%rbx),%xmm4
movapd %xmm3,%xmm1
xorpd  0xdf9(%rip),%xmm1
unpcklpd %xmm1,%xmm3
movapd %xmm0,%xmm1
movapd 0x10(%rsp),%xmm6
unpcklpd %xmm6,%xmm1
mulpd  %xmm3,%xmm1
addpd  %xmm2,%xmm1
unpcklpd %xmm4,%xmm4
movapd %xmm6,%xmm2
unpcklpd %xmm0,%xmm2
mulpd  %xmm4,%xmm2
subpd  %xmm2,%xmm1
movapd %xmm1,%xmm0
unpckhpd %xmm1,%xmm0
divsd  %xmm0,%xmm1
movapd %xmm5,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rsp)
movapd %xmm1,%xmm0
call   1040 <fabs@plt>
test   %eax,%eax
jg     11d0 <EllipticToRectangular+0x80>
movsd  (%rsp),%xmm0
call   1060 <cos@plt>
movsd  %xmm0,0x10(%rsp)
movsd  (%rsp),%xmm0
call   1050 <sin@plt>
movapd %xmm0,%xmm9
movsd  0x10(%rbx),%xmm1
movsd  0x18(%rbx),%xmm2
movsd  0x10(%rsp),%xmm6
mulsd  %xmm2,%xmm6
movapd %xmm1,%xmm0
mulsd  %xmm9,%xmm0
subsd  %xmm0,%xmm6
mulsd  %xmm1,%xmm1
movsd  0xd4c(%rip),%xmm11
movapd %xmm11,%xmm0
subsd  %xmm1,%xmm0
mulsd  %xmm2,%xmm2
subsd  %xmm2,%xmm0
xorpd  %xmm4,%xmm4
ucomisd %xmm4,%xmm0
movsd  %xmm9,(%rsp)
jb     12e3 <EllipticToRectangular+0x193>
sqrtsd %xmm0,%xmm0
jmp    1307 <EllipticToRectangular+0x1b7>
movsd  %xmm6,0x8(%rsp)
call   1070 <sqrt@plt>
xorpd  %xmm4,%xmm4
movsd  0x8(%rsp),%xmm6
movsd  0xd07(%rip),%xmm11
movsd  (%rsp),%xmm9
addsd  %xmm11,%xmm0
movapd %xmm11,%xmm5
divsd  %xmm0,%xmm5
movsd  0x10(%rbx),%xmm0
movsd  0x18(%rbx),%xmm1
movsd  0x10(%rsp),%xmm7
subsd  %xmm0,%xmm7
mulsd  %xmm5,%xmm6
movapd %xmm6,%xmm2
mulsd  %xmm1,%xmm2
subsd  %xmm2,%xmm7
movsd  0x20(%rsp),%xmm3
mulsd  %xmm3,%xmm7
movapd %xmm9,%xmm2
subsd  %xmm1,%xmm2
mulsd  %xmm0,%xmm6
addsd  %xmm2,%xmm6
mulsd  %xmm3,%xmm6
movsd  0x20(%rbx),%xmm2
mulsd  %xmm2,%xmm2
movsd  0x28(%rbx),%xmm3
mulsd  %xmm3,%xmm3
movapd %xmm11,%xmm12
subsd  %xmm2,%xmm12
movapd %xmm12,%xmm0
subsd  %xmm3,%xmm0
ucomisd %xmm4,%xmm0
jb     1389 <EllipticToRectangular+0x239>
sqrtsd %xmm0,%xmm0
jmp    13e7 <EllipticToRectangular+0x297>
movapd %xmm5,0x50(%rsp)
movsd  %xmm6,0x8(%rsp)
movsd  %xmm12,0x40(%rsp)
movsd  %xmm7,0x38(%rsp)
movsd  %xmm2,0x30(%rsp)
movsd  %xmm3,0x28(%rsp)
call   1070 <sqrt@plt>
movsd  0x28(%rsp),%xmm3
movsd  0x30(%rsp),%xmm2
movsd  0x38(%rsp),%xmm7
movsd  0x40(%rsp),%xmm12
movsd  0x8(%rsp),%xmm6
movsd  0xc2d(%rip),%xmm11
movapd 0x50(%rsp),%xmm5
movsd  (%rsp),%xmm9
movapd %xmm11,%xmm1
subsd  %xmm3,%xmm1
subsd  %xmm3,%xmm1
subsd  %xmm2,%xmm12
movsd  0x28(%rbx),%xmm2
addsd  %xmm2,%xmm2
mulsd  0x20(%rbx),%xmm2
movapd %xmm6,%xmm3
mulsd  %xmm2,%xmm3
movapd %xmm7,%xmm4
mulsd  %xmm1,%xmm4
addsd  %xmm3,%xmm4
movsd  %xmm4,(%r14)
movapd %xmm6,%xmm3
mulsd  %xmm12,%xmm3
movapd %xmm7,%xmm4
mulsd  %xmm2,%xmm4
addsd  %xmm3,%xmm4
movsd  %xmm4,0x8(%r14)
mulsd  0x20(%rbx),%xmm6
addsd  %xmm0,%xmm0
mulsd  0x28(%rbx),%xmm7
subsd  %xmm7,%xmm6
mulsd  %xmm0,%xmm6
movsd  %xmm6,0x10(%r14)
movsd  0x10(%rbx),%xmm8
movsd  0x18(%rbx),%xmm4
movapd 0xba6(%rip),%xmm3
movapd %xmm8,%xmm6
xorpd  %xmm3,%xmm6
movapd %xmm9,%xmm7
mulsd  %xmm4,%xmm7
movsd  0x10(%rsp),%xmm10
mulsd  %xmm10,%xmm6
subsd  %xmm7,%xmm6
movsd  0x20(%rsp),%xmm7
mulsd  0x48(%rsp),%xmm7
addsd  %xmm6,%xmm11
divsd  %xmm11,%xmm7
mulsd  %xmm6,%xmm5
xorpd  %xmm5,%xmm3
mulsd  %xmm4,%xmm3
subsd  %xmm9,%xmm3
mulsd  %xmm7,%xmm3
mulsd  %xmm8,%xmm5
addsd  %xmm10,%xmm5
mulsd  %xmm7,%xmm5
movapd %xmm3,%xmm4
mulsd  %xmm2,%xmm4
mulsd  %xmm5,%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,0x18(%r14)
mulsd  %xmm5,%xmm12
addsd  %xmm4,%xmm12
movsd  %xmm12,0x20(%r14)
mulsd  0x20(%rbx),%xmm5
mulsd  0x28(%rbx),%xmm3
subsd  %xmm3,%xmm5
mulsd  %xmm0,%xmm5
movsd  %xmm5,0x28(%r14)
add    $0x68,%rsp
pop    %rbx
pop    %r14
ret

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

LLVM IR (-O2):
define void @EllipticToRectangular(double %0, double %1, double* nocapture %2, double %3, double* nocapture writeonly %4) {
%6 = getelementptr inbounds double, double* %2, i64 1
%7 = load double, double* %6
%8 = tail call double @fmod(double %7, double 0x401921FB54442D18)
%9 = getelementptr inbounds double, double* %2, i64 2
%10 = load double, double* %9
%11 = tail call double @sin(double %8)
%12 = fneg double %10
%13 = tail call double @llvm.fmuladd.f64(double %12, double %11, double %8)
%14 = getelementptr inbounds double, double* %2, i64 3
%15 = load double, double* %14
%16 = tail call double @cos(double %8)
%17 = tail call double @llvm.fmuladd.f64(double %15, double %16, double %13)
br label %18
18: ; preds = %18, %5
%19 = phi double [ %17, %5 ], [ %41, %18 ]
%20 = tail call double @cos(double %19)
%21 = tail call double @sin(double %19)
%22 = fsub double %8, %19
%23 = load double, double* %9
%24 = load double, double* %14
%25 = fneg double %24
%26 = fneg double %23
%27 = insertelement <2 x double> poison, double %23, i64 0
%28 = insertelement <2 x double> %27, double %26, i64 1
%29 = insertelement <2 x double> poison, double %21, i64 0
%30 = insertelement <2 x double> %29, double %20, i64 1
%31 = insertelement <2 x double> <double poison, double 1.000000e+00>, double %22, i64 0
%32 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %28, <2 x double> %30, <2 x double> %31)
%33 = insertelement <2 x double> poison, double %25, i64 0
%34 = shufflevector <2 x double> %33, <2 x double> poison, <2 x i32> zeroinitializer
%35 = insertelement <2 x double> poison, double %20, i64 0
%36 = insertelement <2 x double> %35, double %21, i64 1
%37 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %34, <2 x double> %36, <2 x double> %32)
%38 = extractelement <2 x double> %37, i64 0
%39 = extractelement <2 x double> %37, i64 1
%40 = fdiv double %38, %39
%41 = fadd double %19, %40
%42 = tail call i32 @fabs(double %40)
%43 = icmp slt i32 %42, 1
br i1 %43, label %44, label %18
44: ; preds = %18
%45 = tail call double @cos(double %41)
%46 = tail call double @sin(double %41)
%47 = load double, double* %9
%48 = fneg double %47
%49 = load double, double* %14
%50 = fmul double %45, %49
%51 = tail call double @llvm.fmuladd.f64(double %48, double %46, double %50)
%52 = tail call double @llvm.fmuladd.f64(double %48, double %47, double 1.000000e+00)
%53 = fneg double %49
%54 = tail call double @llvm.fmuladd.f64(double %53, double %49, double %52)
%55 = tail call double @sqrt(double %54)
%56 = fadd double %55, 1.000000e+00
%57 = fdiv double 1.000000e+00, %56
%58 = load double, double* %9
%59 = fsub double %45, %58
%60 = fmul double %51, %57
%61 = load double, double* %14
%62 = fneg double %60
%63 = tail call double @llvm.fmuladd.f64(double %62, double %61, double %59)
%64 = fmul double %63, %0
%65 = fsub double %46, %61
%66 = tail call double @llvm.fmuladd.f64(double %60, double %58, double %65)
%67 = fmul double %66, %0
%68 = getelementptr inbounds double, double* %2, i64 4
%69 = load double, double* %68
%70 = fmul double %69, %69
%71 = getelementptr inbounds double, double* %2, i64 5
%72 = load double, double* %71
%73 = fmul double %72, %72
%74 = fsub double 1.000000e+00, %70
%75 = fsub double %74, %73
%76 = tail call double @sqrt(double %75)
%77 = fmul double %76, 2.000000e+00
%78 = fsub double 1.000000e+00, %73
%79 = fsub double %78, %73
%80 = fsub double %74, %70
%81 = load double, double* %71
%82 = fmul double %81, 2.000000e+00
%83 = load double, double* %68
%84 = fmul double %82, %83
%85 = fmul double %67, %84
%86 = tail call double @llvm.fmuladd.f64(double %64, double %79, double %85)
store double %86, double* %4
%87 = fmul double %67, %80
%88 = tail call double @llvm.fmuladd.f64(double %64, double %84, double %87)
%89 = getelementptr inbounds double, double* %4, i64 1
store double %88, double* %89
%90 = fneg double %64
%91 = load double, double* %71
%92 = load double, double* %68
%93 = fmul double %67, %92
%94 = tail call double @llvm.fmuladd.f64(double %90, double %91, double %93)
%95 = fmul double %77, %94
%96 = getelementptr inbounds double, double* %4, i64 2
store double %95, double* %96
%97 = load double, double* %9
%98 = fneg double %97
%99 = load double, double* %14
%100 = fneg double %99
%101 = fmul double %46, %100
%102 = tail call double @llvm.fmuladd.f64(double %98, double %45, double %101)
%103 = fmul double %0, %1
%104 = fadd double %102, 1.000000e+00
%105 = fdiv double %103, %104
%106 = fneg double %46
%107 = fmul double %57, %102
%108 = fneg double %107
%109 = tail call double @llvm.fmuladd.f64(double %108, double %99, double %106)
%110 = fmul double %105, %109
%111 = tail call double @llvm.fmuladd.f64(double %107, double %97, double %45)
%112 = fmul double %105, %111
%113 = fmul double %84, %112
%114 = tail call double @llvm.fmuladd.f64(double %110, double %79, double %113)
%115 = getelementptr inbounds double, double* %4, i64 3
store double %114, double* %115
%116 = fmul double %80, %112
%117 = tail call double @llvm.fmuladd.f64(double %110, double %84, double %116)
%118 = getelementptr inbounds double, double* %4, i64 4
store double %117, double* %118
%119 = fneg double %110
%120 = load double, double* %71
%121 = load double, double* %68
%122 = fmul double %121, %112
%123 = tail call double @llvm.fmuladd.f64(double %119, double %120, double %122)
%124 = fmul double %77, %123
%125 = getelementptr inbounds double, double* %4, i64 5
store double %124, double* %125
ret void
}