void icdf97(float *out, float *in, int N, int SO, int SI) {
  float a = -1.586134342f, b = -0.05298011854f, c = 0.8829110762f,
        d = 0.4435068522f, e = 1.149604398f;
  int M = N & ~1, K = N + (N & 1);
  for (int i = 0; i < M; i += 2) {
    out[(i + 0) * SO] = in[(i + 0) / 2 * SI] / e;
    out[(i + 1) * SO] = in[(i + K) / 2 * SI] * e;
  }
  if (N & 1)
    out[(N - 1) * SO] = in[(N - 1) / 2 * SI] / e;
  out[0] -= d * (2.f * out[SO]);
  for (int i = 2; i < M; i += 2)
    out[i * SO] -= d * (out[(i - 1) * SO] + out[(i + 1) * SO]);
  for (int i = 1; i < N - 1; i += 2)
    out[i * SO] -= c * (out[(i - 1) * SO] + out[(i + 1) * SO]);
  if (!(N & 1))
    out[(N - 1) * SO] -= c * (2.f * out[(N - 2) * SO]);
  out[0] -= b * (2.f * out[SO]);
  for (int i = 2; i < M; i += 2)
    out[i * SO] -= b * (out[(i - 1) * SO] + out[(i + 1) * SO]);
  for (int i = 1; i < N - 1; i += 2)
    out[i * SO] -= a * (out[(i - 1) * SO] + out[(i + 1) * SO]);
  if (!(N & 1))
    out[(N - 1) * SO] -= a * (2.f * out[(N - 2) * SO]);
}


x86 assembly (clang, -O0):
<icdf97>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
movss  0xef2(%rip),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  0xee1(%rip),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  0xed0(%rip),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  0xebf(%rip),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  0xeae(%rip),%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x14(%rbp),%eax
and    $0xfffffffe,%eax
mov    %eax,-0x34(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x14(%rbp),%ecx
and    $0x1,%ecx
add    %ecx,%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1215 <icdf97+0x115>
mov    -0x10(%rbp),%rax
mov    %rax,-0x60(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x0,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x60(%rbp),%rax
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
divss  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x3c(%rbp),%ecx
add    $0x0,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x3c(%rbp),%eax
add    -0x38(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x58(%rbp),%rax
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mulss  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x3c(%rbp)
jmp    1175 <icdf97+0x75>
mov    -0x14(%rbp),%eax
and    $0x1,%eax
cmp    $0x0,%eax
je     1267 <icdf97+0x167>
mov    -0x10(%rbp),%rax
mov    %rax,-0x68(%rbp)
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x68(%rbp),%rax
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
divss  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x2c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  0xd98(%rip),%xmm2
mulss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm1
movd   %xmm0,%ecx
xor    $0x80000000,%ecx
movd   %ecx,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
movl   $0x2,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1323 <icdf97+0x223>
movss  -0x2c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
add    $0x1,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movd   %xmm0,%edx
xor    $0x80000000,%edx
movd   %edx,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x40(%rbp)
jmp    12aa <icdf97+0x1aa>
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    13a8 <icdf97+0x2a8>
movss  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x44(%rbp),%ecx
add    $0x1,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x44(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movd   %xmm0,%edx
xor    $0x80000000,%edx
movd   %edx,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x44(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x44(%rbp)
jmp    132a <icdf97+0x22a>
mov    -0x14(%rbp),%eax
and    $0x1,%eax
cmp    $0x0,%eax
jne    140b <icdf97+0x30b>
movss  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x2,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  0xc3f(%rip),%xmm2
mulss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movd   %xmm0,%edx
xor    $0x80000000,%edx
movd   %edx,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x24(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  0xbf4(%rip),%xmm2
mulss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm1
movd   %xmm0,%ecx
xor    $0x80000000,%ecx
movd   %ecx,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
movl   $0x2,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    14c7 <icdf97+0x3c7>
movss  -0x24(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    $0x1,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movd   %xmm0,%edx
xor    $0x80000000,%edx
movd   %edx,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x48(%rbp)
jmp    144e <icdf97+0x34e>
movl   $0x1,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    154c <icdf97+0x44c>
movss  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x4c(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    $0x1,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x4c(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movd   %xmm0,%edx
xor    $0x80000000,%edx
movd   %edx,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x4c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x4c(%rbp)
jmp    14ce <icdf97+0x3ce>
mov    -0x14(%rbp),%eax
and    $0x1,%eax
cmp    $0x0,%eax
jne    15af <icdf97+0x4af>
movss  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x2,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  0xa9b(%rip),%xmm2
mulss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movd   %xmm0,%edx
xor    $0x80000000,%edx
movd   %edx,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
pop    %rbp
ret

x86 assembly (clang, -O2):
<icdf97>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %edx,%r11d
and    $0xfffffffe,%r11d
mov    %edx,%eax
and    $0x1,%eax
mov    %rax,-0x8(%rsp)
mov    %ecx,-0xc(%rsp)
movslq %ecx,%r13
mov    %r11d,%r10d
test   %r11d,%r11d
jle    11a4 <icdf97+0xa4>
movslq %r8d,%r15
mov    -0x8(%rsp),%rax
lea    (%rdx,%rax,1),%r14d
lea    0x0(,%r13,8),%r12
shl    $0x2,%r15
xor    %ebx,%ebx
movss  0xeb4(%rip),%xmm0
mov    %rsi,%rax
mov    %rdi,%rbp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%rax),%xmm1
divss  %xmm0,%xmm1
movss  %xmm1,0x0(%rbp)
lea    (%r14,%rbx,1),%ecx
mov    %ecx,%r9d
shr    $0x1f,%r9d
add    %ecx,%r9d
sar    %r9d
imul   %r8d,%r9d
movslq %r9d,%rcx
movss  (%rsi,%rcx,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,0x0(%rbp,%r13,4)
add    $0x2,%rbx
add    %r12,%rbp
add    %r15,%rax
cmp    %r10,%rbx
jb     1160 <icdf97+0x60>
mov    -0x8(%rsp),%r15
test   %r15d,%r15d
mov    -0xc(%rsp),%r9d
je     11e1 <icdf97+0xe1>
lea    -0x1(%rdx),%eax
mov    %eax,%ecx
shr    $0x1f,%ecx
add    %edx,%ecx
add    $0xffffffff,%ecx
sar    %ecx
imul   %r8d,%ecx
movslq %ecx,%rcx
movss  (%rsi,%rcx,4),%xmm0
divss  0xe2a(%rip),%xmm0
imul   %r9d,%eax
cltq
movss  %xmm0,(%rdi,%rax,4)
movss  (%rdi,%r13,4),%xmm0
addss  %xmm0,%xmm0
mulss  0xe11(%rip),%xmm0
addss  (%rdi),%xmm0
movss  %xmm0,(%rdi)
cmp    $0x3,%r11d
jl     1257 <icdf97+0x157>
lea    0x0(,%r13,4),%rax
lea    (%rax,%rax,2),%rax
lea    0x0(,%r13,8),%rsi
mov    $0x2,%ebp
movss  0xde2(%rip),%xmm0
mov    %rdi,%rbx
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%rbx,%r13,4),%xmm1
addss  (%rbx,%rax,1),%xmm1
mulss  %xmm0,%xmm1
addss  (%rbx,%r13,8),%xmm1
movss  %xmm1,(%rbx,%r13,8)
add    $0x2,%rbp
add    %rsi,%rbx
cmp    %r10,%rbp
jb     1230 <icdf97+0x130>
lea    -0x1(%rdx),%r8d
mov    %r8d,%r14d
cmp    $0x3,%edx
jl     12aa <icdf97+0x1aa>
lea    0x0(,%r13,8),%rax
mov    $0x1,%ebx
movss  0xd90(%rip),%xmm0
mov    %rdi,%rbp
nopl   0x0(%rax,%rax,1)
movss  0x0(%rbp),%xmm1
addss  0x0(%rbp,%r13,8),%xmm1
mulss  %xmm0,%xmm1
addss  0x0(%rbp,%r13,4),%xmm1
movss  %xmm1,0x0(%rbp,%r13,4)
add    $0x2,%rbx
add    %rax,%rbp
cmp    %r14,%rbx
jb     1280 <icdf97+0x180>
test   %r15d,%r15d
jne    12dc <icdf97+0x1dc>
lea    -0x2(%rdx),%eax
imul   %r9d,%eax
cltq
movss  (%rdi,%rax,4),%xmm0
addss  %xmm0,%xmm0
mov    %r8d,%eax
imul   %r9d,%eax
mulss  0xd38(%rip),%xmm0
cltq
addss  (%rdi,%rax,4),%xmm0
movss  %xmm0,(%rdi,%rax,4)
movss  (%rdi,%r13,4),%xmm0
addss  %xmm0,%xmm0
mulss  0xd1e(%rip),%xmm0
addss  (%rdi),%xmm0
movss  %xmm0,(%rdi)
cmp    $0x3,%r11d
jl     1347 <icdf97+0x247>
lea    0x0(,%r13,4),%rax
lea    (%rax,%rax,2),%rax
lea    0x0(,%r13,8),%rbp
mov    $0x2,%ebx
movss  0xcef(%rip),%xmm0
mov    %rdi,%rsi
movss  (%rsi,%r13,4),%xmm1
addss  (%rsi,%rax,1),%xmm1
mulss  %xmm0,%xmm1
addss  (%rsi,%r13,8),%xmm1
movss  %xmm1,(%rsi,%r13,8)
add    $0x2,%rbx
add    %rbp,%rsi
cmp    %r10,%rbx
jb     1320 <icdf97+0x220>
cmp    $0x3,%edx
jl     1396 <icdf97+0x296>
lea    0x0(,%r13,8),%rax
mov    $0x1,%ebp
movss  0xcaf(%rip),%xmm0
mov    %rdi,%rsi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rsi),%xmm1
addss  (%rsi,%r13,8),%xmm1
mulss  %xmm0,%xmm1
addss  (%rsi,%r13,4),%xmm1
movss  %xmm1,(%rsi,%r13,4)
add    $0x2,%rbp
add    %rax,%rsi
cmp    %r14,%rbp
jb     1370 <icdf97+0x270>
test   %r15d,%r15d
jne    13c7 <icdf97+0x2c7>
add    $0xfffffffe,%edx
imul   %r9d,%edx
movslq %edx,%rax
movss  (%rdi,%rax,4),%xmm0
addss  %xmm0,%xmm0
imul   %r9d,%r8d
mulss  0xc56(%rip),%xmm0
movslq %r8d,%rax
addss  (%rdi,%rax,4),%xmm0
movss  %xmm0,(%rdi,%rax,4)
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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