void radf2(int ido, int l1, const double cc[], double ch[],
           const double wa1[]) {
  int i, k, ic;
  double ti2, tr2;
  for (k = 0; k < l1; k++) {
    ch[2 * k * ido] = cc[k * ido] + cc[(k + l1) * ido];
    ch[(2 * k + 1) * ido + ido - 1] = cc[k * ido] - cc[(k + l1) * ido];
  }
  if (ido < 2)
    return;
  if (ido != 2) {
    for (k = 0; k < l1; k++) {
      for (i = 2; i < ido; i += 2) {
        ic = ido - i;
        tr2 = wa1[i - 2] * cc[i - 1 + (k + l1) * ido] +
              wa1[i - 1] * cc[i + (k + l1) * ido];
        ti2 = wa1[i - 2] * cc[i + (k + l1) * ido] -
              wa1[i - 1] * cc[i - 1 + (k + l1) * ido];
        ch[i + 2 * k * ido] = cc[i + k * ido] + ti2;
        ch[ic + (2 * k + 1) * ido] = ti2 - cc[i + k * ido];
        ch[i - 1 + 2 * k * ido] = cc[i - 1 + k * ido] + tr2;
        ch[ic - 1 + (2 * k + 1) * ido] = cc[i - 1 + k * ido] - tr2;
      }
    }
    if (ido % 2 == 1)
      return;
  }
  for (k = 0; k < l1; k++) {
    ch[(2 * k + 1) * ido] = -cc[ido - 1 + (k + l1) * ido];
    ch[ido - 1 + 2 * k * ido] = cc[ido - 1 + k * ido];
  }
}


x86 assembly (clang, -O0):
<radf2>:
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,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11be <radf2+0xbe>
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x8(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x8(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
imul   -0x4(%rbp),%ecx
add    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    111d <radf2+0x1d>
cmpl   $0x2,-0x4(%rbp)
jge    11cd <radf2+0xcd>
jmp    14ca <radf2+0x3ca>
cmpl   $0x2,-0x4(%rbp)
je     1422 <radf2+0x322>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1404 <radf2+0x304>
movl   $0x2,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    13f1 <radf2+0x2f1>
mov    -0x4(%rbp),%eax
sub    -0x24(%rbp),%eax
mov    %eax,-0x2c(%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    -0x24(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
add    -0x8(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%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
mov    -0x24(%rbp),%ecx
mov    -0x28(%rbp),%edx
add    -0x8(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%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    -0x24(%rbp),%ecx
mov    -0x28(%rbp),%edx
add    -0x8(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%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
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
add    -0x8(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%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,-0x38(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x28(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x28(%rbp),%edx
shl    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x28(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x28(%rbp),%edx
shl    $0x1,%edx
add    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
shl    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x40(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
shl    $0x1,%edx
add    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x24(%rbp)
jmp    11f1 <radf2+0xf1>
jmp    13f6 <radf2+0x2f6>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11de <radf2+0xde>
mov    -0x4(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x1,%edx
jne    141d <radf2+0x31d>
jmp    14ca <radf2+0x3ca>
jmp    1422 <radf2+0x322>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    14ca <radf2+0x3ca>
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
add    -0x8(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x28(%rbp),%edx
shl    $0x1,%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1429 <radf2+0x329>
pop    %rbp
ret

x86 assembly (clang, -O2):
<radf2>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %esi,%r14d
test   %esi,%esi
jle    116d <radf2+0x6d>
movslq %edi,%rsi
mov    %r14d,%eax
lea    0x0(,%rsi,8),%r9
mov    %esi,%r11d
lea    (%rsi,%rsi,1),%ebp
add    $0xffffffff,%ebp
imul   %rax,%rsi
lea    (%rdx,%rsi,8),%rbx
add    %r11d,%r11d
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rdx,%rsi,1),%xmm0
addsd  (%rbx,%rsi,1),%xmm0
movsd  %xmm0,(%rcx,%rsi,2)
movsd  (%rdx,%rsi,1),%xmm0
subsd  (%rbx,%rsi,1),%xmm0
movslq %ebp,%rbp
movsd  %xmm0,(%rcx,%rbp,8)
add    %r9,%rsi
add    %r11d,%ebp
add    $0xffffffffffffffff,%rax
jne    1140 <radf2+0x40>
cmp    $0x2,%edi
jl     1332 <radf2+0x232>
je     12d2 <radf2+0x1d2>
test   %r14d,%r14d
jle    1332 <radf2+0x232>
mov    %edi,%r10d
mov    %r14d,-0x28(%rsp)
mov    %r14d,%esi
lea    (%rdi,%rdi,1),%eax
mov    %eax,-0x24(%rsp)
lea    0x0(,%r10,8),%rbp
mov    %rbp,-0x10(%rsp)
mov    %rsi,-0x8(%rsp)
imul   %r10,%rsi
lea    (%rdx,%rsi,8),%rbx
lea    (%rdi,%rdi,1),%r15d
add    $0xfffffffd,%r15d
lea    -0x10(%rcx),%rsi
mov    %rsi,-0x18(%rsp)
xor    %r12d,%r12d
mov    %rdi,-0x20(%rsp)
mov    %edi,%r13d
mov    %rdx,%rbp
xor    %r9d,%r9d
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movslq %r12d,%rsi
lea    (%rcx,%rsi,8),%rsi
movslq %r13d,%rdi
add    %r10,%rdi
mov    -0x18(%rsp),%rax
lea    (%rax,%rdi,8),%rdi
mov    $0x2,%r14d
mov    %r15d,%r11d
nop
movsd  -0x10(%r8,%r14,8),%xmm0
movsd  -0x8(%r8,%r14,8),%xmm1
movsd  -0x8(%rbx,%r14,8),%xmm2
movsd  (%rbx,%r14,8),%xmm3
movapd %xmm0,%xmm4
mulsd  %xmm2,%xmm4
mulsd  %xmm1,%xmm2
mulsd  %xmm3,%xmm1
addsd  %xmm1,%xmm4
mulsd  %xmm3,%xmm0
subsd  %xmm2,%xmm0
movsd  0x0(%rbp,%r14,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsi,%r14,8)
subsd  0x0(%rbp,%r14,8),%xmm0
movsd  %xmm0,(%rdi)
movsd  -0x8(%rbp,%r14,8),%xmm0
addsd  %xmm4,%xmm0
movsd  %xmm0,-0x8(%rsi,%r14,8)
movsd  -0x8(%rbp,%r14,8),%xmm0
subsd  %xmm4,%xmm0
movslq %r11d,%r11
movsd  %xmm0,(%rcx,%r11,8)
add    $0x2,%r14
add    $0xfffffffe,%r11d
add    $0xfffffffffffffff0,%rdi
cmp    %r10,%r14
jb     1200 <radf2+0x100>
add    $0x1,%r9
mov    -0x24(%rsp),%eax
add    %eax,%r12d
mov    -0x10(%rsp),%rsi
add    %rsi,%rbp
add    %rsi,%rbx
add    %eax,%r15d
add    %eax,%r13d
cmp    -0x8(%rsp),%r9
jne    11e0 <radf2+0xe0>
mov    -0x20(%rsp),%rdi
mov    %edi,%eax
shr    $0x1f,%eax
add    %edi,%eax
and    $0xfffffffe,%eax
mov    %edi,%esi
sub    %eax,%esi
cmp    $0x1,%esi
mov    -0x28(%rsp),%r14d
je     1332 <radf2+0x232>
test   %r14d,%r14d
jle    1332 <radf2+0x232>
movslq %edi,%rdi
mov    %r14d,%eax
lea    (%rdx,%rdi,8),%rdx
add    $0xfffffffffffffff8,%rdx
lea    0x0(,%rdi,8),%rsi
lea    (%rcx,%rdi,8),%rcx
imul   %rax,%rdi
lea    (%rdx,%rdi,8),%rdi
xor    %ebp,%ebp
movapd 0xcfd(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rdi,%rbp,1),%xmm1
xorpd  %xmm0,%xmm1
movlpd %xmm1,(%rcx,%rbp,2)
movsd  (%rdx,%rbp,1),%xmm1
movsd  %xmm1,-0x8(%rcx,%rbp,2)
add    %rsi,%rbp
add    $0xffffffffffffffff,%rax
jne    1310 <radf2+0x210>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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