void fourn(data, nn, ndim, isign) float data[];
int nn[], ndim, isign;
{
  int i1, i2, i3, i2rev, i3rev, ip1, ip2, ip3, ifp1, ifp2;
  int ibit, idim, k1, k2, n, nprev, nrem, ntot;
  float tempi, tempr;
  double theta, wi, wpi, wpr, wr, wtemp;
  ntot = 1;
  for (idim = 1; idim <= ndim; idim++)
    ntot *= nn[idim];
  nprev = 1;
  for (idim = ndim; idim >= 1; idim--) {
    n = nn[idim];
    nrem = ntot / (n * nprev);
    ip1 = nprev << 1;
    ip2 = ip1 * n;
    ip3 = ip2 * nrem;
    i2rev = 1;
    for (i2 = 1; i2 <= ip2; i2 += ip1) {
      if (i2 < i2rev) {
        for (i1 = i2; i1 <= i2 + ip1 - 2; i1 += 2) {
          for (i3 = i1; i3 <= ip3; i3 += ip2) {
            i3rev = i2rev + i3 - i2;
            tempr = (data[i3]);
            (data[i3]) = (data[i3rev]);
            (data[i3rev]) = tempr;
            tempr = (data[i3 + 1]);
            (data[i3 + 1]) = (data[i3rev + 1]);
            (data[i3rev + 1]) = tempr;
          }
        }
      }
      ibit = ip2 >> 1;
      while (ibit >= ip1 && i2rev > ibit) {
        i2rev -= ibit;
        ibit >>= 1;
      }
      i2rev += ibit;
    }
    ifp1 = ip1;
    while (ifp1 < ip2) {
      ifp2 = ifp1 << 1;
      theta = isign * 6.28318530717959 / (ifp2 / ip1);
      wtemp = sin(0.5 * theta);
      wpr = -2.0 * wtemp * wtemp;
      wpi = sin(theta);
      wr = 1.0;
      wi = 0.0;
      for (i3 = 1; i3 <= ifp1; i3 += ip1) {
        for (i1 = i3; i1 <= i3 + ip1 - 2; i1 += 2) {
          for (i2 = i1; i2 <= ip3; i2 += ifp2) {
            k1 = i2;
            k2 = k1 + ifp1;
            tempr = wr * data[k2] - wi * data[k2 + 1];
            tempi = wr * data[k2 + 1] + wi * data[k2];
            data[k2] = data[k1] - tempr;
            data[k2 + 1] = data[k1 + 1] - tempi;
            data[k1] += tempr;
            data[k1 + 1] += tempi;
          }
        }
        wr = (wtemp = wr) * wpr - wi * wpi + wr;
        wi = wi * wpr + wtemp * wpi + wi;
      }
      ifp1 = ifp2;
    }
    nprev *= n;
  }
}


x86 assembly (clang, -O0):
<fourn>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x1,-0x60(%rbp)
movl   $0x1,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1163 <fourn+0x53>
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
imul   -0x60(%rbp),%eax
mov    %eax,-0x60(%rbp)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1137 <fourn+0x27>
movl   $0x1,-0x58(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x48(%rbp)
cmpl   $0x1,-0x48(%rbp)
jl     161a <fourn+0x50a>
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x54(%rbp)
mov    -0x60(%rbp),%eax
mov    -0x54(%rbp),%ecx
imul   -0x58(%rbp),%ecx
cltd
idiv   %ecx
mov    %eax,-0x5c(%rbp)
mov    -0x58(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
imul   -0x54(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
imul   -0x5c(%rbp),%eax
mov    %eax,-0x38(%rbp)
movl   $0x1,-0x28(%rbp)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x34(%rbp),%eax
jg     1338 <fourn+0x228>
mov    -0x20(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    12c9 <fourn+0x1b9>
mov    -0x20(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x20(%rbp),%ecx
add    -0x30(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jg     12c4 <fourn+0x1b4>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x38(%rbp),%eax
jg     12b1 <fourn+0x1a1>
mov    -0x28(%rbp),%eax
add    -0x24(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x68(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x68(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x68(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x68(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    -0x24(%rbp),%eax
mov    %eax,-0x24(%rbp)
jmp    11fb <fourn+0xeb>
jmp    12b6 <fourn+0x1a6>
mov    -0x1c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x1c(%rbp)
jmp    11e1 <fourn+0xd1>
jmp    12c9 <fourn+0x1b9>
mov    -0x34(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%ecx
xor    %eax,%eax
cmp    -0x30(%rbp),%ecx
mov    %al,-0x99(%rbp)
jl     12f5 <fourn+0x1e5>
mov    -0x28(%rbp),%eax
cmp    -0x44(%rbp),%eax
setg   %al
mov    %al,-0x99(%rbp)
mov    -0x99(%rbp),%al
test   $0x1,%al
jne    1308 <fourn+0x1f8>
jmp    1321 <fourn+0x211>
mov    -0x44(%rbp),%ecx
mov    -0x28(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x28(%rbp)
mov    -0x44(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    12d2 <fourn+0x1c2>
mov    -0x44(%rbp),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x30(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x20(%rbp)
jmp    11c3 <fourn+0xb3>
mov    -0x30(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1602 <fourn+0x4f2>
mov    -0x3c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x40(%rbp)
cvtsi2sdl -0x18(%rbp),%xmm0
movsd  0xcb8(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x40(%rbp),%eax
cltd
idivl  -0x30(%rbp)
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  0xc90(%rip),%xmm0
mulsd  -0x70(%rbp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,-0x98(%rbp)
ff
movsd  0xc6e(%rip),%xmm0
mulsd  -0x98(%rbp),%xmm0
ff
mulsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,-0x80(%rbp)
movsd  0xc37(%rip),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x78(%rbp)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jg     15f7 <fourn+0x4e7>
mov    -0x24(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    -0x30(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jg     156c <fourn+0x45c>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x38(%rbp),%eax
jg     1559 <fourn+0x449>
mov    -0x20(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x50(%rbp)
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
movslq -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  -0x78(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm3,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x68(%rbp)
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  -0x78(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x68(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x50(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x64(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x68(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x4c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x64(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x20(%rbp)
jmp    140c <fourn+0x2fc>
jmp    155e <fourn+0x44e>
mov    -0x1c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x1c(%rbp)
jmp    13f2 <fourn+0x2e2>
movsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x88(%rbp),%xmm2
ff
movsd  -0x78(%rbp),%xmm1
mulsd  -0x80(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
movsd  -0x88(%rbp),%xmm2
ff
movsd  -0x98(%rbp),%xmm1
ff
mulsd  -0x80(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x30(%rbp),%eax
add    -0x24(%rbp),%eax
mov    %eax,-0x24(%rbp)
jmp    13e0 <fourn+0x2d0>
mov    -0x40(%rbp),%eax
mov    %eax,-0x3c(%rbp)
jmp    133e <fourn+0x22e>
mov    -0x54(%rbp),%eax
imul   -0x58(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x48(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x48(%rbp)
jmp    1170 <fourn+0x60>
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<fourn>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
test   %edx,%edx
jle    1630 <fourn+0x520>
mov    %rdi,%rbx
mov    %edx,%r9d
cmp    $0x8,%edx
jae    1141 <fourn+0x31>
mov    $0x1,%eax
mov    $0x1,%r10d
jmp    130e <fourn+0x1fe>
mov    %r9d,%r10d
and    $0xfffffff8,%r10d
lea    -0x8(%r10),%rdi
mov    %rdi,%rbp
shr    $0x3,%rbp
add    $0x1,%rbp
mov    %ebp,%r8d
and    $0x3,%r8d
cmp    $0x18,%rdi
jae    1177 <fourn+0x67>
movdqa 0xe94(%rip),%xmm0
xor    %edi,%edi
movdqa %xmm0,%xmm1
jmp    125c <fourn+0x14c>
and    $0xfffffffffffffffc,%rbp
movdqa 0xe7d(%rip),%xmm0
xor    %edi,%edi
movdqa %xmm0,%xmm1
nopl   0x0(%rax)
movdqu 0x4(%rsi,%rdi,4),%xmm3
movdqu 0x14(%rsi,%rdi,4),%xmm2
movdqu 0x24(%rsi,%rdi,4),%xmm4
movdqu 0x34(%rsi,%rdi,4),%xmm5
pshufd $0xf5,%xmm3,%xmm6
pshufd $0xf5,%xmm0,%xmm7
pmuludq %xmm6,%xmm7
pmuludq %xmm0,%xmm3
pshufd $0xf5,%xmm2,%xmm0
pshufd $0xf5,%xmm1,%xmm6
pmuludq %xmm0,%xmm6
pmuludq %xmm1,%xmm2
pmuludq %xmm4,%xmm3
pshufd $0xf5,%xmm4,%xmm0
pmuludq %xmm7,%xmm0
pmuludq %xmm5,%xmm2
pshufd $0xf5,%xmm5,%xmm1
pmuludq %xmm6,%xmm1
movdqu 0x44(%rsi,%rdi,4),%xmm4
movdqu 0x54(%rsi,%rdi,4),%xmm5
pshufd $0xf5,%xmm4,%xmm6
pmuludq %xmm0,%xmm6
pmuludq %xmm4,%xmm3
pshufd $0xf5,%xmm5,%xmm4
pmuludq %xmm1,%xmm4
pmuludq %xmm5,%xmm2
movdqu 0x64(%rsi,%rdi,4),%xmm1
movdqu 0x74(%rsi,%rdi,4),%xmm5
pmuludq %xmm1,%xmm3
pshufd $0xe8,%xmm3,%xmm0
pshufd $0xf5,%xmm1,%xmm1
pmuludq %xmm6,%xmm1
pshufd $0xe8,%xmm1,%xmm1
punpckldq %xmm1,%xmm0
pmuludq %xmm5,%xmm2
pshufd $0xe8,%xmm2,%xmm1
pshufd $0xf5,%xmm5,%xmm2
pmuludq %xmm4,%xmm2
pshufd $0xe8,%xmm2,%xmm2
punpckldq %xmm2,%xmm1
add    $0x20,%rdi
add    $0xfffffffffffffffc,%rbp
jne    1190 <fourn+0x80>
movdqa %xmm0,%xmm2
movdqa %xmm1,%xmm3
test   %r8,%r8
je     12d9 <fourn+0x1c9>
nopl   0x0(%rax)
lea    0x0(,%rdi,4),%rbp
or     $0x4,%rbp
movdqu (%rsi,%rbp,1),%xmm2
movdqu 0x10(%rsi,%rbp,1),%xmm3
pshufd $0xf5,%xmm2,%xmm4
pmuludq %xmm0,%xmm2
pshufd $0xe8,%xmm2,%xmm2
pshufd $0xf5,%xmm0,%xmm0
pmuludq %xmm4,%xmm0
pshufd $0xe8,%xmm0,%xmm0
punpckldq %xmm0,%xmm2
pshufd $0xf5,%xmm3,%xmm0
pmuludq %xmm1,%xmm3
pshufd $0xe8,%xmm3,%xmm3
pshufd $0xf5,%xmm1,%xmm1
pmuludq %xmm0,%xmm1
pshufd $0xe8,%xmm1,%xmm0
punpckldq %xmm0,%xmm3
add    $0x8,%rdi
movdqa %xmm2,%xmm0
movdqa %xmm3,%xmm1
add    $0xffffffffffffffff,%r8
jne    1270 <fourn+0x160>
pshufd $0xf5,%xmm2,%xmm0
pshufd $0xf5,%xmm3,%xmm1
pmuludq %xmm0,%xmm1
pmuludq %xmm2,%xmm3
pshufd $0xee,%xmm3,%xmm0
pmuludq %xmm3,%xmm0
pshufd $0xaa,%xmm1,%xmm2
pmuludq %xmm1,%xmm2
pmuludq %xmm0,%xmm2
movd   %xmm2,%eax
cmp    %r10,%r9
je     132e <fourn+0x21e>
or     $0x1,%r10
lea    0x1(%r9),%rdi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
imul   (%rsi,%r10,4),%eax
add    $0x1,%r10
cmp    %r10,%rdi
jne    1320 <fourn+0x210>
test   %edx,%edx
jle    1630 <fourn+0x520>
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
mulsd  0xcdb(%rip),%xmm0
movsd  %xmm0,0x28(%rsp)
mov    $0x1,%ecx
mov    %rsi,0x18(%rsp)
mov    %eax,0x8(%rsp)
jmp    1380 <fourn+0x270>
nopl   0x0(%rax,%rax,1)
mov    0x20(%rsp),%rcx
lea    -0x1(%rcx),%r9
cmp    $0x1,%rcx
mov    0x18(%rsp),%rsi
mov    0x8(%rsp),%eax
mov    0xc(%rsp),%ecx
jle    1630 <fourn+0x520>
mov    %ecx,%r11d
mov    %r9,0x20(%rsp)
mov    (%rsi,%r9,4),%r13d
mov    %r13d,%ecx
imul   %r11d,%ecx
cltd
mov    %ecx,0xc(%rsp)
idiv   %ecx
mov    %eax,%r14d
add    %r11d,%r11d
imul   %r11d,%r13d
imul   %r13d,%r14d
movslq %r11d,%rax
mov    %rax,0x10(%rsp)
movslq %r14d,%r12
test   %r13d,%r13d
jle    1490 <fourn+0x380>
mov    %r13d,%eax
mov    $0x1,%edx
mov    $0x1,%r8d
jmp    13e2 <fourn+0x2d2>
nopl   0x0(%rax)
add    %esi,%ecx
add    0x10(%rsp),%r8
mov    %ecx,%edx
cmp    %r8d,%r13d
jl     1490 <fourn+0x380>
cmp    %r8d,%edx
jle    13f4 <fourn+0x2e4>
lea    (%r11,%r8,1),%r9d
add    $0xfffffffe,%r9d
cmp    %r8d,%r9d
jge    1420 <fourn+0x310>
mov    %r13d,%esi
nopw   0x0(%rax,%rax,1)
mov    %edx,%ecx
sar    %esi
cmp    %r11d,%esi
jl     13d0 <fourn+0x2c0>
mov    %ecx,%edx
sub    %esi,%edx
jg     1400 <fourn+0x2f0>
jmp    13d0 <fourn+0x2c0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %edx,%r10d
mov    %r8,%rdi
jmp    143d <fourn+0x32d>
nopl   0x0(%rax,%rax,1)
add    $0x2,%rdi
add    $0x2,%r10d
cmp    %edi,%r9d
jl     13f4 <fourn+0x2e4>
cmp    %edi,%r14d
jl     1430 <fourn+0x320>
mov    %r10d,%esi
mov    %rdi,%rcx
nopl   0x0(%rax,%rax,1)
movss  (%rbx,%rcx,4),%xmm0
movslq %esi,%rsi
movss  (%rbx,%rsi,4),%xmm1
movss  %xmm1,(%rbx,%rcx,4)
movss  %xmm0,(%rbx,%rsi,4)
movss  0x4(%rbx,%rcx,4),%xmm0
lea    0x1(%rsi),%ebp
movslq %ebp,%rbp
movd   (%rbx,%rbp,4),%xmm1
movd   %xmm1,0x4(%rbx,%rcx,4)
movss  %xmm0,(%rbx,%rbp,4)
add    %rax,%rcx
add    %r13d,%esi
cmp    %r12,%rcx
jle    1450 <fourn+0x340>
jmp    1430 <fourn+0x320>
cmp    %r13d,%r11d
jge    1360 <fourn+0x250>
mov    %r11d,%r9d
mov    %r11,0x30(%rsp)
jmp    14be <fourn+0x3ae>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
cmp    %r13d,%r9d
mov    0x30(%rsp),%r11
jge    1360 <fourn+0x250>
mov    %r9d,%r15d
lea    (%r15,%r15,1),%ebp
mov    %ebp,%eax
cltd
idiv   %r11d
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x28(%rsp),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x40(%rsp)
movapd %xmm1,%xmm0
mulsd  0xb3a(%rip),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,0x38(%rsp)
movsd  0x40(%rsp),%xmm0
call   1030 <sin@plt>
mov    %ebp,%r9d
mov    0x10(%rsp),%rbp
movapd 0xafb(%rip),%xmm9
test   %r15d,%r15d
jle    14b0 <fourn+0x3a0>
movsd  0x38(%rsp),%xmm1
movapd %xmm1,%xmm8
mulsd  0xb02(%rip),%xmm8
mulsd  %xmm1,%xmm8
movslq %r9d,%rax
mov    %r15d,%ecx
unpcklpd %xmm8,%xmm8
unpcklpd %xmm0,%xmm0
lea    (%rbx,%rcx,4),%rcx
add    $0x4,%rcx
mov    $0x1,%edx
movsd  0xae1(%rip),%xmm4
jmp    1585 <fourn+0x475>
nopl   0x0(%rax)
unpcklpd %xmm4,%xmm1
mulpd  %xmm0,%xmm1
movapd %xmm4,%xmm2
mulpd  %xmm8,%xmm2
addpd  %xmm1,%xmm2
addpd  %xmm2,%xmm4
mov    %r8,%rdx
cmp    %r8d,%r15d
jl     14b0 <fourn+0x3a0>
lea    (%rdx,%rbp,1),%r8
lea    -0x2(%r8),%edi
movapd %xmm4,%xmm3
unpckhpd %xmm4,%xmm3
movapd %xmm3,%xmm1
xorpd  %xmm9,%xmm1
jmp    15a4 <fourn+0x494>
add    $0x2,%rdx
cmp    %edx,%edi
jl     1560 <fourn+0x450>
cmp    %edx,%r14d
jl     15a0 <fourn+0x490>
mov    %rdx,%rsi
movss  -0x4(%rcx,%rsi,4),%xmm5
movss  (%rcx,%rsi,4),%xmm6
cvtss2sd %xmm5,%xmm5
cvtss2sd %xmm6,%xmm6
movapd %xmm1,%xmm7
mulsd  %xmm6,%xmm7
movapd %xmm4,%xmm2
mulsd  %xmm5,%xmm2
addsd  %xmm7,%xmm2
cvtsd2ss %xmm2,%xmm2
mulsd  %xmm3,%xmm5
mulsd  %xmm4,%xmm6
addsd  %xmm5,%xmm6
xorps  %xmm5,%xmm5
cvtsd2ss %xmm6,%xmm5
movss  (%rbx,%rsi,4),%xmm6
subss  %xmm2,%xmm6
movss  %xmm6,-0x4(%rcx,%rsi,4)
movss  0x4(%rbx,%rsi,4),%xmm6
movaps %xmm6,%xmm7
subss  %xmm5,%xmm7
movss  %xmm7,(%rcx,%rsi,4)
addss  (%rbx,%rsi,4),%xmm2
movss  %xmm2,(%rbx,%rsi,4)
addss  %xmm5,%xmm6
movss  %xmm6,0x4(%rbx,%rsi,4)
add    %rax,%rsi
cmp    %r12,%rsi
jle    15b0 <fourn+0x4a0>
jmp    15a0 <fourn+0x490>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @fourn(float* nocapture %0, i32* nocapture %1, i32 %2, i32 %3) {
%5 = icmp slt i32 %2, 1
br i1 %5, label %263, label %6
6: ; preds = %4
%7 = add nuw i32 %2, 1
%8 = zext i32 %7 to i64
%9 = add nsw i64 %8, -1
%10 = icmp ult i64 %9, 8
br i1 %10, label %95, label %11
11: ; preds = %6
%12 = and i64 %9, -8
%13 = or i64 %12, 1
%14 = add nsw i64 %12, -8
%15 = lshr exact i64 %14, 3
%16 = add nuw nsw i64 %15, 1
%17 = and i64 %16, 3
%18 = icmp ult i64 %14, 24
br i1 %18, label %65, label %19
19: ; preds = %11
%20 = and i64 %16, 4611686018427387900
br label %21
21: ; preds = %21, %19
%22 = phi i64 [ 0, %19 ], [ %62, %21 ]
%23 = phi <4 x i32> [ <i32 1, i32 1, i32 1, i32 1>, %19 ], [ %60, %21 ]
%24 = phi <4 x i32> [ <i32 1, i32 1, i32 1, i32 1>, %19 ], [ %61, %21 ]
%25 = phi i64 [ 0, %19 ], [ %63, %21 ]
%26 = or i64 %22, 1
%27 = getelementptr inbounds i32, i32* %1, i64 %26
%28 = bitcast i32* %27 to <4 x i32>*
%29 = load <4 x i32>, <4 x i32>* %28
%30 = getelementptr inbounds i32, i32* %27, i64 4
%31 = bitcast i32* %30 to <4 x i32>*
%32 = load <4 x i32>, <4 x i32>* %31
%33 = mul <4 x i32> %29, %23
%34 = mul <4 x i32> %32, %24
%35 = or i64 %22, 9
%36 = getelementptr inbounds i32, i32* %1, i64 %35
%37 = bitcast i32* %36 to <4 x i32>*
%38 = load <4 x i32>, <4 x i32>* %37
%39 = getelementptr inbounds i32, i32* %36, i64 4
%40 = bitcast i32* %39 to <4 x i32>*
%41 = load <4 x i32>, <4 x i32>* %40
%42 = mul <4 x i32> %38, %33
%43 = mul <4 x i32> %41, %34
%44 = or i64 %22, 17
%45 = getelementptr inbounds i32, i32* %1, i64 %44
%46 = bitcast i32* %45 to <4 x i32>*
%47 = load <4 x i32>, <4 x i32>* %46
%48 = getelementptr inbounds i32, i32* %45, i64 4
%49 = bitcast i32* %48 to <4 x i32>*
%50 = load <4 x i32>, <4 x i32>* %49
%51 = mul <4 x i32> %47, %42
%52 = mul <4 x i32> %50, %43
%53 = or i64 %22, 25
%54 = getelementptr inbounds i32, i32* %1, i64 %53
%55 = bitcast i32* %54 to <4 x i32>*
%56 = load <4 x i32>, <4 x i32>* %55
%57 = getelementptr inbounds i32, i32* %54, i64 4
%58 = bitcast i32* %57 to <4 x i32>*
%59 = load <4 x i32>, <4 x i32>* %58
%60 = mul <4 x i32> %56, %51
%61 = mul <4 x i32> %59, %52
%62 = add nuw i64 %22, 32
%63 = add i64 %25, 4
%64 = icmp eq i64 %63, %20
br i1 %64, label %65, label %21
65: ; preds = %21, %11
%66 = phi <4 x i32> [ undef, %11 ], [ %60, %21 ]
%67 = phi <4 x i32> [ undef, %11 ], [ %61, %21 ]
%68 = phi i64 [ 0, %11 ], [ %62, %21 ]
%69 = phi <4 x i32> [ <i32 1, i32 1, i32 1, i32 1>, %11 ], [ %60, %21 ]
%70 = phi <4 x i32> [ <i32 1, i32 1, i32 1, i32 1>, %11 ], [ %61, %21 ]
%71 = icmp eq i64 %17, 0
br i1 %71, label %89, label %72
72: ; preds = %65, %72
%73 = phi i64 [ %86, %72 ], [ %68, %65 ]
%74 = phi <4 x i32> [ %84, %72 ], [ %69, %65 ]
%75 = phi <4 x i32> [ %85, %72 ], [ %70, %65 ]
%76 = phi i64 [ %87, %72 ], [ 0, %65 ]
%77 = or i64 %73, 1
%78 = getelementptr inbounds i32, i32* %1, i64 %77
%79 = bitcast i32* %78 to <4 x i32>*
%80 = load <4 x i32>, <4 x i32>* %79
%81 = getelementptr inbounds i32, i32* %78, i64 4
%82 = bitcast i32* %81 to <4 x i32>*
%83 = load <4 x i32>, <4 x i32>* %82
%84 = mul <4 x i32> %80, %74
%85 = mul <4 x i32> %83, %75
%86 = add nuw i64 %73, 8
%87 = add i64 %76, 1
%88 = icmp eq i64 %87, %17
br i1 %88, label %89, label %72
89: ; preds = %72, %65
%90 = phi <4 x i32> [ %66, %65 ], [ %84, %72 ]
%91 = phi <4 x i32> [ %67, %65 ], [ %85, %72 ]
%92 = mul <4 x i32> %91, %90
%93 = call i32 @llvm.vector.reduce.mul.v4i32(<4 x i32> %92)
%94 = icmp eq i64 %9, %12
br i1 %94, label %98, label %95
95: ; preds = %6, %89
%96 = phi i64 [ 1, %6 ], [ %13, %89 ]
%97 = phi i32 [ 1, %6 ], [ %93, %89 ]
br label %105
98: ; preds = %105, %89
%99 = phi i32 [ %93, %89 ], [ %110, %105 ]
%100 = icmp sgt i32 %2, 0
br i1 %100, label %101, label %263
101: ; preds = %98
%102 = sitofp i32 %3 to double
%103 = fmul double %102, 0x401921FB54442D1C
%104 = zext i32 %2 to i64
br label %113
105: ; preds = %95, %105
%106 = phi i64 [ %111, %105 ], [ %96, %95 ]
%107 = phi i32 [ %110, %105 ], [ %97, %95 ]
%108 = getelementptr inbounds i32, i32* %1, i64 %106
%109 = load i32, i32* %108
%110 = mul nsw i32 %109, %107
%111 = add nuw nsw i64 %106, 1
%112 = icmp eq i64 %111, %8
br i1 %112, label %98, label %105
113: ; preds = %101, %260
%114 = phi i64 [ %104, %101 ], [ %261, %260 ]
%115 = phi i32 [ 1, %101 ], [ %118, %260 ]
%116 = getelementptr inbounds i32, i32* %1, i64 %114
%117 = load i32, i32* %116
%118 = mul nsw i32 %117, %115
%119 = sdiv i32 %99, %118
%120 = shl i32 %115, 1
%121 = mul i32 %117, %120
%122 = mul nsw i32 %121, %119
%123 = icmp slt i32 %121, 1
br i1 %123, label %128, label %124
124: ; preds = %113
%125 = sext i32 %120 to i64
%126 = zext i32 %121 to i64
%127 = sext i32 %122 to i64
br label %133
128: ; preds = %179, %113
%129 = icmp slt i32 %120, %121
br i1 %129, label %130, label %260
130: ; preds = %128
%131 = sext i32 %120 to i64
%132 = sext i32 %122 to i64
br label %186
133: ; preds = %124, %179
%134 = phi i64 [ 1, %124 ], [ %181, %179 ]
%135 = phi i32 [ 1, %124 ], [ %180, %179 ]
%136 = trunc i64 %134 to i32
%137 = icmp sgt i32 %135, %136
br i1 %137, label %138, label %170
138: ; preds = %133
%139 = add nsw i32 %120, %136
%140 = add nsw i32 %139, -2
%141 = icmp slt i32 %140, %136
br i1 %141, label %170, label %142
142: ; preds = %138
%143 = sub i32 %135, %136
br label %144
144: ; preds = %142, %166
%145 = phi i64 [ %134, %142 ], [ %167, %166 ]
%146 = trunc i64 %145 to i32
%147 = icmp slt i32 %122, %146
br i1 %147, label %166, label %148
148: ; preds = %144, %148
%149 = phi i64 [ %164, %148 ], [ %145, %144 ]
%150 = trunc i64 %149 to i32
%151 = add i32 %143, %150
%152 = getelementptr inbounds float, float* %0, i64 %149
%153 = load float, float* %152
%154 = sext i32 %151 to i64
%155 = getelementptr inbounds float, float* %0, i64 %154
%156 = load float, float* %155
store float %156, float* %152
store float %153, float* %155
%157 = add nsw i64 %149, 1
%158 = getelementptr inbounds float, float* %0, i64 %157
%159 = load float, float* %158
%160 = add nsw i32 %151, 1
%161 = sext i32 %160 to i64
%162 = getelementptr inbounds float, float* %0, i64 %161
%163 = load float, float* %162
store float %163, float* %158
store float %159, float* %162
%164 = add i64 %149, %126
%165 = icmp sgt i64 %164, %127
br i1 %165, label %166, label %148
166: ; preds = %148, %144
%167 = add nsw i64 %145, 2
%168 = trunc i64 %167 to i32
%169 = icmp slt i32 %140, %168
br i1 %169, label %170, label %144
170: ; preds = %166, %138, %133
br label %171
171: ; preds = %170, %171
%172 = phi i32 [ %174, %171 ], [ %121, %170 ]
%173 = phi i32 [ %178, %171 ], [ %135, %170 ]
%174 = ashr i32 %172, 1
%175 = icmp sge i32 %174, %120
%176 = icmp sgt i32 %173, %174
%177 = select i1 %175, i1 %176, i1 false
%178 = sub nsw i32 %173, %174
br i1 %177, label %171, label %179
179: ; preds = %171
%180 = add nsw i32 %173, %174
%181 = add nsw i64 %134, %125
%182 = trunc i64 %181 to i32
%183 = icmp slt i32 %121, %182
br i1 %183, label %128, label %133
184: ; preds = %253, %186
%185 = icmp slt i32 %188, %121
br i1 %185, label %186, label %260
186: ; preds = %130, %184
%187 = phi i32 [ %188, %184 ], [ %120, %130 ]
%188 = shl i32 %187, 1
%189 = sdiv i32 %188, %120
%190 = sitofp i32 %189 to double
%191 = fdiv double %103, %190
%192 = fmul double %191, 5.000000e-01
%193 = tail call double @sin(double %192)
%194 = tail call double @sin(double %191)
%195 = icmp slt i32 %187, 1
br i1 %195, label %184, label %196
196: ; preds = %186
%197 = fmul double %193, -2.000000e+00
%198 = fmul double %193, %197
%199 = sext i32 %188 to i64
%200 = zext i32 %187 to i64
%201 = insertelement <2 x double> poison, double %198, i64 0
%202 = shufflevector <2 x double> %201, <2 x double> poison, <2 x i32> zeroinitializer
%203 = insertelement <2 x double> poison, double %194, i64 0
%204 = shufflevector <2 x double> %203, <2 x double> poison, <2 x i32> zeroinitializer
br label %205
205: ; preds = %196, %253
%206 = phi i64 [ 1, %196 ], [ %209, %253 ]
%207 = phi <2 x double> [ <double 1.000000e+00, double 0.000000e+00>, %196 ], [ %258, %253 ]
%208 = trunc i64 %206 to i32
%209 = add nsw i64 %206, %131
%210 = trunc i64 %209 to i32
%211 = add nsw i32 %210, -2
%212 = icmp slt i32 %211, %208
%213 = extractelement <2 x double> %207, i64 1
%214 = fneg double %213
br i1 %212, label %253, label %215
215: ; preds = %205
%216 = extractelement <2 x double> %207, i64 0
br label %217
217: ; preds = %215, %249
%218 = phi i64 [ %250, %249 ], [ %206, %215 ]
%219 = trunc i64 %218 to i32
%220 = icmp slt i32 %122, %219
br i1 %220, label %249, label %221
221: ; preds = %217, %221
%222 = phi i64 [ %247, %221 ], [ %218, %217 ]
%223 = add nsw i64 %222, %200
%224 = getelementptr inbounds float, float* %0, i64 %223
%225 = load float, float* %224
%226 = fpext float %225 to double
%227 = add nsw i64 %223, 1
%228 = getelementptr inbounds float, float* %0, i64 %227
%229 = load float, float* %228
%230 = fpext float %229 to double
%231 = fmul double %214, %230
%232 = tail call double @llvm.fmuladd.f64(double %216, double %226, double %231)
%233 = fptrunc double %232 to float
%234 = fmul double %213, %226
%235 = tail call double @llvm.fmuladd.f64(double %216, double %230, double %234)
%236 = fptrunc double %235 to float
%237 = getelementptr inbounds float, float* %0, i64 %222
%238 = load float, float* %237
%239 = fsub float %238, %233
store float %239, float* %224
%240 = add nsw i64 %222, 1
%241 = getelementptr inbounds float, float* %0, i64 %240
%242 = load float, float* %241
%243 = fsub float %242, %236
store float %243, float* %228
%244 = load float, float* %237
%245 = fadd float %244, %233
store float %245, float* %237
%246 = fadd float %242, %236
store float %246, float* %241
%247 = add i64 %222, %199
%248 = icmp sgt i64 %247, %132
br i1 %248, label %249, label %221
249: ; preds = %221, %217
%250 = add nsw i64 %218, 2
%251 = trunc i64 %250 to i32
%252 = icmp slt i32 %211, %251
br i1 %252, label %253, label %217
253: ; preds = %249, %205
%254 = insertelement <2 x double> poison, double %214, i64 0
%255 = shufflevector <2 x double> %254, <2 x double> %207, <2 x i32> <i32 0, i32 2>
%256 = fmul <2 x double> %204, %255
%257 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %207, <2 x double> %202, <2 x double> %256)
%258 = fadd <2 x double> %207, %257
%259 = icmp slt i32 %187, %210
br i1 %259, label %184, label %205
260: ; preds = %184, %128
%261 = add nsw i64 %114, -1
%262 = icmp sgt i64 %114, 1
br i1 %262, label %113, label %263
263: ; preds = %260, %4, %98
ret void
}