long scd_parallel_form_iir_up_kernel(lenx, x, y, iup, nblocks, direct_cof, gain,
                                     b, c, T)
long lenx;
float *x, *y;
long iup, nblocks;
double direct_cof, gain;
float (*b)[3], (*c)[2], (*T)[2];
{
  long kx, ky, n;
  float Ttmp;
  kx = 0;
  for (ky = 0; ky < iup * lenx; ky++) {
    if (ky % iup == 0) {
      y[ky] = direct_cof * x[kx];
      for (n = 0; n < nblocks; n++) {
        Ttmp = 2. * (x[kx] - c[n][0] * T[n][0] - c[n][1] * T[n][1]);
        y[ky] += b[n][2] * Ttmp + b[n][1] * T[n][1] + b[n][0] * T[n][0];
        T[n][0] = T[n][1];
        T[n][1] = Ttmp;
      }
      y[ky] *= gain;
      kx++;
    } else {
      y[ky] = 0.0;
      for (n = 0; n < nblocks; n++) {
        Ttmp = 2. * (0.0 - c[n][0] * T[n][0] - c[n][1] * T[n][1]);
        y[ky] += b[n][2] * Ttmp + b[n][1] * T[n][1] + b[n][0] * T[n][0];
        T[n][0] = T[n][1];
        T[n][1] = Ttmp;
      }
      y[ky] *= gain;
    }
  }
  return ky;
}


x86 assembly (clang, -O0):
<scd_parallel_form_iir_up_kernel>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movsd  %xmm0,-0x30(%rbp)
movsd  %xmm1,-0x38(%rbp)
mov    %r9,-0x40(%rbp)
movq   $0x0,-0x48(%rbp)
movq   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%rax
mov    -0x20(%rbp),%rcx
imul   -0x8(%rbp),%rcx
cmp    %rcx,%rax
jge    1486 <scd_parallel_form_iir_up_kernel+0x386>
mov    -0x50(%rbp),%rax
cqto
idivq  -0x20(%rbp)
cmp    $0x0,%rdx
jne    1316 <scd_parallel_form_iir_up_kernel+0x216>
movsd  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movq   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    12e2 <scd_parallel_form_iir_up_kernel+0x1e2>
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x10(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  (%rax,%rcx,8),%xmm1
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  (%rax,%rcx,8),%xmm2
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    0x10(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  0x4(%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  0x4(%rax,%rcx,8),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdde(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x5c(%rbp)
mov    -0x40(%rbp),%rax
imul   $0xc,-0x58(%rbp),%rcx
add    %rcx,%rax
movss  0x8(%rax),%xmm1
movss  -0x5c(%rbp),%xmm2
mov    -0x40(%rbp),%rax
imul   $0xc,-0x58(%rbp),%rcx
add    %rcx,%rax
movss  0x4(%rax),%xmm0
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
mulss  0x4(%rax,%rcx,8),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x40(%rbp),%rax
imul   $0xc,-0x58(%rbp),%rcx
add    %rcx,%rax
movss  (%rax),%xmm0
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  (%rax,%rcx,8),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  0x4(%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,8)
movss  -0x5c(%rbp),%xmm0
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  %xmm0,0x4(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x58(%rbp)
jmp    119b <scd_parallel_form_iir_up_kernel+0x9b>
movsd  -0x38(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    1470 <scd_parallel_form_iir_up_kernel+0x370>
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movq   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    144d <scd_parallel_form_iir_up_kernel+0x34d>
mov    0x10(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
mulss  (%rax,%rcx,8),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
subsd  %xmm0,%xmm1
mov    0x10(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  0x4(%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
mulss  0x4(%rax,%rcx,8),%xmm0
cvtss2sd %xmm0,%xmm0
subsd  %xmm0,%xmm1
movsd  0xc73(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x5c(%rbp)
mov    -0x40(%rbp),%rax
imul   $0xc,-0x58(%rbp),%rcx
add    %rcx,%rax
movss  0x8(%rax),%xmm1
movss  -0x5c(%rbp),%xmm2
mov    -0x40(%rbp),%rax
imul   $0xc,-0x58(%rbp),%rcx
add    %rcx,%rax
movss  0x4(%rax),%xmm0
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
mulss  0x4(%rax,%rcx,8),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x40(%rbp),%rax
imul   $0xc,-0x58(%rbp),%rcx
add    %rcx,%rax
movss  (%rax),%xmm0
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  (%rax,%rcx,8),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  0x4(%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,8)
movss  -0x5c(%rbp),%xmm0
mov    0x18(%rbp),%rax
mov    -0x58(%rbp),%rcx
movss  %xmm0,0x4(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x58(%rbp)
jmp    132e <scd_parallel_form_iir_up_kernel+0x22e>
movsd  -0x38(%rbp),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    1475 <scd_parallel_form_iir_up_kernel+0x375>
mov    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    113e <scd_parallel_form_iir_up_kernel+0x3e>
mov    -0x50(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<scd_parallel_form_iir_up_kernel>:
mov    %rdi,%r11
imul   %rcx,%r11
test   %r11,%r11
jle    12fe <scd_parallel_form_iir_up_kernel+0x1fe>
push   %r15
push   %r14
push   %rbx
mov    %rdx,%r10
mov    0x28(%rsp),%rbx
mov    0x20(%rsp),%r14
add    $0x8,%r9
xor    %edi,%edi
xor    %r15d,%r15d
jmp    114f <scd_parallel_form_iir_up_kernel+0x4f>
nopl   (%rax)
cvtss2sd %xmm2,%xmm2
mulsd  %xmm1,%xmm2
cvtsd2ss %xmm2,%xmm2
movss  %xmm2,(%r10,%rdi,4)
add    $0x1,%rdi
cmp    %r11,%rdi
je     12f5 <scd_parallel_form_iir_up_kernel+0x1f5>
mov    %rdi,%rax
or     %rcx,%rax
shr    $0x20,%rax
je     1220 <scd_parallel_form_iir_up_kernel+0x120>
mov    %rdi,%rax
cqto
idiv   %rcx
test   %rdx,%rdx
je     122f <scd_parallel_form_iir_up_kernel+0x12f>
movl   $0x0,(%r10,%rdi,4)
xorps  %xmm2,%xmm2
test   %r8,%r8
jle    1130 <scd_parallel_form_iir_up_kernel+0x30>
mov    %r9,%rax
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%rbx,%rdx,8),%xmm3
movss  0x4(%rbx,%rdx,8),%xmm4
movss  (%r14,%rdx,8),%xmm5
mulss  %xmm3,%xmm5
cvtss2sd %xmm5,%xmm5
xorpd  %xmm6,%xmm6
subsd  %xmm5,%xmm6
movss  0x4(%r14,%rdx,8),%xmm5
mulss  %xmm4,%xmm5
cvtss2sd %xmm5,%xmm5
subsd  %xmm5,%xmm6
addsd  %xmm6,%xmm6
xorps  %xmm5,%xmm5
cvtsd2ss %xmm6,%xmm5
mulss  -0x4(%rax),%xmm4
movss  (%rax),%xmm6
mulss  %xmm5,%xmm6
addss  %xmm4,%xmm6
mulss  -0x8(%rax),%xmm3
addss  %xmm6,%xmm3
addss  %xmm2,%xmm3
movss  %xmm3,(%r10,%rdi,4)
movss  0x4(%rbx,%rdx,8),%xmm2
movss  %xmm2,(%rbx,%rdx,8)
movss  %xmm5,0x4(%rbx,%rdx,8)
add    $0x1,%rdx
movss  (%r10,%rdi,4),%xmm2
add    $0xc,%rax
cmp    %rdx,%r8
jne    1190 <scd_parallel_form_iir_up_kernel+0x90>
jmp    1130 <scd_parallel_form_iir_up_kernel+0x30>
mov    %edi,%eax
xor    %edx,%edx
div    %ecx
test   %rdx,%rdx
jne    1170 <scd_parallel_form_iir_up_kernel+0x70>
movss  (%rsi,%r15,4),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm0,%xmm2
cvtsd2ss %xmm2,%xmm2
movss  %xmm2,(%r10,%rdi,4)
test   %r8,%r8
jle    12da <scd_parallel_form_iir_up_kernel+0x1da>
mov    %r9,%rax
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%rsi,%r15,4),%xmm3
movss  (%rbx,%rdx,8),%xmm4
movss  0x4(%rbx,%rdx,8),%xmm5
movss  (%r14,%rdx,8),%xmm6
mulss  %xmm4,%xmm6
subss  %xmm6,%xmm3
movss  0x4(%r14,%rdx,8),%xmm6
mulss  %xmm5,%xmm6
subss  %xmm6,%xmm3
mulss  -0x4(%rax),%xmm5
addss  %xmm3,%xmm3
movss  (%rax),%xmm6
mulss  %xmm3,%xmm6
addss  %xmm5,%xmm6
mulss  -0x8(%rax),%xmm4
addss  %xmm6,%xmm4
addss  %xmm2,%xmm4
movss  %xmm4,(%r10,%rdi,4)
movss  0x4(%rbx,%rdx,8),%xmm2
movss  %xmm2,(%rbx,%rdx,8)
movss  %xmm3,0x4(%rbx,%rdx,8)
add    $0x1,%rdx
movss  (%r10,%rdi,4),%xmm2
add    $0xc,%rax
cmp    %rdx,%r8
jne    1260 <scd_parallel_form_iir_up_kernel+0x160>
cvtss2sd %xmm2,%xmm2
mulsd  %xmm1,%xmm2
cvtsd2ss %xmm2,%xmm2
movss  %xmm2,(%r10,%rdi,4)
add    $0x1,%r15
jmp    1142 <scd_parallel_form_iir_up_kernel+0x42>
pop    %rbx
pop    %r14
pop    %r15
mov    %r11,%rax
ret
xor    %r11d,%r11d
mov    %r11,%rax
ret

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

LLVM IR (-O2):
define i64 @scd_parallel_form_iir_up_kernel(i64 %0, float* nocapture %1, float* nocapture %2, i64 %3, i64 %4, double %5, double %6, [3 x float]* nocapture %7, [2 x float]* nocapture %8, [2 x float]* nocapture %9) {
%11 = mul i64 %3, %0
%12 = icmp sgt i64 %11, 0
br i1 %12, label %13, label %108
13: ; preds = %10
%14 = icmp sgt i64 %4, 0
br label %15
15: ; preds = %13, %104
%16 = phi i64 [ 0, %13 ], [ %106, %104 ]
%17 = phi i64 [ 0, %13 ], [ %105, %104 ]
%18 = srem i64 %16, %3
%19 = icmp eq i64 %18, 0
br i1 %19, label %20, label %64
20: ; preds = %15
%21 = getelementptr inbounds float, float* %1, i64 %17
%22 = load float, float* %21
%23 = fpext float %22 to double
%24 = fmul double %23, %5
%25 = fptrunc double %24 to float
%26 = getelementptr inbounds float, float* %2, i64 %16
store float %25, float* %26
br i1 %14, label %27, label %58
27: ; preds = %20, %27
%28 = phi float [ %57, %27 ], [ %25, %20 ]
%29 = phi i64 [ %55, %27 ], [ 0, %20 ]
%30 = load float, float* %21
%31 = getelementptr inbounds [2 x float], [2 x float]* %8, i64 %29, i64 0
%32 = load float, float* %31
%33 = getelementptr inbounds [2 x float], [2 x float]* %9, i64 %29, i64 0
%34 = load float, float* %33
%35 = fneg float %32
%36 = tail call float @llvm.fmuladd.f32(float %35, float %34, float %30)
%37 = getelementptr inbounds [2 x float], [2 x float]* %8, i64 %29, i64 1
%38 = load float, float* %37
%39 = getelementptr inbounds [2 x float], [2 x float]* %9, i64 %29, i64 1
%40 = load float, float* %39
%41 = fneg float %38
%42 = tail call float @llvm.fmuladd.f32(float %41, float %40, float %36)
%43 = fmul float %42, 2.000000e+00
%44 = getelementptr inbounds [3 x float], [3 x float]* %7, i64 %29, i64 2
%45 = load float, float* %44
%46 = getelementptr inbounds [3 x float], [3 x float]* %7, i64 %29, i64 1
%47 = load float, float* %46
%48 = fmul float %40, %47
%49 = tail call float @llvm.fmuladd.f32(float %45, float %43, float %48)
%50 = getelementptr inbounds [3 x float], [3 x float]* %7, i64 %29, i64 0
%51 = load float, float* %50
%52 = tail call float @llvm.fmuladd.f32(float %51, float %34, float %49)
%53 = fadd float %28, %52
store float %53, float* %26
%54 = load float, float* %39
store float %54, float* %33
store float %43, float* %39
%55 = add nuw nsw i64 %29, 1
%56 = icmp eq i64 %55, %4
%57 = load float, float* %26
br i1 %56, label %58, label %27
58: ; preds = %27, %20
%59 = phi float [ %25, %20 ], [ %57, %27 ]
%60 = fpext float %59 to double
%61 = fmul double %60, %6
%62 = fptrunc double %61 to float
store float %62, float* %26
%63 = add nsw i64 %17, 1
br label %104
64: ; preds = %15
%65 = getelementptr inbounds float, float* %2, i64 %16
store float 0.000000e+00, float* %65
br i1 %14, label %66, label %99
66: ; preds = %64, %66
%67 = phi float [ %98, %66 ], [ 0.000000e+00, %64 ]
%68 = phi i64 [ %96, %66 ], [ 0, %64 ]
%69 = getelementptr inbounds [2 x float], [2 x float]* %8, i64 %68, i64 0
%70 = load float, float* %69
%71 = getelementptr inbounds [2 x float], [2 x float]* %9, i64 %68, i64 0
%72 = load float, float* %71
%73 = fmul float %70, %72
%74 = fpext float %73 to double
%75 = fsub double 0.000000e+00, %74
%76 = getelementptr inbounds [2 x float], [2 x float]* %8, i64 %68, i64 1
%77 = load float, float* %76
%78 = getelementptr inbounds [2 x float], [2 x float]* %9, i64 %68, i64 1
%79 = load float, float* %78
%80 = fmul float %77, %79
%81 = fpext float %80 to double
%82 = fsub double %75, %81
%83 = fmul double %82, 2.000000e+00
%84 = fptrunc double %83 to float
%85 = getelementptr inbounds [3 x float], [3 x float]* %7, i64 %68, i64 2
%86 = load float, float* %85
%87 = getelementptr inbounds [3 x float], [3 x float]* %7, i64 %68, i64 1
%88 = load float, float* %87
%89 = fmul float %79, %88
%90 = tail call float @llvm.fmuladd.f32(float %86, float %84, float %89)
%91 = getelementptr inbounds [3 x float], [3 x float]* %7, i64 %68, i64 0
%92 = load float, float* %91
%93 = tail call float @llvm.fmuladd.f32(float %92, float %72, float %90)
%94 = fadd float %67, %93
store float %94, float* %65
%95 = load float, float* %78
store float %95, float* %71
store float %84, float* %78
%96 = add nuw nsw i64 %68, 1
%97 = icmp eq i64 %96, %4
%98 = load float, float* %65
br i1 %97, label %99, label %66
99: ; preds = %66, %64
%100 = phi float [ 0.000000e+00, %64 ], [ %98, %66 ]
%101 = fpext float %100 to double
%102 = fmul double %101, %6
%103 = fptrunc double %102 to float
store float %103, float* %65
br label %104
104: ; preds = %58, %99
%105 = phi i64 [ %63, %58 ], [ %17, %99 ]
%106 = add nuw nsw i64 %16, 1
%107 = icmp eq i64 %106, %11
br i1 %107, label %108, label %15
108: ; preds = %104, %10
%109 = phi i64 [ 0, %10 ], [ %11, %104 ]
ret i64 %109
}