void dradb4(int ido, int l1, double *cc, double *ch, double *wa1, double *wa2,
            double *wa3) {
  static double sqrt2 = 1.4142135623730950488016887242097;
  int i, k, t0, t1, t2, t3, t4, t5, t6, t7, t8;
  double ci2, ci3, ci4, cr2, cr3, cr4, ti1, ti2, ti3, ti4, tr1, tr2, tr3, tr4;
  t0 = l1 * ido;
  t1 = 0;
  t2 = ido << 2;
  t3 = 0;
  t6 = ido << 1;
  for (k = 0; k < l1; k++) {
    t4 = t3 + t6;
    t5 = t1;
    tr3 = cc[t4 - 1] + cc[t4 - 1];
    tr4 = cc[t4] + cc[t4];
    tr1 = cc[t3] - cc[(t4 += t6) - 1];
    tr2 = cc[t3] + cc[t4 - 1];
    ch[t5] = tr2 + tr3;
    ch[t5 += t0] = tr1 - tr4;
    ch[t5 += t0] = tr2 - tr3;
    ch[t5 += t0] = tr1 + tr4;
    t1 += ido;
    t3 += t2;
  }
  if (ido < 2)
    return;
  if (ido == 2)
    goto L105;
  t1 = 0;
  for (k = 0; k < l1; k++) {
    t5 = (t4 = (t3 = (t2 = t1 << 2) + t6)) + t6;
    t7 = t1;
    for (i = 2; i < ido; i += 2) {
      t2 += 2;
      t3 += 2;
      t4 -= 2;
      t5 -= 2;
      t7 += 2;
      ti1 = cc[t2] + cc[t5];
      ti2 = cc[t2] - cc[t5];
      ti3 = cc[t3] - cc[t4];
      tr4 = cc[t3] + cc[t4];
      tr1 = cc[t2 - 1] - cc[t5 - 1];
      tr2 = cc[t2 - 1] + cc[t5 - 1];
      ti4 = cc[t3 - 1] - cc[t4 - 1];
      tr3 = cc[t3 - 1] + cc[t4 - 1];
      ch[t7 - 1] = tr2 + tr3;
      cr3 = tr2 - tr3;
      ch[t7] = ti2 + ti3;
      ci3 = ti2 - ti3;
      cr2 = tr1 - tr4;
      cr4 = tr1 + tr4;
      ci2 = ti1 + ti4;
      ci4 = ti1 - ti4;
      ch[(t8 = t7 + t0) - 1] = wa1[i - 2] * cr2 - wa1[i - 1] * ci2;
      ch[t8] = wa1[i - 2] * ci2 + wa1[i - 1] * cr2;
      ch[(t8 += t0) - 1] = wa2[i - 2] * cr3 - wa2[i - 1] * ci3;
      ch[t8] = wa2[i - 2] * ci3 + wa2[i - 1] * cr3;
      ch[(t8 += t0) - 1] = wa3[i - 2] * cr4 - wa3[i - 1] * ci4;
      ch[t8] = wa3[i - 2] * ci4 + wa3[i - 1] * cr4;
    }
    t1 += ido;
  }
  if (ido % 2 == 1)
    return;
L105:
  t1 = ido;
  t2 = ido << 2;
  t3 = ido - 1;
  t4 = ido + (ido << 1);
  for (k = 0; k < l1; k++) {
    t5 = t3;
    ti1 = cc[t1] + cc[t4];
    ti2 = cc[t4] - cc[t1];
    tr1 = cc[t1 - 1] - cc[t4 - 1];
    tr2 = cc[t1 - 1] + cc[t4 - 1];
    ch[t5] = tr2 + tr2;
    ch[t5 += t0] = sqrt2 * (tr1 - ti1);
    ch[t5 += t0] = ti2 + ti2;
    ch[t5 += t0] = -sqrt2 * (tr1 + ti1);
    t3 += ido;
    t1 += t2;
    t4 += t2;
  }
}


x86 assembly (clang, -O0):
<dradb4>:
push   %rbp
mov    %rsp,%rbp
sub    $0x48,%rsp
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    -0x8(%rbp),%eax
imul   -0x4(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x4(%rbp),%eax
shl    $0x2,%eax
mov    %eax,-0x3c(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x4(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x4c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    12bc <dradb4+0x1bc>
mov    -0x40(%rbp),%eax
add    -0x4c(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x38(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    -0x44(%rbp),%ecx
mov    %ecx,-0x44(%rbp)
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
ff
addsd  -0xc0(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xb0(%rbp),%xmm0
ff
subsd  -0xc8(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x48(%rbp),%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xb8(%rbp),%xmm0
ff
subsd  -0xc0(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x48(%rbp),%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xb0(%rbp),%xmm0
ff
addsd  -0xc8(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x48(%rbp),%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x3c(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1153 <dradb4+0x53>
cmpl   $0x2,-0x4(%rbp)
jge    12cb <dradb4+0x1cb>
jmp    194a <dradb4+0x84a>
cmpl   $0x2,-0x4(%rbp)
jne    12da <dradb4+0x1da>
jmp    178b <dradb4+0x68b>
movl   $0x0,-0x38(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    176d <dradb4+0x66d>
mov    -0x38(%rbp),%eax
shl    $0x2,%eax
mov    %eax,-0x3c(%rbp)
add    -0x4c(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    %eax,-0x44(%rbp)
add    -0x4c(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x38(%rbp),%eax
mov    %eax,-0x50(%rbp)
movl   $0x2,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1756 <dradb4+0x656>
mov    -0x3c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x40(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x40(%rbp)
mov    -0x44(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x44(%rbp)
mov    -0x48(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x48(%rbp)
mov    -0x50(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
ff
addsd  -0xc0(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xb8(%rbp),%xmm0
ff
subsd  -0xc0(%rbp),%xmm0
ff
movsd  %xmm0,-0x80(%rbp)
movsd  -0x98(%rbp),%xmm0
ff
addsd  -0xa0(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x98(%rbp),%xmm0
ff
subsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0x68(%rbp)
movsd  -0xb0(%rbp),%xmm0
ff
subsd  -0xc8(%rbp),%xmm0
ff
movsd  %xmm0,-0x78(%rbp)
movsd  -0xb0(%rbp),%xmm0
ff
addsd  -0xc8(%rbp),%xmm0
ff
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
addsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0x60(%rbp)
movsd  -0x90(%rbp),%xmm0
ff
subsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0x70(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x78(%rbp),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x60(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%ecx
add    -0x34(%rbp),%ecx
mov    %ecx,-0x54(%rbp)
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x60(%rbp),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x78(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x80(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x68(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x54(%rbp),%ecx
mov    %ecx,-0x54(%rbp)
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x68(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x80(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x88(%rbp),%xmm2
ff
mov    0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x70(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x54(%rbp),%ecx
mov    %ecx,-0x54(%rbp)
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x70(%rbp),%xmm2
mov    0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x2c(%rbp)
jmp    1319 <dradb4+0x219>
mov    -0x4(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12e8 <dradb4+0x1e8>
mov    -0x4(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x1,%edx
jne    1786 <dradb4+0x686>
jmp    194a <dradb4+0x84a>
jmp    178b <dradb4+0x68b>
mov    -0x4(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x4(%rbp),%eax
shl    $0x2,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x4(%rbp),%eax
mov    -0x4(%rbp),%ecx
shl    $0x1,%ecx
add    %ecx,%eax
mov    %eax,-0x44(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    194a <dradb4+0x84a>
mov    -0x40(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0xb8(%rbp),%xmm0
ff
addsd  -0xb8(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  0x2795(%rip),%xmm0
movsd  -0xb0(%rbp),%xmm1
ff
subsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x48(%rbp),%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x98(%rbp),%xmm0
ff
addsd  -0x98(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x48(%rbp),%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  0x273f(%rip),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0xb0(%rbp),%xmm1
ff
addsd  -0x90(%rbp),%xmm1
ff
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x48(%rbp),%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x3c(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x3c(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    17b8 <dradb4+0x6b8>
add    $0x48,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dradb4>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rcx,-0x38(%rsp)
mov    %esi,%eax
imul   %edi,%eax
mov    %eax,-0x44(%rsp)
lea    0x0(,%rdi,4),%ecx
mov    %rdi,-0x40(%rsp)
lea    (%rdi,%rdi,1),%eax
mov    %eax,-0x48(%rsp)
test   %esi,%esi
jle    11d8 <dradb4+0xd8>
movslq -0x40(%rsp),%r10
movslq -0x44(%rsp),%r11
movslq %ecx,%r14
movslq -0x48(%rsp),%r13
lea    -0x1(%r14),%ebp
shl    $0x3,%r14
shl    $0x3,%r10
lea    0x0(,%r11,8),%rax
lea    (%rax,%rax,2),%r15
mov    %r11,%r12
shl    $0x4,%r12
mov    %esi,%eax
mov    -0x38(%rsp),%rdi
mov    %rdx,%rbx
xchg   %ax,%ax
movsd  -0x8(%rbx,%r13,8),%xmm0
movsd  (%rbx,%r13,8),%xmm1
addsd  %xmm0,%xmm0
addsd  %xmm1,%xmm1
movsd  (%rbx),%xmm2
movslq %ebp,%rbp
movsd  (%rdx,%rbp,8),%xmm3
movapd %xmm2,%xmm4
subsd  %xmm3,%xmm4
addsd  %xmm3,%xmm2
movapd %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%rdi)
movapd %xmm4,%xmm3
subsd  %xmm1,%xmm3
movsd  %xmm3,(%rdi,%r11,8)
subsd  %xmm0,%xmm2
movsd  %xmm2,(%rdi,%r12,1)
addsd  %xmm1,%xmm4
movsd  %xmm4,(%rdi,%r15,1)
add    %ecx,%ebp
add    %r14,%rbx
add    %r10,%rdi
add    $0xffffffff,%eax
jne    1170 <dradb4+0x70>
cmpl   $0x2,-0x40(%rsp)
jl     154b <dradb4+0x44b>
mov    %ecx,%r15d
jne    11f5 <dradb4+0xf5>
test   %esi,%esi
jg     1474 <dradb4+0x374>
jmp    154b <dradb4+0x44b>
mov    %rdx,-0x30(%rsp)
test   %esi,%esi
mov    -0x48(%rsp),%r11d
jle    144d <dradb4+0x34d>
mov    0x38(%rsp),%r13
mov    -0x40(%rsp),%ecx
movslq -0x44(%rsp),%rax
mov    %rcx,-0x8(%rsp)
lea    0x0(,%rcx,8),%rcx
mov    %rcx,-0x10(%rsp)
mov    -0x38(%rsp),%r14
lea    (%r14,%rax,8),%rbp
lea    (%rax,%rax,2),%rcx
shl    $0x4,%rax
add    %r14,%rax
lea    (%r14,%rcx,8),%rbx
xor    %r10d,%r10d
mov    %r15d,%edx
xor    %edi,%edi
mov    %esi,-0x20(%rsp)
mov    %r15d,-0x24(%rsp)
mov    %edi,-0x1c(%rsp)
mov    %edx,-0x18(%rsp)
movslq %edx,%rcx
mov    -0x30(%rsp),%rdx
lea    (%rdx,%rcx,8),%r15
mov    %r11d,-0x48(%rsp)
movslq %r11d,%r12
lea    (%rdx,%r12,8),%rsi
shl    $0x20,%r12
mov    %r10d,-0x14(%rsp)
movslq %r10d,%rcx
lea    (%rdx,%rcx,8),%r11
mov    $0x2,%edi
mov    $0xfffffffffffffffe,%rcx
nopl   0x0(%rax)
mov    %r12,%r10
sar    $0x1d,%r10
or     $0x8,%r10
mov    -0x30(%rsp),%rdx
movsd  (%rdx,%r10,1),%xmm4
movsd  -0x8(%r11,%rdi,8),%xmm1
movsd  (%r11,%rdi,8),%xmm2
movsd  -0x8(%r15,%rcx,8),%xmm8
movsd  (%r15,%rcx,8),%xmm3
movapd %xmm2,%xmm0
addsd  %xmm3,%xmm0
subsd  %xmm3,%xmm2
movsd  (%rsi,%rdi,8),%xmm6
movsd  -0x8(%rsi,%rcx,8),%xmm9
movsd  (%rsi,%rcx,8),%xmm3
movapd %xmm6,%xmm5
subsd  %xmm3,%xmm5
addsd  %xmm3,%xmm6
movapd %xmm1,%xmm10
subsd  %xmm8,%xmm10
addsd  %xmm8,%xmm1
movapd %xmm4,%xmm7
subsd  %xmm9,%xmm7
addsd  %xmm9,%xmm4
movapd %xmm1,%xmm3
addsd  %xmm4,%xmm3
movsd  %xmm3,-0x8(%r14,%rdi,8)
subsd  %xmm4,%xmm1
movapd %xmm2,%xmm3
addsd  %xmm5,%xmm3
movsd  %xmm3,(%r14,%rdi,8)
subsd  %xmm5,%xmm2
movapd %xmm10,%xmm3
subsd  %xmm6,%xmm3
addsd  %xmm6,%xmm10
movapd %xmm0,%xmm4
addsd  %xmm7,%xmm4
subsd  %xmm7,%xmm0
movsd  -0x8(%r8,%rdi,8),%xmm5
mulsd  %xmm4,%xmm5
movsd  -0x10(%r8,%rdi,8),%xmm6
mulsd  %xmm3,%xmm6
subsd  %xmm5,%xmm6
movsd  %xmm6,-0x8(%rbp,%rdi,8)
mulsd  -0x8(%r8,%rdi,8),%xmm3
mulsd  -0x10(%r8,%rdi,8),%xmm4
addsd  %xmm3,%xmm4
movsd  %xmm4,0x0(%rbp,%rdi,8)
movsd  -0x8(%r9,%rdi,8),%xmm3
mulsd  %xmm2,%xmm3
movsd  -0x10(%r9,%rdi,8),%xmm4
mulsd  %xmm1,%xmm4
subsd  %xmm3,%xmm4
movsd  %xmm4,-0x8(%rax,%rdi,8)
mulsd  -0x8(%r9,%rdi,8),%xmm1
mulsd  -0x10(%r9,%rdi,8),%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%rax,%rdi,8)
movsd  -0x8(%r13,%rdi,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  -0x10(%r13,%rdi,8),%xmm2
mulsd  %xmm10,%xmm2
subsd  %xmm1,%xmm2
movsd  %xmm2,-0x8(%rbx,%rdi,8)
mulsd  -0x8(%r13,%rdi,8),%xmm10
mulsd  -0x10(%r13,%rdi,8),%xmm0
addsd  %xmm10,%xmm0
movsd  %xmm0,(%rbx,%rdi,8)
add    $0x2,%rdi
add    $0xfffffffffffffffe,%rcx
movabs $0x200000000,%rdx
add    %rdx,%r12
cmp    -0x8(%rsp),%rdi
jb     1290 <dradb4+0x190>
mov    -0x1c(%rsp),%edi
add    $0x1,%edi
mov    -0x24(%rsp),%r15d
mov    -0x14(%rsp),%r10d
add    %r15d,%r10d
mov    -0x48(%rsp),%r11d
add    %r15d,%r11d
mov    -0x18(%rsp),%edx
add    %r15d,%edx
mov    -0x10(%rsp),%rcx
add    %rcx,%r14
add    %rcx,%rbp
add    %rcx,%rax
add    %rcx,%rbx
mov    -0x20(%rsp),%esi
cmp    %esi,%edi
jne    1250 <dradb4+0x150>
test   %esi,%esi
mov    -0x30(%rsp),%rdx
jle    154b <dradb4+0x44b>
mov    -0x40(%rsp),%rdi
mov    %edi,%eax
shr    $0x1f,%eax
add    %edi,%eax
and    $0xfffffffe,%eax
sub    %eax,%edi
cmp    $0x1,%edi
je     154b <dradb4+0x44b>
mov    -0x40(%rsp),%rbp
lea    0x0(%rbp,%rbp,2),%edi
movslq %ebp,%rax
add    $0xffffffff,%ebp
movslq %r15d,%r8
movslq %ebp,%r9
movslq -0x44(%rsp),%r11
movslq %edi,%rbx
shl    $0x3,%r8
mov    -0x38(%rsp),%rcx
lea    (%rcx,%r9,8),%rdi
shl    $0x3,%rax
lea    0x0(,%r11,8),%rbp
lea    0x0(,%rbp,2),%r10
add    %rbp,%r10
mov    %r11,%rbp
shl    $0x4,%rbp
add    %r9,%r11
lea    (%rcx,%r11,8),%rcx
movsd  0xb35(%rip),%xmm0
movsd  0xb35(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  -0x8(%rdx,%rax,1),%xmm2
movsd  (%rdx,%rax,1),%xmm3
movsd  -0x8(%rdx,%rbx,8),%xmm4
movsd  (%rdx,%rbx,8),%xmm5
movapd %xmm3,%xmm6
addsd  %xmm5,%xmm6
subsd  %xmm3,%xmm5
movapd %xmm2,%xmm3
subsd  %xmm4,%xmm3
addsd  %xmm4,%xmm2
addsd  %xmm2,%xmm2
movsd  %xmm2,(%rdi)
movapd %xmm3,%xmm2
subsd  %xmm6,%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,(%rcx)
addsd  %xmm5,%xmm5
movsd  %xmm5,(%rdi,%rbp,1)
addsd  %xmm6,%xmm3
mulsd  %xmm1,%xmm3
movsd  %xmm3,(%rdi,%r10,1)
add    %r8,%rdx
add    %rax,%rdi
add    %rax,%rcx
add    $0xffffffff,%esi
jne    14e0 <dradb4+0x3e0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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