void dradb4(int ido, int l1, float *cc, float *ch, float *wa1, float *wa2,
            float *wa3) {
  static float sqrt2 = 1.4142135623730950488016887242097;
  int i, k, t0, t1, t2, t3, t4, t5, t6, t7, t8;
  float 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    $0x10,%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    12b3 <dradb4+0x1b3>
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
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x8c(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    -0x44(%rbp),%ecx
mov    %ecx,-0x44(%rbp)
sub    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x84(%rbp)
ff
movss  -0x84(%rbp),%xmm0
ff
addss  -0x88(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x80(%rbp),%xmm0
subss  -0x8c(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x48(%rbp),%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x84(%rbp),%xmm0
ff
subss  -0x88(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x48(%rbp),%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x80(%rbp),%xmm0
addss  -0x8c(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x48(%rbp),%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
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    12c2 <dradb4+0x1c2>
jmp    18c8 <dradb4+0x7c8>
cmpl   $0x2,-0x4(%rbp)
jne    12d1 <dradb4+0x1d1>
jmp    172e <dradb4+0x62e>
movl   $0x0,-0x38(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1710 <dradb4+0x610>
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    16f9 <dradb4+0x5f9>
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
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x74(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x8c(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x84(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x7c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x88(%rbp)
ff
movss  -0x84(%rbp),%xmm0
ff
addss  -0x88(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x84(%rbp),%xmm0
ff
subss  -0x88(%rbp),%xmm0
ff
movss  %xmm0,-0x68(%rbp)
movss  -0x74(%rbp),%xmm0
addss  -0x78(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x50(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x74(%rbp),%xmm0
subss  -0x78(%rbp),%xmm0
movss  %xmm0,-0x5c(%rbp)
movss  -0x80(%rbp),%xmm0
subss  -0x8c(%rbp),%xmm0
ff
movss  %xmm0,-0x64(%rbp)
movss  -0x80(%rbp),%xmm0
addss  -0x8c(%rbp),%xmm0
ff
movss  %xmm0,-0x6c(%rbp)
movss  -0x70(%rbp),%xmm0
addss  -0x7c(%rbp),%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x70(%rbp),%xmm0
subss  -0x7c(%rbp),%xmm0
movss  %xmm0,-0x60(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x64(%rbp),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x58(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%ecx
add    -0x34(%rbp),%ecx
mov    %ecx,-0x54(%rbp)
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x58(%rbp),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x64(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x68(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x5c(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x54(%rbp),%ecx
mov    %ecx,-0x54(%rbp)
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x5c(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x68(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x6c(%rbp),%xmm2
mov    0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x60(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x54(%rbp),%ecx
mov    %ecx,-0x54(%rbp)
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x60(%rbp),%xmm2
mov    0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x6c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x2c(%rbp)
jmp    1310 <dradb4+0x210>
mov    -0x4(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12df <dradb4+0x1df>
mov    -0x4(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x1,%edx
jne    1729 <dradb4+0x629>
jmp    18c8 <dradb4+0x7c8>
jmp    172e <dradb4+0x62e>
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    18c8 <dradb4+0x7c8>
mov    -0x40(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x74(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x84(%rbp)
ff
movss  -0x84(%rbp),%xmm0
ff
addss  -0x84(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  0x27fb(%rip),%xmm0
movss  -0x80(%rbp),%xmm1
subss  -0x70(%rbp),%xmm1
mulss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x48(%rbp),%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x74(%rbp),%xmm0
addss  -0x74(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x48(%rbp),%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  0x27b1(%rip),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  -0x80(%rbp),%xmm1
addss  -0x70(%rbp),%xmm1
mulss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x48(%rbp),%ecx
mov    %ecx,-0x48(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
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    175b <dradb4+0x65b>
add    $0x10,%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,%ecx
imul   %edi,%ecx
lea    0x0(,%rdi,4),%r13d
mov    %rdi,-0x48(%rsp)
lea    (%rdi,%rdi,1),%eax
mov    %eax,-0x4c(%rsp)
test   %esi,%esi
jle    11d6 <dradb4+0xd6>
movslq -0x48(%rsp),%r10
movslq %ecx,%r15
movslq %r13d,%r11
movslq -0x4c(%rsp),%r12
lea    -0x1(%r11),%ebp
shl    $0x2,%r11
shl    $0x2,%r10
lea    0x0(,%r15,4),%rax
lea    (%rax,%rax,2),%r14
mov    %esi,%eax
mov    -0x38(%rsp),%rdi
mov    %rdx,%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  -0x4(%rbx,%r12,4),%xmm0
movss  (%rbx,%r12,4),%xmm1
addss  %xmm0,%xmm0
addss  %xmm1,%xmm1
movss  (%rbx),%xmm2
movslq %ebp,%rbp
movss  (%rdx,%rbp,4),%xmm3
movaps %xmm2,%xmm4
subss  %xmm3,%xmm4
addss  %xmm3,%xmm2
movaps %xmm0,%xmm3
addss  %xmm2,%xmm3
movss  %xmm3,(%rdi)
movaps %xmm4,%xmm3
subss  %xmm1,%xmm3
movss  %xmm3,(%rdi,%r15,4)
subss  %xmm0,%xmm2
movss  %xmm2,(%rdi,%r15,8)
addss  %xmm1,%xmm4
movss  %xmm4,(%rdi,%r14,1)
add    %r13d,%ebp
add    %r11,%rbx
add    %r10,%rdi
add    $0xffffffff,%eax
jne    1170 <dradb4+0x70>
cmpl   $0x2,-0x48(%rsp)
jl     152b <dradb4+0x42b>
mov    %r13d,%r15d
jne    11f3 <dradb4+0xf3>
test   %esi,%esi
jg     1470 <dradb4+0x370>
jmp    152b <dradb4+0x42b>
mov    %ecx,-0x3c(%rsp)
mov    %rdx,-0x30(%rsp)
test   %esi,%esi
mov    -0x4c(%rsp),%r11d
jle    1445 <dradb4+0x345>
mov    0x38(%rsp),%r13
mov    -0x48(%rsp),%eax
movslq -0x3c(%rsp),%rcx
mov    %rax,-0x8(%rsp)
lea    0x0(,%rax,4),%rax
mov    %rax,-0x10(%rsp)
mov    -0x38(%rsp),%r14
lea    (%r14,%rcx,4),%rbp
lea    (%r14,%rcx,8),%rax
lea    (%rcx,%rcx,2),%rcx
lea    (%r14,%rcx,4),%rbx
xor    %r10d,%r10d
mov    %r15d,%edx
xor    %edi,%edi
mov    %esi,-0x20(%rsp)
mov    %r15d,-0x24(%rsp)
nop
mov    %edi,-0x1c(%rsp)
mov    %edx,-0x18(%rsp)
movslq %edx,%rcx
mov    -0x30(%rsp),%rdx
lea    (%rdx,%rcx,4),%r15
mov    %r11d,-0x4c(%rsp)
movslq %r11d,%r12
lea    (%rdx,%r12,4),%rsi
shl    $0x20,%r12
mov    %r10d,-0x14(%rsp)
movslq %r10d,%rcx
lea    (%rdx,%rcx,4),%r11
mov    $0x2,%edi
mov    $0xfffffffffffffffe,%rcx
nopl   0x0(%rax)
mov    %r12,%r10
sar    $0x1e,%r10
or     $0x4,%r10
mov    -0x30(%rsp),%rdx
movss  (%rdx,%r10,1),%xmm4
movss  -0x4(%r11,%rdi,4),%xmm1
movss  (%r11,%rdi,4),%xmm2
movss  -0x4(%r15,%rcx,4),%xmm8
movss  (%r15,%rcx,4),%xmm3
movaps %xmm2,%xmm0
addss  %xmm3,%xmm0
subss  %xmm3,%xmm2
movss  (%rsi,%rdi,4),%xmm6
movss  -0x4(%rsi,%rcx,4),%xmm9
movss  (%rsi,%rcx,4),%xmm3
movaps %xmm6,%xmm5
subss  %xmm3,%xmm5
addss  %xmm3,%xmm6
movaps %xmm1,%xmm10
subss  %xmm8,%xmm10
addss  %xmm8,%xmm1
movaps %xmm4,%xmm7
subss  %xmm9,%xmm7
addss  %xmm9,%xmm4
movaps %xmm1,%xmm3
addss  %xmm4,%xmm3
movss  %xmm3,-0x4(%r14,%rdi,4)
subss  %xmm4,%xmm1
movaps %xmm2,%xmm3
addss  %xmm5,%xmm3
movss  %xmm3,(%r14,%rdi,4)
subss  %xmm5,%xmm2
movaps %xmm10,%xmm3
subss  %xmm6,%xmm3
addss  %xmm6,%xmm10
movaps %xmm0,%xmm4
addss  %xmm7,%xmm4
subss  %xmm7,%xmm0
movss  -0x4(%r8,%rdi,4),%xmm5
mulss  %xmm4,%xmm5
movss  -0x8(%r8,%rdi,4),%xmm6
mulss  %xmm3,%xmm6
subss  %xmm5,%xmm6
movss  %xmm6,-0x4(%rbp,%rdi,4)
mulss  -0x4(%r8,%rdi,4),%xmm3
mulss  -0x8(%r8,%rdi,4),%xmm4
addss  %xmm3,%xmm4
movss  %xmm4,0x0(%rbp,%rdi,4)
movss  -0x4(%r9,%rdi,4),%xmm3
mulss  %xmm2,%xmm3
movss  -0x8(%r9,%rdi,4),%xmm4
mulss  %xmm1,%xmm4
subss  %xmm3,%xmm4
movss  %xmm4,-0x4(%rax,%rdi,4)
mulss  -0x4(%r9,%rdi,4),%xmm1
mulss  -0x8(%r9,%rdi,4),%xmm2
addss  %xmm1,%xmm2
movss  %xmm2,(%rax,%rdi,4)
movss  -0x4(%r13,%rdi,4),%xmm1
mulss  %xmm0,%xmm1
movss  -0x8(%r13,%rdi,4),%xmm2
mulss  %xmm10,%xmm2
subss  %xmm1,%xmm2
movss  %xmm2,-0x4(%rbx,%rdi,4)
mulss  -0x4(%r13,%rdi,4),%xmm10
mulss  -0x8(%r13,%rdi,4),%xmm0
addss  %xmm10,%xmm0
movss  %xmm0,(%rbx,%rdi,4)
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    -0x4c(%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
mov    -0x3c(%rsp),%ecx
jle    152b <dradb4+0x42b>
mov    -0x48(%rsp),%rdi
mov    %edi,%eax
shr    $0x1f,%eax
add    %edi,%eax
and    $0xfffffffe,%eax
sub    %eax,%edi
cmp    $0x1,%edi
je     152b <dradb4+0x42b>
mov    -0x48(%rsp),%rbp
lea    0x0(%rbp,%rbp,2),%edi
movslq %ebp,%rax
add    $0xffffffff,%ebp
movslq %r15d,%r8
movslq %ebp,%r9
movslq %ecx,%r10
movslq %edi,%rdi
shl    $0x2,%r8
mov    -0x38(%rsp),%rcx
lea    (%rcx,%r9,4),%rbp
shl    $0x2,%rax
lea    0x0(,%r10,4),%rbx
lea    (%rbx,%rbx,2),%rbx
add    %r10,%r9
lea    (%rcx,%r9,4),%rcx
movss  0xb49(%rip),%xmm0
movss  0xb45(%rip),%xmm1
nop
movss  -0x4(%rdx,%rax,1),%xmm2
movss  (%rdx,%rax,1),%xmm3
movss  -0x4(%rdx,%rdi,4),%xmm4
movss  (%rdx,%rdi,4),%xmm5
movaps %xmm3,%xmm6
addss  %xmm5,%xmm6
subss  %xmm3,%xmm5
movaps %xmm2,%xmm3
subss  %xmm4,%xmm3
addss  %xmm4,%xmm2
addss  %xmm2,%xmm2
movss  %xmm2,0x0(%rbp)
movaps %xmm3,%xmm2
subss  %xmm6,%xmm2
mulss  %xmm0,%xmm2
movss  %xmm2,(%rcx)
addss  %xmm5,%xmm5
movss  %xmm5,0x0(%rbp,%r10,8)
addss  %xmm6,%xmm3
mulss  %xmm1,%xmm3
movss  %xmm3,0x0(%rbp,%rbx,1)
add    %r8,%rdx
add    %rax,%rbp
add    %rax,%rcx
add    $0xffffffff,%esi
jne    14c0 <dradb4+0x3c0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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