void dradf2(int ido, int l1, double *cc, double *ch, double *wa1) {
  int i, k;
  double ti2, tr2;
  int t0, t1, t2, t3, t4, t5, t6;
  t1 = 0;
  t0 = (t2 = l1 * ido);
  t3 = ido << 1;
  for (k = 0; k < l1; k++) {
    ch[t1 << 1] = cc[t1] + cc[t2];
    ch[(t1 << 1) + t3 - 1] = cc[t1] - cc[t2];
    t1 += ido;
    t2 += ido;
  }
  if (ido < 2)
    return;
  if (ido == 2)
    goto L105;
  t1 = 0;
  t2 = t0;
  for (k = 0; k < l1; k++) {
    t3 = t2;
    t4 = (t1 << 1) + (ido << 1);
    t5 = t1;
    t6 = t1 + t1;
    for (i = 2; i < ido; i += 2) {
      t3 += 2;
      t4 -= 2;
      t5 += 2;
      t6 += 2;
      tr2 = wa1[i - 2] * cc[t3 - 1] + wa1[i - 1] * cc[t3];
      ti2 = wa1[i - 2] * cc[t3] - wa1[i - 1] * cc[t3 - 1];
      ch[t6] = cc[t5] + ti2;
      ch[t4] = ti2 - cc[t5];
      ch[t6 - 1] = cc[t5 - 1] + tr2;
      ch[t4 - 1] = cc[t5 - 1] - tr2;
    }
    t1 += ido;
    t2 += ido;
  }
  if (ido % 2 == 1)
    return;
L105:
  t3 = (t2 = (t1 = ido) - 1);
  t2 += t0;
  for (k = 0; k < l1; k++) {
    ch[t1] = -cc[t2];
    ch[t1 - 1] = cc[t3];
    t1 += ido << 1;
    t2 += ido;
    t3 += ido;
  }
}


x86 assembly (clang, -O0):
<dradf2>:
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)
movl   $0x0,-0x40(%rbp)
mov    -0x8(%rbp),%eax
imul   -0x4(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    %eax,-0x3c(%rbp)
mov    -0x4(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x48(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11c4 <dradf2+0xc4>
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
shl    $0x1,%ecx
add    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x4(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    113a <dradf2+0x3a>
cmpl   $0x2,-0x4(%rbp)
jge    11d3 <dradf2+0xd3>
jmp    1499 <dradf2+0x399>
cmpl   $0x2,-0x4(%rbp)
jne    11e2 <dradf2+0xe2>
jmp    13f2 <dradf2+0x2f2>
movl   $0x0,-0x40(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x44(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    13d4 <dradf2+0x2d4>
mov    -0x44(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x40(%rbp),%eax
shl    $0x1,%eax
mov    -0x4(%rbp),%ecx
shl    $0x1,%ecx
add    %ecx,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x40(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x40(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x54(%rbp)
movl   $0x2,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    13b4 <dradf2+0x2b4>
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    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x54(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x4c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x24(%rbp)
jmp    122f <dradf2+0x12f>
mov    -0x4(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x4(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11f6 <dradf2+0xf6>
mov    -0x4(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x1,%edx
jne    13ed <dradf2+0x2ed>
jmp    1499 <dradf2+0x399>
jmp    13f2 <dradf2+0x2f2>
mov    -0x4(%rbp),%eax
mov    %eax,-0x40(%rbp)
sub    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    %eax,-0x48(%rbp)
mov    -0x3c(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x44(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1499 <dradf2+0x399>
mov    -0x10(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4(%rbp),%eax
shl    $0x1,%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x4(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x4(%rbp),%eax
add    -0x48(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1411 <dradf2+0x311>
pop    %rbp
ret

x86 assembly (clang, -O2):
<dradf2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rdx,-0x18(%rsp)
mov    %esi,%r15d
imul   %edi,%r15d
lea    (%rdi,%rdi,1),%r10d
test   %esi,%esi
jle    1173 <dradf2+0x73>
lea    -0x1(%r10),%r9d
movslq %r15d,%r14
movslq %edi,%r11
shl    $0x3,%r11
xor    %ebp,%ebp
mov    %esi,%ebx
mov    -0x18(%rsp),%rax
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rax),%xmm0
addsd  (%rax,%r14,8),%xmm0
movslq %ebp,%rbp
movsd  %xmm0,(%rcx,%rbp,8)
movsd  (%rax),%xmm0
subsd  (%rax,%r14,8),%xmm0
lea    (%r9,%rbp,1),%edx
movslq %edx,%rdx
movsd  %xmm0,(%rcx,%rdx,8)
add    %r11,%rax
add    %r10d,%ebp
add    $0xffffffff,%ebx
jne    1140 <dradf2+0x40>
cmp    $0x2,%edi
jl     1398 <dradf2+0x298>
mov    %esi,-0xc(%rsp)
je     12a4 <dradf2+0x1a4>
mov    %r15d,-0x10(%rsp)
test   %esi,%esi
jle    1288 <dradf2+0x188>
movslq -0x10(%rsp),%rax
mov    %edi,%r11d
mov    -0x18(%rsp),%rdx
lea    (%rdx,%rax,8),%rax
lea    0x0(,%r11,8),%rdx
mov    %rdx,-0x8(%rsp)
lea    -0x10(%rcx),%r15
xor    %r13d,%r13d
mov    %r10d,%ebx
xor    %r12d,%r12d
xor    %ebp,%ebp
movslq %ebx,%rdx
lea    (%r15,%rdx,8),%r9
movslq %r12d,%r12
mov    -0x18(%rsp),%rdx
lea    (%rdx,%r12,8),%r14
add    %edi,%r12d
movslq %r13d,%rdx
lea    (%rcx,%rdx,8),%rdx
mov    $0x2,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  -0x10(%r8,%rsi,8),%xmm0
movsd  -0x8(%r8,%rsi,8),%xmm1
movsd  -0x8(%rax,%rsi,8),%xmm2
movsd  (%rax,%rsi,8),%xmm3
movapd %xmm0,%xmm4
mulsd  %xmm2,%xmm4
mulsd  %xmm1,%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm1,%xmm4
mulsd  %xmm3,%xmm0
subsd  %xmm2,%xmm0
movsd  (%r14,%rsi,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rdx,%rsi,8)
subsd  (%r14,%rsi,8),%xmm0
movsd  %xmm0,(%r9)
movsd  -0x8(%r14,%rsi,8),%xmm0
addsd  %xmm4,%xmm0
movsd  %xmm0,-0x8(%rdx,%rsi,8)
movsd  -0x8(%r14,%rsi,8),%xmm0
subsd  %xmm4,%xmm0
movsd  %xmm0,-0x8(%r9)
add    $0x2,%rsi
add    $0xfffffffffffffff0,%r9
cmp    %r11,%rsi
jb     11f0 <dradf2+0xf0>
add    $0x1,%ebp
add    %r10d,%r13d
add    %r10d,%ebx
add    -0x8(%rsp),%rax
mov    -0xc(%rsp),%esi
cmp    %esi,%ebp
jne    11c0 <dradf2+0xc0>
mov    %edi,%eax
shr    $0x1f,%eax
add    %edi,%eax
and    $0xfffffffe,%eax
mov    %edi,%edx
sub    %eax,%edx
cmp    $0x1,%edx
mov    -0x10(%rsp),%r15d
je     1398 <dradf2+0x298>
test   %esi,%esi
jle    1398 <dradf2+0x298>
lea    -0x1(%rdi),%eax
add    %eax,%r15d
movslq %eax,%r8
movslq %edi,%rdi
movslq %r15d,%r9
cmp    $0x1,%esi
je     1370 <dradf2+0x270>
movslq %r10d,%rdx
mov    %esi,%ebx
and    $0xfffffffe,%ebx
lea    (%rcx,%rdi,8),%rax
mov    %rdx,%r14
shl    $0x4,%r14
lea    (%rcx,%rdx,8),%rsi
mov    -0x18(%rsp),%r13
lea    0x0(,%r9,8),%r10
add    %r13,%r10
lea    (%rdi,%rdi,1),%r11
lea    (%rdi,%r9,1),%rbp
lea    0x0(%r13,%rbp,8),%r15
lea    0x0(,%r8,8),%r12
add    %r13,%r12
lea    (%rdi,%r8,1),%rbp
lea    0x0(%r13,%rbp,8),%r13
xor    %ebp,%ebp
movapd 0xcea(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
movsd  (%r10,%rbp,8),%xmm1
xorpd  %xmm0,%xmm1
movlpd %xmm1,(%rax)
movsd  (%r12,%rbp,8),%xmm1
movsd  %xmm1,-0x8(%rax)
movsd  (%r15,%rbp,8),%xmm1
xorpd  %xmm0,%xmm1
movlpd %xmm1,(%rsi,%rdi,8)
movsd  0x0(%r13,%rbp,8),%xmm1
movsd  %xmm1,-0x8(%rsi,%rdi,8)
add    %rdx,%rdi
add    %rdx,%rdi
add    %r14,%rax
add    %r11,%rbp
add    $0xfffffffe,%ebx
jne    1320 <dradf2+0x220>
add    %rbp,%r9
add    %rbp,%r8
mov    -0xc(%rsp),%esi
test   $0x1,%sil
je     1398 <dradf2+0x298>
mov    -0x18(%rsp),%rax
movsd  (%rax,%r9,8),%xmm0
xorps  0xc78(%rip),%xmm0
movlps %xmm0,(%rcx,%rdi,8)
movsd  (%rax,%r8,8),%xmm0
movsd  %xmm0,-0x8(%rcx,%rdi,8)
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @dradf2(i32 %0, i32 %1, double* %2, double* %3, double* %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca double*
%9 = alloca double*
%10 = alloca double*
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
%14 = alloca double
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
store i32 %0, i32* %6
store i32 %1, i32* %7
store double* %2, double** %8
store double* %3, double** %9
store double* %4, double** %10
store i32 0, i32* %16
%22 = load i32, i32* %7
%23 = load i32, i32* %6
%24 = mul nsw i32 %22, %23
store i32 %24, i32* %17
store i32 %24, i32* %15
%25 = load i32, i32* %6
%26 = shl i32 %25, 1
store i32 %26, i32* %18
store i32 0, i32* %12
br label %27
27: ; preds = %73, %5
%28 = load i32, i32* %12
%29 = load i32, i32* %7
%30 = icmp slt i32 %28, %29
br i1 %30, label %31, label %76
31: ; preds = %27
%32 = load double*, double** %8
%33 = load i32, i32* %16
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds double, double* %32, i64 %34
%36 = load double, double* %35
%37 = load double*, double** %8
%38 = load i32, i32* %17
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds double, double* %37, i64 %39
%41 = load double, double* %40
%42 = fadd double %36, %41
%43 = load double*, double** %9
%44 = load i32, i32* %16
%45 = shl i32 %44, 1
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds double, double* %43, i64 %46
store double %42, double* %47
%48 = load double*, double** %8
%49 = load i32, i32* %16
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double, double* %48, i64 %50
%52 = load double, double* %51
%53 = load double*, double** %8
%54 = load i32, i32* %17
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %53, i64 %55
%57 = load double, double* %56
%58 = fsub double %52, %57
%59 = load double*, double** %9
%60 = load i32, i32* %16
%61 = shl i32 %60, 1
%62 = load i32, i32* %18
%63 = add nsw i32 %61, %62
%64 = sub nsw i32 %63, 1
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %59, i64 %65
store double %58, double* %66
%67 = load i32, i32* %6
%68 = load i32, i32* %16
%69 = add nsw i32 %68, %67
store i32 %69, i32* %16
%70 = load i32, i32* %6
%71 = load i32, i32* %17
%72 = add nsw i32 %71, %70
store i32 %72, i32* %17
br label %73
73: ; preds = %31
%74 = load i32, i32* %12
%75 = add nsw i32 %74, 1
store i32 %75, i32* %12
br label %27
76: ; preds = %27
%77 = load i32, i32* %6
%78 = icmp slt i32 %77, 2
br i1 %78, label %79, label %80
79: ; preds = %76
br label %276
80: ; preds = %76
%81 = load i32, i32* %6
%82 = icmp eq i32 %81, 2
br i1 %82, label %83, label %84
83: ; preds = %80
br label %232
84: ; preds = %80
store i32 0, i32* %16
%85 = load i32, i32* %15
store i32 %85, i32* %17
store i32 0, i32* %12
br label %86
86: ; preds = %223, %84
%87 = load i32, i32* %12
%88 = load i32, i32* %7
%89 = icmp slt i32 %87, %88
br i1 %89, label %90, label %226
90: ; preds = %86
%91 = load i32, i32* %17
store i32 %91, i32* %18
%92 = load i32, i32* %16
%93 = shl i32 %92, 1
%94 = load i32, i32* %6
%95 = shl i32 %94, 1
%96 = add nsw i32 %93, %95
store i32 %96, i32* %19
%97 = load i32, i32* %16
store i32 %97, i32* %20
%98 = load i32, i32* %16
%99 = load i32, i32* %16
%100 = add nsw i32 %98, %99
store i32 %100, i32* %21
store i32 2, i32* %11
br label %101
101: ; preds = %213, %90
%102 = load i32, i32* %11
%103 = load i32, i32* %6
%104 = icmp slt i32 %102, %103
br i1 %104, label %105, label %216
105: ; preds = %101
%106 = load i32, i32* %18
%107 = add nsw i32 %106, 2
store i32 %107, i32* %18
%108 = load i32, i32* %19
%109 = sub nsw i32 %108, 2
store i32 %109, i32* %19
%110 = load i32, i32* %20
%111 = add nsw i32 %110, 2
store i32 %111, i32* %20
%112 = load i32, i32* %21
%113 = add nsw i32 %112, 2
store i32 %113, i32* %21
%114 = load double*, double** %10
%115 = load i32, i32* %11
%116 = sub nsw i32 %115, 2
%117 = sext i32 %116 to i64
%118 = getelementptr inbounds double, double* %114, i64 %117
%119 = load double, double* %118
%120 = load double*, double** %8
%121 = load i32, i32* %18
%122 = sub nsw i32 %121, 1
%123 = sext i32 %122 to i64
%124 = getelementptr inbounds double, double* %120, i64 %123
%125 = load double, double* %124
%126 = load double*, double** %10
%127 = load i32, i32* %11
%128 = sub nsw i32 %127, 1
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds double, double* %126, i64 %129
%131 = load double, double* %130
%132 = load double*, double** %8
%133 = load i32, i32* %18
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds double, double* %132, i64 %134
%136 = load double, double* %135
%137 = fmul double %131, %136
%138 = call double @llvm.fmuladd.f64(double %119, double %125, double %137)
store double %138, double* %14
%139 = load double*, double** %10
%140 = load i32, i32* %11
%141 = sub nsw i32 %140, 2
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds double, double* %139, i64 %142
%144 = load double, double* %143
%145 = load double*, double** %8
%146 = load i32, i32* %18
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds double, double* %145, i64 %147
%149 = load double, double* %148
%150 = load double*, double** %10
%151 = load i32, i32* %11
%152 = sub nsw i32 %151, 1
%153 = sext i32 %152 to i64
%154 = getelementptr inbounds double, double* %150, i64 %153
%155 = load double, double* %154
%156 = load double*, double** %8
%157 = load i32, i32* %18
%158 = sub nsw i32 %157, 1
%159 = sext i32 %158 to i64
%160 = getelementptr inbounds double, double* %156, i64 %159
%161 = load double, double* %160
%162 = fmul double %155, %161
%163 = fneg double %162
%164 = call double @llvm.fmuladd.f64(double %144, double %149, double %163)
store double %164, double* %13
%165 = load double*, double** %8
%166 = load i32, i32* %20
%167 = sext i32 %166 to i64
%168 = getelementptr inbounds double, double* %165, i64 %167
%169 = load double, double* %168
%170 = load double, double* %13
%171 = fadd double %169, %170
%172 = load double*, double** %9
%173 = load i32, i32* %21
%174 = sext i32 %173 to i64
%175 = getelementptr inbounds double, double* %172, i64 %174
store double %171, double* %175
%176 = load double, double* %13
%177 = load double*, double** %8
%178 = load i32, i32* %20
%179 = sext i32 %178 to i64
%180 = getelementptr inbounds double, double* %177, i64 %179
%181 = load double, double* %180
%182 = fsub double %176, %181
%183 = load double*, double** %9
%184 = load i32, i32* %19
%185 = sext i32 %184 to i64
%186 = getelementptr inbounds double, double* %183, i64 %185
store double %182, double* %186
%187 = load double*, double** %8
%188 = load i32, i32* %20
%189 = sub nsw i32 %188, 1
%190 = sext i32 %189 to i64
%191 = getelementptr inbounds double, double* %187, i64 %190
%192 = load double, double* %191
%193 = load double, double* %14
%194 = fadd double %192, %193
%195 = load double*, double** %9
%196 = load i32, i32* %21
%197 = sub nsw i32 %196, 1
%198 = sext i32 %197 to i64
%199 = getelementptr inbounds double, double* %195, i64 %198
store double %194, double* %199
%200 = load double*, double** %8
%201 = load i32, i32* %20
%202 = sub nsw i32 %201, 1
%203 = sext i32 %202 to i64
%204 = getelementptr inbounds double, double* %200, i64 %203
%205 = load double, double* %204
%206 = load double, double* %14
%207 = fsub double %205, %206
%208 = load double*, double** %9
%209 = load i32, i32* %19
%210 = sub nsw i32 %209, 1
%211 = sext i32 %210 to i64
%212 = getelementptr inbounds double, double* %208, i64 %211
store double %207, double* %212
br label %213
213: ; preds = %105
%214 = load i32, i32* %11
%215 = add nsw i32 %214, 2
store i32 %215, i32* %11
br label %101
216: ; preds = %101
%217 = load i32, i32* %6
%218 = load i32, i32* %16
%219 = add nsw i32 %218, %217
store i32 %219, i32* %16
%220 = load i32, i32* %6
%221 = load i32, i32* %17
%222 = add nsw i32 %221, %220
store i32 %222, i32* %17
br label %223
223: ; preds = %216
%224 = load i32, i32* %12
%225 = add nsw i32 %224, 1
store i32 %225, i32* %12
br label %86
226: ; preds = %86
%227 = load i32, i32* %6
%228 = srem i32 %227, 2
%229 = icmp eq i32 %228, 1
br i1 %229, label %230, label %231
230: ; preds = %226
br label %276
231: ; preds = %226
br label %232
232: ; preds = %231, %83
%233 = load i32, i32* %6
store i32 %233, i32* %16
%234 = sub nsw i32 %233, 1
store i32 %234, i32* %17
store i32 %234, i32* %18
%235 = load i32, i32* %15
%236 = load i32, i32* %17
%237 = add nsw i32 %236, %235
store i32 %237, i32* %17
store i32 0, i32* %12
br label %238
238: ; preds = %273, %232
%239 = load i32, i32* %12
%240 = load i32, i32* %7
%241 = icmp slt i32 %239, %240
br i1 %241, label %242, label %276
242: ; preds = %238
%243 = load double*, double** %8
%244 = load i32, i32* %17
%245 = sext i32 %244 to i64
%246 = getelementptr inbounds double, double* %243, i64 %245
%247 = load double, double* %246
%248 = fneg double %247
%249 = load double*, double** %9
%250 = load i32, i32* %16
%251 = sext i32 %250 to i64
%252 = getelementptr inbounds double, double* %249, i64 %251
store double %248, double* %252
%253 = load double*, double** %8
%254 = load i32, i32* %18
%255 = sext i32 %254 to i64
%256 = getelementptr inbounds double, double* %253, i64 %255
%257 = load double, double* %256
%258 = load double*, double** %9
%259 = load i32, i32* %16
%260 = sub nsw i32 %259, 1
%261 = sext i32 %260 to i64
%262 = getelementptr inbounds double, double* %258, i64 %261
store double %257, double* %262
%263 = load i32, i32* %6
%264 = shl i32 %263, 1
%265 = load i32, i32* %16
%266 = add nsw i32 %265, %264
store i32 %266, i32* %16
%267 = load i32, i32* %6
%268 = load i32, i32* %17
%269 = add nsw i32 %268, %267
store i32 %269, i32* %17
%270 = load i32, i32* %6
%271 = load i32, i32* %18
%272 = add nsw i32 %271, %270
store i32 %272, i32* %18
br label %273
273: ; preds = %242
%274 = load i32, i32* %12
%275 = add nsw i32 %274, 1
store i32 %275, i32* %12
br label %238
276: ; preds = %79, %230, %238
ret void
}

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