int is_local_maximum(int start, int size, int type, int p_nb_genes, int p_max,
                     double ***store) {
  int i, j, result = 1, hood = 5;
  for (i = start - hood; i <= start + hood; i++) {
    for (j = size - hood; j <= size + hood; j++) {
      if (i >= 0 && j >= 0 && i < p_nb_genes && j <= p_max) {
        if (store[i][j][type] > store[start][size][type]) {
          result = 0;
        }
      }
    }
  }
  return result;
}


x86 assembly (clang, -O0):
<is_local_maximum>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8d,-0x14(%rbp)
mov    %r9,-0x20(%rbp)
movl   $0x1,-0x2c(%rbp)
movl   $0x5,-0x30(%rbp)
mov    -0x4(%rbp),%eax
sub    -0x30(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x4(%rbp),%ecx
add    -0x30(%rbp),%ecx
cmp    %ecx,%eax
jg     11f8 <is_local_maximum+0xf8>
mov    -0x8(%rbp),%eax
sub    -0x30(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x8(%rbp),%ecx
add    -0x30(%rbp),%ecx
cmp    %ecx,%eax
jg     11e5 <is_local_maximum+0xe5>
cmpl   $0x0,-0x24(%rbp)
jl     11d2 <is_local_maximum+0xd2>
cmpl   $0x0,-0x28(%rbp)
jl     11d2 <is_local_maximum+0xd2>
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    11d2 <is_local_maximum+0xd2>
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     11d2 <is_local_maximum+0xd2>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x8(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    11cd <is_local_maximum+0xcd>
movl   $0x0,-0x2c(%rbp)
jmp    11d2 <is_local_maximum+0xd2>
jmp    11d7 <is_local_maximum+0xd7>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1149 <is_local_maximum+0x49>
jmp    11ea <is_local_maximum+0xea>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    112f <is_local_maximum+0x2f>
mov    -0x2c(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<is_local_maximum>:
push   %r15
push   %r14
push   %r12
push   %rbx
lea    -0x5(%rdi),%eax
lea    -0x5(%rsi),%ebx
movslq %edx,%r12
movslq %edi,%r10
add    $0x5,%edi
movslq %esi,%r11
add    $0x5,%esi
movslq %ebx,%r14
movslq %r8d,%r8
movslq %esi,%rsi
movslq %eax,%rbx
movslq %ecx,%rcx
movslq %edi,%r15
mov    $0x1,%eax
jmp    114c <is_local_maximum+0x4c>
cs nopw 0x0(%rax,%rax,1)
nop
lea    0x1(%rbx),%rdx
cmp    %r15,%rbx
mov    %rdx,%rbx
jge    11a2 <is_local_maximum+0xa2>
mov    %r14,%rdi
jmp    116c <is_local_maximum+0x6c>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    0x1(%rdi),%rdx
cmp    %rsi,%rdi
mov    %rdx,%rdi
jge    1140 <is_local_maximum+0x40>
cmp    %rcx,%rbx
jge    1160 <is_local_maximum+0x60>
test   %rbx,%rbx
js     1160 <is_local_maximum+0x60>
test   %rdi,%rdi
js     1160 <is_local_maximum+0x60>
cmp    %r8,%rdi
jg     1160 <is_local_maximum+0x60>
mov    (%r9,%rbx,8),%rdx
mov    (%rdx,%rdi,8),%rdx
movsd  (%rdx,%r12,8),%xmm0
mov    (%r9,%r10,8),%rdx
mov    (%rdx,%r11,8),%rdx
ucomisd (%rdx,%r12,8),%xmm0
jbe    1160 <is_local_maximum+0x60>
xor    %eax,%eax
jmp    1160 <is_local_maximum+0x60>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

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

LLVM IR (-O2):
define i32 @is_local_maximum(i32 %0, i32 %1, i32 %2, i32 %3, i32 %4, double*** nocapture %5) {
%7 = add i32 %0, -5
%8 = add nsw i32 %0, 5
%9 = add i32 %1, -5
%10 = add nsw i32 %1, 5
%11 = sext i32 %2 to i64
%12 = sext i32 %0 to i64
%13 = getelementptr inbounds double**, double*** %5, i64 %12
%14 = sext i32 %1 to i64
%15 = sext i32 %9 to i64
%16 = sext i32 %4 to i64
%17 = sext i32 %10 to i64
%18 = sext i32 %7 to i64
%19 = sext i32 %3 to i64
%20 = sext i32 %8 to i64
br label %21
21: ; preds = %6, %53
%22 = phi i64 [ %18, %6 ], [ %54, %53 ]
%23 = phi i32 [ 1, %6 ], [ %50, %53 ]
%24 = icmp sgt i64 %22, -1
%25 = icmp slt i64 %22, %19
%26 = getelementptr inbounds double**, double*** %5, i64 %22
br label %27
27: ; preds = %21, %49
%28 = phi i64 [ %15, %21 ], [ %51, %49 ]
%29 = phi i32 [ %23, %21 ], [ %50, %49 ]
%30 = icmp sgt i64 %28, -1
%31 = select i1 %24, i1 %30, i1 false
%32 = and i1 %25, %31
%33 = xor i1 %32, true
%34 = icmp sgt i64 %28, %16
%35 = select i1 %33, i1 true, i1 %34
br i1 %35, label %49, label %36
36: ; preds = %27
%37 = load double**, double*** %26
%38 = getelementptr inbounds double*, double** %37, i64 %28
%39 = load double*, double** %38
%40 = getelementptr inbounds double, double* %39, i64 %11
%41 = load double, double* %40
%42 = load double**, double*** %13
%43 = getelementptr inbounds double*, double** %42, i64 %14
%44 = load double*, double** %43
%45 = getelementptr inbounds double, double* %44, i64 %11
%46 = load double, double* %45
%47 = fcmp ogt double %41, %46
br i1 %47, label %48, label %49
48: ; preds = %36
br label %49
49: ; preds = %27, %48, %36
%50 = phi i32 [ 0, %48 ], [ %29, %36 ], [ %29, %27 ]
%51 = add nsw i64 %28, 1
%52 = icmp slt i64 %28, %17
br i1 %52, label %27, label %53
53: ; preds = %49
%54 = add nsw i64 %22, 1
%55 = icmp slt i64 %22, %20
br i1 %55, label %21, label %56
56: ; preds = %53
ret i32 %50
}