void rstage(int n, int n2, int n4, double *x1, double *x2, double *x3,
            double *x4) {
  int n8, is, id, i1, i2, j, jn, i;
  double t1, t2, e, t3, t4, t5;
  double ss1, sd1, ss3, sd3;
  double cc1, cd1, cc3, cd3;
  n8 = n2 / 8;
  is = 0;
  id = n2 * 2;
  do {
    for (i1 = is + 1; i1 <= n; i1 += id) {
      t1 = x4[i1] + x3[i1];
      x4[i1] = x4[i1] - x3[i1];
      x3[i1] = x1[i1] - t1;
      x1[i1] = x1[i1] + t1;
    }
    is = 2 * id - n2;
    id = 4 * id;
  } while (is < n);
  if (n4 - 1 <= 0)
    return;
  is = 0;
  id = n2 * 2;
  do {
    for (i2 = is + 1 + n8; i2 <= n; i2 += id) {
      t1 = (x3[i2] + x4[i2]) * .7071067811865475;
      t2 = (x3[i2] - x4[i2]) * .7071067811865475;
      x4[i2] = x2[i2] - t1;
      x3[i2] = -x2[i2] - t1;
      x2[i2] = x1[i2] - t2;
      x1[i2] = x1[i2] + t2;
    }
    is = 2 * id - n2;
    id = 4 * id;
  } while (is < n);
  if (n8 - 1 <= 0)
    return;
  e = 2 * 3.14159265358979323 / n2;
  ss1 = sin(e);
  sd1 = ss1;
  sd3 = 3 * sd1 - 4 * (sd1 * sd1 * sd1);
  ss3 = sd3;
  cc1 = cos(e);
  cd1 = cc1;
  cd3 = 4 * (cd1 * cd1 * cd1) - 3. * cd1;
  cc3 = cd3;
  for (j = 2; j <= n8; j++) {
    is = 0;
    id = 2 * n2;
    jn = n4 - 2 * j + 2;
    do {
      for (i1 = is + j; i1 <= n; i1 += id) {
        i2 = i1 + jn;
        t1 = x3[i1] * cc1 + x3[i2] * ss1;
        t2 = x3[i2] * cc1 - x3[i1] * ss1;
        t3 = x4[i1] * cc3 + x4[i2] * ss3;
        t4 = x4[i2] * cc3 - x4[i1] * ss3;
        t5 = t1 + t3;
        t3 = t1 - t3;
        t1 = t2 + t4;
        t4 = t2 - t4;
        x3[i1] = t1 - x2[i2];
        x4[i2] = t1 + x2[i2];
        x3[i2] = -x2[i1] - t3;
        x4[i1] = x2[i1] - t3;
        x2[i2] = x1[i1] - t5;
        x1[i1] = x1[i1] + t5;
        x2[i1] = x1[i2] + t4;
        x1[i2] = x1[i2] - t4;
      }
      is = 2 * id - n2;
      id = 4 * id;
    } while (is < n);
    t1 = cc1 * cd1 - ss1 * sd1;
    ss1 = cc1 * sd1 + ss1 * cd1;
    cc1 = t1;
    t3 = cc3 * cd3 - ss3 * sd3;
    ss3 = cc3 * sd3 + ss3 * cd3;
    cc3 = t3;
  }
  return;
}


x86 assembly (clang, -O0):
<rstage>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    -0x8(%rbp),%eax
mov    $0x8,%ecx
cltd
idiv   %ecx
mov    %eax,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x8(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1209 <rstage+0xe9>
mov    0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x50(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x50(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    116b <rstage+0x4b>
mov    -0x34(%rbp),%eax
shl    $0x1,%eax
sub    -0x8(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x34(%rbp),%eax
shl    $0x2,%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jl     1162 <rstage+0x42>
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
cmp    $0x0,%eax
jg     123e <rstage+0x11e>
jmp    1894 <rstage+0x774>
movl   $0x0,-0x30(%rbp)
mov    -0x8(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     135d <rstage+0x23d>
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  0xd90(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  0xd65(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x50(%rbp),%xmm0
mov    0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
subsd  -0x50(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x58(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x58(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
jmp    125a <rstage+0x13a>
mov    -0x34(%rbp),%eax
shl    $0x1,%eax
sub    -0x8(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x34(%rbp),%eax
shl    $0x2,%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jl     124e <rstage+0x12e>
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
cmp    $0x0,%eax
jg     1392 <rstage+0x272>
jmp    1894 <rstage+0x774>
cvtsi2sdl -0x8(%rbp),%xmm1
movsd  0xc71(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0x88(%rbp),%xmm2
ff
mulsd  -0x88(%rbp),%xmm2
ff
mulsd  -0x88(%rbp),%xmm2
ff
movsd  0xc14(%rip),%xmm1
mulsd  %xmm2,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xbf9(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
call   1040 <cos@plt>
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
mulsd  -0xa8(%rbp),%xmm0
ff
mulsd  -0xa8(%rbp),%xmm0
ff
movsd  0xb97(%rip),%xmm1
mulsd  -0xa8(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xb68(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
ff
movsd  %xmm0,-0xb0(%rbp)
ff
movl   $0x2,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     188f <rstage+0x76f>
movl   $0x0,-0x30(%rbp)
mov    -0x8(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0xc(%rbp),%eax
mov    -0x40(%rbp),%ecx
shl    $0x1,%ecx
sub    %ecx,%eax
add    $0x2,%eax
mov    %eax,-0x44(%rbp)
mov    -0x30(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1767 <rstage+0x647>
mov    -0x38(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0xa0(%rbp),%xmm2
ff
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x80(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0xa0(%rbp),%xmm2
ff
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x80(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0xb0(%rbp),%xmm2
ff
mov    0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0xb0(%rbp),%xmm2
ff
mov    0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x90(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x50(%rbp),%xmm0
addsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x50(%rbp),%xmm0
subsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x58(%rbp),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x58(%rbp),%xmm0
subsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
subsd  -0x68(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x68(%rbp),%xmm0
mov    0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x78(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x78(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x70(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x70(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    14f5 <rstage+0x3d5>
mov    -0x34(%rbp),%eax
shl    $0x1,%eax
sub    -0x8(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x34(%rbp),%eax
shl    $0x2,%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jl     14ec <rstage+0x3cc>
movsd  -0xa0(%rbp),%xmm0
ff
movsd  -0xa8(%rbp),%xmm2
ff
movsd  -0x80(%rbp),%xmm1
mulsd  -0x88(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0xa0(%rbp),%xmm0
ff
movsd  -0x88(%rbp),%xmm2
ff
movsd  -0x80(%rbp),%xmm1
mulsd  -0xa8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
movsd  -0xb8(%rbp),%xmm2
ff
movsd  -0x90(%rbp),%xmm1
ff
mulsd  -0x98(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0xb0(%rbp),%xmm0
ff
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x90(%rbp),%xmm1
ff
mulsd  -0xb8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    14bf <rstage+0x39f>
jmp    1894 <rstage+0x774>
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<rstage>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r9,%r14
mov    %r8,%r15
mov    %rcx,%rbx
mov    %edi,%ebp
mov    0x80(%rsp),%r12
lea    (%rsi,%rsi,1),%r9d
movslq %edi,%r13
xor    %ecx,%ecx
mov    %r9d,%eax
jmp    115c <rstage+0x3c>
nop
lea    (%rax,%rax,1),%ecx
sub    %esi,%ecx
shl    $0x2,%eax
cmp    %ebp,%ecx
jge    11b1 <rstage+0x91>
cmp    %ebp,%ecx
jge    1150 <rstage+0x30>
add    $0x1,%ecx
movslq %ecx,%rcx
movslq %eax,%rdi
nopl   0x0(%rax)
movsd  (%r12,%rcx,8),%xmm0
movsd  (%r14,%rcx,8),%xmm1
movapd %xmm0,%xmm2
addsd  %xmm1,%xmm2
subsd  %xmm1,%xmm0
movsd  %xmm0,(%r12,%rcx,8)
movsd  (%rbx,%rcx,8),%xmm0
subsd  %xmm2,%xmm0
movsd  %xmm0,(%r14,%rcx,8)
addsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,(%rbx,%rcx,8)
add    %rdi,%rcx
cmp    %r13,%rcx
jle    1170 <rstage+0x50>
jmp    1150 <rstage+0x30>
lea    0x7(%rsi),%r10d
test   %esi,%esi
cmovns %esi,%r10d
cmp    $0x2,%edx
jl     1273 <rstage+0x153>
sar    $0x3,%r10d
lea    0x1(%r10),%r8d
xor    %eax,%eax
movsd  0xe2a(%rip),%xmm0
mov    %r9d,%ecx
jmp    11f0 <rstage+0xd0>
nopl   0x0(%rax,%rax,1)
lea    (%rcx,%rcx,1),%eax
sub    %esi,%eax
shl    $0x2,%ecx
cmp    %ebp,%eax
jge    126e <rstage+0x14e>
add    %r8d,%eax
cmp    %ebp,%eax
jg     11e0 <rstage+0xc0>
movslq %eax,%rdi
movslq %ecx,%rax
nopl   (%rax)
movsd  (%r14,%rdi,8),%xmm1
movsd  (%r12,%rdi,8),%xmm2
movapd %xmm1,%xmm3
addsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm3
subsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
movsd  (%r15,%rdi,8),%xmm2
subsd  %xmm3,%xmm2
movsd  %xmm2,(%r12,%rdi,8)
movsd  (%r15,%rdi,8),%xmm2
xorpd  0xdf2(%rip),%xmm2
subsd  %xmm3,%xmm2
movsd  %xmm2,(%r14,%rdi,8)
movsd  (%rbx,%rdi,8),%xmm2
subsd  %xmm1,%xmm2
movsd  %xmm2,(%r15,%rdi,8)
addsd  (%rbx,%rdi,8),%xmm1
movsd  %xmm1,(%rbx,%rdi,8)
add    %rax,%rdi
cmp    %r13,%rdi
jle    1200 <rstage+0xe0>
jmp    11e0 <rstage+0xc0>
cmp    $0x10,%esi
jge    1282 <rstage+0x162>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xorps  %xmm1,%xmm1
cvtsi2sd %esi,%xmm1
movsd  0xd77(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rsp)
mov    %edx,0x4(%rsp)
mov    %rsi,0x18(%rsp)
mov    %r9d,(%rsp)
mov    %r10,0x10(%rsp)
call   1030 <sin@plt>
movapd %xmm0,0x30(%rsp)
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  0xd44(%rip),%xmm1
movsd  0xd44(%rip),%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
movapd %xmm2,0x20(%rsp)
movsd  0x8(%rsp),%xmm0
call   1040 <cos@plt>
movapd 0x20(%rsp),%xmm14
movapd 0xd33(%rip),%xmm13
mov    0x10(%rsp),%r11
mov    (%rsp),%r10d
mov    0x18(%rsp),%r9
mov    0x4(%rsp),%r8d
movapd %xmm0,%xmm9
mulsd  %xmm0,%xmm9
mulsd  %xmm0,%xmm9
movsd  0xcf9(%rip),%xmm2
mulsd  %xmm0,%xmm2
mulsd  0xcf4(%rip),%xmm9
cmp    $0x3,%r11d
mov    $0x2,%eax
cmovl  %eax,%r11d
addsd  %xmm2,%xmm9
unpcklpd 0x30(%rsp),%xmm14
unpcklpd %xmm0,%xmm9
movapd 0xcd5(%rip),%xmm8
movapd %xmm9,%xmm10
movapd %xmm14,%xmm11
jmp    13bd <rstage+0x29d>
nopw   0x0(%rax,%rax,1)
xorpd  %xmm8,%xmm12
movapd %xmm9,%xmm0
mulpd  %xmm11,%xmm0
xorpd  %xmm8,%xmm11
unpcklpd %xmm12,%xmm11
mulpd  %xmm14,%xmm11
movapd %xmm10,%xmm1
mulpd  %xmm14,%xmm1
addpd  %xmm0,%xmm1
mulpd  %xmm9,%xmm10
addpd  %xmm11,%xmm10
lea    0x1(%rax),%ecx
add    $0xffffffff,%r8d
movapd %xmm1,%xmm11
cmp    %r11d,%eax
mov    %ecx,%eax
je     1273 <rstage+0x153>
movapd %xmm11,%xmm12
unpckhpd %xmm11,%xmm12
movapd %xmm10,%xmm5
unpckhpd %xmm10,%xmm5
xor    %edx,%edx
mov    %r10d,%ecx
jmp    13ed <rstage+0x2cd>
nopl   0x0(%rax,%rax,1)
lea    (%rcx,%rcx,1),%edx
sub    %r9d,%edx
shl    $0x2,%ecx
cmp    %ebp,%edx
jge    1370 <rstage+0x250>
lea    (%rdx,%rax,1),%esi
cmp    %ebp,%esi
jg     13e0 <rstage+0x2c0>
movslq %esi,%rsi
movslq %ecx,%rdi
add    %r8d,%edx
nopl   (%rax)
movsd  (%r14,%rsi,8),%xmm7
movslq %edx,%rdx
movsd  (%r14,%rdx,8),%xmm6
movapd %xmm12,%xmm4
mulsd  %xmm6,%xmm4
movapd %xmm12,%xmm3
mulsd  %xmm7,%xmm3
mulsd  %xmm5,%xmm7
addsd  %xmm4,%xmm7
mulsd  %xmm5,%xmm6
subsd  %xmm3,%xmm6
movsd  (%r12,%rsi,8),%xmm3
movsd  (%r12,%rdx,8),%xmm4
movapd %xmm11,%xmm1
mulsd  %xmm4,%xmm1
movapd %xmm11,%xmm2
mulsd  %xmm3,%xmm2
mulsd  %xmm10,%xmm3
addsd  %xmm1,%xmm3
mulsd  %xmm10,%xmm4
subsd  %xmm2,%xmm4
movapd %xmm7,%xmm1
addsd  %xmm3,%xmm1
movapd %xmm6,%xmm2
addsd  %xmm4,%xmm2
movapd %xmm2,%xmm0
subsd  (%r15,%rdx,8),%xmm0
subsd  %xmm3,%xmm7
movsd  %xmm0,(%r14,%rsi,8)
addsd  (%r15,%rdx,8),%xmm2
subsd  %xmm4,%xmm6
movsd  %xmm2,(%r12,%rdx,8)
movsd  (%r15,%rsi,8),%xmm0
xorpd  %xmm13,%xmm0
subsd  %xmm7,%xmm0
movsd  %xmm0,(%r14,%rdx,8)
movsd  (%r15,%rsi,8),%xmm0
subsd  %xmm7,%xmm0
movsd  %xmm0,(%r12,%rsi,8)
movsd  (%rbx,%rsi,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%r15,%rdx,8)
addsd  (%rbx,%rsi,8),%xmm1
movsd  %xmm1,(%rbx,%rsi,8)
movsd  (%rbx,%rdx,8),%xmm0
addsd  %xmm6,%xmm0
movsd  %xmm0,(%r15,%rsi,8)
movsd  (%rbx,%rdx,8),%xmm0
subsd  %xmm6,%xmm0
movsd  %xmm0,(%rbx,%rdx,8)
add    %rdi,%rsi
add    %ecx,%edx
cmp    %r13,%rsi
jle    1400 <rstage+0x2e0>
jmp    13e0 <rstage+0x2c0>

LLVM IR (-O0):
define void @rstage(i32 %0, i32 %1, i32 %2, double* %3, double* %4, double* %5, double* %6) {
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca double*
%12 = alloca double*
%13 = alloca double*
%14 = alloca double*
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca i32
%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
store i32 %0, i32* %8
store i32 %1, i32* %9
store i32 %2, i32* %10
store double* %3, double** %11
store double* %4, double** %12
store double* %5, double** %13
store double* %6, double** %14
%37 = load i32, i32* %9
%38 = sdiv i32 %37, 8
store i32 %38, i32* %15
store i32 0, i32* %16
%39 = load i32, i32* %9
%40 = mul nsw i32 %39, 2
store i32 %40, i32* %17
br label %41
41: ; preds = %108, %7
%42 = load i32, i32* %16
%43 = add nsw i32 %42, 1
store i32 %43, i32* %18
br label %44
44: ; preds = %97, %41
%45 = load i32, i32* %18
%46 = load i32, i32* %8
%47 = icmp sle i32 %45, %46
br i1 %47, label %48, label %101
48: ; preds = %44
%49 = load double*, double** %14
%50 = load i32, i32* %18
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %49, i64 %51
%53 = load double, double* %52
%54 = load double*, double** %13
%55 = load i32, i32* %18
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds double, double* %54, i64 %56
%58 = load double, double* %57
%59 = fadd double %53, %58
store double %59, double* %23
%60 = load double*, double** %14
%61 = load i32, i32* %18
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %60, i64 %62
%64 = load double, double* %63
%65 = load double*, double** %13
%66 = load i32, i32* %18
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds double, double* %65, i64 %67
%69 = load double, double* %68
%70 = fsub double %64, %69
%71 = load double*, double** %14
%72 = load i32, i32* %18
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %71, i64 %73
store double %70, double* %74
%75 = load double*, double** %11
%76 = load i32, i32* %18
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %75, i64 %77
%79 = load double, double* %78
%80 = load double, double* %23
%81 = fsub double %79, %80
%82 = load double*, double** %13
%83 = load i32, i32* %18
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %82, i64 %84
store double %81, double* %85
%86 = load double*, double** %11
%87 = load i32, i32* %18
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds double, double* %86, i64 %88
%90 = load double, double* %89
%91 = load double, double* %23
%92 = fadd double %90, %91
%93 = load double*, double** %11
%94 = load i32, i32* %18
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds double, double* %93, i64 %95
store double %92, double* %96
br label %97
97: ; preds = %48
%98 = load i32, i32* %17
%99 = load i32, i32* %18
%100 = add nsw i32 %99, %98
store i32 %100, i32* %18
br label %44
101: ; preds = %44
%102 = load i32, i32* %17
%103 = mul nsw i32 2, %102
%104 = load i32, i32* %9
%105 = sub nsw i32 %103, %104
store i32 %105, i32* %16
%106 = load i32, i32* %17
%107 = mul nsw i32 4, %106
store i32 %107, i32* %17
br label %108
108: ; preds = %101
%109 = load i32, i32* %16
%110 = load i32, i32* %8
%111 = icmp slt i32 %109, %110
br i1 %111, label %41, label %112
112: ; preds = %108
%113 = load i32, i32* %10
%114 = sub nsw i32 %113, 1
%115 = icmp sle i32 %114, 0
br i1 %115, label %116, label %117
116: ; preds = %112
br label %480
117: ; preds = %112
store i32 0, i32* %16
%118 = load i32, i32* %9
%119 = mul nsw i32 %118, 2
store i32 %119, i32* %17
br label %120
120: ; preds = %210, %117
%121 = load i32, i32* %16
%122 = add nsw i32 %121, 1
%123 = load i32, i32* %15
%124 = add nsw i32 %122, %123
store i32 %124, i32* %19
br label %125
125: ; preds = %199, %120
%126 = load i32, i32* %19
%127 = load i32, i32* %8
%128 = icmp sle i32 %126, %127
br i1 %128, label %129, label %203
129: ; preds = %125
%130 = load double*, double** %13
%131 = load i32, i32* %19
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds double, double* %130, i64 %132
%134 = load double, double* %133
%135 = load double*, double** %14
%136 = load i32, i32* %19
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds double, double* %135, i64 %137
%139 = load double, double* %138
%140 = fadd double %134, %139
%141 = fmul double %140, 0x3FE6A09E667F3BCC
store double %141, double* %23
%142 = load double*, double** %13
%143 = load i32, i32* %19
%144 = sext i32 %143 to i64
%145 = getelementptr inbounds double, double* %142, i64 %144
%146 = load double, double* %145
%147 = load double*, double** %14
%148 = load i32, i32* %19
%149 = sext i32 %148 to i64
%150 = getelementptr inbounds double, double* %147, i64 %149
%151 = load double, double* %150
%152 = fsub double %146, %151
%153 = fmul double %152, 0x3FE6A09E667F3BCC
store double %153, double* %24
%154 = load double*, double** %12
%155 = load i32, i32* %19
%156 = sext i32 %155 to i64
%157 = getelementptr inbounds double, double* %154, i64 %156
%158 = load double, double* %157
%159 = load double, double* %23
%160 = fsub double %158, %159
%161 = load double*, double** %14
%162 = load i32, i32* %19
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds double, double* %161, i64 %163
store double %160, double* %164
%165 = load double*, double** %12
%166 = load i32, i32* %19
%167 = sext i32 %166 to i64
%168 = getelementptr inbounds double, double* %165, i64 %167
%169 = load double, double* %168
%170 = fneg double %169
%171 = load double, double* %23
%172 = fsub double %170, %171
%173 = load double*, double** %13
%174 = load i32, i32* %19
%175 = sext i32 %174 to i64
%176 = getelementptr inbounds double, double* %173, i64 %175
store double %172, double* %176
%177 = load double*, double** %11
%178 = load i32, i32* %19
%179 = sext i32 %178 to i64
%180 = getelementptr inbounds double, double* %177, i64 %179
%181 = load double, double* %180
%182 = load double, double* %24
%183 = fsub double %181, %182
%184 = load double*, double** %12
%185 = load i32, i32* %19
%186 = sext i32 %185 to i64
%187 = getelementptr inbounds double, double* %184, i64 %186
store double %183, double* %187
%188 = load double*, double** %11
%189 = load i32, i32* %19
%190 = sext i32 %189 to i64
%191 = getelementptr inbounds double, double* %188, i64 %190
%192 = load double, double* %191
%193 = load double, double* %24
%194 = fadd double %192, %193
%195 = load double*, double** %11
%196 = load i32, i32* %19
%197 = sext i32 %196 to i64
%198 = getelementptr inbounds double, double* %195, i64 %197
store double %194, double* %198
br label %199
199: ; preds = %129
%200 = load i32, i32* %17
%201 = load i32, i32* %19
%202 = add nsw i32 %201, %200
store i32 %202, i32* %19
br label %125
203: ; preds = %125
%204 = load i32, i32* %17
%205 = mul nsw i32 2, %204
%206 = load i32, i32* %9
%207 = sub nsw i32 %205, %206
store i32 %207, i32* %16
%208 = load i32, i32* %17
%209 = mul nsw i32 4, %208
store i32 %209, i32* %17
br label %210
210: ; preds = %203
%211 = load i32, i32* %16
%212 = load i32, i32* %8
%213 = icmp slt i32 %211, %212
br i1 %213, label %120, label %214
214: ; preds = %210
%215 = load i32, i32* %15
%216 = sub nsw i32 %215, 1
%217 = icmp sle i32 %216, 0
br i1 %217, label %218, label %219
218: ; preds = %214
br label %480
219: ; preds = %214
%220 = load i32, i32* %9
%221 = sitofp i32 %220 to double
%222 = fdiv double 0x401921FB54442D18, %221
store double %222, double* %25
%223 = load double, double* %25
%224 = call double @sin(double %223)
store double %224, double* %29
%225 = load double, double* %29
store double %225, double* %30
%226 = load double, double* %30
%227 = load double, double* %30
%228 = load double, double* %30
%229 = fmul double %227, %228
%230 = load double, double* %30
%231 = fmul double %229, %230
%232 = fmul double 4.000000e+00, %231
%233 = fneg double %232
%234 = call double @llvm.fmuladd.f64(double 3.000000e+00, double %226, double %233)
store double %234, double* %32
%235 = load double, double* %32
store double %235, double* %31
%236 = load double, double* %25
%237 = call double @cos(double %236)
store double %237, double* %33
%238 = load double, double* %33
store double %238, double* %34
%239 = load double, double* %34
%240 = load double, double* %34
%241 = fmul double %239, %240
%242 = load double, double* %34
%243 = fmul double %241, %242
%244 = load double, double* %34
%245 = fmul double 3.000000e+00, %244
%246 = fneg double %245
%247 = call double @llvm.fmuladd.f64(double 4.000000e+00, double %243, double %246)
store double %247, double* %36
%248 = load double, double* %36
store double %248, double* %35
store i32 2, i32* %20
br label %249
249: ; preds = %476, %219
%250 = load i32, i32* %20
%251 = load i32, i32* %15
%252 = icmp sle i32 %250, %251
br i1 %252, label %253, label %479
253: ; preds = %249
store i32 0, i32* %16
%254 = load i32, i32* %9
%255 = mul nsw i32 2, %254
store i32 %255, i32* %17
%256 = load i32, i32* %10
%257 = load i32, i32* %20
%258 = mul nsw i32 2, %257
%259 = sub nsw i32 %256, %258
%260 = add nsw i32 %259, 2
store i32 %260, i32* %21
br label %261
261: ; preds = %443, %253
%262 = load i32, i32* %16
%263 = load i32, i32* %20
%264 = add nsw i32 %262, %263
store i32 %264, i32* %18
br label %265
265: ; preds = %432, %261
%266 = load i32, i32* %18
%267 = load i32, i32* %8
%268 = icmp sle i32 %266, %267
br i1 %268, label %269, label %436
269: ; preds = %265
%270 = load i32, i32* %18
%271 = load i32, i32* %21
%272 = add nsw i32 %270, %271
store i32 %272, i32* %19
%273 = load double*, double** %13
%274 = load i32, i32* %18
%275 = sext i32 %274 to i64
%276 = getelementptr inbounds double, double* %273, i64 %275
%277 = load double, double* %276
%278 = load double, double* %33
%279 = load double*, double** %13
%280 = load i32, i32* %19
%281 = sext i32 %280 to i64
%282 = getelementptr inbounds double, double* %279, i64 %281
%283 = load double, double* %282
%284 = load double, double* %29
%285 = fmul double %283, %284
%286 = call double @llvm.fmuladd.f64(double %277, double %278, double %285)
store double %286, double* %23
%287 = load double*, double** %13
%288 = load i32, i32* %19
%289 = sext i32 %288 to i64
%290 = getelementptr inbounds double, double* %287, i64 %289
%291 = load double, double* %290
%292 = load double, double* %33
%293 = load double*, double** %13
%294 = load i32, i32* %18
%295 = sext i32 %294 to i64
%296 = getelementptr inbounds double, double* %293, i64 %295
%297 = load double, double* %296
%298 = load double, double* %29
%299 = fmul double %297, %298
%300 = fneg double %299
%301 = call double @llvm.fmuladd.f64(double %291, double %292, double %300)
store double %301, double* %24
%302 = load double*, double** %14
%303 = load i32, i32* %18
%304 = sext i32 %303 to i64
%305 = getelementptr inbounds double, double* %302, i64 %304
%306 = load double, double* %305
%307 = load double, double* %35
%308 = load double*, double** %14
%309 = load i32, i32* %19
%310 = sext i32 %309 to i64
%311 = getelementptr inbounds double, double* %308, i64 %310
%312 = load double, double* %311
%313 = load double, double* %31
%314 = fmul double %312, %313
%315 = call double @llvm.fmuladd.f64(double %306, double %307, double %314)
store double %315, double* %26
%316 = load double*, double** %14
%317 = load i32, i32* %19
%318 = sext i32 %317 to i64
%319 = getelementptr inbounds double, double* %316, i64 %318
%320 = load double, double* %319
%321 = load double, double* %35
%322 = load double*, double** %14
%323 = load i32, i32* %18
%324 = sext i32 %323 to i64
%325 = getelementptr inbounds double, double* %322, i64 %324
%326 = load double, double* %325
%327 = load double, double* %31
%328 = fmul double %326, %327
%329 = fneg double %328
%330 = call double @llvm.fmuladd.f64(double %320, double %321, double %329)
store double %330, double* %27
%331 = load double, double* %23
%332 = load double, double* %26
%333 = fadd double %331, %332
store double %333, double* %28
%334 = load double, double* %23
%335 = load double, double* %26
%336 = fsub double %334, %335
store double %336, double* %26
%337 = load double, double* %24
%338 = load double, double* %27
%339 = fadd double %337, %338
store double %339, double* %23
%340 = load double, double* %24
%341 = load double, double* %27
%342 = fsub double %340, %341
store double %342, double* %27
%343 = load double, double* %23
%344 = load double*, double** %12
%345 = load i32, i32* %19
%346 = sext i32 %345 to i64
%347 = getelementptr inbounds double, double* %344, i64 %346
%348 = load double, double* %347
%349 = fsub double %343, %348
%350 = load double*, double** %13
%351 = load i32, i32* %18
%352 = sext i32 %351 to i64
%353 = getelementptr inbounds double, double* %350, i64 %352
store double %349, double* %353
%354 = load double, double* %23
%355 = load double*, double** %12
%356 = load i32, i32* %19
%357 = sext i32 %356 to i64
%358 = getelementptr inbounds double, double* %355, i64 %357
%359 = load double, double* %358
%360 = fadd double %354, %359
%361 = load double*, double** %14
%362 = load i32, i32* %19
%363 = sext i32 %362 to i64
%364 = getelementptr inbounds double, double* %361, i64 %363
store double %360, double* %364
%365 = load double*, double** %12
%366 = load i32, i32* %18
%367 = sext i32 %366 to i64
%368 = getelementptr inbounds double, double* %365, i64 %367
%369 = load double, double* %368
%370 = fneg double %369
%371 = load double, double* %26
%372 = fsub double %370, %371
%373 = load double*, double** %13
%374 = load i32, i32* %19
%375 = sext i32 %374 to i64
%376 = getelementptr inbounds double, double* %373, i64 %375
store double %372, double* %376
%377 = load double*, double** %12
%378 = load i32, i32* %18
%379 = sext i32 %378 to i64
%380 = getelementptr inbounds double, double* %377, i64 %379
%381 = load double, double* %380
%382 = load double, double* %26
%383 = fsub double %381, %382
%384 = load double*, double** %14
%385 = load i32, i32* %18
%386 = sext i32 %385 to i64
%387 = getelementptr inbounds double, double* %384, i64 %386
store double %383, double* %387
%388 = load double*, double** %11
%389 = load i32, i32* %18
%390 = sext i32 %389 to i64
%391 = getelementptr inbounds double, double* %388, i64 %390
%392 = load double, double* %391
%393 = load double, double* %28
%394 = fsub double %392, %393
%395 = load double*, double** %12
%396 = load i32, i32* %19
%397 = sext i32 %396 to i64
%398 = getelementptr inbounds double, double* %395, i64 %397
store double %394, double* %398
%399 = load double*, double** %11
%400 = load i32, i32* %18
%401 = sext i32 %400 to i64
%402 = getelementptr inbounds double, double* %399, i64 %401
%403 = load double, double* %402
%404 = load double, double* %28
%405 = fadd double %403, %404
%406 = load double*, double** %11
%407 = load i32, i32* %18
%408 = sext i32 %407 to i64
%409 = getelementptr inbounds double, double* %406, i64 %408
store double %405, double* %409
%410 = load double*, double** %11
%411 = load i32, i32* %19
%412 = sext i32 %411 to i64
%413 = getelementptr inbounds double, double* %410, i64 %412
%414 = load double, double* %413
%415 = load double, double* %27
%416 = fadd double %414, %415
%417 = load double*, double** %12
%418 = load i32, i32* %18
%419 = sext i32 %418 to i64
%420 = getelementptr inbounds double, double* %417, i64 %419
store double %416, double* %420
%421 = load double*, double** %11
%422 = load i32, i32* %19
%423 = sext i32 %422 to i64
%424 = getelementptr inbounds double, double* %421, i64 %423
%425 = load double, double* %424
%426 = load double, double* %27
%427 = fsub double %425, %426
%428 = load double*, double** %11
%429 = load i32, i32* %19
%430 = sext i32 %429 to i64
%431 = getelementptr inbounds double, double* %428, i64 %430
store double %427, double* %431
br label %432
432: ; preds = %269
%433 = load i32, i32* %17
%434 = load i32, i32* %18
%435 = add nsw i32 %434, %433
store i32 %435, i32* %18
br label %265
436: ; preds = %265
%437 = load i32, i32* %17
%438 = mul nsw i32 2, %437
%439 = load i32, i32* %9
%440 = sub nsw i32 %438, %439
store i32 %440, i32* %16
%441 = load i32, i32* %17
%442 = mul nsw i32 4, %441
store i32 %442, i32* %17
br label %443
443: ; preds = %436
%444 = load i32, i32* %16
%445 = load i32, i32* %8
%446 = icmp slt i32 %444, %445
br i1 %446, label %261, label %447
447: ; preds = %443
%448 = load double, double* %33
%449 = load double, double* %34
%450 = load double, double* %29
%451 = load double, double* %30
%452 = fmul double %450, %451
%453 = fneg double %452
%454 = call double @llvm.fmuladd.f64(double %448, double %449, double %453)
store double %454, double* %23
%455 = load double, double* %33
%456 = load double, double* %30
%457 = load double, double* %29
%458 = load double, double* %34
%459 = fmul double %457, %458
%460 = call double @llvm.fmuladd.f64(double %455, double %456, double %459)
store double %460, double* %29
%461 = load double, double* %23
store double %461, double* %33
%462 = load double, double* %35
%463 = load double, double* %36
%464 = load double, double* %31
%465 = load double, double* %32
%466 = fmul double %464, %465
%467 = fneg double %466
%468 = call double @llvm.fmuladd.f64(double %462, double %463, double %467)
store double %468, double* %26
%469 = load double, double* %35
%470 = load double, double* %32
%471 = load double, double* %31
%472 = load double, double* %36
%473 = fmul double %471, %472
%474 = call double @llvm.fmuladd.f64(double %469, double %470, double %473)
store double %474, double* %31
%475 = load double, double* %26
store double %475, double* %35
br label %476
476: ; preds = %447
%477 = load i32, i32* %20
%478 = add nsw i32 %477, 1
store i32 %478, i32* %20
br label %249
479: ; preds = %249
br label %480
480: ; preds = %479, %218, %116
ret void
}

LLVM IR (-O2):
define void @rstage(i32 %0, i32 %1, i32 %2, double* nocapture %3, double* nocapture %4, double* nocapture %5, double* nocapture %6) {
%8 = shl nsw i32 %1, 1
%9 = sext i32 %0 to i64
br label %10
10: ; preds = %33, %7
%11 = phi i32 [ %8, %7 ], [ %36, %33 ]
%12 = phi i32 [ 0, %7 ], [ %35, %33 ]
%13 = icmp slt i32 %12, %0
br i1 %13, label %14, label %33
14: ; preds = %10
%15 = add nsw i32 %12, 1
%16 = sext i32 %15 to i64
%17 = sext i32 %11 to i64
br label %18
18: ; preds = %14, %18
%19 = phi i64 [ %16, %14 ], [ %31, %18 ]
%20 = getelementptr inbounds double, double* %6, i64 %19
%21 = load double, double* %20
%22 = getelementptr inbounds double, double* %5, i64 %19
%23 = load double, double* %22
%24 = fadd double %21, %23
%25 = fsub double %21, %23
store double %25, double* %20
%26 = getelementptr inbounds double, double* %3, i64 %19
%27 = load double, double* %26
%28 = fsub double %27, %24
store double %28, double* %22
%29 = load double, double* %26
%30 = fadd double %24, %29
store double %30, double* %26
%31 = add i64 %19, %17
%32 = icmp sgt i64 %31, %9
br i1 %32, label %33, label %18
33: ; preds = %18, %10
%34 = shl nsw i32 %11, 1
%35 = sub nsw i32 %34, %1
%36 = shl nsw i32 %11, 2
%37 = icmp slt i32 %35, %0
br i1 %37, label %10, label %38
38: ; preds = %33
%39 = sdiv i32 %1, 8
%40 = icmp slt i32 %2, 2
br i1 %40, label %189, label %41
41: ; preds = %38
%42 = add nsw i32 %39, 1
br label %43
43: ; preds = %41, %74
%44 = phi i32 [ %77, %74 ], [ %8, %41 ]
%45 = phi i32 [ %76, %74 ], [ 0, %41 ]
%46 = add i32 %42, %45
%47 = icmp sgt i32 %46, %0
br i1 %47, label %74, label %48
48: ; preds = %43
%49 = sext i32 %46 to i64
%50 = sext i32 %44 to i64
br label %51
51: ; preds = %48, %51
%52 = phi i64 [ %49, %48 ], [ %72, %51 ]
%53 = getelementptr inbounds double, double* %5, i64 %52
%54 = load double, double* %53
%55 = getelementptr inbounds double, double* %6, i64 %52
%56 = load double, double* %55
%57 = fadd double %54, %56
%58 = fmul double %57, 0x3FE6A09E667F3BCC
%59 = fsub double %54, %56
%60 = fmul double %59, 0x3FE6A09E667F3BCC
%61 = getelementptr inbounds double, double* %4, i64 %52
%62 = load double, double* %61
%63 = fsub double %62, %58
store double %63, double* %55
%64 = load double, double* %61
%65 = fneg double %64
%66 = fsub double %65, %58
store double %66, double* %53
%67 = getelementptr inbounds double, double* %3, i64 %52
%68 = load double, double* %67
%69 = fsub double %68, %60
store double %69, double* %61
%70 = load double, double* %67
%71 = fadd double %60, %70
store double %71, double* %67
%72 = add i64 %52, %50
%73 = icmp sgt i64 %72, %9
br i1 %73, label %74, label %51
74: ; preds = %51, %43
%75 = shl nsw i32 %44, 1
%76 = sub nsw i32 %75, %1
%77 = shl nsw i32 %44, 2
%78 = icmp slt i32 %76, %0
br i1 %78, label %43, label %79
79: ; preds = %74
%80 = icmp slt i32 %1, 16
br i1 %80, label %189, label %81
81: ; preds = %79
%82 = sitofp i32 %1 to double
%83 = fdiv double 0x401921FB54442D18, %82
%84 = tail call double @sin(double %83)
%85 = fmul double %84, %84
%86 = fmul double %84, %85
%87 = fmul double %86, -4.000000e+00
%88 = tail call double @llvm.fmuladd.f64(double %84, double 3.000000e+00, double %87)
%89 = tail call double @cos(double %83)
%90 = fmul double %89, %89
%91 = fmul double %89, %90
%92 = fmul double %89, -3.000000e+00
%93 = tail call double @llvm.fmuladd.f64(double %91, double 4.000000e+00, double %92)
%94 = add i32 %2, 2
%95 = call i32 @llvm.smax.i32(i32 %39, i32 2)
%96 = insertelement <4 x double> poison, double %93, i64 0
%97 = insertelement <4 x double> %96, double %89, i64 1
%98 = insertelement <4 x double> %97, double %88, i64 2
%99 = insertelement <4 x double> %98, double %84, i64 3
%100 = insertelement <4 x double> poison, double %88, i64 0
%101 = insertelement <4 x double> %100, double %84, i64 1
%102 = insertelement <4 x double> %101, double %93, i64 2
%103 = insertelement <4 x double> %102, double %89, i64 3
br label %104
104: ; preds = %81, %176
%105 = phi i32 [ 2, %81 ], [ %187, %176 ]
%106 = phi <4 x double> [ %99, %81 ], [ %186, %176 ]
%107 = mul nsw i32 %105, -2
%108 = add i32 %94, %107
%109 = extractelement <4 x double> %106, i64 3
%110 = extractelement <4 x double> %106, i64 1
%111 = extractelement <4 x double> %106, i64 2
%112 = extractelement <4 x double> %106, i64 0
br label %113
113: ; preds = %171, %104
%114 = phi i32 [ %8, %104 ], [ %174, %171 ]
%115 = phi i32 [ 0, %104 ], [ %173, %171 ]
%116 = add i32 %115, %105
%117 = icmp sgt i32 %116, %0
br i1 %117, label %171, label %118
118: ; preds = %113
%119 = sext i32 %116 to i64
%120 = sext i32 %114 to i64
br label %121
121: ; preds = %118, %121
%122 = phi i64 [ %119, %118 ], [ %169, %121 ]
%123 = trunc i64 %122 to i32
%124 = add i32 %108, %123
%125 = getelementptr inbounds double, double* %5, i64 %122
%126 = load double, double* %125
%127 = sext i32 %124 to i64
%128 = getelementptr inbounds double, double* %5, i64 %127
%129 = load double, double* %128
%130 = fmul double %109, %129
%131 = tail call double @llvm.fmuladd.f64(double %126, double %110, double %130)
%132 = fneg double %126
%133 = fmul double %109, %132
%134 = tail call double @llvm.fmuladd.f64(double %129, double %110, double %133)
%135 = getelementptr inbounds double, double* %6, i64 %122
%136 = load double, double* %135
%137 = getelementptr inbounds double, double* %6, i64 %127
%138 = load double, double* %137
%139 = fmul double %111, %138
%140 = tail call double @llvm.fmuladd.f64(double %136, double %112, double %139)
%141 = fneg double %136
%142 = fmul double %111, %141
%143 = tail call double @llvm.fmuladd.f64(double %138, double %112, double %142)
%144 = fadd double %131, %140
%145 = fsub double %131, %140
%146 = fadd double %134, %143
%147 = fsub double %134, %143
%148 = getelementptr inbounds double, double* %4, i64 %127
%149 = load double, double* %148
%150 = fsub double %146, %149
store double %150, double* %125
%151 = load double, double* %148
%152 = fadd double %146, %151
store double %152, double* %137
%153 = getelementptr inbounds double, double* %4, i64 %122
%154 = load double, double* %153
%155 = fneg double %154
%156 = fsub double %155, %145
store double %156, double* %128
%157 = load double, double* %153
%158 = fsub double %157, %145
store double %158, double* %135
%159 = getelementptr inbounds double, double* %3, i64 %122
%160 = load double, double* %159
%161 = fsub double %160, %144
store double %161, double* %148
%162 = load double, double* %159
%163 = fadd double %144, %162
store double %163, double* %159
%164 = getelementptr inbounds double, double* %3, i64 %127
%165 = load double, double* %164
%166 = fadd double %147, %165
store double %166, double* %153
%167 = load double, double* %164
%168 = fsub double %167, %147
store double %168, double* %164
%169 = add i64 %122, %120
%170 = icmp sgt i64 %169, %9
br i1 %170, label %171, label %121
171: ; preds = %121, %113
%172 = shl nsw i32 %114, 1
%173 = sub nsw i32 %172, %1
%174 = shl nsw i32 %114, 2
%175 = icmp slt i32 %173, %0
br i1 %175, label %113, label %176
176: ; preds = %171
%177 = extractelement <4 x double> %106, i64 3
%178 = fneg double %177
%179 = extractelement <4 x double> %106, i64 2
%180 = fneg double %179
%181 = insertelement <4 x double> poison, double %180, i64 0
%182 = insertelement <4 x double> %181, double %178, i64 1
%183 = shufflevector <4 x double> %182, <4 x double> %106, <4 x i32> <i32 0, i32 1, i32 6, i32 7>
%184 = fmul <4 x double> %103, %183
%185 = shufflevector <4 x double> %106, <4 x double> poison, <4 x i32> <i32 0, i32 1, i32 0, i32 1>
%186 = call <4 x double> @llvm.fmuladd.v4f64(<4 x double> %185, <4 x double> %99, <4 x double> %184)
%187 = add nuw nsw i32 %105, 1
%188 = icmp eq i32 %105, %95
br i1 %188, label %189, label %104
189: ; preds = %176, %79, %38
ret void
}