void dradb3(int ido, int l1, double *cc, double *ch, double *wa1, double *wa2) {
  static double taur = -.5;
  static double taui = .86602540378443864676372317075293618;
  int i, k, t0, t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;
  double 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
sub    $0x30,%rsp
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    1276 <dradb3+0x176>
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  0x2e76(%rip),%xmm0
movsd  -0xb0(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0xb0(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  0x2e3c(%rip),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
subsd  -0x70(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x88(%rbp),%xmm0
ff
addsd  -0x70(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
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    115d <dradb3+0x5d>
cmpl   $0x1,-0x4(%rbp)
jne    1285 <dradb3+0x185>
jmp    163e <dradb3+0x53e>
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    163e <dradb3+0x53e>
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    1627 <dradb3+0x527>
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
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  0x2cba(%rip),%xmm0
movsd  -0xb0(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0xb0(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  0x2c3f(%rip),%xmm0
movsd  -0xa8(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0xa8(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  0x2c08(%rip),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  0x2bd0(%rip),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
subsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x88(%rbp),%xmm0
ff
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x68(%rbp),%xmm0
addsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x78(%rbp)
movsd  -0x68(%rbp),%xmm0
subsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x80(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x98(%rbp),%xmm2
ff
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x78(%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    -0x58(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x78(%rbp),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x98(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0xa0(%rbp),%xmm2
ff
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x80(%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    -0x5c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x80(%rbp),%xmm2
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0xa0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x2c(%rbp)
jmp    12de <dradb3+0x1de>
mov    -0x4(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    129c <dradb3+0x19c>
add    $0x30,%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    13ae <dradb3+0x2ae>
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    $0x3,%rsi
lea    0x0(,%rbx,8),%rdi
movsd  0xe97(%rip),%xmm0
movsd  0xe97(%rip),%xmm1
mov    -0x30(%rsp),%rax
mov    -0x38(%rsp),%rcx
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%rcx,%rdx,8),%xmm2
addsd  %xmm2,%xmm2
movsd  (%rcx),%xmm3
movapd %xmm2,%xmm4
mulsd  %xmm0,%xmm4
addsd  %xmm3,%xmm4
addsd  %xmm3,%xmm2
movsd  %xmm2,(%rax)
movsd  (%rcx,%rdx,8),%xmm2
addsd  %xmm2,%xmm2
mulsd  %xmm1,%xmm2
movapd %xmm4,%xmm3
subsd  %xmm2,%xmm3
movsd  %xmm3,(%rax,%r14,8)
addsd  %xmm4,%xmm2
movsd  %xmm2,(%rax,%r10,8)
add    %rsi,%rcx
add    %rdi,%rax
add    $0xffffffff,%ebp
jne    1180 <dradb3+0x80>
mov    -0x18(%rsp),%rcx
cmp    $0x1,%ecx
je     13ae <dradb3+0x2ae>
cmpl   $0x0,-0x3c(%rsp)
jle    13ae <dradb3+0x2ae>
mov    %r11d,%ebp
mov    -0x38(%rsp),%rax
add    $0xfffffffffffffff0,%rax
mov    %rax,-0x10(%rsp)
xor    %r11d,%r11d
mov    -0x30(%rsp),%rax
xor    %r15d,%r15d
mov    %rdi,-0x8(%rsp)
jmp    1242 <dradb3+0x142>
xchg   %ax,%ax
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     13ae <dradb3+0x2ae>
cmp    $0x3,%ecx
mov    %rcx,%rdi
jl     1210 <dradb3+0x110>
movslq -0x28(%rsp),%rcx
mov    -0x30(%rsp),%rdx
lea    (%rdx,%rcx,8),%r14
movslq %ebp,%rcx
lea    (%rdx,%rcx,8),%rcx
movslq %r11d,%rdx
mov    %ebp,%r10d
mov    -0x38(%rsp),%rbp
lea    0x0(%rbp,%rdx,8),%rdx
movslq %r13d,%rdi
mov    -0x10(%rsp),%rsi
lea    (%rsi,%rdi,8),%rsi
lea    0x0(,%rdi,8),%r12
add    %rbp,%r12
mov    %r10d,%ebp
mov    -0x18(%rsp),%rdi
mov    $0x2,%r10d
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  -0x8(%r12,%r10,8),%xmm2
addsd  -0x8(%rsi),%xmm2
movsd  -0x8(%rdx,%r10,8),%xmm3
movapd %xmm2,%xmm4
mulsd  %xmm0,%xmm4
addsd  %xmm3,%xmm4
addsd  %xmm3,%xmm2
movsd  %xmm2,-0x8(%rax,%r10,8)
movsd  (%r12,%r10,8),%xmm3
subsd  (%rsi),%xmm3
movsd  (%rdx,%r10,8),%xmm5
movapd %xmm3,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm5,%xmm2
addsd  %xmm5,%xmm3
movsd  %xmm3,(%rax,%r10,8)
movsd  -0x8(%r12,%r10,8),%xmm5
movsd  (%r12,%r10,8),%xmm3
subsd  -0x8(%rsi),%xmm5
mulsd  %xmm1,%xmm5
addsd  (%rsi),%xmm3
mulsd  %xmm1,%xmm3
movapd %xmm4,%xmm6
subsd  %xmm3,%xmm6
addsd  %xmm4,%xmm3
movapd %xmm2,%xmm4
addsd  %xmm5,%xmm4
subsd  %xmm5,%xmm2
movsd  -0x8(%r8,%r10,8),%xmm5
mulsd  %xmm4,%xmm5
movsd  -0x10(%r8,%r10,8),%xmm7
mulsd  %xmm6,%xmm7
subsd  %xmm5,%xmm7
movsd  %xmm7,-0x8(%r14,%r10,8)
mulsd  -0x8(%r8,%r10,8),%xmm6
mulsd  -0x10(%r8,%r10,8),%xmm4
addsd  %xmm6,%xmm4
movsd  %xmm4,(%r14,%r10,8)
movsd  -0x8(%r9,%r10,8),%xmm4
mulsd  %xmm2,%xmm4
movsd  -0x10(%r9,%r10,8),%xmm5
mulsd  %xmm3,%xmm5
subsd  %xmm4,%xmm5
movsd  %xmm5,-0x8(%rcx,%r10,8)
mulsd  -0x8(%r9,%r10,8),%xmm3
mulsd  -0x10(%r9,%r10,8),%xmm2
addsd  %xmm3,%xmm2
movsd  %xmm2,(%rcx,%r10,8)
add    $0x2,%r10
add    $0xfffffffffffffff0,%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 double -5.000000e-01
@dradb3.taui = internal global double 0x3FEBB67AE8584CAA
define void @dradb3(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 i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca i32
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca double
%27 = alloca double
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca double
%32 = alloca double
%33 = alloca double
%34 = alloca double
%35 = alloca double
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
%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 double*, double** %9
%53 = load i32, i32* %18
%54 = sub nsw i32 %53, 1
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %52, i64 %55
%57 = load double, double* %56
%58 = load double*, double** %9
%59 = load i32, i32* %18
%60 = sub nsw i32 %59, 1
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double, double* %58, i64 %61
%63 = load double, double* %62
%64 = fadd double %57, %63
store double %64, double* %35
%65 = load double*, double** %9
%66 = load i32, i32* %20
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds double, double* %65, i64 %67
%69 = load double, double* %68
%70 = load double, double* @dradb3.taur
%71 = load double, double* %35
%72 = call double @llvm.fmuladd.f64(double %70, double %71, double %69)
store double %72, double* %30
%73 = load double*, double** %9
%74 = load i32, i32* %20
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds double, double* %73, i64 %75
%77 = load double, double* %76
%78 = load double, double* %35
%79 = fadd double %77, %78
%80 = load double*, double** %10
%81 = load i32, i32* %16
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds double, double* %80, i64 %82
store double %79, double* %83
%84 = load double, double* @dradb3.taui
%85 = load double*, double** %9
%86 = load i32, i32* %18
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double, double* %85, i64 %87
%89 = load double, double* %88
%90 = load double*, double** %9
%91 = load i32, i32* %18
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds double, double* %90, i64 %92
%94 = load double, double* %93
%95 = fadd double %89, %94
%96 = fmul double %84, %95
store double %96, double* %27
%97 = load double, double* %30
%98 = load double, double* %27
%99 = fsub double %97, %98
%100 = load double*, double** %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 double, double* %100, i64 %104
store double %99, double* %105
%106 = load double, double* %30
%107 = load double, double* %27
%108 = fadd double %106, %107
%109 = load double*, double** %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 double, double* %109, i64 %113
store double %108, double* %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 double*, double** %9
%170 = load i32, i32* %20
%171 = sub nsw i32 %170, 1
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds double, double* %169, i64 %172
%174 = load double, double* %173
%175 = load double*, double** %9
%176 = load i32, i32* %21
%177 = sub nsw i32 %176, 1
%178 = sext i32 %177 to i64
%179 = getelementptr inbounds double, double* %175, i64 %178
%180 = load double, double* %179
%181 = fadd double %174, %180
store double %181, double* %35
%182 = load double*, double** %9
%183 = load i32, i32* %22
%184 = sub nsw i32 %183, 1
%185 = sext i32 %184 to i64
%186 = getelementptr inbounds double, double* %182, i64 %185
%187 = load double, double* %186
%188 = load double, double* @dradb3.taur
%189 = load double, double* %35
%190 = call double @llvm.fmuladd.f64(double %188, double %189, double %187)
store double %190, double* %30
%191 = load double*, double** %9
%192 = load i32, i32* %22
%193 = sub nsw i32 %192, 1
%194 = sext i32 %193 to i64
%195 = getelementptr inbounds double, double* %191, i64 %194
%196 = load double, double* %195
%197 = load double, double* %35
%198 = fadd double %196, %197
%199 = load double*, double** %10
%200 = load i32, i32* %23
%201 = sub nsw i32 %200, 1
%202 = sext i32 %201 to i64
%203 = getelementptr inbounds double, double* %199, i64 %202
store double %198, double* %203
%204 = load double*, double** %9
%205 = load i32, i32* %20
%206 = sext i32 %205 to i64
%207 = getelementptr inbounds double, double* %204, i64 %206
%208 = load double, double* %207
%209 = load double*, double** %9
%210 = load i32, i32* %21
%211 = sext i32 %210 to i64
%212 = getelementptr inbounds double, double* %209, i64 %211
%213 = load double, double* %212
%214 = fsub double %208, %213
store double %214, double* %34
%215 = load double*, double** %9
%216 = load i32, i32* %22
%217 = sext i32 %216 to i64
%218 = getelementptr inbounds double, double* %215, i64 %217
%219 = load double, double* %218
%220 = load double, double* @dradb3.taur
%221 = load double, double* %34
%222 = call double @llvm.fmuladd.f64(double %220, double %221, double %219)
store double %222, double* %26
%223 = load double*, double** %9
%224 = load i32, i32* %22
%225 = sext i32 %224 to i64
%226 = getelementptr inbounds double, double* %223, i64 %225
%227 = load double, double* %226
%228 = load double, double* %34
%229 = fadd double %227, %228
%230 = load double*, double** %10
%231 = load i32, i32* %23
%232 = sext i32 %231 to i64
%233 = getelementptr inbounds double, double* %230, i64 %232
store double %229, double* %233
%234 = load double, double* @dradb3.taui
%235 = load double*, double** %9
%236 = load i32, i32* %20
%237 = sub nsw i32 %236, 1
%238 = sext i32 %237 to i64
%239 = getelementptr inbounds double, double* %235, i64 %238
%240 = load double, double* %239
%241 = load double*, double** %9
%242 = load i32, i32* %21
%243 = sub nsw i32 %242, 1
%244 = sext i32 %243 to i64
%245 = getelementptr inbounds double, double* %241, i64 %244
%246 = load double, double* %245
%247 = fsub double %240, %246
%248 = fmul double %234, %247
store double %248, double* %31
%249 = load double, double* @dradb3.taui
%250 = load double*, double** %9
%251 = load i32, i32* %20
%252 = sext i32 %251 to i64
%253 = getelementptr inbounds double, double* %250, i64 %252
%254 = load double, double* %253
%255 = load double*, double** %9
%256 = load i32, i32* %21
%257 = sext i32 %256 to i64
%258 = getelementptr inbounds double, double* %255, i64 %257
%259 = load double, double* %258
%260 = fadd double %254, %259
%261 = fmul double %249, %260
store double %261, double* %27
%262 = load double, double* %30
%263 = load double, double* %27
%264 = fsub double %262, %263
store double %264, double* %32
%265 = load double, double* %30
%266 = load double, double* %27
%267 = fadd double %265, %266
store double %267, double* %33
%268 = load double, double* %26
%269 = load double, double* %31
%270 = fadd double %268, %269
store double %270, double* %28
%271 = load double, double* %26
%272 = load double, double* %31
%273 = fsub double %271, %272
store double %273, double* %29
%274 = load double*, double** %11
%275 = load i32, i32* %13
%276 = sub nsw i32 %275, 2
%277 = sext i32 %276 to i64
%278 = getelementptr inbounds double, double* %274, i64 %277
%279 = load double, double* %278
%280 = load double, double* %32
%281 = load double*, double** %11
%282 = load i32, i32* %13
%283 = sub nsw i32 %282, 1
%284 = sext i32 %283 to i64
%285 = getelementptr inbounds double, double* %281, i64 %284
%286 = load double, double* %285
%287 = load double, double* %28
%288 = fmul double %286, %287
%289 = fneg double %288
%290 = call double @llvm.fmuladd.f64(double %279, double %280, double %289)
%291 = load double*, double** %10
%292 = load i32, i32* %24
%293 = sub nsw i32 %292, 1
%294 = sext i32 %293 to i64
%295 = getelementptr inbounds double, double* %291, i64 %294
store double %290, double* %295
%296 = load double*, double** %11
%297 = load i32, i32* %13
%298 = sub nsw i32 %297, 2
%299 = sext i32 %298 to i64
%300 = getelementptr inbounds double, double* %296, i64 %299
%301 = load double, double* %300
%302 = load double, double* %28
%303 = load double*, double** %11
%304 = load i32, i32* %13
%305 = sub nsw i32 %304, 1
%306 = sext i32 %305 to i64
%307 = getelementptr inbounds double, double* %303, i64 %306
%308 = load double, double* %307
%309 = load double, double* %32
%310 = fmul double %308, %309
%311 = call double @llvm.fmuladd.f64(double %301, double %302, double %310)
%312 = load double*, double** %10
%313 = load i32, i32* %24
%314 = sext i32 %313 to i64
%315 = getelementptr inbounds double, double* %312, i64 %314
store double %311, double* %315
%316 = load double*, double** %12
%317 = load i32, i32* %13
%318 = sub nsw i32 %317, 2
%319 = sext i32 %318 to i64
%320 = getelementptr inbounds double, double* %316, i64 %319
%321 = load double, double* %320
%322 = load double, double* %33
%323 = load double*, double** %12
%324 = load i32, i32* %13
%325 = sub nsw i32 %324, 1
%326 = sext i32 %325 to i64
%327 = getelementptr inbounds double, double* %323, i64 %326
%328 = load double, double* %327
%329 = load double, double* %29
%330 = fmul double %328, %329
%331 = fneg double %330
%332 = call double @llvm.fmuladd.f64(double %321, double %322, double %331)
%333 = load double*, double** %10
%334 = load i32, i32* %25
%335 = sub nsw i32 %334, 1
%336 = sext i32 %335 to i64
%337 = getelementptr inbounds double, double* %333, i64 %336
store double %332, double* %337
%338 = load double*, double** %12
%339 = load i32, i32* %13
%340 = sub nsw i32 %339, 2
%341 = sext i32 %340 to i64
%342 = getelementptr inbounds double, double* %338, i64 %341
%343 = load double, double* %342
%344 = load double, double* %29
%345 = load double*, double** %12
%346 = load i32, i32* %13
%347 = sub nsw i32 %346, 1
%348 = sext i32 %347 to i64
%349 = getelementptr inbounds double, double* %345, i64 %348
%350 = load double, double* %349
%351 = load double, double* %33
%352 = fmul double %350, %351
%353 = call double @llvm.fmuladd.f64(double %343, double %344, double %352)
%354 = load double*, double** %10
%355 = load i32, i32* %25
%356 = sext i32 %355 to i64
%357 = getelementptr inbounds double, double* %354, i64 %356
store double %353, double* %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, double* nocapture %2, double* nocapture writeonly %3, double* nocapture %4, double* 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 double, double* %2, i64 %23
%25 = load double, double* %24
%26 = fadd double %25, %25
%27 = getelementptr inbounds double, double* %2, i64 %19
%28 = load double, double* %27
%29 = tail call double @llvm.fmuladd.f64(double %26, double -5.000000e-01, double %28)
%30 = fadd double %28, %26
%31 = getelementptr inbounds double, double* %3, i64 %21
store double %30, double* %31
%32 = getelementptr inbounds double, double* %2, i64 %20
%33 = load double, double* %32
%34 = fadd double %33, %33
%35 = fmul double %34, 0x3FEBB67AE8584CAA
%36 = fsub double %29, %35
%37 = add nsw i64 %21, %15
%38 = getelementptr inbounds double, double* %3, i64 %37
store double %36, double* %38
%39 = fadd double %29, %35
%40 = add nsw i64 %21, %14
%41 = getelementptr inbounds double, double* %3, i64 %40
store double %39, double* %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 double, double* %2, i64 %79
%81 = load double, double* %80
%82 = add nsw i64 %70, -3
%83 = getelementptr inbounds double, double* %2, i64 %82
%84 = load double, double* %83
%85 = fadd double %81, %84
%86 = add nsw i64 %69, 1
%87 = getelementptr inbounds double, double* %2, i64 %86
%88 = load double, double* %87
%89 = tail call double @llvm.fmuladd.f64(double %85, double -5.000000e-01, double %88)
%90 = fadd double %88, %85
%91 = add nsw i64 %66, 1
%92 = getelementptr inbounds double, double* %3, i64 %91
store double %90, double* %92
%93 = getelementptr inbounds double, double* %2, i64 %73
%94 = load double, double* %93
%95 = getelementptr inbounds double, double* %2, i64 %74
%96 = load double, double* %95
%97 = fsub double %94, %96
%98 = getelementptr inbounds double, double* %2, i64 %75
%99 = load double, double* %98
%100 = tail call double @llvm.fmuladd.f64(double %97, double -5.000000e-01, double %99)
%101 = fadd double %99, %97
%102 = getelementptr inbounds double, double* %3, i64 %76
store double %101, double* %102
%103 = load double, double* %80
%104 = load double, double* %83
%105 = fsub double %103, %104
%106 = fmul double %105, 0x3FEBB67AE8584CAA
%107 = load double, double* %93
%108 = load double, double* %95
%109 = fadd double %107, %108
%110 = fmul double %109, 0x3FEBB67AE8584CAA
%111 = fsub double %89, %110
%112 = fadd double %89, %110
%113 = fadd double %100, %106
%114 = fsub double %100, %106
%115 = add nsw i64 %72, -2
%116 = getelementptr inbounds double, double* %4, i64 %115
%117 = load double, double* %116
%118 = add nsw i64 %72, -1
%119 = getelementptr inbounds double, double* %4, i64 %118
%120 = load double, double* %119
%121 = fneg double %120
%122 = fmul double %113, %121
%123 = tail call double @llvm.fmuladd.f64(double %117, double %111, double %122)
%124 = add nsw i64 %67, 1
%125 = getelementptr inbounds double, double* %3, i64 %124
store double %123, double* %125
%126 = load double, double* %116
%127 = load double, double* %119
%128 = fmul double %111, %127
%129 = tail call double @llvm.fmuladd.f64(double %126, double %113, double %128)
%130 = getelementptr inbounds double, double* %3, i64 %77
store double %129, double* %130
%131 = getelementptr inbounds double, double* %5, i64 %115
%132 = load double, double* %131
%133 = getelementptr inbounds double, double* %5, i64 %118
%134 = load double, double* %133
%135 = fneg double %134
%136 = fmul double %114, %135
%137 = tail call double @llvm.fmuladd.f64(double %132, double %112, double %136)
%138 = add nsw i64 %68, 1
%139 = getelementptr inbounds double, double* %3, i64 %138
store double %137, double* %139
%140 = load double, double* %131
%141 = load double, double* %133
%142 = fmul double %112, %141
%143 = tail call double @llvm.fmuladd.f64(double %140, double %114, double %142)
%144 = getelementptr inbounds double, double* %3, i64 %78
store double %143, double* %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
}