int threshold_data(float *data, int *sizes, float threshold) {
  int i, j, k, index;
  for (i = 0; i < sizes[0]; i++) {
    for (j = 0; j < sizes[1]; j++) {
      for (k = 0; k < sizes[2]; k++) {
        index = i * sizes[2] * sizes[1] + j * sizes[2] + k;
        if (data[index] > threshold)
          data[index] = 1;
        else
          data[index] = 0;
      }
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<threshold_data>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x10(%rbp),%rcx
cmp    (%rcx),%eax
jge    11f5 <threshold_data+0xf5>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x10(%rbp),%rcx
cmp    0x4(%rcx),%eax
jge    11e2 <threshold_data+0xe2>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x10(%rbp),%rcx
cmp    0x8(%rcx),%eax
jge    11cf <threshold_data+0xcf>
mov    -0x18(%rbp),%eax
mov    -0x10(%rbp),%rcx
imul   0x8(%rcx),%eax
mov    -0x10(%rbp),%rcx
imul   0x4(%rcx),%eax
mov    -0x1c(%rbp),%ecx
mov    -0x10(%rbp),%rdx
imul   0x8(%rdx),%ecx
add    %ecx,%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x14(%rbp),%xmm0
jbe    11ac <threshold_data+0xac>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  0xe5e(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    11bc <threshold_data+0xbc>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    11c1 <threshold_data+0xc1>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1145 <threshold_data+0x45>
jmp    11d4 <threshold_data+0xd4>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    112e <threshold_data+0x2e>
jmp    11e7 <threshold_data+0xe7>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1118 <threshold_data+0x18>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<threshold_data>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rsi,-0x30(%rsp)
mov    (%rsi),%eax
mov    %eax,-0x34(%rsp)
test   %eax,%eax
jle    139a <threshold_data+0x29a>
mov    -0x30(%rsp),%rax
mov    0x4(%rax),%r11d
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
xor    %eax,%eax
mov    %rax,-0x40(%rsp)
movss  0xed4(%rip),%xmm2
movaps 0xebd(%rip),%xmm3
movl   $0x0,-0x44(%rsp)
jmp    1172 <threshold_data+0x72>
nopl   (%rax)
mov    -0x44(%rsp),%ecx
add    $0x1,%ecx
mov    -0x40(%rsp),%rax
add    %r11d,%eax
mov    %rax,-0x40(%rsp)
mov    %ecx,-0x44(%rsp)
cmp    -0x34(%rsp),%ecx
je     139a <threshold_data+0x29a>
test   %r11d,%r11d
jle    1150 <threshold_data+0x50>
mov    %r11d,%r10d
imul   -0x44(%rsp),%r10d
mov    -0x30(%rsp),%rax
mov    0x8(%rax),%r13d
lea    -0x1(%r13),%r9
mov    %r13d,%eax
and    $0xfffffff8,%eax
mov    %rax,-0x10(%rsp)
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x18(%rsp)
mov    %rax,%rcx
shr    $0x3,%rcx
add    $0x1,%rcx
mov    %rcx,%rax
mov    %rcx,-0x20(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,-0x28(%rsp)
mov    %r13,%r15
neg    %r15
mov    %r13d,%eax
imul   -0x40(%rsp),%eax
lea    0x1(%rax),%r8d
xor    %esi,%esi
mov    %r9,-0x8(%rsp)
jmp    11f2 <threshold_data+0xf2>
nopl   0x0(%rax,%rax,1)
add    $0x1,%esi
add    %r13d,%eax
add    %r13d,%r8d
cmp    %r11d,%esi
je     1150 <threshold_data+0x50>
test   %r13d,%r13d
jle    11e0 <threshold_data+0xe0>
lea    (%r10,%rsi,1),%ecx
mov    %ecx,%r12d
imul   %r13d,%r12d
cmp    $0x8,%r13d
jb     12c0 <threshold_data+0x1c0>
imul   %r13d,%ecx
lea    (%rcx,%r9,1),%edx
cmp    %ecx,%edx
jl     12c0 <threshold_data+0x1c0>
mov    %r9,%rcx
shr    $0x20,%rcx
mov    $0x0,%edx
jne    12c2 <threshold_data+0x1c2>
cmpq   $0x0,-0x18(%rsp)
je     1349 <threshold_data+0x249>
mov    -0x28(%rsp),%r9
mov    %eax,%ebp
xor    %r14d,%r14d
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movslq %ebp,%rbp
movups (%rdi,%rbp,4),%xmm4
movups 0x10(%rdi,%rbp,4),%xmm5
movaps %xmm1,%xmm6
cmpltps %xmm4,%xmm6
movaps %xmm1,%xmm4
cmpltps %xmm5,%xmm4
andps  %xmm3,%xmm6
andps  %xmm3,%xmm4
movups %xmm6,(%rdi,%rbp,4)
movups %xmm4,0x10(%rdi,%rbp,4)
lea    0x8(%rbp),%ecx
movslq %ecx,%rcx
movups (%rdi,%rcx,4),%xmm4
movups 0x10(%rdi,%rcx,4),%xmm5
movaps %xmm1,%xmm6
cmpltps %xmm4,%xmm6
movaps %xmm1,%xmm4
cmpltps %xmm5,%xmm4
andps  %xmm3,%xmm6
andps  %xmm3,%xmm4
movups %xmm6,(%rdi,%rcx,4)
movups %xmm4,0x10(%rdi,%rcx,4)
add    $0x10,%r14
add    $0x10,%ebp
add    $0xfffffffffffffffe,%r9
jne    1250 <threshold_data+0x150>
jmp    134c <threshold_data+0x24c>
nopw   0x0(%rax,%rax,1)
xor    %edx,%edx
mov    %rdx,%rcx
test   $0x1,%r13b
je     12e9 <threshold_data+0x1e9>
add    %edx,%r12d
movslq %r12d,%rcx
movaps %xmm0,%xmm4
cmpltss (%rdi,%rcx,4),%xmm4
andps  %xmm2,%xmm4
movss  %xmm4,(%rdi,%rcx,4)
mov    %rdx,%rcx
or     $0x1,%rcx
not    %rdx
cmp    %r15,%rdx
je     11e0 <threshold_data+0xe0>
lea    (%r8,%rcx,1),%ebp
mov    %r13,%rbx
sub    %rcx,%rbx
lea    (%rax,%rcx,1),%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movslq %edx,%rdx
movaps %xmm0,%xmm4
cmpltss (%rdi,%rdx,4),%xmm4
andps  %xmm2,%xmm4
movss  %xmm4,(%rdi,%rdx,4)
movslq %ebp,%rbp
movaps %xmm0,%xmm4
cmpltss (%rdi,%rbp,4),%xmm4
andps  %xmm2,%xmm4
movss  %xmm4,(%rdi,%rbp,4)
add    $0x2,%ebp
add    $0x2,%edx
add    $0xfffffffffffffffe,%rbx
jne    1310 <threshold_data+0x210>
jmp    11e0 <threshold_data+0xe0>
xor    %r14d,%r14d
testb  $0x1,-0x20(%rsp)
mov    -0x8(%rsp),%r9
je     1384 <threshold_data+0x284>
add    %r12d,%r14d
movslq %r14d,%rcx
movups (%rdi,%rcx,4),%xmm4
movups 0x10(%rdi,%rcx,4),%xmm5
movaps %xmm1,%xmm6
cmpltps %xmm4,%xmm6
movaps %xmm1,%xmm4
cmpltps %xmm5,%xmm4
andps  %xmm3,%xmm6
andps  %xmm3,%xmm4
movups %xmm6,(%rdi,%rcx,4)
movups %xmm4,0x10(%rdi,%rcx,4)
mov    -0x10(%rsp),%rcx
mov    %rcx,%rdx
cmp    %r13,%rcx
je     11e0 <threshold_data+0xe0>
jmp    12c2 <threshold_data+0x1c2>
xor    %eax,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @threshold_data(float* %0, i32* %1, float %2) {
%4 = alloca float*
%5 = alloca i32*
%6 = alloca float
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store float* %0, float** %4
store i32* %1, i32** %5
store float %2, float* %6
store i32 0, i32* %7
br label %11
11: ; preds = %75, %3
%12 = load i32, i32* %7
%13 = load i32*, i32** %5
%14 = getelementptr inbounds i32, i32* %13, i64 0
%15 = load i32, i32* %14
%16 = icmp slt i32 %12, %15
br i1 %16, label %17, label %78
17: ; preds = %11
store i32 0, i32* %8
br label %18
18: ; preds = %71, %17
%19 = load i32, i32* %8
%20 = load i32*, i32** %5
%21 = getelementptr inbounds i32, i32* %20, i64 1
%22 = load i32, i32* %21
%23 = icmp slt i32 %19, %22
br i1 %23, label %24, label %74
24: ; preds = %18
store i32 0, i32* %9
br label %25
25: ; preds = %67, %24
%26 = load i32, i32* %9
%27 = load i32*, i32** %5
%28 = getelementptr inbounds i32, i32* %27, i64 2
%29 = load i32, i32* %28
%30 = icmp slt i32 %26, %29
br i1 %30, label %31, label %70
31: ; preds = %25
%32 = load i32, i32* %7
%33 = load i32*, i32** %5
%34 = getelementptr inbounds i32, i32* %33, i64 2
%35 = load i32, i32* %34
%36 = mul nsw i32 %32, %35
%37 = load i32*, i32** %5
%38 = getelementptr inbounds i32, i32* %37, i64 1
%39 = load i32, i32* %38
%40 = mul nsw i32 %36, %39
%41 = load i32, i32* %8
%42 = load i32*, i32** %5
%43 = getelementptr inbounds i32, i32* %42, i64 2
%44 = load i32, i32* %43
%45 = mul nsw i32 %41, %44
%46 = add nsw i32 %40, %45
%47 = load i32, i32* %9
%48 = add nsw i32 %46, %47
store i32 %48, i32* %10
%49 = load float*, float** %4
%50 = load i32, i32* %10
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds float, float* %49, i64 %51
%53 = load float, float* %52
%54 = load float, float* %6
%55 = fcmp ogt float %53, %54
br i1 %55, label %56, label %61
56: ; preds = %31
%57 = load float*, float** %4
%58 = load i32, i32* %10
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds float, float* %57, i64 %59
store float 1.000000e+00, float* %60
br label %66
61: ; preds = %31
%62 = load float*, float** %4
%63 = load i32, i32* %10
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds float, float* %62, i64 %64
store float 0.000000e+00, float* %65
br label %66
66: ; preds = %61, %56
br label %67
67: ; preds = %66
%68 = load i32, i32* %9
%69 = add nsw i32 %68, 1
store i32 %69, i32* %9
br label %25
70: ; preds = %25
br label %71
71: ; preds = %70
%72 = load i32, i32* %8
%73 = add nsw i32 %72, 1
store i32 %73, i32* %8
br label %18
74: ; preds = %18
br label %75
75: ; preds = %74
%76 = load i32, i32* %7
%77 = add nsw i32 %76, 1
store i32 %77, i32* %7
br label %11
78: ; preds = %11
ret i32 0
}

LLVM IR (-O2):
define i32 @threshold_data(float* nocapture %0, i32* nocapture %1, float %2) {
%4 = load i32, i32* %1
%5 = icmp sgt i32 %4, 0
br i1 %5, label %6, label %148
6: ; preds = %3
%7 = getelementptr inbounds i32, i32* %1, i64 1
%8 = load i32, i32* %7
%9 = icmp sgt i32 %8, 0
%10 = getelementptr inbounds i32, i32* %1, i64 2
%11 = insertelement <4 x float> poison, float %2, i64 0
%12 = shufflevector <4 x float> %11, <4 x float> poison, <4 x i32> zeroinitializer
%13 = insertelement <4 x float> poison, float %2, i64 0
%14 = shufflevector <4 x float> %13, <4 x float> poison, <4 x i32> zeroinitializer
br label %15
15: ; preds = %6, %145
%16 = phi i32 [ 0, %6 ], [ %146, %145 ]
br i1 %9, label %17, label %145
17: ; preds = %15
%18 = mul i32 %8, %16
%19 = load i32, i32* %10
%20 = icmp sgt i32 %19, 0
%21 = mul i32 %8, %16
%22 = zext i32 %19 to i64
%23 = add nsw i64 %22, -1
%24 = and i64 %22, 4294967288
%25 = add nsw i64 %24, -8
%26 = lshr exact i64 %25, 3
%27 = add nuw nsw i64 %26, 1
%28 = icmp ult i32 %19, 8
%29 = trunc i64 %23 to i32
%30 = icmp ugt i64 %23, 4294967295
%31 = and i64 %22, 4294967288
%32 = and i64 %27, 1
%33 = icmp eq i64 %25, 0
%34 = and i64 %27, 4611686018427387902
%35 = icmp eq i64 %32, 0
%36 = icmp eq i64 %31, %22
%37 = and i64 %22, 1
%38 = icmp eq i64 %37, 0
%39 = sub nsw i64 0, %22
br label %40
40: ; preds = %17, %142
%41 = phi i32 [ 0, %17 ], [ %143, %142 ]
%42 = add i32 %18, %41
%43 = mul i32 %19, %42
br i1 %20, label %44, label %142
44: ; preds = %40
%45 = add i32 %21, %41
%46 = mul i32 %45, %19
br i1 %28, label %108, label %47
47: ; preds = %44
%48 = add i32 %43, %29
%49 = icmp slt i32 %48, %43
%50 = or i1 %49, %30
br i1 %50, label %108, label %51
51: ; preds = %47
br i1 %33, label %89, label %52
52: ; preds = %51, %52
%53 = phi i64 [ %86, %52 ], [ 0, %51 ]
%54 = phi i64 [ %87, %52 ], [ 0, %51 ]
%55 = trunc i64 %53 to i32
%56 = add i32 %46, %55
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds float, float* %0, i64 %57
%59 = bitcast float* %58 to <4 x float>*
%60 = load <4 x float>, <4 x float>* %59
%61 = getelementptr inbounds float, float* %58, i64 4
%62 = bitcast float* %61 to <4 x float>*
%63 = load <4 x float>, <4 x float>* %62
%64 = fcmp ogt <4 x float> %60, %12
%65 = fcmp ogt <4 x float> %63, %14
%66 = select <4 x i1> %64, <4 x float> <float 1.000000e+00, float 1.000000e+00, float 1.000000e+00, float 1.000000e+00>, <4 x float> zeroinitializer
%67 = select <4 x i1> %65, <4 x float> <float 1.000000e+00, float 1.000000e+00, float 1.000000e+00, float 1.000000e+00>, <4 x float> zeroinitializer
%68 = bitcast float* %58 to <4 x float>*
store <4 x float> %66, <4 x float>* %68
%69 = bitcast float* %61 to <4 x float>*
store <4 x float> %67, <4 x float>* %69
%70 = trunc i64 %53 to i32
%71 = or i32 %70, 8
%72 = add i32 %46, %71
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds float, float* %0, i64 %73
%75 = bitcast float* %74 to <4 x float>*
%76 = load <4 x float>, <4 x float>* %75
%77 = getelementptr inbounds float, float* %74, i64 4
%78 = bitcast float* %77 to <4 x float>*
%79 = load <4 x float>, <4 x float>* %78
%80 = fcmp ogt <4 x float> %76, %12
%81 = fcmp ogt <4 x float> %79, %14
%82 = select <4 x i1> %80, <4 x float> <float 1.000000e+00, float 1.000000e+00, float 1.000000e+00, float 1.000000e+00>, <4 x float> zeroinitializer
%83 = select <4 x i1> %81, <4 x float> <float 1.000000e+00, float 1.000000e+00, float 1.000000e+00, float 1.000000e+00>, <4 x float> zeroinitializer
%84 = bitcast float* %74 to <4 x float>*
store <4 x float> %82, <4 x float>* %84
%85 = bitcast float* %77 to <4 x float>*
store <4 x float> %83, <4 x float>* %85
%86 = add nuw i64 %53, 16
%87 = add i64 %54, 2
%88 = icmp eq i64 %87, %34
br i1 %88, label %89, label %52
89: ; preds = %52, %51
%90 = phi i64 [ 0, %51 ], [ %86, %52 ]
br i1 %35, label %107, label %91
91: ; preds = %89
%92 = trunc i64 %90 to i32
%93 = add i32 %46, %92
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds float, float* %0, i64 %94
%96 = bitcast float* %95 to <4 x float>*
%97 = load <4 x float>, <4 x float>* %96
%98 = getelementptr inbounds float, float* %95, i64 4
%99 = bitcast float* %98 to <4 x float>*
%100 = load <4 x float>, <4 x float>* %99
%101 = fcmp ogt <4 x float> %97, %12
%102 = fcmp ogt <4 x float> %100, %14
%103 = select <4 x i1> %101, <4 x float> <float 1.000000e+00, float 1.000000e+00, float 1.000000e+00, float 1.000000e+00>, <4 x float> zeroinitializer
%104 = select <4 x i1> %102, <4 x float> <float 1.000000e+00, float 1.000000e+00, float 1.000000e+00, float 1.000000e+00>, <4 x float> zeroinitializer
%105 = bitcast float* %95 to <4 x float>*
store <4 x float> %103, <4 x float>* %105
%106 = bitcast float* %98 to <4 x float>*
store <4 x float> %104, <4 x float>* %106
br label %107
107: ; preds = %89, %91
br i1 %36, label %142, label %108
108: ; preds = %47, %44, %107
%109 = phi i64 [ 0, %47 ], [ 0, %44 ], [ %31, %107 ]
%110 = xor i64 %109, -1
br i1 %38, label %120, label %111
111: ; preds = %108
%112 = trunc i64 %109 to i32
%113 = add i32 %46, %112
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds float, float* %0, i64 %114
%116 = load float, float* %115
%117 = fcmp ogt float %116, %2
%118 = select i1 %117, float 1.000000e+00, float 0.000000e+00
store float %118, float* %115
%119 = or i64 %109, 1
br label %120
120: ; preds = %111, %108
%121 = phi i64 [ %109, %108 ], [ %119, %111 ]
%122 = icmp eq i64 %110, %39
br i1 %122, label %142, label %123
123: ; preds = %120, %123
%124 = phi i64 [ %140, %123 ], [ %121, %120 ]
%125 = trunc i64 %124 to i32
%126 = add i32 %46, %125
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds float, float* %0, i64 %127
%129 = load float, float* %128
%130 = fcmp ogt float %129, %2
%131 = select i1 %130, float 1.000000e+00, float 0.000000e+00
store float %131, float* %128
%132 = trunc i64 %124 to i32
%133 = add i32 %132, 1
%134 = add i32 %46, %133
%135 = sext i32 %134 to i64
%136 = getelementptr inbounds float, float* %0, i64 %135
%137 = load float, float* %136
%138 = fcmp ogt float %137, %2
%139 = select i1 %138, float 1.000000e+00, float 0.000000e+00
store float %139, float* %136
%140 = add nuw nsw i64 %124, 2
%141 = icmp eq i64 %140, %22
br i1 %141, label %142, label %123
142: ; preds = %120, %123, %107, %40
%143 = add nuw nsw i32 %41, 1
%144 = icmp eq i32 %143, %8
br i1 %144, label %145, label %40
145: ; preds = %142, %15
%146 = add nuw nsw i32 %16, 1
%147 = icmp eq i32 %146, %4
br i1 %147, label %148, label %15
148: ; preds = %145, %3
ret i32 0
}