void radf2(int ido, int l1, const float cc[], float ch[], const float wa1[]) {
  int i, k, ic;
  float 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
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x8(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x8(%rbp),%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%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
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    111d <radf2+0x1d>
cmpl   $0x2,-0x4(%rbp)
jge    11cd <radf2+0xcd>
jmp    14b6 <radf2+0x3b6>
cmpl   $0x2,-0x4(%rbp)
je     1418 <radf2+0x318>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    13fa <radf2+0x2fa>
movl   $0x2,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    13e7 <radf2+0x2e7>
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
movss  (%rax,%rcx,4),%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
movss  (%rax,%rcx,4),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%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
mulss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%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
movss  (%rax,%rcx,4),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%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
mulss  (%rax,%rcx,4),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x28(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x30(%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
movss  %xmm0,(%rax,%rcx,4)
movss  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x28(%rbp),%edx
imul   -0x4(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%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
movss  %xmm0,(%rax,%rcx,4)
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
movss  (%rax,%rcx,4),%xmm0
addss  -0x34(%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
movss  %xmm0,(%rax,%rcx,4)
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
movss  (%rax,%rcx,4),%xmm0
subss  -0x34(%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
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x24(%rbp)
jmp    11f1 <radf2+0xf1>
jmp    13ec <radf2+0x2ec>
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    1413 <radf2+0x313>
jmp    14b6 <radf2+0x3b6>
jmp    1418 <radf2+0x318>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    14b6 <radf2+0x3b6>
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
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
imul   -0x4(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
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
movss  (%rax,%rcx,4),%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
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    141f <radf2+0x31f>
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,4),%r9
mov    %esi,%r11d
lea    (%rsi,%rsi,1),%ebp
add    $0xffffffff,%ebp
imul   %rax,%rsi
lea    (%rdx,%rsi,4),%rbx
add    %r11d,%r11d
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%rdx,%rsi,1),%xmm0
addss  (%rbx,%rsi,1),%xmm0
movss  %xmm0,(%rcx,%rsi,2)
movss  (%rdx,%rsi,1),%xmm0
subss  (%rbx,%rsi,1),%xmm0
movslq %ebp,%rbp
movss  %xmm0,(%rcx,%rbp,4)
add    %r9,%rsi
add    %r11d,%ebp
add    $0xffffffffffffffff,%rax
jne    1140 <radf2+0x40>
cmp    $0x2,%edi
jl     1331 <radf2+0x231>
je     12d1 <radf2+0x1d1>
test   %r14d,%r14d
jle    1331 <radf2+0x231>
mov    %edi,%r10d
mov    %r14d,-0x28(%rsp)
mov    %r14d,%esi
lea    (%rdi,%rdi,1),%eax
mov    %eax,-0x24(%rsp)
lea    0x0(,%r10,4),%rbp
mov    %rbp,-0x10(%rsp)
mov    %rsi,-0x8(%rsp)
imul   %r10,%rsi
lea    (%rdx,%rsi,4),%rbx
lea    (%rdi,%rdi,1),%r15d
add    $0xfffffffd,%r15d
lea    -0x8(%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,4),%rsi
movslq %r13d,%rdi
add    %r10,%rdi
mov    -0x18(%rsp),%rax
lea    (%rax,%rdi,4),%rdi
mov    $0x2,%r14d
mov    %r15d,%r11d
nop
movss  -0x8(%r8,%r14,4),%xmm0
movss  -0x4(%r8,%r14,4),%xmm1
movss  -0x4(%rbx,%r14,4),%xmm2
movss  (%rbx,%r14,4),%xmm3
movaps %xmm0,%xmm4
mulss  %xmm2,%xmm4
mulss  %xmm1,%xmm2
mulss  %xmm3,%xmm1
addss  %xmm1,%xmm4
mulss  %xmm3,%xmm0
subss  %xmm2,%xmm0
movss  0x0(%rbp,%r14,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rsi,%r14,4)
subss  0x0(%rbp,%r14,4),%xmm0
movss  %xmm0,(%rdi)
movss  -0x4(%rbp,%r14,4),%xmm0
addss  %xmm4,%xmm0
movss  %xmm0,-0x4(%rsi,%r14,4)
movss  -0x4(%rbp,%r14,4),%xmm0
subss  %xmm4,%xmm0
movslq %r11d,%r11
movss  %xmm0,(%rcx,%r11,4)
add    $0x2,%r14
add    $0xfffffffe,%r11d
add    $0xfffffffffffffff8,%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     1331 <radf2+0x231>
test   %r14d,%r14d
jle    1331 <radf2+0x231>
movslq %edi,%rdi
mov    %r14d,%eax
lea    (%rdx,%rdi,4),%rdx
add    $0xfffffffffffffffc,%rdx
lea    0x0(,%rdi,4),%rsi
lea    (%rcx,%rdi,4),%rcx
imul   %rax,%rdi
lea    (%rdx,%rdi,4),%rdi
xor    %ebp,%ebp
movaps 0xcff(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  (%rdi,%rbp,1),%xmm1
xorps  %xmm0,%xmm1
movss  %xmm1,(%rcx,%rbp,2)
movss  (%rdx,%rbp,1),%xmm1
movss  %xmm1,-0x4(%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, float* %2, float* %3, float* %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca float*
%9 = alloca float*
%10 = alloca float*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca float
%15 = alloca float
store i32 %0, i32* %6
store i32 %1, i32* %7
store float* %2, float** %8
store float* %3, float** %9
store float* %4, float** %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 float*, float** %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 float, float* %21, i64 %25
%27 = load float, float* %26
%28 = load float*, float** %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 float, float* %28, i64 %34
%36 = load float, float* %35
%37 = fadd float %27, %36
%38 = load float*, float** %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 float, float* %38, i64 %43
store float %37, float* %44
%45 = load float*, float** %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 float, float* %45, i64 %49
%51 = load float, float* %50
%52 = load float*, float** %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 float, float* %52, i64 %58
%60 = load float, float* %59
%61 = fsub float %51, %60
%62 = load float*, float** %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 float, float* %62, i64 %71
store float %61, float* %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 float*, float** %10
%98 = load i32, i32* %11
%99 = sub nsw i32 %98, 2
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds float, float* %97, i64 %100
%102 = load float, float* %101
%103 = load float*, float** %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 float, float* %103, i64 %112
%114 = load float, float* %113
%115 = load float*, float** %10
%116 = load i32, i32* %11
%117 = sub nsw i32 %116, 1
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds float, float* %115, i64 %118
%120 = load float, float* %119
%121 = load float*, float** %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 float, float* %121, i64 %129
%131 = load float, float* %130
%132 = fmul float %120, %131
%133 = call float @llvm.fmuladd.f32(float %102, float %114, float %132)
store float %133, float* %15
%134 = load float*, float** %10
%135 = load i32, i32* %11
%136 = sub nsw i32 %135, 2
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds float, float* %134, i64 %137
%139 = load float, float* %138
%140 = load float*, float** %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 float, float* %140, i64 %148
%150 = load float, float* %149
%151 = load float*, float** %10
%152 = load i32, i32* %11
%153 = sub nsw i32 %152, 1
%154 = sext i32 %153 to i64
%155 = getelementptr inbounds float, float* %151, i64 %154
%156 = load float, float* %155
%157 = load float*, float** %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 float, float* %157, i64 %166
%168 = load float, float* %167
%169 = fmul float %156, %168
%170 = fneg float %169
%171 = call float @llvm.fmuladd.f32(float %139, float %150, float %170)
store float %171, float* %14
%172 = load float*, float** %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 float, float* %172, i64 %178
%180 = load float, float* %179
%181 = load float, float* %14
%182 = fadd float %180, %181
%183 = load float*, float** %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 float, float* %183, i64 %190
store float %182, float* %191
%192 = load float, float* %14
%193 = load float*, float** %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 float, float* %193, i64 %199
%201 = load float, float* %200
%202 = fsub float %192, %201
%203 = load float*, float** %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 float, float* %203, i64 %211
store float %202, float* %212
%213 = load float*, float** %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 float, float* %213, i64 %220
%222 = load float, float* %221
%223 = load float, float* %15
%224 = fadd float %222, %223
%225 = load float*, float** %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 float, float* %225, i64 %233
store float %224, float* %234
%235 = load float*, float** %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 float, float* %235, i64 %242
%244 = load float, float* %243
%245 = load float, float* %15
%246 = fsub float %244, %245
%247 = load float*, float** %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 float, float* %247, i64 %256
store float %246, float* %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 float*, float** %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 float, float* %277, i64 %286
%288 = load float, float* %287
%289 = fneg float %288
%290 = load float*, float** %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 float, float* %290, i64 %296
store float %289, float* %297
%298 = load float*, float** %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 float, float* %298, i64 %305
%307 = load float, float* %306
%308 = load float*, float** %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 float, float* %308, i64 %316
store float %307, float* %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, float* nocapture %2, float* nocapture writeonly %3, float* 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 float, float* %2, i64 %14
%16 = load float, float* %15
%17 = add nuw nsw i64 %13, %10
%18 = mul nsw i64 %17, %9
%19 = getelementptr inbounds float, float* %2, i64 %18
%20 = load float, float* %19
%21 = fadd float %16, %20
%22 = shl nuw nsw i64 %13, 1
%23 = mul nsw i64 %22, %9
%24 = getelementptr inbounds float, float* %3, i64 %23
store float %21, float* %24
%25 = load float, float* %15
%26 = load float, float* %19
%27 = fsub float %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 float, float* %3, i64 %32
store float %27, float* %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 float, float* %4, i64 %61
%63 = load float, float* %62
%64 = add nsw i64 %59, -1
%65 = add nsw i64 %64, %48
%66 = getelementptr inbounds float, float* %2, i64 %65
%67 = load float, float* %66
%68 = getelementptr inbounds float, float* %4, i64 %64
%69 = load float, float* %68
%70 = add nuw nsw i64 %59, %48
%71 = getelementptr inbounds float, float* %2, i64 %70
%72 = load float, float* %71
%73 = fmul float %69, %72
%74 = tail call float @llvm.fmuladd.f32(float %63, float %67, float %73)
%75 = fneg float %69
%76 = fmul float %67, %75
%77 = tail call float @llvm.fmuladd.f32(float %63, float %72, float %76)
%78 = add nuw nsw i64 %59, %49
%79 = getelementptr inbounds float, float* %2, i64 %78
%80 = load float, float* %79
%81 = fadd float %80, %77
%82 = add nsw i64 %59, %56
%83 = getelementptr inbounds float, float* %3, i64 %82
store float %81, float* %83
%84 = load float, float* %79
%85 = fsub float %77, %84
%86 = add nsw i64 %60, %57
%87 = getelementptr inbounds float, float* %3, i64 %86
store float %85, float* %87
%88 = add nsw i64 %64, %49
%89 = getelementptr inbounds float, float* %2, i64 %88
%90 = load float, float* %89
%91 = fadd float %74, %90
%92 = add nsw i64 %64, %56
%93 = getelementptr inbounds float, float* %3, i64 %92
store float %91, float* %93
%94 = load float, float* %89
%95 = fsub float %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 float, float* %3, i64 %99
store float %95, float* %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 float, float* %2, i64 %121
%123 = load float, float* %122
%124 = fneg float %123
%125 = shl nuw nsw i64 %118, 1
%126 = or i64 %125, 1
%127 = mul nsw i64 %126, %113
%128 = getelementptr inbounds float, float* %3, i64 %127
store float %124, float* %128
%129 = mul nsw i64 %118, %113
%130 = add nsw i64 %129, %114
%131 = getelementptr inbounds float, float* %2, i64 %130
%132 = load float, float* %131
%133 = mul nsw i64 %125, %113
%134 = add nsw i64 %133, %114
%135 = getelementptr inbounds float, float* %3, i64 %134
store float %132, float* %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
}