#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ printf(char *);
int pack2bitNchan(float **in, int nchan, int off, char *out, float mean,
                  float stddev, int len) {
  int i, j, n, k, ch[4];
  float maxposThresh, maxnegThresh;
  if (len * nchan % 4 != 0) {
    printf("Can only pack multiple of 4 samples*channels!\n");
    return (1);
  }
  maxposThresh = mean + stddev * 0.95;
  maxnegThresh = mean - stddev * 0.95;
  j = 0;
  k = 0;
  for (i = off; i < len + off; i++) {
    for (n = 0; n < nchan; n++) {
      {
        if (in[n][i] >= maxposThresh)
          ch[k] = 3;
        else if (in[n][i] <= maxnegThresh)
          ch[k] = 0;
        else if (in[n][i] > mean)
          ch[k] = 2;
        else
          ch[k] = 1;
      };
      k++;
      if (k == 4) {
        out[j] = (ch[0]) | ((ch[1] << 2)) | ((ch[2] << 4)) | ((ch[3] << 6));
        j++;
        k = 0;
      }
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<pack2bitNchan>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movss  %xmm0,-0x24(%rbp)
movss  %xmm1,-0x28(%rbp)
mov    %r8d,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
je     1164 <pack2bitNchan+0x54>
lea    0xeb5(%rip),%rdi
call   1030 <printf@plt>
movl   $0x1,-0x4(%rbp)
jmp    132b <pack2bitNchan+0x21b>
movss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movss  -0x28(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe82(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x24(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movss  -0x28(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0xe40(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x58(%rbp)
movl   $0x0,-0x34(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x2c(%rbp),%ecx
add    -0x18(%rbp),%ecx
cmp    %ecx,%eax
jge    1324 <pack2bitNchan+0x214>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1311 <pack2bitNchan+0x201>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x54(%rbp),%xmm0
jb     1239 <pack2bitNchan+0x129>
movslq -0x3c(%rbp),%rax
movl   $0x3,-0x50(%rbp,%rax,4)
jmp    12b3 <pack2bitNchan+0x1a3>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movss  -0x58(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jb     126d <pack2bitNchan+0x15d>
movslq -0x3c(%rbp),%rax
movl   $0x0,-0x50(%rbp,%rax,4)
jmp    12ae <pack2bitNchan+0x19e>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x24(%rbp),%xmm0
jbe    129d <pack2bitNchan+0x18d>
movslq -0x3c(%rbp),%rax
movl   $0x2,-0x50(%rbp,%rax,4)
jmp    12a9 <pack2bitNchan+0x199>
movslq -0x3c(%rbp),%rax
movl   $0x1,-0x50(%rbp,%rax,4)
jmp    12ae <pack2bitNchan+0x19e>
jmp    12b3 <pack2bitNchan+0x1a3>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
cmpl   $0x4,-0x3c(%rbp)
jne    12fe <pack2bitNchan+0x1ee>
mov    -0x50(%rbp),%eax
mov    -0x4c(%rbp),%ecx
shl    $0x2,%ecx
or     %ecx,%eax
mov    -0x48(%rbp),%ecx
shl    $0x4,%ecx
or     %ecx,%eax
mov    -0x44(%rbp),%ecx
shl    $0x6,%ecx
or     %ecx,%eax
mov    %al,%dl
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x3c(%rbp)
jmp    1303 <pack2bitNchan+0x1f3>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11fd <pack2bitNchan+0xed>
jmp    1316 <pack2bitNchan+0x206>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11e5 <pack2bitNchan+0xd5>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<pack2bitNchan>:
push   %rbp
push   %rbx
sub    $0x18,%rsp
mov    %r8d,%eax
imul   %esi,%eax
test   $0x3,%al
je     1136 <pack2bitNchan+0x26>
lea    0xee1(%rip),%rdi
call   1030 <puts@plt>
mov    $0x1,%eax
jmp    1238 <pack2bitNchan+0x128>
test   %r8d,%r8d
jle    1236 <pack2bitNchan+0x126>
cvtss2sd %xmm0,%xmm2
cvtss2sd %xmm1,%xmm3
mulsd  0xeb1(%rip),%xmm3
movapd %xmm3,%xmm1
addsd  %xmm2,%xmm1
cvtsd2ss %xmm1,%xmm1
subsd  %xmm3,%xmm2
cvtsd2ss %xmm2,%xmm2
add    %edx,%r8d
movslq %edx,%r10
movslq %r8d,%r8
mov    %esi,%r9d
xor    %eax,%eax
xor    %r11d,%r11d
jmp    118d <pack2bitNchan+0x7d>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r10
cmp    %r8,%r10
jge    1236 <pack2bitNchan+0x126>
test   %esi,%esi
jle    1180 <pack2bitNchan+0x70>
xor    %edx,%edx
jmp    11a9 <pack2bitNchan+0x99>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rdx
cmp    %rdx,%r9
je     1180 <pack2bitNchan+0x70>
mov    (%rdi,%rdx,8),%rbx
movss  (%rbx,%r10,4),%xmm3
ucomiss %xmm1,%xmm3
jae    11d0 <pack2bitNchan+0xc0>
ucomiss %xmm3,%xmm2
jae    11e0 <pack2bitNchan+0xd0>
ucomiss %xmm0,%xmm3
movslq %eax,%rbx
jbe    11ec <pack2bitNchan+0xdc>
movl   $0x2,(%rsp,%rbx,4)
jmp    1200 <pack2bitNchan+0xf0>
xchg   %ax,%ax
movslq %eax,%rbx
movl   $0x3,(%rsp,%rbx,4)
jmp    1200 <pack2bitNchan+0xf0>
nopl   0x0(%rax)
movslq %eax,%rbx
movl   $0x0,(%rsp,%rbx,4)
jmp    1200 <pack2bitNchan+0xf0>
movl   $0x1,(%rsp,%rbx,4)
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%eax
cmp    $0x4,%eax
jne    11a0 <pack2bitNchan+0x90>
mov    0x4(%rsp),%eax
mov    0x8(%rsp),%ebx
shl    $0x2,%eax
or     (%rsp),%eax
shl    $0x4,%ebx
mov    0xc(%rsp),%ebp
shl    $0x6,%ebp
or     %ebx,%ebp
or     %eax,%ebp
movslq %r11d,%r11
mov    %bpl,(%rcx,%r11,1)
add    $0x1,%r11d
xor    %eax,%eax
jmp    11a0 <pack2bitNchan+0x90>
xor    %eax,%eax
add    $0x18,%rsp
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
@.str = [47 x i8] c"Can only pack multiple of 4 samples*channels!\0A\00"
define i32 @pack2bitNchan(float** %0, i32 %1, i32 %2, i8* %3, float %4, float %5, i32 %6) {
%8 = alloca i32
%9 = alloca float**
%10 = alloca i32
%11 = alloca i32
%12 = alloca i8*
%13 = alloca float
%14 = alloca float
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca [4 x i32]
%21 = alloca float
%22 = alloca float
store float** %0, float*** %9
store i32 %1, i32* %10
store i32 %2, i32* %11
store i8* %3, i8** %12
store float %4, float* %13
store float %5, float* %14
store i32 %6, i32* %15
%23 = load i32, i32* %15
%24 = load i32, i32* %10
%25 = mul nsw i32 %23, %24
%26 = srem i32 %25, 4
%27 = icmp ne i32 %26, 0
br i1 %27, label %28, label %30
28: ; preds = %7
%29 = call i32 @printf(i8* getelementptr inbounds ([47 x i8], [47 x i8]* @.str, i64 0, i64 0))
store i32 1, i32* %8
br label %146
30: ; preds = %7
%31 = load float, float* %13
%32 = fpext float %31 to double
%33 = load float, float* %14
%34 = fpext float %33 to double
%35 = call double @llvm.fmuladd.f64(double %34, double 0x3FEE666666666666, double %32)
%36 = fptrunc double %35 to float
store float %36, float* %21
%37 = load float, float* %13
%38 = fpext float %37 to double
%39 = load float, float* %14
%40 = fpext float %39 to double
%41 = fneg double %40
%42 = call double @llvm.fmuladd.f64(double %41, double 0x3FEE666666666666, double %38)
%43 = fptrunc double %42 to float
store float %43, float* %22
store i32 0, i32* %17
store i32 0, i32* %19
%44 = load i32, i32* %11
store i32 %44, i32* %16
br label %45
45: ; preds = %142, %30
%46 = load i32, i32* %16
%47 = load i32, i32* %15
%48 = load i32, i32* %11
%49 = add nsw i32 %47, %48
%50 = icmp slt i32 %46, %49
br i1 %50, label %51, label %145
51: ; preds = %45
store i32 0, i32* %18
br label %52
52: ; preds = %138, %51
%53 = load i32, i32* %18
%54 = load i32, i32* %10
%55 = icmp slt i32 %53, %54
br i1 %55, label %56, label %141
56: ; preds = %52
%57 = load float**, float*** %9
%58 = load i32, i32* %18
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds float*, float** %57, i64 %59
%61 = load float*, float** %60
%62 = load i32, i32* %16
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds float, float* %61, i64 %63
%65 = load float, float* %64
%66 = load float, float* %21
%67 = fcmp oge float %65, %66
br i1 %67, label %68, label %72
68: ; preds = %56
%69 = load i32, i32* %19
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds [4 x i32], [4 x i32]* %20, i64 0, i64 %70
store i32 3, i32* %71
br label %110
72: ; preds = %56
%73 = load float**, float*** %9
%74 = load i32, i32* %18
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds float*, float** %73, i64 %75
%77 = load float*, float** %76
%78 = load i32, i32* %16
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds float, float* %77, i64 %79
%81 = load float, float* %80
%82 = load float, float* %22
%83 = fcmp ole float %81, %82
br i1 %83, label %84, label %88
84: ; preds = %72
%85 = load i32, i32* %19
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds [4 x i32], [4 x i32]* %20, i64 0, i64 %86
store i32 0, i32* %87
br label %109
88: ; preds = %72
%89 = load float**, float*** %9
%90 = load i32, i32* %18
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds float*, float** %89, i64 %91
%93 = load float*, float** %92
%94 = load i32, i32* %16
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds float, float* %93, i64 %95
%97 = load float, float* %96
%98 = load float, float* %13
%99 = fcmp ogt float %97, %98
br i1 %99, label %100, label %104
100: ; preds = %88
%101 = load i32, i32* %19
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds [4 x i32], [4 x i32]* %20, i64 0, i64 %102
store i32 2, i32* %103
br label %108
104: ; preds = %88
%105 = load i32, i32* %19
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds [4 x i32], [4 x i32]* %20, i64 0, i64 %106
store i32 1, i32* %107
br label %108
108: ; preds = %104, %100
br label %109
109: ; preds = %108, %84
br label %110
110: ; preds = %109, %68
%111 = load i32, i32* %19
%112 = add nsw i32 %111, 1
store i32 %112, i32* %19
%113 = load i32, i32* %19
%114 = icmp eq i32 %113, 4
br i1 %114, label %115, label %137
115: ; preds = %110
%116 = getelementptr inbounds [4 x i32], [4 x i32]* %20, i64 0, i64 0
%117 = load i32, i32* %116
%118 = getelementptr inbounds [4 x i32], [4 x i32]* %20, i64 0, i64 1
%119 = load i32, i32* %118
%120 = shl i32 %119, 2
%121 = or i32 %117, %120
%122 = getelementptr inbounds [4 x i32], [4 x i32]* %20, i64 0, i64 2
%123 = load i32, i32* %122
%124 = shl i32 %123, 4
%125 = or i32 %121, %124
%126 = getelementptr inbounds [4 x i32], [4 x i32]* %20, i64 0, i64 3
%127 = load i32, i32* %126
%128 = shl i32 %127, 6
%129 = or i32 %125, %128
%130 = trunc i32 %129 to i8
%131 = load i8*, i8** %12
%132 = load i32, i32* %17
%133 = sext i32 %132 to i64
%134 = getelementptr inbounds i8, i8* %131, i64 %133
store i8 %130, i8* %134
%135 = load i32, i32* %17
%136 = add nsw i32 %135, 1
store i32 %136, i32* %17
store i32 0, i32* %19
br label %137
137: ; preds = %115, %110
br label %138
138: ; preds = %137
%139 = load i32, i32* %18
%140 = add nsw i32 %139, 1
store i32 %140, i32* %18
br label %52
141: ; preds = %52
br label %142
142: ; preds = %141
%143 = load i32, i32* %16
%144 = add nsw i32 %143, 1
store i32 %144, i32* %16
br label %45
145: ; preds = %45
store i32 0, i32* %8
br label %146
146: ; preds = %145, %28
%147 = load i32, i32* %8
ret i32 %147
}

LLVM IR (-O2):
@str = [46 x i8] c"Can only pack multiple of 4 samples*channels!\00"
define i32 @pack2bitNchan(float** nocapture %0, i32 %1, i32 %2, i8* nocapture writeonly %3, float %4, float %5, i32 %6) {
%8 = alloca [4 x i32]
%9 = bitcast [4 x i32]* %8 to i8*
%10 = mul nsw i32 %6, %1
%11 = and i32 %10, 3
%12 = icmp eq i32 %11, 0
br i1 %12, label %15, label %13
13: ; preds = %7
%14 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([46 x i8], [46 x i8]* @str, i64 0, i64 0))
br label %89
15: ; preds = %7
%16 = fpext float %4 to double
%17 = fpext float %5 to double
%18 = tail call double @llvm.fmuladd.f64(double %17, double 0x3FEE666666666666, double %16)
%19 = fptrunc double %18 to float
%20 = fneg double %17
%21 = tail call double @llvm.fmuladd.f64(double %20, double 0x3FEE666666666666, double %16)
%22 = fptrunc double %21 to float
%23 = icmp sgt i32 %6, 0
br i1 %23, label %24, label %89
24: ; preds = %15
%25 = add nsw i32 %6, %2
%26 = icmp sgt i32 %1, 0
%27 = getelementptr inbounds [4 x i32], [4 x i32]* %8, i64 0, i64 0
%28 = getelementptr inbounds [4 x i32], [4 x i32]* %8, i64 0, i64 1
%29 = getelementptr inbounds [4 x i32], [4 x i32]* %8, i64 0, i64 2
%30 = getelementptr inbounds [4 x i32], [4 x i32]* %8, i64 0, i64 3
%31 = sext i32 %2 to i64
%32 = sext i32 %25 to i64
%33 = zext i32 %1 to i64
br label %34
34: ; preds = %24, %84
%35 = phi i64 [ %31, %24 ], [ %87, %84 ]
%36 = phi i32 [ 0, %24 ], [ %86, %84 ]
%37 = phi i32 [ 0, %24 ], [ %85, %84 ]
br i1 %26, label %38, label %84
38: ; preds = %34, %79
%39 = phi i64 [ %82, %79 ], [ 0, %34 ]
%40 = phi i32 [ %81, %79 ], [ %36, %34 ]
%41 = phi i32 [ %80, %79 ], [ %37, %34 ]
%42 = getelementptr inbounds float*, float** %0, i64 %39
%43 = load float*, float** %42
%44 = getelementptr inbounds float, float* %43, i64 %35
%45 = load float, float* %44
%46 = fcmp ult float %45, %19
br i1 %46, label %50, label %47
47: ; preds = %38
%48 = sext i32 %40 to i64
%49 = getelementptr inbounds [4 x i32], [4 x i32]* %8, i64 0, i64 %48
store i32 3, i32* %49
br label %61
50: ; preds = %38
%51 = fcmp ugt float %45, %22
br i1 %51, label %55, label %52
52: ; preds = %50
%53 = sext i32 %40 to i64
%54 = getelementptr inbounds [4 x i32], [4 x i32]* %8, i64 0, i64 %53
store i32 0, i32* %54
br label %61
55: ; preds = %50
%56 = fcmp ogt float %45, %4
%57 = sext i32 %40 to i64
%58 = getelementptr inbounds [4 x i32], [4 x i32]* %8, i64 0, i64 %57
br i1 %56, label %59, label %60
59: ; preds = %55
store i32 2, i32* %58
br label %61
60: ; preds = %55
store i32 1, i32* %58
br label %61
61: ; preds = %52, %60, %59, %47
%62 = add nsw i32 %40, 1
%63 = icmp eq i32 %62, 4
br i1 %63, label %64, label %79
64: ; preds = %61
%65 = load i32, i32* %27
%66 = load i32, i32* %28
%67 = shl i32 %66, 2
%68 = or i32 %67, %65
%69 = load i32, i32* %29
%70 = shl i32 %69, 4
%71 = or i32 %68, %70
%72 = load i32, i32* %30
%73 = shl i32 %72, 6
%74 = or i32 %71, %73
%75 = trunc i32 %74 to i8
%76 = sext i32 %41 to i64
%77 = getelementptr inbounds i8, i8* %3, i64 %76
store i8 %75, i8* %77
%78 = add nsw i32 %41, 1
br label %79
79: ; preds = %61, %64
%80 = phi i32 [ %78, %64 ], [ %41, %61 ]
%81 = phi i32 [ 0, %64 ], [ %62, %61 ]
%82 = add nuw nsw i64 %39, 1
%83 = icmp eq i64 %82, %33
br i1 %83, label %84, label %38
84: ; preds = %79, %34
%85 = phi i32 [ %37, %34 ], [ %80, %79 ]
%86 = phi i32 [ %36, %34 ], [ %81, %79 ]
%87 = add nsw i64 %35, 1
%88 = icmp slt i64 %87, %32
br i1 %88, label %34, label %89
89: ; preds = %84, %15, %13
%90 = phi i32 [ 1, %13 ], [ 0, %15 ], [ 0, %84 ]
ret i32 %90
}