void radb3(int ido, int l1, const float cc[], float ch[], const float wa1[],
           const float wa2[]) {
  static const float taur = -0.5;
  static const float taui = 0.866025403784439;
  int i, k, ic;
  float ci2, ci3, di2, di3, cr2, cr3, dr2, dr3, ti2, tr2;
  for (k = 0; k < l1; k++) {
    tr2 = 2 * cc[ido - 1 + (3 * k + 1) * ido];
    cr2 = cc[3 * k * ido] + taur * tr2;
    ch[k * ido] = cc[3 * k * ido] + tr2;
    ci3 = 2 * taui * cc[(3 * k + 2) * ido];
    ch[(k + l1) * ido] = cr2 - ci3;
    ch[(k + 2 * l1) * ido] = cr2 + ci3;
  }
  if (ido == 1)
    return;
  for (k = 0; k < l1; k++) {
    for (i = 2; i < ido; i += 2) {
      ic = ido - i;
      tr2 = cc[i - 1 + (3 * k + 2) * ido] + cc[ic - 1 + (3 * k + 1) * ido];
      cr2 = cc[i - 1 + 3 * k * ido] + taur * tr2;
      ch[i - 1 + k * ido] = cc[i - 1 + 3 * k * ido] + tr2;
      ti2 = cc[i + (3 * k + 2) * ido] - cc[ic + (3 * k + 1) * ido];
      ci2 = cc[i + 3 * k * ido] + taur * ti2;
      ch[i + k * ido] = cc[i + 3 * k * ido] + ti2;
      cr3 = taui *
            (cc[i - 1 + (3 * k + 2) * ido] - cc[ic - 1 + (3 * k + 1) * ido]);
      ci3 = taui * (cc[i + (3 * k + 2) * ido] + cc[ic + (3 * k + 1) * ido]);
      dr2 = cr2 - ci3;
      dr3 = cr2 + ci3;
      di2 = ci2 + cr3;
      di3 = ci2 - cr3;
      ch[i - 1 + (k + l1) * ido] = wa1[i - 2] * dr2 - wa1[i - 1] * di2;
      ch[i + (k + l1) * ido] = wa1[i - 2] * di2 + wa1[i - 1] * dr2;
      ch[i - 1 + (k + 2 * l1) * ido] = wa2[i - 2] * dr3 - wa2[i - 1] * di3;
      ch[i + (k + 2 * l1) * ido] = wa2[i - 2] * di3 + wa2[i - 1] * dr3;
    }
  }
}


x86 assembly (clang, -O0):
<radb3>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    122a <radb3+0x12a>
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
imul   $0x3,-0x30(%rbp),%edx
add    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  0xebd(%rip),%xmm0
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x5c(%rbp)
mov    -0x10(%rbp),%rax
imul   $0x3,-0x30(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movss  -0x5c(%rbp),%xmm0
movss  0xe8a(%rip),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
imul   $0x3,-0x30(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x5c(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
imul   $0x3,-0x30(%rbp),%ecx
add    $0x2,%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movss  0xe3b(%rip),%xmm0
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x48(%rbp),%xmm0
subss  -0x3c(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    -0x8(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x48(%rbp),%xmm0
addss  -0x3c(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
mov    -0x8(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1121 <radb3+0x21>
cmpl   $0x1,-0x4(%rbp)
jne    1239 <radb3+0x139>
jmp    1626 <radb3+0x526>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1626 <radb3+0x526>
movl   $0x2,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1613 <radb3+0x513>
mov    -0x4(%rbp),%eax
sub    -0x2c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
imul   $0x3,-0x30(%rbp),%edx
add    $0x2,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
imul   $0x3,-0x30(%rbp),%edx
add    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x5c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
imul   $0x3,-0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movss  -0x5c(%rbp),%xmm0
movss  0xd30(%rip),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
imul   $0x3,-0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x5c(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   $0x3,-0x30(%rbp),%edx
add    $0x2,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   $0x3,-0x30(%rbp),%edx
add    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   $0x3,-0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movss  -0x58(%rbp),%xmm0
movss  0xc84(%rip),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   $0x3,-0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x58(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
imul   $0x3,-0x30(%rbp),%edx
add    $0x2,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
imul   $0x3,-0x30(%rbp),%edx
add    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm1
movss  0xbf7(%rip),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,-0x4c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   $0x3,-0x30(%rbp),%edx
add    $0x2,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   $0x3,-0x30(%rbp),%edx
add    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
movss  0xbae(%rip),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x48(%rbp),%xmm0
subss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x48(%rbp),%xmm0
addss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x38(%rbp),%xmm0
addss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x38(%rbp),%xmm0
subss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  -0x50(%rbp),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x40(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x30(%rbp),%edx
add    -0x8(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%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  -0x40(%rbp),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x50(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
add    -0x8(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%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  -0x54(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x44(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x30(%rbp),%edx
mov    -0x8(%rbp),%esi
shl    $0x1,%esi
add    %esi,%edx
imul   -0x4(%rbp),%edx
add    %edx,%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  -0x44(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  -0x54(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
mov    -0x8(%rbp),%esi
shl    $0x1,%esi
add    %esi,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x2c(%rbp)
jmp    1253 <radb3+0x153>
jmp    1618 <radb3+0x518>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1240 <radb3+0x140>
pop    %rbp
ret

x86 assembly (clang, -O2):
<radb3>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %edi,-0x44(%rsp)
test   %esi,%esi
jle    13f1 <radb3+0x2f1>
lea    (%rsi,%rsi,1),%eax
movslq -0x44(%rsp),%rbp
movslq %eax,%rdi
mov    %esi,%ebx
mov    %rbp,%r11
mov    %rdi,%rax
mov    %rdi,-0x30(%rsp)
imul   %rdi,%r11
lea    0x0(,%rbp,4),%r10
mov    %rbp,%r12
imul   %rbx,%r12
lea    (%r10,%r10,2),%r14
lea    -0x1(,%rbp,1),%r15
ff
add    %rbp,%r15
movss  0xea8(%rip),%xmm0
movss  0xea4(%rip),%xmm1
mov    %rbx,%rax
mov    %rbx,-0x38(%rsp)
mov    %rbx,%rdi
mov    %rdx,%rax
mov    %rcx,%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  (%rax,%r15,4),%xmm2
addss  %xmm2,%xmm2
movss  (%rax),%xmm3
movaps %xmm2,%xmm4
mulss  %xmm0,%xmm4
addss  %xmm3,%xmm4
addss  %xmm3,%xmm2
movss  %xmm2,(%rbx)
movss  (%rax,%rbp,8),%xmm2
mulss  %xmm1,%xmm2
movaps %xmm4,%xmm3
subss  %xmm2,%xmm3
movss  %xmm3,(%rbx,%r12,4)
addss  %xmm4,%xmm2
movss  %xmm2,(%rbx,%r11,4)
add    %r10,%rbx
add    %r14,%rax
add    $0xffffffffffffffff,%rdi
jne    1180 <radb3+0x80>
cmpl   $0x1,-0x44(%rsp)
je     13f1 <radb3+0x2f1>
test   %esi,%esi
jle    13f1 <radb3+0x2f1>
xor    %eax,%eax
mov    %rax,-0x40(%rsp)
movabs $0xffffffff00000000,%rbx
ff ff ff
movss  0xe0f(%rip),%xmm1
jmp    121c <radb3+0x11c>
nopl   0x0(%rax,%rax,1)
mov    -0x40(%rsp),%rsi
add    $0x1,%rsi
mov    %rsi,%rax
mov    %rsi,-0x40(%rsp)
cmp    -0x38(%rsp),%rsi
je     13f1 <radb3+0x2f1>
cmpl   $0x3,-0x44(%rsp)
jl     1200 <radb3+0x100>
mov    -0x40(%rsp),%rdi
lea    (%rdi,%rdi,2),%rax
lea    (%rdi,%rdi,2),%rsi
add    $0x2,%rsi
imul   %rbp,%rsi
mov    %rsi,-0x10(%rsp)
lea    (%rdi,%rdi,2),%rsi
add    $0x1,%rsi
imul   %rbp,%rsi
mov    %rsi,-0x18(%rsp)
imul   %rbp,%rax
mov    %rax,-0x8(%rsp)
mov    %rdi,%rax
imul   %rbp,%rax
mov    %rax,-0x20(%rsp)
mov    -0x38(%rsp),%rsi
lea    (%rdi,%rsi,1),%rax
imul   %rbp,%rax
mov    %rax,-0x28(%rsp)
mov    -0x30(%rsp),%rsi
lea    (%rdi,%rsi,1),%r10
imul   %rbp,%r10
mov    $0x2,%r11d
nopl   0x0(%rax,%rax,1)
mov    %rbp,%rdi
sub    %r11,%rdi
mov    -0x10(%rsp),%r12
lea    (%r11,%r12,1),%r14
add    $0xffffffffffffffff,%r14
movss  (%rdx,%r14,4),%xmm2
add    -0x18(%rsp),%rdi
mov    %rdi,%r13
shl    $0x20,%r13
add    %rbx,%r13
sar    $0x1e,%r13
addss  (%rdx,%r13,1),%xmm2
mov    -0x8(%rsp),%rsi
lea    (%r11,%rsi,1),%r15
add    $0xffffffffffffffff,%r15
movss  (%rdx,%r15,4),%xmm5
movaps %xmm2,%xmm4
addss  %xmm5,%xmm2
mov    %rbx,%r15
mov    -0x20(%rsp),%rax
lea    -0x1(%r11,%rax,1),%rbx
movss  %xmm2,(%rcx,%rbx,4)
lea    (%r11,%r12,1),%rbx
movss  (%rdx,%rbx,4),%xmm3
subss  (%rdx,%rdi,4),%xmm3
lea    (%r11,%rsi,1),%rsi
movss  (%rdx,%rsi,4),%xmm6
movaps %xmm3,%xmm2
addss  %xmm6,%xmm3
lea    (%r11,%rax,1),%rsi
movss  %xmm3,(%rcx,%rsi,4)
movss  (%rdx,%r14,4),%xmm7
subss  (%rdx,%r13,1),%xmm7
movss  (%rdx,%rbx,4),%xmm3
mov    %r15,%rbx
mulss  %xmm0,%xmm4
addss  %xmm5,%xmm4
mulss  %xmm0,%xmm2
addss  %xmm6,%xmm2
mulss  %xmm1,%xmm7
addss  (%rdx,%rdi,4),%xmm3
mulss  %xmm1,%xmm3
movaps %xmm4,%xmm5
subss  %xmm3,%xmm5
addss  %xmm4,%xmm3
movaps %xmm2,%xmm4
addss  %xmm7,%xmm4
subss  %xmm7,%xmm2
movss  -0x4(%r8,%r11,4),%xmm6
mulss  %xmm4,%xmm6
movss  -0x8(%r8,%r11,4),%xmm7
mulss  %xmm5,%xmm7
subss  %xmm6,%xmm7
mov    -0x28(%rsp),%rax
lea    -0x1(%r11,%rax,1),%rsi
movss  %xmm7,(%rcx,%rsi,4)
mulss  -0x4(%r8,%r11,4),%xmm5
mulss  -0x8(%r8,%r11,4),%xmm4
addss  %xmm5,%xmm4
lea    (%r11,%rax,1),%rsi
movss  %xmm4,(%rcx,%rsi,4)
movss  -0x4(%r9,%r11,4),%xmm4
mulss  %xmm2,%xmm4
movss  -0x8(%r9,%r11,4),%xmm5
mulss  %xmm3,%xmm5
subss  %xmm4,%xmm5
lea    -0x1(%r11,%r10,1),%rsi
movss  %xmm5,(%rcx,%rsi,4)
mulss  -0x4(%r9,%r11,4),%xmm3
mulss  -0x8(%r9,%r11,4),%xmm2
addss  %xmm3,%xmm2
lea    (%r11,%r10,1),%rsi
movss  %xmm2,(%rcx,%rsi,4)
add    $0x2,%r11
cmp    %rbp,%r11
jl     1290 <radb3+0x190>
jmp    1200 <radb3+0x100>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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