#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 map_index(int, int, int, int, int, int);
bool check_inclusion(int i, int j, int k, int m, int n, int p, double *img,
                     int t) {
  double sum;
  int imin, imax, jmin, jmax, kmin, kmax;
  int count, l;
  imin = i - t;
  imin = imin < 0 ? 0 : imin;
  jmin = j - t;
  jmin = jmin < 0 ? 0 : jmin;
  kmin = k - t;
  kmin = kmin < 0 ? 0 : kmin;
  imax = i + t;
  imax = imax > m - 1 ? m - 1 : imax;
  jmax = j + t;
  jmax = jmax > n - 1 ? n - 1 : jmax;
  kmax = k + t;
  kmax = kmax > p - 1 ? p - 1 : kmax;
  sum = 0.0;
  count = 0;
  for (int ii = imin; ii <= imax; ii++) {
    for (int jj = jmin; jj <= jmax; jj++) {
      for (int kk = kmin; kk <= kmax; kk++) {
        l = map_index(ii, jj, kk, m, n, p);
        sum += img[l];
        count += 1;
      }
    }
  }
  l = map_index(i, j, k, m, n, p);
  return (img[l] >= (sum / count));
}


x86 assembly (clang, -O0):
<check_inclusion>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8d,-0x14(%rbp)
mov    %r9d,-0x18(%rbp)
mov    -0x4(%rbp),%eax
sub    0x18(%rbp),%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jge    1150 <check_inclusion+0x40>
xor    %eax,%eax
mov    %eax,-0x50(%rbp)
jmp    1156 <check_inclusion+0x46>
mov    -0x24(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x50(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x8(%rbp),%eax
sub    0x18(%rbp),%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jge    1179 <check_inclusion+0x69>
xor    %eax,%eax
mov    %eax,-0x54(%rbp)
jmp    117f <check_inclusion+0x6f>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0xc(%rbp),%eax
sub    0x18(%rbp),%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x34(%rbp)
jge    11a2 <check_inclusion+0x92>
xor    %eax,%eax
mov    %eax,-0x58(%rbp)
jmp    11a8 <check_inclusion+0x98>
mov    -0x34(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x58(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x4(%rbp),%eax
add    0x18(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jle    11d6 <check_inclusion+0xc6>
mov    -0x10(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    11dc <check_inclusion+0xcc>
mov    -0x28(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x8(%rbp),%eax
add    0x18(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jle    120a <check_inclusion+0xfa>
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    1210 <check_inclusion+0x100>
mov    -0x30(%rbp),%eax
mov    %eax,-0x60(%rbp)
mov    -0x60(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0xc(%rbp),%eax
add    0x18(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jle    123e <check_inclusion+0x12e>
mov    -0x18(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    1244 <check_inclusion+0x134>
mov    -0x38(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
mov    %eax,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x28(%rbp),%eax
jg     12ff <check_inclusion+0x1ef>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x30(%rbp),%eax
jg     12ec <check_inclusion+0x1dc>
mov    -0x34(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jg     12d9 <check_inclusion+0x1c9>
mov    -0x44(%rbp),%edi
mov    -0x48(%rbp),%esi
mov    -0x4c(%rbp),%edx
mov    -0x10(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    -0x18(%rbp),%r9d
call   1030 <map_index@plt>
mov    %eax,-0x40(%rbp)
mov    0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1283 <check_inclusion+0x173>
jmp    12de <check_inclusion+0x1ce>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1271 <check_inclusion+0x161>
jmp    12f1 <check_inclusion+0x1e1>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    125f <check_inclusion+0x14f>
mov    -0x4(%rbp),%edi
mov    -0x8(%rbp),%esi
mov    -0xc(%rbp),%edx
mov    -0x10(%rbp),%ecx
mov    -0x14(%rbp),%r8d
mov    -0x18(%rbp),%r9d
call   1030 <map_index@plt>
mov    %eax,-0x40(%rbp)
mov    0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x20(%rbp),%xmm1
cvtsi2sdl -0x3c(%rbp),%xmm2
divsd  %xmm2,%xmm1
ucomisd %xmm1,%xmm0
setae  %al
and    $0x1,%al
movzbl %al,%eax
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<check_inclusion>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r9d,%r15d
mov    0x88(%rsp),%ebp
xor    %ebx,%ebx
mov    %edi,%r10d
sub    %ebp,%r10d
cmovs  %ebx,%r10d
mov    %esi,%eax
sub    %ebp,%eax
cmovs  %ebx,%eax
mov    %eax,0x10(%rsp)
mov    %edx,%r12d
sub    %ebp,%r12d
cmovs  %ebx,%r12d
mov    %rdi,0x28(%rsp)
lea    (%rdi,%rbp,1),%eax
lea    -0x1(%rcx),%edi
cmp    %edi,%eax
cmovg  %edi,%eax
mov    %rsi,0x30(%rsp)
lea    (%rsi,%rbp,1),%edi
lea    -0x1(%r8),%esi
cmp    %esi,%edi
cmovg  %esi,%edi
mov    0x80(%rsp),%r13
mov    %edx,0x8(%rsp)
add    %edx,%ebp
lea    -0x1(%r15),%edx
cmp    %edx,%ebp
cmovg  %edx,%ebp
xorpd  %xmm0,%xmm0
mov    %eax,0xc(%rsp)
cmp    %eax,%r10d
jle    11d7 <check_inclusion+0xc7>
movsd  %xmm0,(%rsp)
mov    0x28(%rsp),%rdi
mov    0x30(%rsp),%rsi
mov    0x8(%rsp),%edx
mov    %r15d,%r9d
call   1030 <map_index@plt>
cltq
movsd  0x0(%r13,%rax,8),%xmm0
cvtsi2sd %ebx,%xmm1
movsd  (%rsp),%xmm2
divsd  %xmm1,%xmm2
xor    %eax,%eax
ucomisd %xmm2,%xmm0
setae  %al
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    0x1(%rbp),%edx
mov    %edx,%eax
sub    %r12d,%eax
mov    %eax,0x14(%rsp)
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
mov    %rbp,0x38(%rsp)
mov    %r12d,0x1c(%rsp)
mov    %edi,0x18(%rsp)
mov    %edx,0x24(%rsp)
jmp    120e <check_inclusion+0xfe>
nopl   (%rax)
lea    0x1(%r10),%eax
cmp    0xc(%rsp),%r10d
mov    %eax,%r10d
je     118e <check_inclusion+0x7e>
mov    0x10(%rsp),%eax
mov    %eax,%esi
cmp    %edi,%eax
jle    1229 <check_inclusion+0x119>
jmp    1200 <check_inclusion+0xf0>
nopw   0x0(%rax,%rax,1)
lea    0x1(%rsi),%eax
cmp    %edi,%esi
mov    %eax,%esi
je     1200 <check_inclusion+0xf0>
cmp    %ebp,%r12d
jg     1220 <check_inclusion+0x110>
mov    %ebx,0x20(%rsp)
mov    %rsi,0x40(%rsp)
nopw   0x0(%rax,%rax,1)
movsd  %xmm0,(%rsp)
mov    %r10d,%edi
mov    %r12d,%edx
mov    %rcx,%rbp
mov    %r8,%rbx
mov    %r15d,%r9d
mov    %r10,%r14
call   1030 <map_index@plt>
mov    0x40(%rsp),%rsi
mov    0x24(%rsp),%edx
movsd  (%rsp),%xmm0
mov    %r14,%r10
mov    %rbp,%rcx
mov    %rbx,%r8
cltq
addsd  0x0(%r13,%rax,8),%xmm0
add    $0x1,%r12d
cmp    %r12d,%edx
jne    1240 <check_inclusion+0x130>
mov    0x20(%rsp),%ebx
add    0x14(%rsp),%ebx
mov    0x38(%rsp),%rbp
mov    0x1c(%rsp),%r12d
mov    0x18(%rsp),%edi
jmp    1220 <check_inclusion+0x110>

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

LLVM IR (-O2):
define i32 @check_inclusion(i32 %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5, double* nocapture %6, i32 %7) {
%9 = sub nsw i32 %0, %7
%10 = icmp sgt i32 %9, 0
%11 = select i1 %10, i32 %9, i32 0
%12 = sub nsw i32 %1, %7
%13 = icmp sgt i32 %12, 0
%14 = select i1 %13, i32 %12, i32 0
%15 = sub nsw i32 %2, %7
%16 = icmp sgt i32 %15, 0
%17 = select i1 %16, i32 %15, i32 0
%18 = add nsw i32 %7, %0
%19 = add nsw i32 %3, -1
%20 = icmp sgt i32 %18, %19
%21 = select i1 %20, i32 %19, i32 %18
%22 = add nsw i32 %7, %1
%23 = add nsw i32 %4, -1
%24 = icmp sgt i32 %22, %23
%25 = select i1 %24, i32 %23, i32 %22
%26 = add nsw i32 %7, %2
%27 = add nsw i32 %5, -1
%28 = icmp sgt i32 %26, %27
%29 = select i1 %28, i32 %27, i32 %26
%30 = icmp sgt i32 %11, %21
br i1 %30, label %40, label %31
31: ; preds = %8
%32 = icmp sgt i32 %14, %25
%33 = icmp sgt i32 %17, %29
%34 = add nsw i32 %29, 1
%35 = sub i32 %34, %17
br label %36
36: ; preds = %31, %55
%37 = phi i32 [ %11, %31 ], [ %58, %55 ]
%38 = phi i32 [ 0, %31 ], [ %57, %55 ]
%39 = phi double [ 0.000000e+00, %31 ], [ %56, %55 ]
br i1 %32, label %55, label %51
40: ; preds = %55, %8
%41 = phi double [ 0.000000e+00, %8 ], [ %56, %55 ]
%42 = phi i32 [ 0, %8 ], [ %57, %55 ]
%43 = tail call i32 @map_index(i32 %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5)
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds double, double* %6, i64 %44
%46 = load double, double* %45
%47 = sitofp i32 %42 to double
%48 = fdiv double %41, %47
%49 = fcmp oge double %46, %48
%50 = zext i1 %49 to i32
ret i32 %50
51: ; preds = %36, %62
%52 = phi i32 [ %65, %62 ], [ %14, %36 ]
%53 = phi i32 [ %64, %62 ], [ %38, %36 ]
%54 = phi double [ %63, %62 ], [ %39, %36 ]
br i1 %33, label %62, label %67
55: ; preds = %62, %36
%56 = phi double [ %39, %36 ], [ %63, %62 ]
%57 = phi i32 [ %38, %36 ], [ %64, %62 ]
%58 = add nuw i32 %37, 1
%59 = icmp eq i32 %37, %21
br i1 %59, label %40, label %36
60: ; preds = %67
%61 = add i32 %35, %53
br label %62
62: ; preds = %60, %51
%63 = phi double [ %54, %51 ], [ %74, %60 ]
%64 = phi i32 [ %53, %51 ], [ %61, %60 ]
%65 = add nuw i32 %52, 1
%66 = icmp eq i32 %52, %25
br i1 %66, label %55, label %51
67: ; preds = %51, %67
%68 = phi i32 [ %75, %67 ], [ %17, %51 ]
%69 = phi double [ %74, %67 ], [ %54, %51 ]
%70 = tail call i32 @map_index(i32 %37, i32 %52, i32 %68, i32 %3, i32 %4, i32 %5)
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds double, double* %6, i64 %71
%73 = load double, double* %72
%74 = fadd double %69, %73
%75 = add nuw i32 %68, 1
%76 = icmp eq i32 %68, %29
br i1 %76, label %60, label %67
}