void radb3(int ido, int l1, const double cc[], double ch[], const double wa1[],
           const double wa2[]) {
  static const double taur = -0.5;
  static const double taui = 0.866025403784439;
  int i, k, ic;
  double 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
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)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1234 <radb3+0x134>
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
movsd  0xec8(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
imul   $0x3,-0x30(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x88(%rbp),%xmm0
ff
movsd  0xe87(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x10(%rbp),%rax
imul   $0x3,-0x30(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
imul   $0x3,-0x30(%rbp),%ecx
add    $0x2,%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movsd  0xe39(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x60(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    -0x8(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
addsd  -0x48(%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
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1122 <radb3+0x22>
cmpl   $0x1,-0x4(%rbp)
jne    1243 <radb3+0x143>
jmp    164d <radb3+0x54d>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    164d <radb3+0x54d>
movl   $0x2,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    163a <radb3+0x53a>
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
movsd  (%rax,%rcx,8),%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
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
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
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x88(%rbp),%xmm0
ff
movsd  0xd24(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%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
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
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
movsd  %xmm0,(%rax,%rcx,8)
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
movsd  (%rax,%rcx,8),%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
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   $0x3,-0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x80(%rbp),%xmm0
movsd  0xc75(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   $0x3,-0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x80(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
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
movsd  (%rax,%rcx,8),%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
subsd  (%rax,%rcx,8),%xmm1
movsd  0xbe4(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%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
movsd  (%rax,%rcx,8),%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
addsd  (%rax,%rcx,8),%xmm1
movsd  0xb9b(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x60(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x60(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x40(%rbp),%xmm0
addsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x40(%rbp),%xmm0
subsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x70(%rbp),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x50(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x30(%rbp),%edx
add    -0x8(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x50(%rbp),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x70(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %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
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x78(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x58(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -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
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x58(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x78(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
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
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x2c(%rbp)
jmp    125d <radb3+0x15d>
jmp    163f <radb3+0x53f>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    124a <radb3+0x14a>
add    $0x8,%rsp
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    1401 <radb3+0x301>
lea    (%rsi,%rsi,1),%eax
movslq -0x44(%rsp),%rbx
movslq %eax,%rdi
mov    %esi,%ebp
mov    %rbx,%r11
mov    %rdi,%rax
mov    %rdi,-0x28(%rsp)
imul   %rdi,%r11
lea    0x0(,%rbx,8),%r10
mov    %rbx,%rdi
imul   %rbp,%rdi
mov    %rbx,%r13
shl    $0x4,%r13
lea    (%r10,%r10,2),%r14
mov    %rbx,-0x38(%rsp)
lea    (%rbx,%rbx,1),%r12
add    $0xffffffffffffffff,%r12
movsd  0xe9f(%rip),%xmm0
movsd  0xe9f(%rip),%xmm1
mov    %rbp,%rax
mov    %rbp,-0x30(%rsp)
mov    %rbp,%rbx
mov    %rdx,%rax
mov    %rcx,%rbp
nopw   0x0(%rax,%rax,1)
movsd  (%rax,%r12,8),%xmm2
addsd  %xmm2,%xmm2
movsd  (%rax),%xmm3
movapd %xmm2,%xmm4
mulsd  %xmm0,%xmm4
addsd  %xmm3,%xmm4
addsd  %xmm3,%xmm2
movsd  %xmm2,0x0(%rbp)
movsd  (%rax,%r13,1),%xmm2
mulsd  %xmm1,%xmm2
movapd %xmm4,%xmm3
subsd  %xmm2,%xmm3
movsd  %xmm3,0x0(%rbp,%rdi,8)
addsd  %xmm4,%xmm2
movsd  %xmm2,0x0(%rbp,%r11,8)
add    %r10,%rbp
add    %r14,%rax
add    $0xffffffffffffffff,%rbx
jne    1180 <radb3+0x80>
cmpl   $0x1,-0x44(%rsp)
mov    -0x38(%rsp),%rdi
je     1401 <radb3+0x301>
test   %esi,%esi
jle    1401 <radb3+0x301>
xor    %eax,%eax
mov    %rax,-0x40(%rsp)
movabs $0xffffffff00000000,%rbp
ff ff ff
movsd  0xe0d(%rip),%xmm1
jmp    122c <radb3+0x12c>
cs nopw 0x0(%rax,%rax,1)
nop
mov    -0x40(%rsp),%rsi
add    $0x1,%rsi
mov    %rsi,%rax
mov    %rsi,-0x40(%rsp)
cmp    -0x30(%rsp),%rsi
je     1401 <radb3+0x301>
cmpl   $0x3,-0x44(%rsp)
jl     1210 <radb3+0x110>
mov    -0x40(%rsp),%rbx
lea    (%rbx,%rbx,2),%rax
lea    (%rbx,%rbx,2),%rsi
add    $0x2,%rsi
imul   %rdi,%rsi
mov    %rsi,-0x10(%rsp)
lea    (%rbx,%rbx,2),%rsi
add    $0x1,%rsi
imul   %rdi,%rsi
mov    %rsi,-0x18(%rsp)
imul   %rdi,%rax
mov    %rax,-0x8(%rsp)
mov    %rbx,%rax
imul   %rdi,%rax
mov    %rax,-0x20(%rsp)
mov    -0x30(%rsp),%rsi
lea    (%rbx,%rsi,1),%r15
imul   %rdi,%r15
mov    -0x28(%rsp),%rsi
lea    (%rbx,%rsi,1),%r10
imul   %rdi,%r10
mov    $0x2,%r11d
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    -0x38(%rsp),%rdi
sub    %r11,%rdi
mov    -0x10(%rsp),%rax
lea    (%r11,%rax,1),%r14
add    $0xffffffffffffffff,%r14
movsd  (%rdx,%r14,8),%xmm2
add    -0x18(%rsp),%rdi
mov    %rdi,%r13
shl    $0x20,%r13
add    %rbp,%r13
sar    $0x1d,%r13
addsd  (%rdx,%r13,1),%xmm2
mov    -0x8(%rsp),%rsi
lea    (%r11,%rsi,1),%r12
add    $0xffffffffffffffff,%r12
movsd  (%rdx,%r12,8),%xmm5
movapd %xmm2,%xmm4
addsd  %xmm5,%xmm2
mov    -0x20(%rsp),%r12
lea    -0x1(%r11,%r12,1),%rbx
movsd  %xmm2,(%rcx,%rbx,8)
lea    (%r11,%rax,1),%rbx
movsd  (%rdx,%rbx,8),%xmm3
subsd  (%rdx,%rdi,8),%xmm3
lea    (%r11,%rsi,1),%rsi
movsd  (%rdx,%rsi,8),%xmm6
movapd %xmm3,%xmm2
addsd  %xmm6,%xmm3
lea    (%r11,%r12,1),%rsi
movsd  %xmm3,(%rcx,%rsi,8)
movsd  (%rdx,%r14,8),%xmm7
subsd  (%rdx,%r13,1),%xmm7
movsd  (%rdx,%rbx,8),%xmm3
mulsd  %xmm0,%xmm4
addsd  %xmm5,%xmm4
mulsd  %xmm0,%xmm2
addsd  %xmm6,%xmm2
mulsd  %xmm1,%xmm7
addsd  (%rdx,%rdi,8),%xmm3
mov    -0x38(%rsp),%rdi
mulsd  %xmm1,%xmm3
movapd %xmm4,%xmm5
subsd  %xmm3,%xmm5
addsd  %xmm4,%xmm3
movapd %xmm2,%xmm4
addsd  %xmm7,%xmm4
subsd  %xmm7,%xmm2
movsd  -0x8(%r8,%r11,8),%xmm6
mulsd  %xmm4,%xmm6
movsd  -0x10(%r8,%r11,8),%xmm7
mulsd  %xmm5,%xmm7
subsd  %xmm6,%xmm7
lea    -0x1(%r11,%r15,1),%rsi
movsd  %xmm7,(%rcx,%rsi,8)
mulsd  -0x8(%r8,%r11,8),%xmm5
mulsd  -0x10(%r8,%r11,8),%xmm4
addsd  %xmm5,%xmm4
lea    (%r11,%r15,1),%rsi
movsd  %xmm4,(%rcx,%rsi,8)
movsd  -0x8(%r9,%r11,8),%xmm4
mulsd  %xmm2,%xmm4
movsd  -0x10(%r9,%r11,8),%xmm5
mulsd  %xmm3,%xmm5
subsd  %xmm4,%xmm5
lea    -0x1(%r11,%r10,1),%rsi
movsd  %xmm5,(%rcx,%rsi,8)
mulsd  -0x8(%r9,%r11,8),%xmm3
mulsd  -0x10(%r9,%r11,8),%xmm2
addsd  %xmm3,%xmm2
lea    (%r11,%r10,1),%rsi
movsd  %xmm2,(%rcx,%rsi,8)
add    $0x2,%r11
cmp    %rdi,%r11
jl     12a0 <radb3+0x1a0>
jmp    1210 <radb3+0x110>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@radb3.taur = internal constant double -5.000000e-01
@radb3.taui = internal constant double 0x3FEBB67AE8584CAE
define void @radb3(i32 %0, i32 %1, double* %2, double* %3, double* %4, double* %5) {
%7 = alloca i32
%8 = alloca i32
%9 = alloca double*
%10 = alloca double*
%11 = alloca double*
%12 = alloca double*
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double
store i32 %0, i32* %7
store i32 %1, i32* %8
store double* %2, double** %9
store double* %3, double** %10
store double* %4, double** %11
store double* %5, double** %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 double*, double** %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 double, double* %31, i64 %40
%42 = load double, double* %41
%43 = fmul double 2.000000e+00, %42
store double %43, double* %25
%44 = load double*, double** %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 double, double* %44, i64 %49
%51 = load double, double* %50
%52 = load double, double* %25
%53 = call double @llvm.fmuladd.f64(double -5.000000e-01, double %52, double %51)
store double %53, double* %20
%54 = load double*, double** %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 double, double* %54, i64 %59
%61 = load double, double* %60
%62 = load double, double* %25
%63 = fadd double %61, %62
%64 = load double*, double** %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 double, double* %64, i64 %68
store double %63, double* %69
%70 = load double*, double** %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 double, double* %70, i64 %76
%78 = load double, double* %77
%79 = fmul double 0x3FFBB67AE8584CAE, %78
store double %79, double* %17
%80 = load double, double* %20
%81 = load double, double* %17
%82 = fsub double %80, %81
%83 = load double*, double** %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 double, double* %83, i64 %89
store double %82, double* %90
%91 = load double, double* %20
%92 = load double, double* %17
%93 = fadd double %91, %92
%94 = load double*, double** %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 double, double* %94, i64 %101
store double %93, double* %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 double*, double** %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 double, double* %124, i64 %133
%135 = load double, double* %134
%136 = load double*, double** %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 double, double* %136, i64 %145
%147 = load double, double* %146
%148 = fadd double %135, %147
store double %148, double* %25
%149 = load double*, double** %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 double, double* %149, i64 %157
%159 = load double, double* %158
%160 = load double, double* %25
%161 = call double @llvm.fmuladd.f64(double -5.000000e-01, double %160, double %159)
store double %161, double* %20
%162 = load double*, double** %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 double, double* %162, i64 %170
%172 = load double, double* %171
%173 = load double, double* %25
%174 = fadd double %172, %173
%175 = load double*, double** %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 double, double* %175, i64 %182
store double %174, double* %183
%184 = load double*, double** %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 double, double* %184, i64 %192
%194 = load double, double* %193
%195 = load double*, double** %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 double, double* %195, i64 %203
%205 = load double, double* %204
%206 = fsub double %194, %205
store double %206, double* %24
%207 = load double*, double** %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 double, double* %207, i64 %214
%216 = load double, double* %215
%217 = load double, double* %24
%218 = call double @llvm.fmuladd.f64(double -5.000000e-01, double %217, double %216)
store double %218, double* %16
%219 = load double*, double** %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 double, double* %219, i64 %226
%228 = load double, double* %227
%229 = load double, double* %24
%230 = fadd double %228, %229
%231 = load double*, double** %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 double, double* %231, i64 %237
store double %230, double* %238
%239 = load double*, double** %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 double, double* %239, i64 %248
%250 = load double, double* %249
%251 = load double*, double** %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 double, double* %251, i64 %260
%262 = load double, double* %261
%263 = fsub double %250, %262
%264 = fmul double 0x3FEBB67AE8584CAE, %263
store double %264, double* %21
%265 = load double*, double** %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 double, double* %265, i64 %273
%275 = load double, double* %274
%276 = load double*, double** %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 double, double* %276, i64 %284
%286 = load double, double* %285
%287 = fadd double %275, %286
%288 = fmul double 0x3FEBB67AE8584CAE, %287
store double %288, double* %17
%289 = load double, double* %20
%290 = load double, double* %17
%291 = fsub double %289, %290
store double %291, double* %22
%292 = load double, double* %20
%293 = load double, double* %17
%294 = fadd double %292, %293
store double %294, double* %23
%295 = load double, double* %16
%296 = load double, double* %21
%297 = fadd double %295, %296
store double %297, double* %18
%298 = load double, double* %16
%299 = load double, double* %21
%300 = fsub double %298, %299
store double %300, double* %19
%301 = load double*, double** %11
%302 = load i32, i32* %13
%303 = sub nsw i32 %302, 2
%304 = sext i32 %303 to i64
%305 = getelementptr inbounds double, double* %301, i64 %304
%306 = load double, double* %305
%307 = load double, double* %22
%308 = load double*, double** %11
%309 = load i32, i32* %13
%310 = sub nsw i32 %309, 1
%311 = sext i32 %310 to i64
%312 = getelementptr inbounds double, double* %308, i64 %311
%313 = load double, double* %312
%314 = load double, double* %18
%315 = fmul double %313, %314
%316 = fneg double %315
%317 = call double @llvm.fmuladd.f64(double %306, double %307, double %316)
%318 = load double*, double** %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 double, double* %318, i64 %327
store double %317, double* %328
%329 = load double*, double** %11
%330 = load i32, i32* %13
%331 = sub nsw i32 %330, 2
%332 = sext i32 %331 to i64
%333 = getelementptr inbounds double, double* %329, i64 %332
%334 = load double, double* %333
%335 = load double, double* %18
%336 = load double*, double** %11
%337 = load i32, i32* %13
%338 = sub nsw i32 %337, 1
%339 = sext i32 %338 to i64
%340 = getelementptr inbounds double, double* %336, i64 %339
%341 = load double, double* %340
%342 = load double, double* %22
%343 = fmul double %341, %342
%344 = call double @llvm.fmuladd.f64(double %334, double %335, double %343)
%345 = load double*, double** %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 double, double* %345, i64 %353
store double %344, double* %354
%355 = load double*, double** %12
%356 = load i32, i32* %13
%357 = sub nsw i32 %356, 2
%358 = sext i32 %357 to i64
%359 = getelementptr inbounds double, double* %355, i64 %358
%360 = load double, double* %359
%361 = load double, double* %23
%362 = load double*, double** %12
%363 = load i32, i32* %13
%364 = sub nsw i32 %363, 1
%365 = sext i32 %364 to i64
%366 = getelementptr inbounds double, double* %362, i64 %365
%367 = load double, double* %366
%368 = load double, double* %19
%369 = fmul double %367, %368
%370 = fneg double %369
%371 = call double @llvm.fmuladd.f64(double %360, double %361, double %370)
%372 = load double*, double** %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 double, double* %372, i64 %382
store double %371, double* %383
%384 = load double*, double** %12
%385 = load i32, i32* %13
%386 = sub nsw i32 %385, 2
%387 = sext i32 %386 to i64
%388 = getelementptr inbounds double, double* %384, i64 %387
%389 = load double, double* %388
%390 = load double, double* %19
%391 = load double*, double** %12
%392 = load i32, i32* %13
%393 = sub nsw i32 %392, 1
%394 = sext i32 %393 to i64
%395 = getelementptr inbounds double, double* %391, i64 %394
%396 = load double, double* %395
%397 = load double, double* %23
%398 = fmul double %396, %397
%399 = call double @llvm.fmuladd.f64(double %389, double %390, double %398)
%400 = load double*, double** %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 double, double* %400, i64 %409
store double %399, double* %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, double* nocapture %2, double* nocapture writeonly %3, double* nocapture %4, double* 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 double, double* %2, i64 %21
%23 = load double, double* %22
%24 = fmul double %23, 2.000000e+00
%25 = mul nsw i64 %18, %11
%26 = getelementptr inbounds double, double* %2, i64 %25
%27 = load double, double* %26
%28 = tail call double @llvm.fmuladd.f64(double %24, double -5.000000e-01, double %27)
%29 = fadd double %27, %24
%30 = mul nsw i64 %17, %11
%31 = getelementptr inbounds double, double* %3, i64 %30
store double %29, double* %31
%32 = add nuw nsw i64 %18, 2
%33 = mul nsw i64 %32, %11
%34 = getelementptr inbounds double, double* %2, i64 %33
%35 = load double, double* %34
%36 = fmul double %35, 0x3FFBB67AE8584CAE
%37 = fsub double %28, %36
%38 = add nuw nsw i64 %17, %13
%39 = mul nsw i64 %38, %11
%40 = getelementptr inbounds double, double* %3, i64 %39
store double %37, double* %40
%41 = fadd double %28, %36
%42 = add nsw i64 %17, %14
%43 = mul nsw i64 %42, %11
%44 = getelementptr inbounds double, double* %3, i64 %43
store double %41, double* %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 double, double* %2, i64 %75
%77 = load double, double* %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 double, double* %2, i64 %81
%83 = load double, double* %82
%84 = fadd double %77, %83
%85 = add nsw i64 %74, %65
%86 = getelementptr inbounds double, double* %2, i64 %85
%87 = load double, double* %86
%88 = tail call double @llvm.fmuladd.f64(double %84, double -5.000000e-01, double %87)
%89 = fadd double %87, %84
%90 = add nsw i64 %74, %66
%91 = getelementptr inbounds double, double* %3, i64 %90
store double %89, double* %91
%92 = add nsw i64 %72, %62
%93 = getelementptr inbounds double, double* %2, i64 %92
%94 = load double, double* %93
%95 = add nsw i64 %73, %64
%96 = getelementptr inbounds double, double* %2, i64 %95
%97 = load double, double* %96
%98 = fsub double %94, %97
%99 = add nsw i64 %72, %65
%100 = getelementptr inbounds double, double* %2, i64 %99
%101 = load double, double* %100
%102 = tail call double @llvm.fmuladd.f64(double %98, double -5.000000e-01, double %101)
%103 = fadd double %101, %98
%104 = add nsw i64 %72, %66
%105 = getelementptr inbounds double, double* %3, i64 %104
store double %103, double* %105
%106 = load double, double* %76
%107 = load double, double* %82
%108 = fsub double %106, %107
%109 = fmul double %108, 0x3FEBB67AE8584CAE
%110 = load double, double* %93
%111 = load double, double* %96
%112 = fadd double %110, %111
%113 = fmul double %112, 0x3FEBB67AE8584CAE
%114 = fsub double %88, %113
%115 = fadd double %88, %113
%116 = fadd double %102, %109
%117 = fsub double %102, %109
%118 = add nsw i64 %72, -2
%119 = getelementptr inbounds double, double* %4, i64 %118
%120 = load double, double* %119
%121 = getelementptr inbounds double, double* %4, i64 %74
%122 = load double, double* %121
%123 = fneg double %122
%124 = fmul double %116, %123
%125 = tail call double @llvm.fmuladd.f64(double %120, double %114, double %124)
%126 = add nsw i64 %74, %68
%127 = getelementptr inbounds double, double* %3, i64 %126
store double %125, double* %127
%128 = load double, double* %119
%129 = load double, double* %121
%130 = fmul double %114, %129
%131 = tail call double @llvm.fmuladd.f64(double %128, double %116, double %130)
%132 = add nsw i64 %72, %68
%133 = getelementptr inbounds double, double* %3, i64 %132
store double %131, double* %133
%134 = getelementptr inbounds double, double* %5, i64 %118
%135 = load double, double* %134
%136 = getelementptr inbounds double, double* %5, i64 %74
%137 = load double, double* %136
%138 = fneg double %137
%139 = fmul double %117, %138
%140 = tail call double @llvm.fmuladd.f64(double %135, double %115, double %139)
%141 = add nsw i64 %74, %70
%142 = getelementptr inbounds double, double* %3, i64 %141
store double %140, double* %142
%143 = load double, double* %134
%144 = load double, double* %136
%145 = fmul double %115, %144
%146 = tail call double @llvm.fmuladd.f64(double %143, double %117, double %145)
%147 = add nsw i64 %72, %70
%148 = getelementptr inbounds double, double* %3, i64 %147
store double %146, double* %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
}