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    197c <cftmdl1+0x87c>
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
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     149e <cftmdl1+0x39e>
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
mov    %rcx,-0x10(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rbp,-0x8(%rsp)
movsd  -0x18(%rdx),%xmm11
movsd  -0x10(%rdx),%xmm12
movsd  -0x8(%rdx),%xmm10
movsd  (%rdx),%xmm9
movsd  (%rsi,%r11,1),%xmm1
movupd (%rsi,%r8,8),%xmm5
movsd  (%rsi,%r12,8),%xmm2
movapd %xmm5,%xmm7
subsd  %xmm1,%xmm7
unpcklpd %xmm2,%xmm1
movsd  (%rsi,%r9,1),%xmm3
addpd  %xmm5,%xmm1
movupd (%rsi,%rbp,8),%xmm6
mov    -0x20(%rsp),%rcx
movsd  (%rsi,%rcx,8),%xmm4
movapd %xmm6,%xmm0
subsd  %xmm3,%xmm0
unpckhpd %xmm5,%xmm5
unpcklpd %xmm4,%xmm3
addpd  %xmm6,%xmm3
unpckhpd %xmm6,%xmm6
subsd  %xmm2,%xmm5
subsd  %xmm4,%xmm6
movapd %xmm1,%xmm2
addpd  %xmm3,%xmm2
subpd  %xmm3,%xmm1
movupd %xmm2,(%rsi,%r8,8)
movupd %xmm1,(%rsi,%rbp,8)
movapd %xmm7,%xmm1
subsd  %xmm6,%xmm1
movapd %xmm5,%xmm2
addsd  %xmm0,%xmm2
movapd %xmm11,%xmm3
mulsd  %xmm2,%xmm3
mulsd  %xmm12,%xmm2
movapd %xmm11,%xmm4
mulsd  %xmm1,%xmm4
mulsd  %xmm12,%xmm1
subsd  %xmm2,%xmm4
addsd  %xmm3,%xmm1
addsd  %xmm7,%xmm6
subsd  %xmm0,%xmm5
movapd %xmm9,%xmm0
movsd  %xmm4,(%rsi,%r11,1)
mulsd  %xmm5,%xmm0
movapd %xmm10,%xmm2
mulsd  %xmm6,%xmm2
addsd  %xmm0,%xmm2
movsd  %xmm1,(%rsi,%r12,8)
mulsd  %xmm9,%xmm6
mulsd  %xmm10,%xmm5
subsd  %xmm6,%xmm5
movsd  %xmm2,(%rsi,%r9,1)
movsd  %xmm5,(%rsi,%rcx,8)
movsd  (%rsi,%r15,1),%xmm1
movupd (%rsi,%rax,8),%xmm5
movsd  (%rsi,%rdi,8),%xmm2
movapd %xmm5,%xmm13
subsd  %xmm1,%xmm13
unpcklpd %xmm2,%xmm1
movsd  (%rsi,%r10,1),%xmm7
addpd  %xmm5,%xmm1
movsd  (%rsi,%r14,1),%xmm3
mov    %r11,%rbp
mov    %r15,%r11
mov    %r10,%r15
mov    -0x18(%rsp),%r10
movsd  (%rsi,%r10,8),%xmm6
mov    %rbx,%r13
sar    $0x1d,%r13
unpckhpd %xmm5,%xmm5
or     $0x8,%r13
movsd  (%rsi,%r13,1),%xmm4
movapd %xmm7,%xmm0
unpcklpd %xmm6,%xmm0
subsd  %xmm2,%xmm5
subsd  %xmm3,%xmm7
unpcklpd %xmm4,%xmm3
subsd  %xmm4,%xmm6
addpd  %xmm0,%xmm3
movapd %xmm1,%xmm0
addpd  %xmm3,%xmm0
movupd %xmm0,(%rsi,%rax,8)
subpd  %xmm3,%xmm1
movlpd %xmm1,(%rsi,%r15,1)
movhpd %xmm1,(%rsi,%r10,8)
movapd %xmm13,%xmm0
subsd  %xmm6,%xmm0
movapd %xmm5,%xmm1
addsd  %xmm7,%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm11,%xmm2
mulsd  %xmm12,%xmm1
mulsd  %xmm0,%xmm12
mulsd  %xmm11,%xmm0
subsd  %xmm2,%xmm12
addsd  %xmm0,%xmm1
addsd  %xmm13,%xmm6
subsd  %xmm7,%xmm5
movapd %xmm9,%xmm0
movsd  %xmm12,(%rsi,%r11,1)
mulsd  %xmm6,%xmm0
mulsd  %xmm10,%xmm6
mulsd  %xmm5,%xmm10
movsd  %xmm1,(%rsi,%rdi,8)
addsd  %xmm10,%xmm0
mulsd  %xmm9,%xmm5
subsd  %xmm6,%xmm5
movsd  %xmm0,(%rsi,%r14,1)
movsd  %xmm5,(%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  0xab2(%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 = %598, %3
%184 = load i32, i32* %7
%185 = load i32, i32* %14
%186 = icmp slt i32 %184, %185
br i1 %186, label %187, label %601
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
store double %212, double* %19
%213 = load i32, i32* %7
%214 = load i32, i32* %13
%215 = add nsw i32 %213, %214
store i32 %215, i32* %9
%216 = load i32, i32* %9
%217 = load i32, i32* %13
%218 = add nsw i32 %216, %217
store i32 %218, i32* %10
%219 = load i32, i32* %10
%220 = load i32, i32* %13
%221 = add nsw i32 %219, %220
store i32 %221, i32* %11
%222 = load double*, double** %5
%223 = load i32, i32* %7
%224 = sext i32 %223 to i64
%225 = getelementptr inbounds double, double* %222, i64 %224
%226 = load double, double* %225
%227 = load double*, double** %5
%228 = load i32, i32* %10
%229 = sext i32 %228 to i64
%230 = getelementptr inbounds double, double* %227, i64 %229
%231 = load double, double* %230
%232 = fadd double %226, %231
store double %232, double* %20
%233 = load double*, double** %5
%234 = load i32, i32* %7
%235 = add nsw i32 %234, 1
%236 = sext i32 %235 to i64
%237 = getelementptr inbounds double, double* %233, i64 %236
%238 = load double, double* %237
%239 = load double*, double** %5
%240 = load i32, i32* %10
%241 = add nsw i32 %240, 1
%242 = sext i32 %241 to i64
%243 = getelementptr inbounds double, double* %239, i64 %242
%244 = load double, double* %243
%245 = fadd double %238, %244
store double %245, double* %21
%246 = load double*, double** %5
%247 = load i32, i32* %7
%248 = sext i32 %247 to i64
%249 = getelementptr inbounds double, double* %246, i64 %248
%250 = load double, double* %249
%251 = load double*, double** %5
%252 = load i32, i32* %10
%253 = sext i32 %252 to i64
%254 = getelementptr inbounds double, double* %251, i64 %253
%255 = load double, double* %254
%256 = fsub double %250, %255
store double %256, double* %22
%257 = load double*, double** %5
%258 = load i32, i32* %7
%259 = add nsw i32 %258, 1
%260 = sext i32 %259 to i64
%261 = getelementptr inbounds double, double* %257, i64 %260
%262 = load double, double* %261
%263 = load double*, double** %5
%264 = load i32, i32* %10
%265 = add nsw i32 %264, 1
%266 = sext i32 %265 to i64
%267 = getelementptr inbounds double, double* %263, i64 %266
%268 = load double, double* %267
%269 = fsub double %262, %268
store double %269, double* %23
%270 = load double*, double** %5
%271 = load i32, i32* %9
%272 = sext i32 %271 to i64
%273 = getelementptr inbounds double, double* %270, i64 %272
%274 = load double, double* %273
%275 = load double*, double** %5
%276 = load i32, i32* %11
%277 = sext i32 %276 to i64
%278 = getelementptr inbounds double, double* %275, i64 %277
%279 = load double, double* %278
%280 = fadd double %274, %279
store double %280, double* %24
%281 = load double*, double** %5
%282 = load i32, i32* %9
%283 = add nsw i32 %282, 1
%284 = sext i32 %283 to i64
%285 = getelementptr inbounds double, double* %281, i64 %284
%286 = load double, double* %285
%287 = load double*, double** %5
%288 = load i32, i32* %11
%289 = add nsw i32 %288, 1
%290 = sext i32 %289 to i64
%291 = getelementptr inbounds double, double* %287, i64 %290
%292 = load double, double* %291
%293 = fadd double %286, %292
store double %293, double* %25
%294 = load double*, double** %5
%295 = load i32, i32* %9
%296 = sext i32 %295 to i64
%297 = getelementptr inbounds double, double* %294, i64 %296
%298 = load double, double* %297
%299 = load double*, double** %5
%300 = load i32, i32* %11
%301 = sext i32 %300 to i64
%302 = getelementptr inbounds double, double* %299, i64 %301
%303 = load double, double* %302
%304 = fsub double %298, %303
store double %304, double* %26
%305 = load double*, double** %5
%306 = load i32, i32* %9
%307 = add nsw i32 %306, 1
%308 = sext i32 %307 to i64
%309 = getelementptr inbounds double, double* %305, i64 %308
%310 = load double, double* %309
%311 = load double*, double** %5
%312 = load i32, i32* %11
%313 = add nsw i32 %312, 1
%314 = sext i32 %313 to i64
%315 = getelementptr inbounds double, double* %311, i64 %314
%316 = load double, double* %315
%317 = fsub double %310, %316
store double %317, double* %27
%318 = load double, double* %20
%319 = load double, double* %24
%320 = fadd double %318, %319
%321 = load double*, double** %5
%322 = load i32, i32* %7
%323 = sext i32 %322 to i64
%324 = getelementptr inbounds double, double* %321, i64 %323
store double %320, double* %324
%325 = load double, double* %21
%326 = load double, double* %25
%327 = fadd double %325, %326
%328 = load double*, double** %5
%329 = load i32, i32* %7
%330 = add nsw i32 %329, 1
%331 = sext i32 %330 to i64
%332 = getelementptr inbounds double, double* %328, i64 %331
store double %327, double* %332
%333 = load double, double* %20
%334 = load double, double* %24
%335 = fsub double %333, %334
%336 = load double*, double** %5
%337 = load i32, i32* %9
%338 = sext i32 %337 to i64
%339 = getelementptr inbounds double, double* %336, i64 %338
store double %335, double* %339
%340 = load double, double* %21
%341 = load double, double* %25
%342 = fsub double %340, %341
%343 = load double*, double** %5
%344 = load i32, i32* %9
%345 = add nsw i32 %344, 1
%346 = sext i32 %345 to i64
%347 = getelementptr inbounds double, double* %343, i64 %346
store double %342, double* %347
%348 = load double, double* %22
%349 = load double, double* %27
%350 = fsub double %348, %349
store double %350, double* %20
%351 = load double, double* %23
%352 = load double, double* %26
%353 = fadd double %351, %352
store double %353, double* %21
%354 = load double, double* %16
%355 = load double, double* %20
%356 = load double, double* %17
%357 = load double, double* %21
%358 = fmul double %356, %357
%359 = fneg double %358
%360 = call double @llvm.fmuladd.f64(double %354, double %355, double %359)
%361 = load double*, double** %5
%362 = load i32, i32* %10
%363 = sext i32 %362 to i64
%364 = getelementptr inbounds double, double* %361, i64 %363
store double %360, double* %364
%365 = load double, double* %16
%366 = load double, double* %21
%367 = load double, double* %17
%368 = load double, double* %20
%369 = fmul double %367, %368
%370 = call double @llvm.fmuladd.f64(double %365, double %366, double %369)
%371 = load double*, double** %5
%372 = load i32, i32* %10
%373 = add nsw i32 %372, 1
%374 = sext i32 %373 to i64
%375 = getelementptr inbounds double, double* %371, i64 %374
store double %370, double* %375
%376 = load double, double* %22
%377 = load double, double* %27
%378 = fadd double %376, %377
store double %378, double* %20
%379 = load double, double* %23
%380 = load double, double* %26
%381 = fsub double %379, %380
store double %381, double* %21
%382 = load double, double* %18
%383 = load double, double* %20
%384 = load double, double* %19
%385 = load double, double* %21
%386 = fmul double %384, %385
%387 = call double @llvm.fmuladd.f64(double %382, double %383, double %386)
%388 = load double*, double** %5
%389 = load i32, i32* %11
%390 = sext i32 %389 to i64
%391 = getelementptr inbounds double, double* %388, i64 %390
store double %387, double* %391
%392 = load double, double* %18
%393 = load double, double* %21
%394 = load double, double* %19
%395 = load double, double* %20
%396 = fmul double %394, %395
%397 = fneg double %396
%398 = call double @llvm.fmuladd.f64(double %392, double %393, double %397)
%399 = load double*, double** %5
%400 = load i32, i32* %11
%401 = add nsw i32 %400, 1
%402 = sext i32 %401 to i64
%403 = getelementptr inbounds double, double* %399, i64 %402
store double %398, double* %403
%404 = load i32, i32* %13
%405 = load i32, i32* %7
%406 = sub nsw i32 %404, %405
store i32 %406, i32* %8
%407 = load i32, i32* %8
%408 = load i32, i32* %13
%409 = add nsw i32 %407, %408
store i32 %409, i32* %9
%410 = load i32, i32* %9
%411 = load i32, i32* %13
%412 = add nsw i32 %410, %411
store i32 %412, i32* %10
%413 = load i32, i32* %10
%414 = load i32, i32* %13
%415 = add nsw i32 %413, %414
store i32 %415, i32* %11
%416 = load double*, double** %5
%417 = load i32, i32* %8
%418 = sext i32 %417 to i64
%419 = getelementptr inbounds double, double* %416, i64 %418
%420 = load double, double* %419
%421 = load double*, double** %5
%422 = load i32, i32* %10
%423 = sext i32 %422 to i64
%424 = getelementptr inbounds double, double* %421, i64 %423
%425 = load double, double* %424
%426 = fadd double %420, %425
store double %426, double* %20
%427 = load double*, double** %5
%428 = load i32, i32* %8
%429 = add nsw i32 %428, 1
%430 = sext i32 %429 to i64
%431 = getelementptr inbounds double, double* %427, i64 %430
%432 = load double, double* %431
%433 = load double*, double** %5
%434 = load i32, i32* %10
%435 = add nsw i32 %434, 1
%436 = sext i32 %435 to i64
%437 = getelementptr inbounds double, double* %433, i64 %436
%438 = load double, double* %437
%439 = fadd double %432, %438
store double %439, double* %21
%440 = load double*, double** %5
%441 = load i32, i32* %8
%442 = sext i32 %441 to i64
%443 = getelementptr inbounds double, double* %440, i64 %442
%444 = load double, double* %443
%445 = load double*, double** %5
%446 = load i32, i32* %10
%447 = sext i32 %446 to i64
%448 = getelementptr inbounds double, double* %445, i64 %447
%449 = load double, double* %448
%450 = fsub double %444, %449
store double %450, double* %22
%451 = load double*, double** %5
%452 = load i32, i32* %8
%453 = add nsw i32 %452, 1
%454 = sext i32 %453 to i64
%455 = getelementptr inbounds double, double* %451, i64 %454
%456 = load double, double* %455
%457 = load double*, double** %5
%458 = load i32, i32* %10
%459 = add nsw i32 %458, 1
%460 = sext i32 %459 to i64
%461 = getelementptr inbounds double, double* %457, i64 %460
%462 = load double, double* %461
%463 = fsub double %456, %462
store double %463, double* %23
%464 = load double*, double** %5
%465 = load i32, i32* %9
%466 = sext i32 %465 to i64
%467 = getelementptr inbounds double, double* %464, i64 %466
%468 = load double, double* %467
%469 = load double*, double** %5
%470 = load i32, i32* %11
%471 = sext i32 %470 to i64
%472 = getelementptr inbounds double, double* %469, i64 %471
%473 = load double, double* %472
%474 = fadd double %468, %473
store double %474, double* %24
%475 = load double*, double** %5
%476 = load i32, i32* %9
%477 = add nsw i32 %476, 1
%478 = sext i32 %477 to i64
%479 = getelementptr inbounds double, double* %475, i64 %478
%480 = load double, double* %479
%481 = load double*, double** %5
%482 = load i32, i32* %11
%483 = add nsw i32 %482, 1
%484 = sext i32 %483 to i64
%485 = getelementptr inbounds double, double* %481, i64 %484
%486 = load double, double* %485
%487 = fadd double %480, %486
store double %487, double* %25
%488 = load double*, double** %5
%489 = load i32, i32* %9
%490 = sext i32 %489 to i64
%491 = getelementptr inbounds double, double* %488, i64 %490
%492 = load double, double* %491
%493 = load double*, double** %5
%494 = load i32, i32* %11
%495 = sext i32 %494 to i64
%496 = getelementptr inbounds double, double* %493, i64 %495
%497 = load double, double* %496
%498 = fsub double %492, %497
store double %498, double* %26
%499 = load double*, double** %5
%500 = load i32, i32* %9
%501 = add nsw i32 %500, 1
%502 = sext i32 %501 to i64
%503 = getelementptr inbounds double, double* %499, i64 %502
%504 = load double, double* %503
%505 = load double*, double** %5
%506 = load i32, i32* %11
%507 = add nsw i32 %506, 1
%508 = sext i32 %507 to i64
%509 = getelementptr inbounds double, double* %505, i64 %508
%510 = load double, double* %509
%511 = fsub double %504, %510
store double %511, double* %27
%512 = load double, double* %20
%513 = load double, double* %24
%514 = fadd double %512, %513
%515 = load double*, double** %5
%516 = load i32, i32* %8
%517 = sext i32 %516 to i64
%518 = getelementptr inbounds double, double* %515, i64 %517
store double %514, double* %518
%519 = load double, double* %21
%520 = load double, double* %25
%521 = fadd double %519, %520
%522 = load double*, double** %5
%523 = load i32, i32* %8
%524 = add nsw i32 %523, 1
%525 = sext i32 %524 to i64
%526 = getelementptr inbounds double, double* %522, i64 %525
store double %521, double* %526
%527 = load double, double* %20
%528 = load double, double* %24
%529 = fsub double %527, %528
%530 = load double*, double** %5
%531 = load i32, i32* %9
%532 = sext i32 %531 to i64
%533 = getelementptr inbounds double, double* %530, i64 %532
store double %529, double* %533
%534 = load double, double* %21
%535 = load double, double* %25
%536 = fsub double %534, %535
%537 = load double*, double** %5
%538 = load i32, i32* %9
%539 = add nsw i32 %538, 1
%540 = sext i32 %539 to i64
%541 = getelementptr inbounds double, double* %537, i64 %540
store double %536, double* %541
%542 = load double, double* %22
%543 = load double, double* %27
%544 = fsub double %542, %543
store double %544, double* %20
%545 = load double, double* %23
%546 = load double, double* %26
%547 = fadd double %545, %546
store double %547, double* %21
%548 = load double, double* %17
%549 = load double, double* %20
%550 = load double, double* %16
%551 = load double, double* %21
%552 = fmul double %550, %551
%553 = fneg double %552
%554 = call double @llvm.fmuladd.f64(double %548, double %549, double %553)
%555 = load double*, double** %5
%556 = load i32, i32* %10
%557 = sext i32 %556 to i64
%558 = getelementptr inbounds double, double* %555, i64 %557
store double %554, double* %558
%559 = load double, double* %17
%560 = load double, double* %21
%561 = load double, double* %16
%562 = load double, double* %20
%563 = fmul double %561, %562
%564 = call double @llvm.fmuladd.f64(double %559, double %560, double %563)
%565 = load double*, double** %5
%566 = load i32, i32* %10
%567 = add nsw i32 %566, 1
%568 = sext i32 %567 to i64
%569 = getelementptr inbounds double, double* %565, i64 %568
store double %564, double* %569
%570 = load double, double* %22
%571 = load double, double* %27
%572 = fadd double %570, %571
store double %572, double* %20
%573 = load double, double* %23
%574 = load double, double* %26
%575 = fsub double %573, %574
store double %575, double* %21
%576 = load double, double* %19
%577 = load double, double* %20
%578 = load double, double* %18
%579 = load double, double* %21
%580 = fmul double %578, %579
%581 = call double @llvm.fmuladd.f64(double %576, double %577, double %580)
%582 = load double*, double** %5
%583 = load i32, i32* %11
%584 = sext i32 %583 to i64
%585 = getelementptr inbounds double, double* %582, i64 %584
store double %581, double* %585
%586 = load double, double* %19
%587 = load double, double* %21
%588 = load double, double* %18
%589 = load double, double* %20
%590 = fmul double %588, %589
%591 = fneg double %590
%592 = call double @llvm.fmuladd.f64(double %586, double %587, double %591)
%593 = load double*, double** %5
%594 = load i32, i32* %11
%595 = add nsw i32 %594, 1
%596 = sext i32 %595 to i64
%597 = getelementptr inbounds double, double* %593, i64 %596
store double %592, double* %597
br label %598
598: ; preds = %187
%599 = load i32, i32* %7
%600 = add nsw i32 %599, 2
store i32 %600, i32* %7
br label %183
601: ; preds = %183
%602 = load i32, i32* %14
store i32 %602, i32* %8
%603 = load i32, i32* %8
%604 = load i32, i32* %13
%605 = add nsw i32 %603, %604
store i32 %605, i32* %9
%606 = load i32, i32* %9
%607 = load i32, i32* %13
%608 = add nsw i32 %606, %607
store i32 %608, i32* %10
%609 = load i32, i32* %10
%610 = load i32, i32* %13
%611 = add nsw i32 %609, %610
store i32 %611, i32* %11
%612 = load double*, double** %5
%613 = load i32, i32* %8
%614 = sext i32 %613 to i64
%615 = getelementptr inbounds double, double* %612, i64 %614
%616 = load double, double* %615
%617 = load double*, double** %5
%618 = load i32, i32* %10
%619 = sext i32 %618 to i64
%620 = getelementptr inbounds double, double* %617, i64 %619
%621 = load double, double* %620
%622 = fadd double %616, %621
store double %622, double* %20
%623 = load double*, double** %5
%624 = load i32, i32* %8
%625 = add nsw i32 %624, 1
%626 = sext i32 %625 to i64
%627 = getelementptr inbounds double, double* %623, i64 %626
%628 = load double, double* %627
%629 = load double*, double** %5
%630 = load i32, i32* %10
%631 = add nsw i32 %630, 1
%632 = sext i32 %631 to i64
%633 = getelementptr inbounds double, double* %629, i64 %632
%634 = load double, double* %633
%635 = fadd double %628, %634
store double %635, double* %21
%636 = load double*, double** %5
%637 = load i32, i32* %8
%638 = sext i32 %637 to i64
%639 = getelementptr inbounds double, double* %636, i64 %638
%640 = load double, double* %639
%641 = load double*, double** %5
%642 = load i32, i32* %10
%643 = sext i32 %642 to i64
%644 = getelementptr inbounds double, double* %641, i64 %643
%645 = load double, double* %644
%646 = fsub double %640, %645
store double %646, double* %22
%647 = load double*, double** %5
%648 = load i32, i32* %8
%649 = add nsw i32 %648, 1
%650 = sext i32 %649 to i64
%651 = getelementptr inbounds double, double* %647, i64 %650
%652 = load double, double* %651
%653 = load double*, double** %5
%654 = load i32, i32* %10
%655 = add nsw i32 %654, 1
%656 = sext i32 %655 to i64
%657 = getelementptr inbounds double, double* %653, i64 %656
%658 = load double, double* %657
%659 = fsub double %652, %658
store double %659, double* %23
%660 = load double*, double** %5
%661 = load i32, i32* %9
%662 = sext i32 %661 to i64
%663 = getelementptr inbounds double, double* %660, i64 %662
%664 = load double, double* %663
%665 = load double*, double** %5
%666 = load i32, i32* %11
%667 = sext i32 %666 to i64
%668 = getelementptr inbounds double, double* %665, i64 %667
%669 = load double, double* %668
%670 = fadd double %664, %669
store double %670, double* %24
%671 = load double*, double** %5
%672 = load i32, i32* %9
%673 = add nsw i32 %672, 1
%674 = sext i32 %673 to i64
%675 = getelementptr inbounds double, double* %671, i64 %674
%676 = load double, double* %675
%677 = load double*, double** %5
%678 = load i32, i32* %11
%679 = add nsw i32 %678, 1
%680 = sext i32 %679 to i64
%681 = getelementptr inbounds double, double* %677, i64 %680
%682 = load double, double* %681
%683 = fadd double %676, %682
store double %683, double* %25
%684 = load double*, double** %5
%685 = load i32, i32* %9
%686 = sext i32 %685 to i64
%687 = getelementptr inbounds double, double* %684, i64 %686
%688 = load double, double* %687
%689 = load double*, double** %5
%690 = load i32, i32* %11
%691 = sext i32 %690 to i64
%692 = getelementptr inbounds double, double* %689, i64 %691
%693 = load double, double* %692
%694 = fsub double %688, %693
store double %694, double* %26
%695 = load double*, double** %5
%696 = load i32, i32* %9
%697 = add nsw i32 %696, 1
%698 = sext i32 %697 to i64
%699 = getelementptr inbounds double, double* %695, i64 %698
%700 = load double, double* %699
%701 = load double*, double** %5
%702 = load i32, i32* %11
%703 = add nsw i32 %702, 1
%704 = sext i32 %703 to i64
%705 = getelementptr inbounds double, double* %701, i64 %704
%706 = load double, double* %705
%707 = fsub double %700, %706
store double %707, double* %27
%708 = load double, double* %20
%709 = load double, double* %24
%710 = fadd double %708, %709
%711 = load double*, double** %5
%712 = load i32, i32* %8
%713 = sext i32 %712 to i64
%714 = getelementptr inbounds double, double* %711, i64 %713
store double %710, double* %714
%715 = load double, double* %21
%716 = load double, double* %25
%717 = fadd double %715, %716
%718 = load double*, double** %5
%719 = load i32, i32* %8
%720 = add nsw i32 %719, 1
%721 = sext i32 %720 to i64
%722 = getelementptr inbounds double, double* %718, i64 %721
store double %717, double* %722
%723 = load double, double* %20
%724 = load double, double* %24
%725 = fsub double %723, %724
%726 = load double*, double** %5
%727 = load i32, i32* %9
%728 = sext i32 %727 to i64
%729 = getelementptr inbounds double, double* %726, i64 %728
store double %725, double* %729
%730 = load double, double* %21
%731 = load double, double* %25
%732 = fsub double %730, %731
%733 = load double*, double** %5
%734 = load i32, i32* %9
%735 = add nsw i32 %734, 1
%736 = sext i32 %735 to i64
%737 = getelementptr inbounds double, double* %733, i64 %736
store double %732, double* %737
%738 = load double, double* %22
%739 = load double, double* %27
%740 = fsub double %738, %739
store double %740, double* %20
%741 = load double, double* %23
%742 = load double, double* %26
%743 = fadd double %741, %742
store double %743, double* %21
%744 = load double, double* %15
%745 = load double, double* %20
%746 = load double, double* %21
%747 = fsub double %745, %746
%748 = fmul double %744, %747
%749 = load double*, double** %5
%750 = load i32, i32* %10
%751 = sext i32 %750 to i64
%752 = getelementptr inbounds double, double* %749, i64 %751
store double %748, double* %752
%753 = load double, double* %15
%754 = load double, double* %21
%755 = load double, double* %20
%756 = fadd double %754, %755
%757 = fmul double %753, %756
%758 = load double*, double** %5
%759 = load i32, i32* %10
%760 = add nsw i32 %759, 1
%761 = sext i32 %760 to i64
%762 = getelementptr inbounds double, double* %758, i64 %761
store double %757, double* %762
%763 = load double, double* %22
%764 = load double, double* %27
%765 = fadd double %763, %764
store double %765, double* %20
%766 = load double, double* %23
%767 = load double, double* %26
%768 = fsub double %766, %767
store double %768, double* %21
%769 = load double, double* %15
%770 = fneg double %769
%771 = load double, double* %20
%772 = load double, double* %21
%773 = fadd double %771, %772
%774 = fmul double %770, %773
%775 = load double*, double** %5
%776 = load i32, i32* %11
%777 = sext i32 %776 to i64
%778 = getelementptr inbounds double, double* %775, i64 %777
store double %774, double* %778
%779 = load double, double* %15
%780 = fneg double %779
%781 = load double, double* %21
%782 = load double, double* %20
%783 = fsub double %781, %782
%784 = fmul double %780, %783
%785 = load double*, double** %5
%786 = load i32, i32* %11
%787 = add nsw i32 %786, 1
%788 = sext i32 %787 to i64
%789 = getelementptr inbounds double, double* %785, i64 %788
store double %784, double* %789
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 = add nsw i64 %45, %16
%60 = add nsw i64 %59, %16
%61 = add i32 %42, %46
%62 = add nsw i64 %60, %16
%63 = add nsw i32 %61, %5
%64 = getelementptr inbounds double, double* %1, i64 %45
%65 = getelementptr inbounds double, double* %1, i64 %60
%66 = load double, double* %65
%67 = bitcast double* %64 to <2 x double>*
%68 = load <2 x double>, <2 x double>* %67
%69 = or i32 %61, 1
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds double, double* %1, i64 %70
%72 = load double, double* %71
%73 = insertelement <2 x double> poison, double %66, i64 0
%74 = insertelement <2 x double> %73, double %72, i64 1
%75 = fadd <2 x double> %68, %74
%76 = extractelement <2 x double> %68, i64 0
%77 = fsub double %76, %66
%78 = extractelement <2 x double> %68, i64 1
%79 = fsub double %78, %72
%80 = getelementptr inbounds double, double* %1, i64 %59
%81 = getelementptr inbounds double, double* %1, i64 %62
%82 = load double, double* %81
%83 = bitcast double* %80 to <2 x double>*
%84 = load <2 x double>, <2 x double>* %83
%85 = or i32 %63, 1
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %1, i64 %86
%88 = load double, double* %87
%89 = insertelement <2 x double> poison, double %82, i64 0
%90 = insertelement <2 x double> %89, double %88, i64 1
%91 = fadd <2 x double> %84, %90
%92 = extractelement <2 x double> %84, i64 0
%93 = fsub double %92, %82
%94 = extractelement <2 x double> %84, i64 1
%95 = fsub double %94, %88
%96 = fadd <2 x double> %75, %91
%97 = bitcast double* %64 to <2 x double>*
store <2 x double> %96, <2 x double>* %97
%98 = fsub <2 x double> %75, %91
%99 = bitcast double* %80 to <2 x double>*
store <2 x double> %98, <2 x double>* %99
%100 = fsub double %77, %95
%101 = fadd double %79, %93
%102 = fneg double %52
%103 = fmul double %101, %102
%104 = tail call double @llvm.fmuladd.f64(double %49, double %100, double %103)
store double %104, double* %65
%105 = fmul double %52, %100
%106 = tail call double @llvm.fmuladd.f64(double %49, double %101, double %105)
store double %106, double* %71
%107 = fadd double %77, %95
%108 = fsub double %79, %93
%109 = fmul double %58, %108
%110 = tail call double @llvm.fmuladd.f64(double %55, double %107, double %109)
store double %110, double* %81
%111 = fneg double %58
%112 = fmul double %107, %111
%113 = tail call double @llvm.fmuladd.f64(double %55, double %108, double %112)
store double %113, double* %87
%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 %58, 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 %58, 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
}