int design_pfir(unsigned int n, unsigned int k, float *w, float **pw, float g,
                unsigned int flags) {
  int l = (int)n / k;
  int i;
  int j;
  float t;
  if (l < 1 || k < 1 || !w || !pw)
    return -1;
  if (flags & 0x00000002) {
    for (j = l - 1; j > -1; j--) {
      for (i = 0; i < (int)k; i++) {
        t = g * *w++;
        pw[i][j] = t * ((flags & 0x00000010) ? ((j & 1) ? -1 : 1) : 1);
      }
    }
  } else {
    for (j = 0; j < l; j++) {
      for (i = 0; i < (int)k; i++) {
        t = g * *w++;
        pw[i][j] = t * ((flags & 0x00000010) ? ((j & 1) ? 1 : -1) : 1);
      }
    }
  }
  return -1;
}


x86 assembly (clang, -O0):
<design_pfir>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movss  %xmm0,-0x24(%rbp)
mov    %r8d,-0x28(%rbp)
mov    -0x8(%rbp),%eax
xor    %edx,%edx
divl   -0xc(%rbp)
mov    %eax,-0x2c(%rbp)
cmpl   $0x1,-0x2c(%rbp)
jl     1150 <design_pfir+0x50>
cmpl   $0x1,-0xc(%rbp)
jb     1150 <design_pfir+0x50>
cmpq   $0x0,-0x18(%rbp)
je     1150 <design_pfir+0x50>
cmpq   $0x0,-0x20(%rbp)
jne    115c <design_pfir+0x5c>
movl   $0xffffffff,-0x4(%rbp)
jmp    1316 <design_pfir+0x216>
mov    -0x28(%rbp),%eax
and    $0x2,%eax
cmp    $0x0,%eax
je     123d <design_pfir+0x13d>
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x34(%rbp)
cmpl   $0xffffffff,-0x34(%rbp)
jle    1238 <design_pfir+0x138>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1225 <design_pfir+0x125>
movss  -0x24(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
mulss  (%rax),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
mov    -0x28(%rbp),%eax
and    $0x10,%eax
cmp    $0x0,%eax
je     11e5 <design_pfir+0xe5>
mov    -0x34(%rbp),%edx
and    $0x1,%edx
mov    $0x1,%eax
mov    $0xffffffff,%ecx
cmp    $0x0,%edx
cmovne %ecx,%eax
mov    %eax,-0x40(%rbp)
jmp    11f2 <design_pfir+0xf2>
mov    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11f2 <design_pfir+0xf2>
movss  -0x3c(%rbp),%xmm0
mov    -0x40(%rbp),%eax
cvtsi2ss %eax,%xmm1
mulss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1185 <design_pfir+0x85>
jmp    122a <design_pfir+0x12a>
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    1174 <design_pfir+0x74>
jmp    130f <design_pfir+0x20f>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    130a <design_pfir+0x20a>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12f7 <design_pfir+0x1f7>
movss  -0x24(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
mulss  (%rax),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x38(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x28(%rbp),%eax
and    $0x10,%eax
cmp    $0x0,%eax
je     12b7 <design_pfir+0x1b7>
mov    -0x34(%rbp),%edx
and    $0x1,%edx
mov    $0xffffffff,%eax
mov    $0x1,%ecx
cmp    $0x0,%edx
cmovne %ecx,%eax
mov    %eax,-0x48(%rbp)
jmp    12c4 <design_pfir+0x1c4>
mov    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    12c4 <design_pfir+0x1c4>
movss  -0x44(%rbp),%xmm0
mov    -0x48(%rbp),%eax
cvtsi2ss %eax,%xmm1
mulss  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1257 <design_pfir+0x157>
jmp    12fc <design_pfir+0x1fc>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1244 <design_pfir+0x144>
jmp    130f <design_pfir+0x20f>
movl   $0xffffffff,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<design_pfir>:
push   %r15
push   %r14
push   %rbx
mov    %rdx,%r9
mov    %edi,%eax
xor    %edx,%edx
div    %esi
test   %rcx,%rcx
je     12d6 <design_pfir+0x1d6>
test   %r9,%r9
je     12d6 <design_pfir+0x1d6>
test   %esi,%esi
je     12d6 <design_pfir+0x1d6>
test   %eax,%eax
jle    12d6 <design_pfir+0x1d6>
mov    %eax,%r15d
mov    %esi,%r10d
test   $0x2,%r8b
jne    11f4 <design_pfir+0xf4>
mov    %r10d,%edx
and    $0xfffffffe,%edx
xor    %edi,%edi
mov    $0x1,%r11d
jmp    115d <design_pfir+0x5d>
add    $0x1,%rdi
cmp    %r15,%rdi
je     12d6 <design_pfir+0x1d6>
test   %esi,%esi
jle    1150 <design_pfir+0x50>
mov    %edi,%eax
and    $0x1,%eax
cmp    $0x1,%rax
sbb    %eax,%eax
or     $0x1,%eax
test   $0x10,%r8b
cmove  %r11d,%eax
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
xor    %ebx,%ebx
cmp    $0x1,%esi
je     11cb <design_pfir+0xcb>
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%r9),%xmm2
mulss  %xmm0,%xmm2
mulss  %xmm1,%xmm2
mov    (%rcx,%rbx,8),%rax
movss  %xmm2,(%rax,%rdi,4)
movss  0x4(%r9),%xmm2
mulss  %xmm0,%xmm2
add    $0x8,%r9
mulss  %xmm1,%xmm2
mov    0x8(%rcx,%rbx,8),%rax
movss  %xmm2,(%rax,%rdi,4)
add    $0x2,%rbx
cmp    %rbx,%rdx
jne    1190 <design_pfir+0x90>
test   $0x1,%r10b
je     1150 <design_pfir+0x50>
movss  (%r9),%xmm2
mulss  %xmm0,%xmm2
add    $0x4,%r9
mulss  %xmm1,%xmm2
mov    (%rcx,%rbx,8),%rax
movss  %xmm2,(%rax,%rdi,4)
jmp    1150 <design_pfir+0x50>
mov    %r10d,%r14d
and    $0xfffffffe,%r14d
movaps 0xdfe(%rip),%xmm1
jmp    1228 <design_pfir+0x128>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x4,%r9
mov    (%rcx,%rdx,8),%rdx
movss  %xmm2,(%rdx,%r15,4)
cmp    $0x2,%r11
jl     12d6 <design_pfir+0x1d6>
mov    %r15,%r11
add    $0xffffffffffffffff,%r15
test   %esi,%esi
jle    121e <design_pfir+0x11e>
test   $0x10,%r8b
sete   %dl
test   $0x1,%r15b
sete   %bl
or     %dl,%bl
xor    %edx,%edx
cmp    $0x1,%esi
jne    1277 <design_pfir+0x177>
test   $0x1,%r10b
je     121e <design_pfir+0x11e>
jmp    12bd <design_pfir+0x1bd>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x8(%rcx,%rdx,8),%rax
movss  %xmm2,(%rax,%r15,4)
add    $0x2,%rdx
mov    %rdi,%r9
cmp    %rdx,%r14
je     12b0 <design_pfir+0x1b0>
movss  (%r9),%xmm2
mulss  %xmm0,%xmm2
lea    0x8(%r9),%rdi
test   %bl,%bl
jne    128b <design_pfir+0x18b>
xorps  %xmm1,%xmm2
mov    (%rcx,%rdx,8),%rax
movss  %xmm2,(%rax,%r15,4)
movss  0x4(%r9),%xmm2
mulss  %xmm0,%xmm2
jne    1260 <design_pfir+0x160>
xorps  %xmm1,%xmm2
jmp    1260 <design_pfir+0x160>
cs nopw 0x0(%rax,%rax,1)
mov    %rdi,%r9
test   $0x1,%r10b
je     121e <design_pfir+0x11e>
movss  (%r9),%xmm2
mulss  %xmm0,%xmm2
test   %bl,%bl
jne    1210 <design_pfir+0x110>
xorps  %xmm1,%xmm2
jmp    1210 <design_pfir+0x110>
mov    $0xffffffff,%eax
pop    %rbx
pop    %r14
pop    %r15
ret

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

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