#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] + n * dt, 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),%xmm1
movsd  -0x10(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  0xe6b(%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)
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  0xd38(%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  0xcd1(%rip),%xmm0
ucomisd %xmm1,%xmm0
jb     1346 <EllipticToRectangular+0x1f6>
jmp    134b <EllipticToRectangular+0x1fb>
jmp    123c <EllipticToRectangular+0xec>
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  0xc32(%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  0xbeb(%rip),%xmm1
addsd  -0x70(%rbp),%xmm1
movsd  0xbde(%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  0xb11(%rip),%xmm0
subsd  -0x90(%rbp),%xmm0
ff
subsd  -0x98(%rbp),%xmm0
ff
call   1070 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xaf9(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  0xadd(%rip),%xmm0
subsd  -0x98(%rbp),%xmm0
ff
subsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  0xabd(%rip),%xmm0
subsd  -0x90(%rbp),%xmm0
ff
subsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x18(%rbp),%rax
movsd  0xaa1(%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  0x96d(%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  %xmm0,0x20(%rsp)
movsd  %xmm1,0x48(%rsp)
mulsd  %xmm1,%xmm2
addsd  0x8(%rdi),%xmm2
movsd  0xe86(%rip),%xmm1
movapd %xmm2,%xmm0
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)
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
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,0xe0f(%rip),%xmm2
movsd  0x10(%rbx),%xmm3
movsd  0x18(%rbx),%xmm4
movapd %xmm3,%xmm1
xorpd  0xde9(%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     11e0 <EllipticToRectangular+0x90>
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  0xd3c(%rip),%xmm11
movapd %xmm11,%xmm0
subsd  %xmm1,%xmm0
mulsd  %xmm2,%xmm2
subsd  %xmm2,%xmm0
xorpd  %xmm4,%xmm4
ucomisd %xmm4,%xmm0
movsd  %xmm9,(%rsp)
jb     12f3 <EllipticToRectangular+0x1a3>
sqrtsd %xmm0,%xmm0
jmp    1317 <EllipticToRectangular+0x1c7>
movsd  %xmm6,0x8(%rsp)
call   1070 <sqrt@plt>
xorpd  %xmm4,%xmm4
movsd  0x8(%rsp),%xmm6
movsd  0xcf7(%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     1399 <EllipticToRectangular+0x249>
sqrtsd %xmm0,%xmm0
jmp    13f7 <EllipticToRectangular+0x2a7>
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  0xc1d(%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 0xb96(%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 = load double, double* %7
%37 = load double, double* %9
%38 = call double @llvm.fmuladd.f64(double %36, double %37, double %35)
%39 = call double @fmod(double %38, double 0x401921FB54442D18)
store double %39, double* %11
%40 = load double, double* %11
%41 = load double*, double** %8
%42 = getelementptr inbounds double, double* %41, i64 2
%43 = load double, double* %42
%44 = load double, double* %11
%45 = call double @sin(double %44)
%46 = fneg double %43
%47 = call double @llvm.fmuladd.f64(double %46, double %45, double %40)
%48 = load double*, double** %8
%49 = getelementptr inbounds double, double* %48, i64 3
%50 = load double, double* %49
%51 = load double, double* %11
%52 = call double @cos(double %51)
%53 = call double @llvm.fmuladd.f64(double %50, double %52, double %47)
store double %53, double* %12
br label %54
54: ; preds = %94, %5
%55 = load double, double* %12
%56 = call double @cos(double %55)
store double %56, double* %13
%57 = load double, double* %12
%58 = call double @sin(double %57)
store double %58, double* %14
%59 = load double, double* %11
%60 = load double, double* %12
%61 = fsub double %59, %60
%62 = load double*, double** %8
%63 = getelementptr inbounds double, double* %62, i64 2
%64 = load double, double* %63
%65 = load double, double* %14
%66 = call double @llvm.fmuladd.f64(double %64, double %65, double %61)
%67 = load double*, double** %8
%68 = getelementptr inbounds double, double* %67, i64 3
%69 = load double, double* %68
%70 = load double, double* %13
%71 = fneg double %69
%72 = call double @llvm.fmuladd.f64(double %71, double %70, double %66)
%73 = load double*, double** %8
%74 = getelementptr inbounds double, double* %73, i64 2
%75 = load double, double* %74
%76 = load double, double* %13
%77 = fneg double %75
%78 = call double @llvm.fmuladd.f64(double %77, double %76, double 1.000000e+00)
%79 = load double*, double** %8
%80 = getelementptr inbounds double, double* %79, i64 3
%81 = load double, double* %80
%82 = load double, double* %14
%83 = fneg double %81
%84 = call double @llvm.fmuladd.f64(double %83, double %82, double %78)
%85 = fdiv double %72, %84
store double %85, double* %15
%86 = load double, double* %15
%87 = load double, double* %12
%88 = fadd double %87, %86
store double %88, double* %12
%89 = load double, double* %15
%90 = call i32 @fabs(double %89)
%91 = sitofp i32 %90 to double
%92 = fcmp ole double %91, 0x3D06849B86A12B9B
br i1 %92, label %93, label %94
93: ; preds = %54
br label %95
94: ; preds = %54
br label %54
95: ; preds = %93
%96 = load double, double* %12
%97 = call double @cos(double %96)
store double %97, double* %16
%98 = load double, double* %12
%99 = call double @sin(double %98)
store double %99, double* %17
%100 = load double*, double** %8
%101 = getelementptr inbounds double, double* %100, i64 2
%102 = load double, double* %101
%103 = fneg double %102
%104 = load double, double* %17
%105 = load double*, double** %8
%106 = getelementptr inbounds double, double* %105, i64 3
%107 = load double, double* %106
%108 = load double, double* %16
%109 = fmul double %107, %108
%110 = call double @llvm.fmuladd.f64(double %103, double %104, double %109)
store double %110, double* %18
%111 = load double*, double** %8
%112 = getelementptr inbounds double, double* %111, i64 2
%113 = load double, double* %112
%114 = load double*, double** %8
%115 = getelementptr inbounds double, double* %114, i64 2
%116 = load double, double* %115
%117 = fneg double %113
%118 = call double @llvm.fmuladd.f64(double %117, double %116, double 1.000000e+00)
%119 = load double*, double** %8
%120 = getelementptr inbounds double, double* %119, i64 3
%121 = load double, double* %120
%122 = load double*, double** %8
%123 = getelementptr inbounds double, double* %122, i64 3
%124 = load double, double* %123
%125 = fneg double %121
%126 = call double @llvm.fmuladd.f64(double %125, double %124, double %118)
%127 = call double @sqrt(double %126)
store double %127, double* %19
%128 = load double, double* %19
%129 = fadd double 1.000000e+00, %128
%130 = fdiv double 1.000000e+00, %129
store double %130, double* %20
%131 = load double, double* %6
%132 = load double, double* %16
%133 = load double*, double** %8
%134 = getelementptr inbounds double, double* %133, i64 2
%135 = load double, double* %134
%136 = fsub double %132, %135
%137 = load double, double* %20
%138 = load double, double* %18
%139 = fmul double %137, %138
%140 = load double*, double** %8
%141 = getelementptr inbounds double, double* %140, i64 3
%142 = load double, double* %141
%143 = fneg double %139
%144 = call double @llvm.fmuladd.f64(double %143, double %142, double %136)
%145 = fmul double %131, %144
store double %145, double* %21
%146 = load double, double* %6
%147 = load double, double* %17
%148 = load double*, double** %8
%149 = getelementptr inbounds double, double* %148, i64 3
%150 = load double, double* %149
%151 = fsub double %147, %150
%152 = load double, double* %20
%153 = load double, double* %18
%154 = fmul double %152, %153
%155 = load double*, double** %8
%156 = getelementptr inbounds double, double* %155, i64 2
%157 = load double, double* %156
%158 = call double @llvm.fmuladd.f64(double %154, double %157, double %151)
%159 = fmul double %146, %158
store double %159, double* %22
%160 = load double*, double** %8
%161 = getelementptr inbounds double, double* %160, i64 4
%162 = load double, double* %161
%163 = load double*, double** %8
%164 = getelementptr inbounds double, double* %163, i64 4
%165 = load double, double* %164
%166 = fmul double %162, %165
store double %166, double* %23
%167 = load double*, double** %8
%168 = getelementptr inbounds double, double* %167, i64 5
%169 = load double, double* %168
%170 = load double*, double** %8
%171 = getelementptr inbounds double, double* %170, i64 5
%172 = load double, double* %171
%173 = fmul double %169, %172
store double %173, double* %24
%174 = load double, double* %23
%175 = fsub double 1.000000e+00, %174
%176 = load double, double* %24
%177 = fsub double %175, %176
%178 = call double @sqrt(double %177)
%179 = fmul double 2.000000e+00, %178
store double %179, double* %25
%180 = load double, double* %24
%181 = fsub double 1.000000e+00, %180
%182 = load double, double* %24
%183 = fsub double %181, %182
store double %183, double* %26
%184 = load double, double* %23
%185 = fsub double 1.000000e+00, %184
%186 = load double, double* %23
%187 = fsub double %185, %186
store double %187, double* %27
%188 = load double*, double** %8
%189 = getelementptr inbounds double, double* %188, i64 5
%190 = load double, double* %189
%191 = fmul double 2.000000e+00, %190
%192 = load double*, double** %8
%193 = getelementptr inbounds double, double* %192, i64 4
%194 = load double, double* %193
%195 = fmul double %191, %194
store double %195, double* %28
%196 = load double, double* %21
%197 = load double, double* %26
%198 = load double, double* %22
%199 = load double, double* %28
%200 = fmul double %198, %199
%201 = call double @llvm.fmuladd.f64(double %196, double %197, double %200)
%202 = load double*, double** %10
%203 = getelementptr inbounds double, double* %202, i64 0
store double %201, double* %203
%204 = load double, double* %21
%205 = load double, double* %28
%206 = load double, double* %22
%207 = load double, double* %27
%208 = fmul double %206, %207
%209 = call double @llvm.fmuladd.f64(double %204, double %205, double %208)
%210 = load double*, double** %10
%211 = getelementptr inbounds double, double* %210, i64 1
store double %209, double* %211
%212 = load double, double* %21
%213 = fneg double %212
%214 = load double*, double** %8
%215 = getelementptr inbounds double, double* %214, i64 5
%216 = load double, double* %215
%217 = load double, double* %22
%218 = load double*, double** %8
%219 = getelementptr inbounds double, double* %218, i64 4
%220 = load double, double* %219
%221 = fmul double %217, %220
%222 = call double @llvm.fmuladd.f64(double %213, double %216, double %221)
%223 = load double, double* %25
%224 = fmul double %222, %223
%225 = load double*, double** %10
%226 = getelementptr inbounds double, double* %225, i64 2
store double %224, double* %226
%227 = load double*, double** %8
%228 = getelementptr inbounds double, double* %227, i64 2
%229 = load double, double* %228
%230 = fneg double %229
%231 = load double, double* %16
%232 = load double*, double** %8
%233 = getelementptr inbounds double, double* %232, i64 3
%234 = load double, double* %233
%235 = load double, double* %17
%236 = fmul double %234, %235
%237 = fneg double %236
%238 = call double @llvm.fmuladd.f64(double %230, double %231, double %237)
store double %238, double* %29
%239 = load double, double* %6
%240 = load double, double* %7
%241 = fmul double %239, %240
%242 = load double, double* %29
%243 = fadd double 1.000000e+00, %242
%244 = fdiv double %241, %243
store double %244, double* %30
%245 = load double, double* %30
%246 = load double, double* %17
%247 = fneg double %246
%248 = load double, double* %20
%249 = load double, double* %29
%250 = fmul double %248, %249
%251 = load double*, double** %8
%252 = getelementptr inbounds double, double* %251, i64 3
%253 = load double, double* %252
%254 = fneg double %250
%255 = call double @llvm.fmuladd.f64(double %254, double %253, double %247)
%256 = fmul double %245, %255
store double %256, double* %31
%257 = load double, double* %30
%258 = load double, double* %16
%259 = load double, double* %20
%260 = load double, double* %29
%261 = fmul double %259, %260
%262 = load double*, double** %8
%263 = getelementptr inbounds double, double* %262, i64 2
%264 = load double, double* %263
%265 = call double @llvm.fmuladd.f64(double %261, double %264, double %258)
%266 = fmul double %257, %265
store double %266, double* %32
%267 = load double, double* %31
%268 = load double, double* %26
%269 = load double, double* %32
%270 = load double, double* %28
%271 = fmul double %269, %270
%272 = call double @llvm.fmuladd.f64(double %267, double %268, double %271)
%273 = load double*, double** %10
%274 = getelementptr inbounds double, double* %273, i64 3
store double %272, double* %274
%275 = load double, double* %31
%276 = load double, double* %28
%277 = load double, double* %32
%278 = load double, double* %27
%279 = fmul double %277, %278
%280 = call double @llvm.fmuladd.f64(double %275, double %276, double %279)
%281 = load double*, double** %10
%282 = getelementptr inbounds double, double* %281, i64 4
store double %280, double* %282
%283 = load double, double* %31
%284 = fneg double %283
%285 = load double*, double** %8
%286 = getelementptr inbounds double, double* %285, i64 5
%287 = load double, double* %286
%288 = load double, double* %32
%289 = load double*, double** %8
%290 = getelementptr inbounds double, double* %289, i64 4
%291 = load double, double* %290
%292 = fmul double %288, %291
%293 = call double @llvm.fmuladd.f64(double %284, double %287, double %292)
%294 = load double, double* %25
%295 = fmul double %293, %294
%296 = load double*, double** %10
%297 = getelementptr inbounds double, double* %296, i64 5
store double %295, double* %297
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 @llvm.fmuladd.f64(double %1, double %3, double %7)
%9 = tail call double @fmod(double %8, double 0x401921FB54442D18)
%10 = getelementptr inbounds double, double* %2, i64 2
%11 = load double, double* %10
%12 = tail call double @sin(double %9)
%13 = fneg double %11
%14 = tail call double @llvm.fmuladd.f64(double %13, double %12, double %9)
%15 = getelementptr inbounds double, double* %2, i64 3
%16 = load double, double* %15
%17 = tail call double @cos(double %9)
%18 = tail call double @llvm.fmuladd.f64(double %16, double %17, double %14)
br label %19
19: ; preds = %19, %5
%20 = phi double [ %18, %5 ], [ %42, %19 ]
%21 = tail call double @cos(double %20)
%22 = tail call double @sin(double %20)
%23 = fsub double %9, %20
%24 = load double, double* %10
%25 = load double, double* %15
%26 = fneg double %25
%27 = fneg double %24
%28 = insertelement <2 x double> poison, double %24, i64 0
%29 = insertelement <2 x double> %28, double %27, i64 1
%30 = insertelement <2 x double> poison, double %22, i64 0
%31 = insertelement <2 x double> %30, double %21, i64 1
%32 = insertelement <2 x double> <double poison, double 1.000000e+00>, double %23, i64 0
%33 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %29, <2 x double> %31, <2 x double> %32)
%34 = insertelement <2 x double> poison, double %26, i64 0
%35 = shufflevector <2 x double> %34, <2 x double> poison, <2 x i32> zeroinitializer
%36 = insertelement <2 x double> poison, double %21, i64 0
%37 = insertelement <2 x double> %36, double %22, i64 1
%38 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %35, <2 x double> %37, <2 x double> %33)
%39 = extractelement <2 x double> %38, i64 0
%40 = extractelement <2 x double> %38, i64 1
%41 = fdiv double %39, %40
%42 = fadd double %20, %41
%43 = tail call i32 @fabs(double %41)
%44 = icmp slt i32 %43, 1
br i1 %44, label %45, label %19
45: ; preds = %19
%46 = tail call double @cos(double %42)
%47 = tail call double @sin(double %42)
%48 = load double, double* %10
%49 = fneg double %48
%50 = load double, double* %15
%51 = fmul double %46, %50
%52 = tail call double @llvm.fmuladd.f64(double %49, double %47, double %51)
%53 = tail call double @llvm.fmuladd.f64(double %49, double %48, double 1.000000e+00)
%54 = fneg double %50
%55 = tail call double @llvm.fmuladd.f64(double %54, double %50, double %53)
%56 = tail call double @sqrt(double %55)
%57 = fadd double %56, 1.000000e+00
%58 = fdiv double 1.000000e+00, %57
%59 = load double, double* %10
%60 = fsub double %46, %59
%61 = fmul double %52, %58
%62 = load double, double* %15
%63 = fneg double %61
%64 = tail call double @llvm.fmuladd.f64(double %63, double %62, double %60)
%65 = fmul double %64, %0
%66 = fsub double %47, %62
%67 = tail call double @llvm.fmuladd.f64(double %61, double %59, double %66)
%68 = fmul double %67, %0
%69 = getelementptr inbounds double, double* %2, i64 4
%70 = load double, double* %69
%71 = fmul double %70, %70
%72 = getelementptr inbounds double, double* %2, i64 5
%73 = load double, double* %72
%74 = fmul double %73, %73
%75 = fsub double 1.000000e+00, %71
%76 = fsub double %75, %74
%77 = tail call double @sqrt(double %76)
%78 = fmul double %77, 2.000000e+00
%79 = fsub double 1.000000e+00, %74
%80 = fsub double %79, %74
%81 = fsub double %75, %71
%82 = load double, double* %72
%83 = fmul double %82, 2.000000e+00
%84 = load double, double* %69
%85 = fmul double %83, %84
%86 = fmul double %68, %85
%87 = tail call double @llvm.fmuladd.f64(double %65, double %80, double %86)
store double %87, double* %4
%88 = fmul double %68, %81
%89 = tail call double @llvm.fmuladd.f64(double %65, double %85, double %88)
%90 = getelementptr inbounds double, double* %4, i64 1
store double %89, double* %90
%91 = fneg double %65
%92 = load double, double* %72
%93 = load double, double* %69
%94 = fmul double %68, %93
%95 = tail call double @llvm.fmuladd.f64(double %91, double %92, double %94)
%96 = fmul double %78, %95
%97 = getelementptr inbounds double, double* %4, i64 2
store double %96, double* %97
%98 = load double, double* %10
%99 = fneg double %98
%100 = load double, double* %15
%101 = fneg double %100
%102 = fmul double %47, %101
%103 = tail call double @llvm.fmuladd.f64(double %99, double %46, double %102)
%104 = fmul double %0, %1
%105 = fadd double %103, 1.000000e+00
%106 = fdiv double %104, %105
%107 = fneg double %47
%108 = fmul double %58, %103
%109 = fneg double %108
%110 = tail call double @llvm.fmuladd.f64(double %109, double %100, double %107)
%111 = fmul double %106, %110
%112 = tail call double @llvm.fmuladd.f64(double %108, double %98, double %46)
%113 = fmul double %106, %112
%114 = fmul double %85, %113
%115 = tail call double @llvm.fmuladd.f64(double %111, double %80, double %114)
%116 = getelementptr inbounds double, double* %4, i64 3
store double %115, double* %116
%117 = fmul double %81, %113
%118 = tail call double @llvm.fmuladd.f64(double %111, double %85, double %117)
%119 = getelementptr inbounds double, double* %4, i64 4
store double %118, double* %119
%120 = fneg double %111
%121 = load double, double* %72
%122 = load double, double* %69
%123 = fmul double %122, %113
%124 = tail call double @llvm.fmuladd.f64(double %120, double %121, double %123)
%125 = fmul double %78, %124
%126 = getelementptr inbounds double, double* %4, i64 5
store double %125, double* %126
ret void
}