void cftmdl1(int n, double *a, double *w) {
  int j, j0, j1, j2, j3, k, m, mh;
  double wn4r, wk1r, wk1i, wk3r, wk3i;
  double x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;
  mh = n >> 3;
  m = 2 * mh;
  j1 = m;
  j2 = j1 + m;
  j3 = j2 + m;
  x0r = a[0] + a[j2];
  x0i = a[1] + a[j2 + 1];
  x1r = a[0] - a[j2];
  x1i = a[1] - a[j2 + 1];
  x2r = a[j1] + a[j3];
  x2i = a[j1 + 1] + a[j3 + 1];
  x3r = a[j1] - a[j3];
  x3i = a[j1 + 1] - a[j3 + 1];
  a[0] = x0r + x2r;
  a[1] = x0i + x2i;
  a[j1] = x0r - x2r;
  a[j1 + 1] = x0i - x2i;
  a[j2] = x1r - x3i;
  a[j2 + 1] = x1i + x3r;
  a[j3] = x1r + x3i;
  a[j3 + 1] = x1i - x3r;
  wn4r = w[1];
  k = 0;
  for (j = 2; j < mh; j += 2) {
    k += 4;
    wk1r = w[k];
    wk1i = w[k + 1];
    wk3r = w[k + 2];
    wk3i = -w[k + 3];
    j1 = j + m;
    j2 = j1 + m;
    j3 = j2 + m;
    x0r = a[j] + a[j2];
    x0i = a[j + 1] + a[j2 + 1];
    x1r = a[j] - a[j2];
    x1i = a[j + 1] - a[j2 + 1];
    x2r = a[j1] + a[j3];
    x2i = a[j1 + 1] + a[j3 + 1];
    x3r = a[j1] - a[j3];
    x3i = a[j1 + 1] - a[j3 + 1];
    a[j] = x0r + x2r;
    a[j + 1] = x0i + x2i;
    a[j1] = x0r - x2r;
    a[j1 + 1] = x0i - x2i;
    x0r = x1r - x3i;
    x0i = x1i + x3r;
    a[j2] = wk1r * x0r - wk1i * x0i;
    a[j2 + 1] = wk1r * x0i + wk1i * x0r;
    x0r = x1r + x3i;
    x0i = x1i - x3r;
    a[j3] = wk3r * x0r + wk3i * x0i;
    a[j3 + 1] = wk3r * x0i - wk3i * x0r;
    j0 = m - j;
    j1 = j0 + m;
    j2 = j1 + m;
    j3 = j2 + m;
    x0r = a[j0] + a[j2];
    x0i = a[j0 + 1] + a[j2 + 1];
    x1r = a[j0] - a[j2];
    x1i = a[j0 + 1] - a[j2 + 1];
    x2r = a[j1] + a[j3];
    x2i = a[j1 + 1] + a[j3 + 1];
    x3r = a[j1] - a[j3];
    x3i = a[j1 + 1] - a[j3 + 1];
    a[j0] = x0r + x2r;
    a[j0 + 1] = x0i + x2i;
    a[j1] = x0r - x2r;
    a[j1 + 1] = x0i - x2i;
    x0r = x1r - x3i;
    x0i = x1i + x3r;
    a[j2] = wk1i * x0r - wk1r * x0i;
    a[j2 + 1] = wk1i * x0i + wk1r * x0r;
    x0r = x1r + x3i;
    x0i = x1i - x3r;
    a[j3] = wk3i * x0r + wk3r * x0i;
    a[j3 + 1] = wk3i * x0i - wk3r * x0r;
  }
  j0 = mh;
  j1 = j0 + m;
  j2 = j1 + m;
  j3 = j2 + m;
  x0r = a[j0] + a[j2];
  x0i = a[j0 + 1] + a[j2 + 1];
  x1r = a[j0] - a[j2];
  x1i = a[j0 + 1] - a[j2 + 1];
  x2r = a[j1] + a[j3];
  x2i = a[j1 + 1] + a[j3 + 1];
  x3r = a[j1] - a[j3];
  x3i = a[j1 + 1] - a[j3 + 1];
  a[j0] = x0r + x2r;
  a[j0 + 1] = x0i + x2i;
  a[j1] = x0r - x2r;
  a[j1 + 1] = x0i - x2i;
  x0r = x1r - x3i;
  x0i = x1i + x3r;
  a[j2] = wn4r * (x0r - x0i);
  a[j2 + 1] = wn4r * (x0i + x0r);
  x0r = x1r + x3i;
  x0i = x1i - x3r;
  a[j3] = -wn4r * (x0r + x0i);
  a[j3 + 1] = -wn4r * (x0i - x0r);
}


x86 assembly (clang, -O0):
<cftmdl1>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x4(%rbp),%eax
sar    $0x3,%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x70(%rbp),%xmm0
addsd  -0x90(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0x68(%rbp),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
subsd  -0x90(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x78(%rbp),%xmm0
subsd  -0xa0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x80(%rbp),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x78(%rbp),%xmm0
addsd  -0xa0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x80(%rbp),%xmm0
subsd  -0x98(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x2,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    1993 <cftmdl1+0x893>
mov    -0x30(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x1c(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
addsd  -0x90(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x68(%rbp),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
subsd  -0x90(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x78(%rbp),%xmm0
subsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0x68(%rbp)
movsd  -0x80(%rbp),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x70(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x70(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x78(%rbp),%xmm0
addsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0x68(%rbp)
movsd  -0x80(%rbp),%xmm0
subsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x70(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm2
movsd  -0x60(%rbp),%xmm1
mulsd  -0x70(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm2
movsd  -0x60(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
sub    -0x1c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
addsd  -0x90(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x68(%rbp),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
subsd  -0x90(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x78(%rbp),%xmm0
subsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0x68(%rbp)
movsd  -0x80(%rbp),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x70(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  -0x70(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x78(%rbp),%xmm0
addsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0x68(%rbp)
movsd  -0x80(%rbp),%xmm0
subsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x70(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  -0x70(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x1c(%rbp)
jmp    133f <cftmdl1+0x23f>
mov    -0x38(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
addsd  -0x90(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x68(%rbp),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
subsd  -0x90(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x78(%rbp),%xmm0
subsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0x68(%rbp)
movsd  -0x80(%rbp),%xmm0
addsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x70(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm1
subsd  -0x70(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm1
addsd  -0x68(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x78(%rbp),%xmm0
addsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0x68(%rbp)
movsd  -0x80(%rbp),%xmm0
subsd  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x70(%rbp)
movsd  -0x40(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x68(%rbp),%xmm1
addsd  -0x70(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x70(%rbp),%xmm1
subsd  -0x68(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cftmdl1>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %edi,%eax
sar    $0x3,%eax
lea    (%rax,%rax,1),%ecx
lea    0x0(,%rax,4),%ebp
lea    (%rcx,%rcx,2),%r8d
movslq %ebp,%rbx
movupd (%rsi),%xmm0
movupd (%rsi,%rbx,8),%xmm1
movapd %xmm0,%xmm2
addpd  %xmm1,%xmm2
subpd  %xmm1,%xmm0
movslq %ecx,%rbp
movslq %r8d,%rcx
movupd (%rsi,%rbp,8),%xmm1
movupd (%rsi,%rcx,8),%xmm3
movapd %xmm1,%xmm4
addpd  %xmm3,%xmm4
subpd  %xmm3,%xmm1
shufpd $0x1,%xmm1,%xmm1
movapd %xmm2,%xmm3
addpd  %xmm4,%xmm3
movupd %xmm3,(%rsi)
subpd  %xmm4,%xmm2
movupd %xmm2,(%rsi,%rbp,8)
movapd %xmm0,%xmm2
subpd  %xmm1,%xmm2
addpd  %xmm0,%xmm1
movapd %xmm1,%xmm0
movsd  %xmm2,%xmm0
movupd %xmm0,(%rsi,%rbx,8)
shufpd $0x2,%xmm2,%xmm1
movupd %xmm1,(%rsi,%rcx,8)
movsd  0x8(%rdx),%xmm8
movslq %eax,%rcx
cmp    $0x18,%edi
jl     14a4 <cftmdl1+0x3a4>
mov    %rbx,%r8
movabs $0xfffffffe00000000,%r9
ff ff ff
add    $0x38,%rdx
lea    0x0(%rbp,%rbp,2),%eax
movslq %eax,%rdi
lea    0x3(%rdi),%rax
mov    %rax,-0x20(%rsp)
lea    0x3(%r8),%r12
add    $0xffffffffffffffff,%rbx
mov    %rbx,-0x18(%rsp)
add    $0xffffffffffffffff,%rdi
lea    -0x2(%rbp),%rax
mov    %rbp,%r11
shl    $0x4,%r11
lea    -0x10(%r11),%r10
lea    0x0(,%rbp,8),%r8
lea    (%r8,%r8,2),%r15
add    $0xfffffffffffffff0,%r15
mov    %rbp,%rbx
shl    $0x22,%rbx
add    %r9,%rbx
mov    %rbp,%r14
shl    $0x5,%r14
add    $0xfffffffffffffff0,%r14
add    $0x2,%rbp
add    $0x10,%r11
lea    (%r8,%r8,2),%r9
add    $0x10,%r9
mov    $0x2,%r8d
movapd 0xddb(%rip),%xmm9
mov    %rcx,-0x10(%rsp)
nopw   0x0(%rax,%rax,1)
mov    %rbp,-0x8(%rsp)
movsd  -0x18(%rdx),%xmm12
movsd  -0x10(%rdx),%xmm13
movsd  -0x8(%rdx),%xmm11
movsd  (%rdx),%xmm14
movapd %xmm14,%xmm10
movsd  (%rsi,%r11,1),%xmm4
movupd (%rsi,%r8,8),%xmm7
movsd  (%rsi,%r12,8),%xmm5
xorpd  %xmm9,%xmm10
movapd %xmm7,%xmm1
subsd  %xmm4,%xmm1
unpcklpd %xmm5,%xmm4
movsd  (%rsi,%r9,1),%xmm6
addpd  %xmm7,%xmm4
movupd (%rsi,%rbp,8),%xmm0
mov    -0x20(%rsp),%rcx
movsd  (%rsi,%rcx,8),%xmm3
movapd %xmm0,%xmm2
subsd  %xmm6,%xmm2
unpckhpd %xmm7,%xmm7
unpcklpd %xmm3,%xmm6
addpd  %xmm0,%xmm6
unpckhpd %xmm0,%xmm0
subsd  %xmm5,%xmm7
subsd  %xmm3,%xmm0
movapd %xmm4,%xmm3
addpd  %xmm6,%xmm3
subpd  %xmm6,%xmm4
movupd %xmm3,(%rsi,%r8,8)
movupd %xmm4,(%rsi,%rbp,8)
movapd %xmm1,%xmm3
subsd  %xmm0,%xmm3
movapd %xmm7,%xmm4
addsd  %xmm2,%xmm4
movapd %xmm12,%xmm5
mulsd  %xmm4,%xmm5
mulsd  %xmm13,%xmm4
movapd %xmm12,%xmm6
mulsd  %xmm3,%xmm6
mulsd  %xmm13,%xmm3
subsd  %xmm4,%xmm6
addsd  %xmm5,%xmm3
addsd  %xmm1,%xmm0
subsd  %xmm2,%xmm7
movapd %xmm11,%xmm1
movsd  %xmm6,(%rsi,%r11,1)
mulsd  %xmm7,%xmm1
mulsd  %xmm14,%xmm7
movapd %xmm11,%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm3,(%rsi,%r12,8)
subsd  %xmm7,%xmm2
mulsd  %xmm14,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm2,(%rsi,%r9,1)
movsd  %xmm0,(%rsi,%rcx,8)
movsd  (%rsi,%r15,1),%xmm4
movupd (%rsi,%rax,8),%xmm7
movsd  (%rsi,%rdi,8),%xmm3
movapd %xmm7,%xmm15
subsd  %xmm4,%xmm15
unpcklpd %xmm3,%xmm4
movsd  (%rsi,%r10,1),%xmm1
addpd  %xmm7,%xmm4
movsd  (%rsi,%r14,1),%xmm5
mov    %r11,%rbp
mov    %r15,%r11
mov    %r10,%r15
mov    -0x18(%rsp),%r10
movsd  (%rsi,%r10,8),%xmm0
mov    %rbx,%r13
sar    $0x1d,%r13
unpckhpd %xmm7,%xmm7
or     $0x8,%r13
movsd  (%rsi,%r13,1),%xmm6
movapd %xmm1,%xmm2
unpcklpd %xmm0,%xmm2
subsd  %xmm3,%xmm7
subsd  %xmm5,%xmm1
unpcklpd %xmm6,%xmm5
subsd  %xmm6,%xmm0
addpd  %xmm2,%xmm5
movapd %xmm4,%xmm2
addpd  %xmm5,%xmm2
movupd %xmm2,(%rsi,%rax,8)
subpd  %xmm5,%xmm4
movlpd %xmm4,(%rsi,%r15,1)
movhpd %xmm4,(%rsi,%r10,8)
movapd %xmm15,%xmm2
subsd  %xmm0,%xmm2
movapd %xmm7,%xmm3
addsd  %xmm1,%xmm3
movapd %xmm3,%xmm4
mulsd  %xmm12,%xmm4
mulsd  %xmm13,%xmm3
mulsd  %xmm2,%xmm13
mulsd  %xmm12,%xmm2
subsd  %xmm4,%xmm13
addsd  %xmm2,%xmm3
addsd  %xmm15,%xmm0
subsd  %xmm1,%xmm7
movsd  %xmm13,(%rsi,%r11,1)
mulsd  %xmm0,%xmm14
mulsd  %xmm11,%xmm0
mulsd  %xmm7,%xmm11
movsd  %xmm3,(%rsi,%rdi,8)
subsd  %xmm14,%xmm11
mulsd  %xmm10,%xmm7
subsd  %xmm0,%xmm7
movsd  %xmm11,(%rsi,%r14,1)
movsd  %xmm7,(%rsi,%r13,1)
mov    -0x10(%rsp),%r13
add    $0x2,%r8
add    $0x20,%rdx
add    $0x2,%rcx
mov    %rcx,-0x20(%rsp)
add    $0x2,%r12
add    $0xfffffffffffffffe,%r10
mov    %r10,-0x18(%rsp)
mov    %r15,%r10
mov    %r11,%r15
mov    %rbp,%r11
mov    -0x8(%rsp),%rbp
add    $0xfffffffffffffffe,%rax
add    $0xfffffffffffffffe,%rdi
add    $0xfffffffffffffff0,%r10
add    $0xfffffffffffffff0,%r15
movabs $0xfffffffe00000000,%rcx
ff ff ff
add    %rcx,%rbx
add    $0xfffffffffffffff0,%r14
add    $0x2,%rbp
add    $0x10,%r11
add    $0x10,%r9
mov    %r13,%rcx
cmp    %r13,%r8
jl     1230 <cftmdl1+0x130>
mov    %rcx,%r13
lea    (%rcx,%rcx,2),%ecx
lea    0x0(%r13,%r13,4),%eax
lea    0x0(,%r13,8),%edi
sub    %r13d,%edi
cltq
movupd (%rsi,%r13,8),%xmm1
movupd (%rsi,%rax,8),%xmm2
movapd %xmm1,%xmm0
addpd  %xmm2,%xmm0
subpd  %xmm2,%xmm1
movslq %ecx,%rdx
movslq %edi,%rcx
movupd (%rsi,%rdx,8),%xmm2
movupd (%rsi,%rcx,8),%xmm4
movapd %xmm2,%xmm5
subpd  %xmm4,%xmm2
movapd %xmm2,%xmm3
unpckhpd %xmm2,%xmm3
movapd %xmm1,%xmm6
subsd  %xmm3,%xmm6
addsd  %xmm1,%xmm3
unpckhpd %xmm1,%xmm1
addpd  %xmm4,%xmm5
movapd %xmm0,%xmm4
addpd  %xmm5,%xmm4
movupd %xmm4,(%rsi,%r13,8)
subpd  %xmm5,%xmm0
movupd %xmm0,(%rsi,%rdx,8)
movapd %xmm1,%xmm0
addsd  %xmm2,%xmm0
movapd %xmm6,%xmm4
subsd  %xmm0,%xmm4
mulsd  %xmm8,%xmm4
movsd  %xmm4,(%rsi,%rax,8)
addsd  %xmm6,%xmm0
mulsd  %xmm8,%xmm0
movsd  %xmm0,0x8(%rsi,%rax,8)
subsd  %xmm2,%xmm1
xorpd  0xaac(%rip),%xmm8
movapd %xmm1,%xmm0
addsd  %xmm3,%xmm0
mulsd  %xmm8,%xmm0
movsd  %xmm0,(%rsi,%rcx,8)
subsd  %xmm3,%xmm1
mulsd  %xmm8,%xmm1
movsd  %xmm1,0x8(%rsi,%rcx,8)
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @cftmdl1(i32 %0, double* %1, double* %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca double*
%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 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
store i32 %0, i32* %4
store double* %1, double** %5
store double* %2, double** %6
%28 = load i32, i32* %4
%29 = ashr i32 %28, 3
store i32 %29, i32* %14
%30 = load i32, i32* %14
%31 = mul nsw i32 2, %30
store i32 %31, i32* %13
%32 = load i32, i32* %13
store i32 %32, i32* %9
%33 = load i32, i32* %9
%34 = load i32, i32* %13
%35 = add nsw i32 %33, %34
store i32 %35, i32* %10
%36 = load i32, i32* %10
%37 = load i32, i32* %13
%38 = add nsw i32 %36, %37
store i32 %38, i32* %11
%39 = load double*, double** %5
%40 = getelementptr inbounds double, double* %39, i64 0
%41 = load double, double* %40
%42 = load double*, double** %5
%43 = load i32, i32* %10
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds double, double* %42, i64 %44
%46 = load double, double* %45
%47 = fadd double %41, %46
store double %47, double* %20
%48 = load double*, double** %5
%49 = getelementptr inbounds double, double* %48, i64 1
%50 = load double, double* %49
%51 = load double*, double** %5
%52 = load i32, i32* %10
%53 = add nsw i32 %52, 1
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds double, double* %51, i64 %54
%56 = load double, double* %55
%57 = fadd double %50, %56
store double %57, double* %21
%58 = load double*, double** %5
%59 = getelementptr inbounds double, double* %58, i64 0
%60 = load double, double* %59
%61 = load double*, double** %5
%62 = load i32, i32* %10
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds double, double* %61, i64 %63
%65 = load double, double* %64
%66 = fsub double %60, %65
store double %66, double* %22
%67 = load double*, double** %5
%68 = getelementptr inbounds double, double* %67, i64 1
%69 = load double, double* %68
%70 = load double*, double** %5
%71 = load i32, i32* %10
%72 = add nsw i32 %71, 1
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %70, i64 %73
%75 = load double, double* %74
%76 = fsub double %69, %75
store double %76, double* %23
%77 = load double*, double** %5
%78 = load i32, i32* %9
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds double, double* %77, i64 %79
%81 = load double, double* %80
%82 = load double*, double** %5
%83 = load i32, i32* %11
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %82, i64 %84
%86 = load double, double* %85
%87 = fadd double %81, %86
store double %87, double* %24
%88 = load double*, double** %5
%89 = load i32, i32* %9
%90 = add nsw i32 %89, 1
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double, double* %88, i64 %91
%93 = load double, double* %92
%94 = load double*, double** %5
%95 = load i32, i32* %11
%96 = add nsw i32 %95, 1
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds double, double* %94, i64 %97
%99 = load double, double* %98
%100 = fadd double %93, %99
store double %100, double* %25
%101 = load double*, double** %5
%102 = load i32, i32* %9
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds double, double* %101, i64 %103
%105 = load double, double* %104
%106 = load double*, double** %5
%107 = load i32, i32* %11
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds double, double* %106, i64 %108
%110 = load double, double* %109
%111 = fsub double %105, %110
store double %111, double* %26
%112 = load double*, double** %5
%113 = load i32, i32* %9
%114 = add nsw i32 %113, 1
%115 = sext i32 %114 to i64
%116 = getelementptr inbounds double, double* %112, i64 %115
%117 = load double, double* %116
%118 = load double*, double** %5
%119 = load i32, i32* %11
%120 = add nsw i32 %119, 1
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds double, double* %118, i64 %121
%123 = load double, double* %122
%124 = fsub double %117, %123
store double %124, double* %27
%125 = load double, double* %20
%126 = load double, double* %24
%127 = fadd double %125, %126
%128 = load double*, double** %5
%129 = getelementptr inbounds double, double* %128, i64 0
store double %127, double* %129
%130 = load double, double* %21
%131 = load double, double* %25
%132 = fadd double %130, %131
%133 = load double*, double** %5
%134 = getelementptr inbounds double, double* %133, i64 1
store double %132, double* %134
%135 = load double, double* %20
%136 = load double, double* %24
%137 = fsub double %135, %136
%138 = load double*, double** %5
%139 = load i32, i32* %9
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds double, double* %138, i64 %140
store double %137, double* %141
%142 = load double, double* %21
%143 = load double, double* %25
%144 = fsub double %142, %143
%145 = load double*, double** %5
%146 = load i32, i32* %9
%147 = add nsw i32 %146, 1
%148 = sext i32 %147 to i64
%149 = getelementptr inbounds double, double* %145, i64 %148
store double %144, double* %149
%150 = load double, double* %22
%151 = load double, double* %27
%152 = fsub double %150, %151
%153 = load double*, double** %5
%154 = load i32, i32* %10
%155 = sext i32 %154 to i64
%156 = getelementptr inbounds double, double* %153, i64 %155
store double %152, double* %156
%157 = load double, double* %23
%158 = load double, double* %26
%159 = fadd double %157, %158
%160 = load double*, double** %5
%161 = load i32, i32* %10
%162 = add nsw i32 %161, 1
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds double, double* %160, i64 %163
store double %159, double* %164
%165 = load double, double* %22
%166 = load double, double* %27
%167 = fadd double %165, %166
%168 = load double*, double** %5
%169 = load i32, i32* %11
%170 = sext i32 %169 to i64
%171 = getelementptr inbounds double, double* %168, i64 %170
store double %167, double* %171
%172 = load double, double* %23
%173 = load double, double* %26
%174 = fsub double %172, %173
%175 = load double*, double** %5
%176 = load i32, i32* %11
%177 = add nsw i32 %176, 1
%178 = sext i32 %177 to i64
%179 = getelementptr inbounds double, double* %175, i64 %178
store double %174, double* %179
%180 = load double*, double** %6
%181 = getelementptr inbounds double, double* %180, i64 1
%182 = load double, double* %181
store double %182, double* %15
store i32 0, i32* %12
store i32 2, i32* %7
br label %183
183: ; preds = %599, %3
%184 = load i32, i32* %7
%185 = load i32, i32* %14
%186 = icmp slt i32 %184, %185
br i1 %186, label %187, label %602
187: ; preds = %183
%188 = load i32, i32* %12
%189 = add nsw i32 %188, 4
store i32 %189, i32* %12
%190 = load double*, double** %6
%191 = load i32, i32* %12
%192 = sext i32 %191 to i64
%193 = getelementptr inbounds double, double* %190, i64 %192
%194 = load double, double* %193
store double %194, double* %16
%195 = load double*, double** %6
%196 = load i32, i32* %12
%197 = add nsw i32 %196, 1
%198 = sext i32 %197 to i64
%199 = getelementptr inbounds double, double* %195, i64 %198
%200 = load double, double* %199
store double %200, double* %17
%201 = load double*, double** %6
%202 = load i32, i32* %12
%203 = add nsw i32 %202, 2
%204 = sext i32 %203 to i64
%205 = getelementptr inbounds double, double* %201, i64 %204
%206 = load double, double* %205
store double %206, double* %18
%207 = load double*, double** %6
%208 = load i32, i32* %12
%209 = add nsw i32 %208, 3
%210 = sext i32 %209 to i64
%211 = getelementptr inbounds double, double* %207, i64 %210
%212 = load double, double* %211
%213 = fneg double %212
store double %213, double* %19
%214 = load i32, i32* %7
%215 = load i32, i32* %13
%216 = add nsw i32 %214, %215
store i32 %216, i32* %9
%217 = load i32, i32* %9
%218 = load i32, i32* %13
%219 = add nsw i32 %217, %218
store i32 %219, i32* %10
%220 = load i32, i32* %10
%221 = load i32, i32* %13
%222 = add nsw i32 %220, %221
store i32 %222, i32* %11
%223 = load double*, double** %5
%224 = load i32, i32* %7
%225 = sext i32 %224 to i64
%226 = getelementptr inbounds double, double* %223, i64 %225
%227 = load double, double* %226
%228 = load double*, double** %5
%229 = load i32, i32* %10
%230 = sext i32 %229 to i64
%231 = getelementptr inbounds double, double* %228, i64 %230
%232 = load double, double* %231
%233 = fadd double %227, %232
store double %233, double* %20
%234 = load double*, double** %5
%235 = load i32, i32* %7
%236 = add nsw i32 %235, 1
%237 = sext i32 %236 to i64
%238 = getelementptr inbounds double, double* %234, i64 %237
%239 = load double, double* %238
%240 = load double*, double** %5
%241 = load i32, i32* %10
%242 = add nsw i32 %241, 1
%243 = sext i32 %242 to i64
%244 = getelementptr inbounds double, double* %240, i64 %243
%245 = load double, double* %244
%246 = fadd double %239, %245
store double %246, double* %21
%247 = load double*, double** %5
%248 = load i32, i32* %7
%249 = sext i32 %248 to i64
%250 = getelementptr inbounds double, double* %247, i64 %249
%251 = load double, double* %250
%252 = load double*, double** %5
%253 = load i32, i32* %10
%254 = sext i32 %253 to i64
%255 = getelementptr inbounds double, double* %252, i64 %254
%256 = load double, double* %255
%257 = fsub double %251, %256
store double %257, double* %22
%258 = load double*, double** %5
%259 = load i32, i32* %7
%260 = add nsw i32 %259, 1
%261 = sext i32 %260 to i64
%262 = getelementptr inbounds double, double* %258, i64 %261
%263 = load double, double* %262
%264 = load double*, double** %5
%265 = load i32, i32* %10
%266 = add nsw i32 %265, 1
%267 = sext i32 %266 to i64
%268 = getelementptr inbounds double, double* %264, i64 %267
%269 = load double, double* %268
%270 = fsub double %263, %269
store double %270, double* %23
%271 = load double*, double** %5
%272 = load i32, i32* %9
%273 = sext i32 %272 to i64
%274 = getelementptr inbounds double, double* %271, i64 %273
%275 = load double, double* %274
%276 = load double*, double** %5
%277 = load i32, i32* %11
%278 = sext i32 %277 to i64
%279 = getelementptr inbounds double, double* %276, i64 %278
%280 = load double, double* %279
%281 = fadd double %275, %280
store double %281, double* %24
%282 = load double*, double** %5
%283 = load i32, i32* %9
%284 = add nsw i32 %283, 1
%285 = sext i32 %284 to i64
%286 = getelementptr inbounds double, double* %282, i64 %285
%287 = load double, double* %286
%288 = load double*, double** %5
%289 = load i32, i32* %11
%290 = add nsw i32 %289, 1
%291 = sext i32 %290 to i64
%292 = getelementptr inbounds double, double* %288, i64 %291
%293 = load double, double* %292
%294 = fadd double %287, %293
store double %294, double* %25
%295 = load double*, double** %5
%296 = load i32, i32* %9
%297 = sext i32 %296 to i64
%298 = getelementptr inbounds double, double* %295, i64 %297
%299 = load double, double* %298
%300 = load double*, double** %5
%301 = load i32, i32* %11
%302 = sext i32 %301 to i64
%303 = getelementptr inbounds double, double* %300, i64 %302
%304 = load double, double* %303
%305 = fsub double %299, %304
store double %305, double* %26
%306 = load double*, double** %5
%307 = load i32, i32* %9
%308 = add nsw i32 %307, 1
%309 = sext i32 %308 to i64
%310 = getelementptr inbounds double, double* %306, i64 %309
%311 = load double, double* %310
%312 = load double*, double** %5
%313 = load i32, i32* %11
%314 = add nsw i32 %313, 1
%315 = sext i32 %314 to i64
%316 = getelementptr inbounds double, double* %312, i64 %315
%317 = load double, double* %316
%318 = fsub double %311, %317
store double %318, double* %27
%319 = load double, double* %20
%320 = load double, double* %24
%321 = fadd double %319, %320
%322 = load double*, double** %5
%323 = load i32, i32* %7
%324 = sext i32 %323 to i64
%325 = getelementptr inbounds double, double* %322, i64 %324
store double %321, double* %325
%326 = load double, double* %21
%327 = load double, double* %25
%328 = fadd double %326, %327
%329 = load double*, double** %5
%330 = load i32, i32* %7
%331 = add nsw i32 %330, 1
%332 = sext i32 %331 to i64
%333 = getelementptr inbounds double, double* %329, i64 %332
store double %328, double* %333
%334 = load double, double* %20
%335 = load double, double* %24
%336 = fsub double %334, %335
%337 = load double*, double** %5
%338 = load i32, i32* %9
%339 = sext i32 %338 to i64
%340 = getelementptr inbounds double, double* %337, i64 %339
store double %336, double* %340
%341 = load double, double* %21
%342 = load double, double* %25
%343 = fsub double %341, %342
%344 = load double*, double** %5
%345 = load i32, i32* %9
%346 = add nsw i32 %345, 1
%347 = sext i32 %346 to i64
%348 = getelementptr inbounds double, double* %344, i64 %347
store double %343, double* %348
%349 = load double, double* %22
%350 = load double, double* %27
%351 = fsub double %349, %350
store double %351, double* %20
%352 = load double, double* %23
%353 = load double, double* %26
%354 = fadd double %352, %353
store double %354, double* %21
%355 = load double, double* %16
%356 = load double, double* %20
%357 = load double, double* %17
%358 = load double, double* %21
%359 = fmul double %357, %358
%360 = fneg double %359
%361 = call double @llvm.fmuladd.f64(double %355, double %356, double %360)
%362 = load double*, double** %5
%363 = load i32, i32* %10
%364 = sext i32 %363 to i64
%365 = getelementptr inbounds double, double* %362, i64 %364
store double %361, double* %365
%366 = load double, double* %16
%367 = load double, double* %21
%368 = load double, double* %17
%369 = load double, double* %20
%370 = fmul double %368, %369
%371 = call double @llvm.fmuladd.f64(double %366, double %367, double %370)
%372 = load double*, double** %5
%373 = load i32, i32* %10
%374 = add nsw i32 %373, 1
%375 = sext i32 %374 to i64
%376 = getelementptr inbounds double, double* %372, i64 %375
store double %371, double* %376
%377 = load double, double* %22
%378 = load double, double* %27
%379 = fadd double %377, %378
store double %379, double* %20
%380 = load double, double* %23
%381 = load double, double* %26
%382 = fsub double %380, %381
store double %382, double* %21
%383 = load double, double* %18
%384 = load double, double* %20
%385 = load double, double* %19
%386 = load double, double* %21
%387 = fmul double %385, %386
%388 = call double @llvm.fmuladd.f64(double %383, double %384, double %387)
%389 = load double*, double** %5
%390 = load i32, i32* %11
%391 = sext i32 %390 to i64
%392 = getelementptr inbounds double, double* %389, i64 %391
store double %388, double* %392
%393 = load double, double* %18
%394 = load double, double* %21
%395 = load double, double* %19
%396 = load double, double* %20
%397 = fmul double %395, %396
%398 = fneg double %397
%399 = call double @llvm.fmuladd.f64(double %393, double %394, double %398)
%400 = load double*, double** %5
%401 = load i32, i32* %11
%402 = add nsw i32 %401, 1
%403 = sext i32 %402 to i64
%404 = getelementptr inbounds double, double* %400, i64 %403
store double %399, double* %404
%405 = load i32, i32* %13
%406 = load i32, i32* %7
%407 = sub nsw i32 %405, %406
store i32 %407, i32* %8
%408 = load i32, i32* %8
%409 = load i32, i32* %13
%410 = add nsw i32 %408, %409
store i32 %410, i32* %9
%411 = load i32, i32* %9
%412 = load i32, i32* %13
%413 = add nsw i32 %411, %412
store i32 %413, i32* %10
%414 = load i32, i32* %10
%415 = load i32, i32* %13
%416 = add nsw i32 %414, %415
store i32 %416, i32* %11
%417 = load double*, double** %5
%418 = load i32, i32* %8
%419 = sext i32 %418 to i64
%420 = getelementptr inbounds double, double* %417, i64 %419
%421 = load double, double* %420
%422 = load double*, double** %5
%423 = load i32, i32* %10
%424 = sext i32 %423 to i64
%425 = getelementptr inbounds double, double* %422, i64 %424
%426 = load double, double* %425
%427 = fadd double %421, %426
store double %427, double* %20
%428 = load double*, double** %5
%429 = load i32, i32* %8
%430 = add nsw i32 %429, 1
%431 = sext i32 %430 to i64
%432 = getelementptr inbounds double, double* %428, i64 %431
%433 = load double, double* %432
%434 = load double*, double** %5
%435 = load i32, i32* %10
%436 = add nsw i32 %435, 1
%437 = sext i32 %436 to i64
%438 = getelementptr inbounds double, double* %434, i64 %437
%439 = load double, double* %438
%440 = fadd double %433, %439
store double %440, double* %21
%441 = load double*, double** %5
%442 = load i32, i32* %8
%443 = sext i32 %442 to i64
%444 = getelementptr inbounds double, double* %441, i64 %443
%445 = load double, double* %444
%446 = load double*, double** %5
%447 = load i32, i32* %10
%448 = sext i32 %447 to i64
%449 = getelementptr inbounds double, double* %446, i64 %448
%450 = load double, double* %449
%451 = fsub double %445, %450
store double %451, double* %22
%452 = load double*, double** %5
%453 = load i32, i32* %8
%454 = add nsw i32 %453, 1
%455 = sext i32 %454 to i64
%456 = getelementptr inbounds double, double* %452, i64 %455
%457 = load double, double* %456
%458 = load double*, double** %5
%459 = load i32, i32* %10
%460 = add nsw i32 %459, 1
%461 = sext i32 %460 to i64
%462 = getelementptr inbounds double, double* %458, i64 %461
%463 = load double, double* %462
%464 = fsub double %457, %463
store double %464, double* %23
%465 = load double*, double** %5
%466 = load i32, i32* %9
%467 = sext i32 %466 to i64
%468 = getelementptr inbounds double, double* %465, i64 %467
%469 = load double, double* %468
%470 = load double*, double** %5
%471 = load i32, i32* %11
%472 = sext i32 %471 to i64
%473 = getelementptr inbounds double, double* %470, i64 %472
%474 = load double, double* %473
%475 = fadd double %469, %474
store double %475, double* %24
%476 = load double*, double** %5
%477 = load i32, i32* %9
%478 = add nsw i32 %477, 1
%479 = sext i32 %478 to i64
%480 = getelementptr inbounds double, double* %476, i64 %479
%481 = load double, double* %480
%482 = load double*, double** %5
%483 = load i32, i32* %11
%484 = add nsw i32 %483, 1
%485 = sext i32 %484 to i64
%486 = getelementptr inbounds double, double* %482, i64 %485
%487 = load double, double* %486
%488 = fadd double %481, %487
store double %488, double* %25
%489 = load double*, double** %5
%490 = load i32, i32* %9
%491 = sext i32 %490 to i64
%492 = getelementptr inbounds double, double* %489, i64 %491
%493 = load double, double* %492
%494 = load double*, double** %5
%495 = load i32, i32* %11
%496 = sext i32 %495 to i64
%497 = getelementptr inbounds double, double* %494, i64 %496
%498 = load double, double* %497
%499 = fsub double %493, %498
store double %499, double* %26
%500 = load double*, double** %5
%501 = load i32, i32* %9
%502 = add nsw i32 %501, 1
%503 = sext i32 %502 to i64
%504 = getelementptr inbounds double, double* %500, i64 %503
%505 = load double, double* %504
%506 = load double*, double** %5
%507 = load i32, i32* %11
%508 = add nsw i32 %507, 1
%509 = sext i32 %508 to i64
%510 = getelementptr inbounds double, double* %506, i64 %509
%511 = load double, double* %510
%512 = fsub double %505, %511
store double %512, double* %27
%513 = load double, double* %20
%514 = load double, double* %24
%515 = fadd double %513, %514
%516 = load double*, double** %5
%517 = load i32, i32* %8
%518 = sext i32 %517 to i64
%519 = getelementptr inbounds double, double* %516, i64 %518
store double %515, double* %519
%520 = load double, double* %21
%521 = load double, double* %25
%522 = fadd double %520, %521
%523 = load double*, double** %5
%524 = load i32, i32* %8
%525 = add nsw i32 %524, 1
%526 = sext i32 %525 to i64
%527 = getelementptr inbounds double, double* %523, i64 %526
store double %522, double* %527
%528 = load double, double* %20
%529 = load double, double* %24
%530 = fsub double %528, %529
%531 = load double*, double** %5
%532 = load i32, i32* %9
%533 = sext i32 %532 to i64
%534 = getelementptr inbounds double, double* %531, i64 %533
store double %530, double* %534
%535 = load double, double* %21
%536 = load double, double* %25
%537 = fsub double %535, %536
%538 = load double*, double** %5
%539 = load i32, i32* %9
%540 = add nsw i32 %539, 1
%541 = sext i32 %540 to i64
%542 = getelementptr inbounds double, double* %538, i64 %541
store double %537, double* %542
%543 = load double, double* %22
%544 = load double, double* %27
%545 = fsub double %543, %544
store double %545, double* %20
%546 = load double, double* %23
%547 = load double, double* %26
%548 = fadd double %546, %547
store double %548, double* %21
%549 = load double, double* %17
%550 = load double, double* %20
%551 = load double, double* %16
%552 = load double, double* %21
%553 = fmul double %551, %552
%554 = fneg double %553
%555 = call double @llvm.fmuladd.f64(double %549, double %550, double %554)
%556 = load double*, double** %5
%557 = load i32, i32* %10
%558 = sext i32 %557 to i64
%559 = getelementptr inbounds double, double* %556, i64 %558
store double %555, double* %559
%560 = load double, double* %17
%561 = load double, double* %21
%562 = load double, double* %16
%563 = load double, double* %20
%564 = fmul double %562, %563
%565 = call double @llvm.fmuladd.f64(double %560, double %561, double %564)
%566 = load double*, double** %5
%567 = load i32, i32* %10
%568 = add nsw i32 %567, 1
%569 = sext i32 %568 to i64
%570 = getelementptr inbounds double, double* %566, i64 %569
store double %565, double* %570
%571 = load double, double* %22
%572 = load double, double* %27
%573 = fadd double %571, %572
store double %573, double* %20
%574 = load double, double* %23
%575 = load double, double* %26
%576 = fsub double %574, %575
store double %576, double* %21
%577 = load double, double* %19
%578 = load double, double* %20
%579 = load double, double* %18
%580 = load double, double* %21
%581 = fmul double %579, %580
%582 = call double @llvm.fmuladd.f64(double %577, double %578, double %581)
%583 = load double*, double** %5
%584 = load i32, i32* %11
%585 = sext i32 %584 to i64
%586 = getelementptr inbounds double, double* %583, i64 %585
store double %582, double* %586
%587 = load double, double* %19
%588 = load double, double* %21
%589 = load double, double* %18
%590 = load double, double* %20
%591 = fmul double %589, %590
%592 = fneg double %591
%593 = call double @llvm.fmuladd.f64(double %587, double %588, double %592)
%594 = load double*, double** %5
%595 = load i32, i32* %11
%596 = add nsw i32 %595, 1
%597 = sext i32 %596 to i64
%598 = getelementptr inbounds double, double* %594, i64 %597
store double %593, double* %598
br label %599
599: ; preds = %187
%600 = load i32, i32* %7
%601 = add nsw i32 %600, 2
store i32 %601, i32* %7
br label %183
602: ; preds = %183
%603 = load i32, i32* %14
store i32 %603, i32* %8
%604 = load i32, i32* %8
%605 = load i32, i32* %13
%606 = add nsw i32 %604, %605
store i32 %606, i32* %9
%607 = load i32, i32* %9
%608 = load i32, i32* %13
%609 = add nsw i32 %607, %608
store i32 %609, i32* %10
%610 = load i32, i32* %10
%611 = load i32, i32* %13
%612 = add nsw i32 %610, %611
store i32 %612, i32* %11
%613 = load double*, double** %5
%614 = load i32, i32* %8
%615 = sext i32 %614 to i64
%616 = getelementptr inbounds double, double* %613, i64 %615
%617 = load double, double* %616
%618 = load double*, double** %5
%619 = load i32, i32* %10
%620 = sext i32 %619 to i64
%621 = getelementptr inbounds double, double* %618, i64 %620
%622 = load double, double* %621
%623 = fadd double %617, %622
store double %623, double* %20
%624 = load double*, double** %5
%625 = load i32, i32* %8
%626 = add nsw i32 %625, 1
%627 = sext i32 %626 to i64
%628 = getelementptr inbounds double, double* %624, i64 %627
%629 = load double, double* %628
%630 = load double*, double** %5
%631 = load i32, i32* %10
%632 = add nsw i32 %631, 1
%633 = sext i32 %632 to i64
%634 = getelementptr inbounds double, double* %630, i64 %633
%635 = load double, double* %634
%636 = fadd double %629, %635
store double %636, double* %21
%637 = load double*, double** %5
%638 = load i32, i32* %8
%639 = sext i32 %638 to i64
%640 = getelementptr inbounds double, double* %637, i64 %639
%641 = load double, double* %640
%642 = load double*, double** %5
%643 = load i32, i32* %10
%644 = sext i32 %643 to i64
%645 = getelementptr inbounds double, double* %642, i64 %644
%646 = load double, double* %645
%647 = fsub double %641, %646
store double %647, double* %22
%648 = load double*, double** %5
%649 = load i32, i32* %8
%650 = add nsw i32 %649, 1
%651 = sext i32 %650 to i64
%652 = getelementptr inbounds double, double* %648, i64 %651
%653 = load double, double* %652
%654 = load double*, double** %5
%655 = load i32, i32* %10
%656 = add nsw i32 %655, 1
%657 = sext i32 %656 to i64
%658 = getelementptr inbounds double, double* %654, i64 %657
%659 = load double, double* %658
%660 = fsub double %653, %659
store double %660, double* %23
%661 = load double*, double** %5
%662 = load i32, i32* %9
%663 = sext i32 %662 to i64
%664 = getelementptr inbounds double, double* %661, i64 %663
%665 = load double, double* %664
%666 = load double*, double** %5
%667 = load i32, i32* %11
%668 = sext i32 %667 to i64
%669 = getelementptr inbounds double, double* %666, i64 %668
%670 = load double, double* %669
%671 = fadd double %665, %670
store double %671, double* %24
%672 = load double*, double** %5
%673 = load i32, i32* %9
%674 = add nsw i32 %673, 1
%675 = sext i32 %674 to i64
%676 = getelementptr inbounds double, double* %672, i64 %675
%677 = load double, double* %676
%678 = load double*, double** %5
%679 = load i32, i32* %11
%680 = add nsw i32 %679, 1
%681 = sext i32 %680 to i64
%682 = getelementptr inbounds double, double* %678, i64 %681
%683 = load double, double* %682
%684 = fadd double %677, %683
store double %684, double* %25
%685 = load double*, double** %5
%686 = load i32, i32* %9
%687 = sext i32 %686 to i64
%688 = getelementptr inbounds double, double* %685, i64 %687
%689 = load double, double* %688
%690 = load double*, double** %5
%691 = load i32, i32* %11
%692 = sext i32 %691 to i64
%693 = getelementptr inbounds double, double* %690, i64 %692
%694 = load double, double* %693
%695 = fsub double %689, %694
store double %695, double* %26
%696 = load double*, double** %5
%697 = load i32, i32* %9
%698 = add nsw i32 %697, 1
%699 = sext i32 %698 to i64
%700 = getelementptr inbounds double, double* %696, i64 %699
%701 = load double, double* %700
%702 = load double*, double** %5
%703 = load i32, i32* %11
%704 = add nsw i32 %703, 1
%705 = sext i32 %704 to i64
%706 = getelementptr inbounds double, double* %702, i64 %705
%707 = load double, double* %706
%708 = fsub double %701, %707
store double %708, double* %27
%709 = load double, double* %20
%710 = load double, double* %24
%711 = fadd double %709, %710
%712 = load double*, double** %5
%713 = load i32, i32* %8
%714 = sext i32 %713 to i64
%715 = getelementptr inbounds double, double* %712, i64 %714
store double %711, double* %715
%716 = load double, double* %21
%717 = load double, double* %25
%718 = fadd double %716, %717
%719 = load double*, double** %5
%720 = load i32, i32* %8
%721 = add nsw i32 %720, 1
%722 = sext i32 %721 to i64
%723 = getelementptr inbounds double, double* %719, i64 %722
store double %718, double* %723
%724 = load double, double* %20
%725 = load double, double* %24
%726 = fsub double %724, %725
%727 = load double*, double** %5
%728 = load i32, i32* %9
%729 = sext i32 %728 to i64
%730 = getelementptr inbounds double, double* %727, i64 %729
store double %726, double* %730
%731 = load double, double* %21
%732 = load double, double* %25
%733 = fsub double %731, %732
%734 = load double*, double** %5
%735 = load i32, i32* %9
%736 = add nsw i32 %735, 1
%737 = sext i32 %736 to i64
%738 = getelementptr inbounds double, double* %734, i64 %737
store double %733, double* %738
%739 = load double, double* %22
%740 = load double, double* %27
%741 = fsub double %739, %740
store double %741, double* %20
%742 = load double, double* %23
%743 = load double, double* %26
%744 = fadd double %742, %743
store double %744, double* %21
%745 = load double, double* %15
%746 = load double, double* %20
%747 = load double, double* %21
%748 = fsub double %746, %747
%749 = fmul double %745, %748
%750 = load double*, double** %5
%751 = load i32, i32* %10
%752 = sext i32 %751 to i64
%753 = getelementptr inbounds double, double* %750, i64 %752
store double %749, double* %753
%754 = load double, double* %15
%755 = load double, double* %21
%756 = load double, double* %20
%757 = fadd double %755, %756
%758 = fmul double %754, %757
%759 = load double*, double** %5
%760 = load i32, i32* %10
%761 = add nsw i32 %760, 1
%762 = sext i32 %761 to i64
%763 = getelementptr inbounds double, double* %759, i64 %762
store double %758, double* %763
%764 = load double, double* %22
%765 = load double, double* %27
%766 = fadd double %764, %765
store double %766, double* %20
%767 = load double, double* %23
%768 = load double, double* %26
%769 = fsub double %767, %768
store double %769, double* %21
%770 = load double, double* %15
%771 = fneg double %770
%772 = load double, double* %20
%773 = load double, double* %21
%774 = fadd double %772, %773
%775 = fmul double %771, %774
%776 = load double*, double** %5
%777 = load i32, i32* %11
%778 = sext i32 %777 to i64
%779 = getelementptr inbounds double, double* %776, i64 %778
store double %775, double* %779
%780 = load double, double* %15
%781 = fneg double %780
%782 = load double, double* %21
%783 = load double, double* %20
%784 = fsub double %782, %783
%785 = fmul double %781, %784
%786 = load double*, double** %5
%787 = load i32, i32* %11
%788 = add nsw i32 %787, 1
%789 = sext i32 %788 to i64
%790 = getelementptr inbounds double, double* %786, i64 %789
store double %785, double* %790
ret void
}

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