void dradb3(int ido, int l1, float *cc, float *ch, float *wa1, float *wa2) {
  static float taur = -.5f;
  static float taui = .8660254037844386f;
  int i, k, t0, t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;
  float ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
  t0 = l1 * ido;
  t1 = 0;
  t2 = t0 << 1;
  t3 = ido << 1;
  t4 = ido + (ido << 1);
  t5 = 0;
  for (k = 0; k < l1; k++) {
    tr2 = cc[t3 - 1] + cc[t3 - 1];
    cr2 = cc[t5] + (taur * tr2);
    ch[t1] = cc[t5] + tr2;
    ci3 = taui * (cc[t3] + cc[t3]);
    ch[t1 + t0] = cr2 - ci3;
    ch[t1 + t2] = cr2 + ci3;
    t1 += ido;
    t3 += t4;
    t5 += t4;
  }
  if (ido == 1)
    return;
  t1 = 0;
  t3 = ido << 1;
  for (k = 0; k < l1; k++) {
    t7 = t1 + (t1 << 1);
    t6 = (t5 = t7 + t3);
    t8 = t1;
    t10 = (t9 = t1 + t0) + t0;
    for (i = 2; i < ido; i += 2) {
      t5 += 2;
      t6 -= 2;
      t7 += 2;
      t8 += 2;
      t9 += 2;
      t10 += 2;
      tr2 = cc[t5 - 1] + cc[t6 - 1];
      cr2 = cc[t7 - 1] + (taur * tr2);
      ch[t8 - 1] = cc[t7 - 1] + tr2;
      ti2 = cc[t5] - cc[t6];
      ci2 = cc[t7] + (taur * ti2);
      ch[t8] = cc[t7] + ti2;
      cr3 = taui * (cc[t5 - 1] - cc[t6 - 1]);
      ci3 = taui * (cc[t5] + cc[t6]);
      dr2 = cr2 - ci3;
      dr3 = cr2 + ci3;
      di2 = ci2 + cr3;
      di3 = ci2 - cr3;
      ch[t9 - 1] = wa1[i - 2] * dr2 - wa1[i - 1] * di2;
      ch[t9] = wa1[i - 2] * di2 + wa1[i - 1] * dr2;
      ch[t10 - 1] = wa2[i - 2] * dr3 - wa2[i - 1] * di3;
      ch[t10] = wa2[i - 2] * di3 + wa2[i - 1] * dr3;
    }
    t1 += ido;
  }
}


x86 assembly (clang, -O0):
<dradb3>:
push   %rbp
mov    %rsp,%rbp
push   %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    -0x34(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x4(%rbp),%eax
shl    $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,-0x48(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    126a <dradb3+0x16a>
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x84(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movss  0x2e79(%rip),%xmm0
movss  -0x84(%rbp),%xmm2
ff
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x84(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  0x2e3e(%rip),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x64(%rbp)
movss  -0x70(%rbp),%xmm0
subss  -0x64(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x70(%rbp),%xmm0
addss  -0x64(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x4(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x44(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x44(%rbp),%eax
add    -0x48(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    115a <dradb3+0x5a>
cmpl   $0x1,-0x4(%rbp)
jne    1279 <dradb3+0x179>
jmp    15f1 <dradb3+0x4f1>
movl   $0x0,-0x38(%rbp)
mov    -0x4(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x40(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    15f1 <dradb3+0x4f1>
mov    -0x38(%rbp),%eax
mov    -0x38(%rbp),%ecx
shl    $0x1,%ecx
add    %ecx,%eax
mov    %eax,-0x50(%rbp)
mov    -0x50(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    %eax,-0x4c(%rbp)
mov    -0x38(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x38(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x58(%rbp)
add    -0x34(%rbp),%eax
mov    %eax,-0x5c(%rbp)
movl   $0x2,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    15da <dradb3+0x4da>
mov    -0x48(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x48(%rbp)
mov    -0x4c(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x50(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x50(%rbp)
mov    -0x54(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x54(%rbp)
mov    -0x58(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x58(%rbp)
mov    -0x5c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x5c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x84(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movss  0x2cc6(%rip),%xmm0
movss  -0x84(%rbp),%xmm2
ff
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x84(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movss  0x2c51(%rip),%xmm0
movss  -0x80(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x60(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x80(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  0x2c1c(%rip),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x74(%rbp)
movss  0x2be7(%rip),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,-0x64(%rbp)
movss  -0x70(%rbp),%xmm0
subss  -0x64(%rbp),%xmm0
movss  %xmm0,-0x78(%rbp)
movss  -0x70(%rbp),%xmm0
addss  -0x64(%rbp),%xmm0
movss  %xmm0,-0x7c(%rbp)
movss  -0x60(%rbp),%xmm0
addss  -0x74(%rbp),%xmm0
movss  %xmm0,-0x68(%rbp)
movss  -0x60(%rbp),%xmm0
subss  -0x74(%rbp),%xmm0
movss  %xmm0,-0x6c(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x78(%rbp),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x68(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x58(%rbp),%ecx
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  -0x68(%rbp),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x78(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x58(%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  -0x7c(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x6c(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x5c(%rbp),%ecx
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  -0x6c(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x7c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x2c(%rbp)
jmp    12d2 <dradb3+0x1d2>
mov    -0x4(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1290 <dradb3+0x190>
add    $0x8,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dradb3>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rcx,-0x30(%rsp)
mov    %rdx,-0x38(%rsp)
mov    %esi,-0x3c(%rsp)
test   %esi,%esi
jle    13aa <dradb3+0x2aa>
mov    -0x3c(%rsp),%ebp
mov    %ebp,%edx
imul   %edi,%edx
lea    (%rdx,%rdx,1),%eax
lea    (%rdi,%rdi,1),%ecx
lea    (%rdi,%rdi,2),%esi
mov    %rdi,-0x18(%rsp)
movslq %edi,%rbx
mov    %eax,%r11d
movslq %eax,%r10
mov    %rdx,-0x28(%rsp)
movslq %edx,%r14
mov    %ecx,%r13d
movslq %ecx,%rdx
mov    %esi,-0x1c(%rsp)
movslq %esi,%rsi
shl    $0x2,%rsi
lea    0x0(,%rbx,4),%rdi
movss  0xe97(%rip),%xmm0
movss  0xe93(%rip),%xmm1
mov    -0x30(%rsp),%rax
mov    -0x38(%rsp),%rcx
nopl   0x0(%rax,%rax,1)
movss  -0x4(%rcx,%rdx,4),%xmm2
addss  %xmm2,%xmm2
movss  (%rcx),%xmm3
movaps %xmm2,%xmm4
mulss  %xmm0,%xmm4
addss  %xmm3,%xmm4
addss  %xmm3,%xmm2
movss  %xmm2,(%rax)
movss  (%rcx,%rdx,4),%xmm2
addss  %xmm2,%xmm2
mulss  %xmm1,%xmm2
movaps %xmm4,%xmm3
subss  %xmm2,%xmm3
movss  %xmm3,(%rax,%r14,4)
addss  %xmm4,%xmm2
movss  %xmm2,(%rax,%r10,4)
add    %rsi,%rcx
add    %rdi,%rax
add    $0xffffffff,%ebp
jne    1180 <dradb3+0x80>
mov    -0x18(%rsp),%rcx
cmp    $0x1,%ecx
je     13aa <dradb3+0x2aa>
cmpl   $0x0,-0x3c(%rsp)
jle    13aa <dradb3+0x2aa>
mov    %r11d,%ebp
mov    -0x38(%rsp),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x10(%rsp)
xor    %r11d,%r11d
mov    -0x30(%rsp),%rax
xor    %r15d,%r15d
mov    %rdi,-0x8(%rsp)
jmp    1242 <dradb3+0x142>
nopl   0x0(%rax)
add    $0x1,%r15d
mov    -0x1c(%rsp),%ecx
add    %ecx,%r13d
add    %ecx,%r11d
add    %edi,%ebp
mov    -0x28(%rsp),%rcx
add    %edi,%ecx
mov    %rcx,-0x28(%rsp)
mov    %rdi,%rcx
mov    -0x8(%rsp),%rdi
add    %rdi,%rax
cmp    -0x3c(%rsp),%r15d
je     13aa <dradb3+0x2aa>
cmp    $0x3,%ecx
mov    %rcx,%rdi
jl     1210 <dradb3+0x110>
movslq -0x28(%rsp),%rcx
mov    -0x30(%rsp),%rdx
lea    (%rdx,%rcx,4),%r14
movslq %ebp,%rcx
lea    (%rdx,%rcx,4),%rcx
movslq %r11d,%rdx
mov    %ebp,%r10d
mov    -0x38(%rsp),%rbp
lea    0x0(%rbp,%rdx,4),%rdx
movslq %r13d,%rdi
mov    -0x10(%rsp),%rsi
lea    (%rsi,%rdi,4),%rsi
lea    0x0(,%rdi,4),%r12
add    %rbp,%r12
mov    %r10d,%ebp
mov    -0x18(%rsp),%rdi
mov    $0x2,%r10d
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  -0x4(%r12,%r10,4),%xmm2
addss  -0x4(%rsi),%xmm2
movss  -0x4(%rdx,%r10,4),%xmm3
movaps %xmm2,%xmm4
mulss  %xmm0,%xmm4
addss  %xmm3,%xmm4
addss  %xmm3,%xmm2
movss  %xmm2,-0x4(%rax,%r10,4)
movss  (%r12,%r10,4),%xmm3
subss  (%rsi),%xmm3
movss  (%rdx,%r10,4),%xmm5
movaps %xmm3,%xmm2
mulss  %xmm0,%xmm2
addss  %xmm5,%xmm2
addss  %xmm5,%xmm3
movss  %xmm3,(%rax,%r10,4)
movss  -0x4(%r12,%r10,4),%xmm5
movss  (%r12,%r10,4),%xmm3
subss  -0x4(%rsi),%xmm5
mulss  %xmm1,%xmm5
addss  (%rsi),%xmm3
mulss  %xmm1,%xmm3
movaps %xmm4,%xmm6
subss  %xmm3,%xmm6
addss  %xmm4,%xmm3
movaps %xmm2,%xmm4
addss  %xmm5,%xmm4
subss  %xmm5,%xmm2
movss  -0x4(%r8,%r10,4),%xmm5
mulss  %xmm4,%xmm5
movss  -0x8(%r8,%r10,4),%xmm7
mulss  %xmm6,%xmm7
subss  %xmm5,%xmm7
movss  %xmm7,-0x4(%r14,%r10,4)
mulss  -0x4(%r8,%r10,4),%xmm6
mulss  -0x8(%r8,%r10,4),%xmm4
addss  %xmm6,%xmm4
movss  %xmm4,(%r14,%r10,4)
movss  -0x4(%r9,%r10,4),%xmm4
mulss  %xmm2,%xmm4
movss  -0x8(%r9,%r10,4),%xmm5
mulss  %xmm3,%xmm5
subss  %xmm4,%xmm5
movss  %xmm5,-0x4(%rcx,%r10,4)
mulss  -0x4(%r9,%r10,4),%xmm3
mulss  -0x8(%r9,%r10,4),%xmm2
addss  %xmm3,%xmm2
movss  %xmm2,(%rcx,%r10,4)
add    $0x2,%r10
add    $0xfffffffffffffff8,%rsi
cmp    %rbx,%r10
jl     12a0 <dradb3+0x1a0>
jmp    1210 <dradb3+0x110>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @dradb3(i32 %0, i32 %1, float* nocapture %2, float* nocapture writeonly %3, float* nocapture %4, float* nocapture %5) {
%7 = mul i32 %1, %0
%8 = shl i32 %7, 1
%9 = shl i32 %0, 1
%10 = mul i32 %0, 3
%11 = icmp sgt i32 %1, 0
br i1 %11, label %12, label %155
12: ; preds = %6
%13 = sext i32 %0 to i64
%14 = sext i32 %8 to i64
%15 = sext i32 %7 to i64
%16 = sext i32 %9 to i64
%17 = sext i32 %10 to i64
br label %18
18: ; preds = %12, %18
%19 = phi i64 [ 0, %12 ], [ %44, %18 ]
%20 = phi i64 [ %16, %12 ], [ %43, %18 ]
%21 = phi i64 [ 0, %12 ], [ %42, %18 ]
%22 = phi i32 [ 0, %12 ], [ %45, %18 ]
%23 = add nsw i64 %20, -1
%24 = getelementptr inbounds float, float* %2, i64 %23
%25 = load float, float* %24
%26 = fadd float %25, %25
%27 = getelementptr inbounds float, float* %2, i64 %19
%28 = load float, float* %27
%29 = tail call float @llvm.fmuladd.f32(float %26, float -5.000000e-01, float %28)
%30 = fadd float %28, %26
%31 = getelementptr inbounds float, float* %3, i64 %21
store float %30, float* %31
%32 = getelementptr inbounds float, float* %2, i64 %20
%33 = load float, float* %32
%34 = fadd float %33, %33
%35 = fmul float %34, 0x3FEBB67AE0000000
%36 = fsub float %29, %35
%37 = add nsw i64 %21, %15
%38 = getelementptr inbounds float, float* %3, i64 %37
store float %36, float* %38
%39 = fadd float %29, %35
%40 = add nsw i64 %21, %14
%41 = getelementptr inbounds float, float* %3, i64 %40
store float %39, float* %41
%42 = add nsw i64 %21, %13
%43 = add i64 %20, %17
%44 = add nsw i64 %19, %17
%45 = add nuw nsw i32 %22, 1
%46 = icmp eq i32 %45, %1
br i1 %46, label %47, label %18
47: ; preds = %18
%48 = icmp ne i32 %0, 1
%49 = and i1 %48, %11
br i1 %49, label %50, label %155
50: ; preds = %47
%51 = icmp sgt i32 %0, 2
%52 = sext i32 %0 to i64
br label %53
53: ; preds = %50, %147
%54 = phi i64 [ 0, %50 ], [ %148, %147 ]
%55 = phi i32 [ %7, %50 ], [ %153, %147 ]
%56 = phi i32 [ %8, %50 ], [ %152, %147 ]
%57 = phi i32 [ 0, %50 ], [ %151, %147 ]
%58 = phi i32 [ %9, %50 ], [ %150, %147 ]
%59 = phi i32 [ 0, %50 ], [ %149, %147 ]
%60 = sext i32 %58 to i64
br i1 %51, label %61, label %147
61: ; preds = %53
%62 = sext i32 %55 to i64
%63 = sext i32 %56 to i64
%64 = sext i32 %57 to i64
br label %65
65: ; preds = %61, %65
%66 = phi i64 [ %54, %61 ], [ %76, %65 ]
%67 = phi i64 [ %62, %61 ], [ %77, %65 ]
%68 = phi i64 [ %63, %61 ], [ %78, %65 ]
%69 = phi i64 [ %64, %61 ], [ %75, %65 ]
%70 = phi i64 [ %60, %61 ], [ %74, %65 ]
%71 = phi i64 [ %60, %61 ], [ %73, %65 ]
%72 = phi i64 [ 2, %61 ], [ %145, %65 ]
%73 = add nsw i64 %71, 2
%74 = add nsw i64 %70, -2
%75 = add nsw i64 %69, 2
%76 = add nsw i64 %66, 2
%77 = add nsw i64 %67, 2
%78 = add nsw i64 %68, 2
%79 = add nsw i64 %71, 1
%80 = getelementptr inbounds float, float* %2, i64 %79
%81 = load float, float* %80
%82 = add nsw i64 %70, -3
%83 = getelementptr inbounds float, float* %2, i64 %82
%84 = load float, float* %83
%85 = fadd float %81, %84
%86 = add nsw i64 %69, 1
%87 = getelementptr inbounds float, float* %2, i64 %86
%88 = load float, float* %87
%89 = tail call float @llvm.fmuladd.f32(float %85, float -5.000000e-01, float %88)
%90 = fadd float %88, %85
%91 = add nsw i64 %66, 1
%92 = getelementptr inbounds float, float* %3, i64 %91
store float %90, float* %92
%93 = getelementptr inbounds float, float* %2, i64 %73
%94 = load float, float* %93
%95 = getelementptr inbounds float, float* %2, i64 %74
%96 = load float, float* %95
%97 = fsub float %94, %96
%98 = getelementptr inbounds float, float* %2, i64 %75
%99 = load float, float* %98
%100 = tail call float @llvm.fmuladd.f32(float %97, float -5.000000e-01, float %99)
%101 = fadd float %99, %97
%102 = getelementptr inbounds float, float* %3, i64 %76
store float %101, float* %102
%103 = load float, float* %80
%104 = load float, float* %83
%105 = fsub float %103, %104
%106 = fmul float %105, 0x3FEBB67AE0000000
%107 = load float, float* %93
%108 = load float, float* %95
%109 = fadd float %107, %108
%110 = fmul float %109, 0x3FEBB67AE0000000
%111 = fsub float %89, %110
%112 = fadd float %89, %110
%113 = fadd float %100, %106
%114 = fsub float %100, %106
%115 = add nsw i64 %72, -2
%116 = getelementptr inbounds float, float* %4, i64 %115
%117 = load float, float* %116
%118 = add nsw i64 %72, -1
%119 = getelementptr inbounds float, float* %4, i64 %118
%120 = load float, float* %119
%121 = fneg float %120
%122 = fmul float %113, %121
%123 = tail call float @llvm.fmuladd.f32(float %117, float %111, float %122)
%124 = add nsw i64 %67, 1
%125 = getelementptr inbounds float, float* %3, i64 %124
store float %123, float* %125
%126 = load float, float* %116
%127 = load float, float* %119
%128 = fmul float %111, %127
%129 = tail call float @llvm.fmuladd.f32(float %126, float %113, float %128)
%130 = getelementptr inbounds float, float* %3, i64 %77
store float %129, float* %130
%131 = getelementptr inbounds float, float* %5, i64 %115
%132 = load float, float* %131
%133 = getelementptr inbounds float, float* %5, i64 %118
%134 = load float, float* %133
%135 = fneg float %134
%136 = fmul float %114, %135
%137 = tail call float @llvm.fmuladd.f32(float %132, float %112, float %136)
%138 = add nsw i64 %68, 1
%139 = getelementptr inbounds float, float* %3, i64 %138
store float %137, float* %139
%140 = load float, float* %131
%141 = load float, float* %133
%142 = fmul float %112, %141
%143 = tail call float @llvm.fmuladd.f32(float %140, float %114, float %142)
%144 = getelementptr inbounds float, float* %3, i64 %78
store float %143, float* %144
%145 = add nuw nsw i64 %72, 2
%146 = icmp slt i64 %145, %52
br i1 %146, label %65, label %147
147: ; preds = %65, %53
%148 = add nsw i64 %54, %52
%149 = add nuw nsw i32 %59, 1
%150 = add i32 %58, %10
%151 = add i32 %57, %10
%152 = add i32 %56, %0
%153 = add i32 %55, %0
%154 = icmp eq i32 %149, %1
br i1 %154, label %155, label %53
155: ; preds = %147, %6, %47
ret void
}