void compress4(short int *src, int *ncol, int *nrow) {
  int i, j, ncolo;
  ncolo = *ncol;
  *ncol /= 2;
  *nrow /= 2;
  for (i = 0; i < *nrow; ++i) {
    for (j = 0; j < *ncol; ++j) {
      *(src + i * (*ncol) + j) = (*(src + (2 * i + 1) * ncolo + 2 * j) +
                                  *(src + (2 * i + 1) * ncolo + 2 * j + 1) +
                                  *(src + (2 * i) * ncolo + 2 * j) +
                                  *(src + (2 * i) * ncolo + 2 * j + 1)) /
                                 4;
    }
  }
}


x86 assembly (clang, -O0):
<compress4>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x24(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    (%rax),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x38(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x18(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    (%rax),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x30(%rbp),%rax
mov    %ecx,(%rax)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x18(%rbp),%rcx
cmp    (%rcx),%eax
jge    1267 <compress4+0x167>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x10(%rbp),%rcx
cmp    (%rcx),%eax
jge    1254 <compress4+0x154>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
imul   -0x24(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x1,%rcx
add    %rcx,%rax
mov    -0x20(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movswl (%rax,%rcx,2),%eax
mov    -0x8(%rbp),%rcx
mov    -0x1c(%rbp),%edx
shl    $0x1,%edx
add    $0x1,%edx
imul   -0x24(%rbp),%edx
movslq %edx,%rdx
shl    $0x1,%rdx
add    %rdx,%rcx
mov    -0x20(%rbp),%edx
shl    $0x1,%edx
movslq %edx,%rdx
movswl 0x2(%rcx,%rdx,2),%ecx
add    %ecx,%eax
mov    -0x8(%rbp),%rcx
mov    -0x1c(%rbp),%edx
shl    $0x1,%edx
imul   -0x24(%rbp),%edx
movslq %edx,%rdx
shl    $0x1,%rdx
add    %rdx,%rcx
mov    -0x20(%rbp),%edx
shl    $0x1,%edx
movslq %edx,%rdx
movswl (%rcx,%rdx,2),%ecx
add    %ecx,%eax
mov    -0x8(%rbp),%rcx
mov    -0x1c(%rbp),%edx
shl    $0x1,%edx
imul   -0x24(%rbp),%edx
movslq %edx,%rdx
shl    $0x1,%rdx
add    %rdx,%rcx
mov    -0x20(%rbp),%edx
shl    $0x1,%edx
movslq %edx,%rdx
movswl 0x2(%rcx,%rdx,2),%ecx
add    %ecx,%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %ax,%dx
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    -0x10(%rbp),%rsi
imul   (%rsi),%ecx
movslq %ecx,%rcx
shl    $0x1,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    116a <compress4+0x6a>
jmp    1259 <compress4+0x159>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1154 <compress4+0x54>
pop    %rbp
ret

x86 assembly (clang, -O2):
<compress4>:
push   %rbp
push   %r15
push   %r14
push   %rbx
movslq (%rsi),%rax
mov    %eax,%ecx
shr    $0x1f,%ecx
add    %eax,%ecx
sar    %ecx
mov    %ecx,(%rsi)
mov    (%rdx),%ecx
mov    %ecx,%ebx
shr    $0x1f,%ebx
add    %ecx,%ebx
sar    %ebx
mov    %ebx,(%rdx)
cmp    $0x2,%ecx
jl     11a5 <compress4+0xa5>
mov    (%rsi),%r15d
movslq %r15d,%r10
mov    %ebx,%r8d
lea    0x2(%rdi),%r11
lea    0x0(,%rax,4),%r9
add    %r10,%r10
xor    %r14d,%r14d
jmp    115f <compress4+0x5f>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r14
add    %r9,%r11
add    %r10,%rdi
cmp    %r8,%r14
je     11a5 <compress4+0xa5>
test   %r15d,%r15d
jle    1150 <compress4+0x50>
mov    %r11,%rsi
xor    %ebx,%ebx
nopl   0x0(%rax)
movswl -0x2(%rsi,%rax,2),%ebp
movswl (%rsi,%rax,2),%edx
add    %ebp,%edx
movswl -0x2(%rsi),%ebp
add    %edx,%ebp
movswl (%rsi),%edx
lea    (%rdx,%rbp,1),%ecx
add    $0x3,%ecx
add    %edx,%ebp
cmovns %ebp,%ecx
shr    $0x2,%ecx
mov    %cx,(%rdi,%rbx,2)
add    $0x1,%rbx
add    $0x4,%rsi
cmp    %rbx,%r15
jne    1170 <compress4+0x70>
jmp    1150 <compress4+0x50>
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @compress4(i16* %0, i32* %1, i32* %2) {
%4 = alloca i16*
%5 = alloca i32*
%6 = alloca i32*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
store i16* %0, i16** %4
store i32* %1, i32** %5
store i32* %2, i32** %6
%10 = load i32*, i32** %5
%11 = load i32, i32* %10
store i32 %11, i32* %9
%12 = load i32*, i32** %5
%13 = load i32, i32* %12
%14 = sdiv i32 %13, 2
store i32 %14, i32* %12
%15 = load i32*, i32** %6
%16 = load i32, i32* %15
%17 = sdiv i32 %16, 2
store i32 %17, i32* %15
store i32 0, i32* %7
br label %18
18: ; preds = %105, %3
%19 = load i32, i32* %7
%20 = load i32*, i32** %6
%21 = load i32, i32* %20
%22 = icmp slt i32 %19, %21
br i1 %22, label %23, label %108
23: ; preds = %18
store i32 0, i32* %8
br label %24
24: ; preds = %101, %23
%25 = load i32, i32* %8
%26 = load i32*, i32** %5
%27 = load i32, i32* %26
%28 = icmp slt i32 %25, %27
br i1 %28, label %29, label %104
29: ; preds = %24
%30 = load i16*, i16** %4
%31 = load i32, i32* %7
%32 = mul nsw i32 2, %31
%33 = add nsw i32 %32, 1
%34 = load i32, i32* %9
%35 = mul nsw i32 %33, %34
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds i16, i16* %30, i64 %36
%38 = load i32, i32* %8
%39 = mul nsw i32 2, %38
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds i16, i16* %37, i64 %40
%42 = load i16, i16* %41
%43 = sext i16 %42 to i32
%44 = load i16*, i16** %4
%45 = load i32, i32* %7
%46 = mul nsw i32 2, %45
%47 = add nsw i32 %46, 1
%48 = load i32, i32* %9
%49 = mul nsw i32 %47, %48
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds i16, i16* %44, i64 %50
%52 = load i32, i32* %8
%53 = mul nsw i32 2, %52
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds i16, i16* %51, i64 %54
%56 = getelementptr inbounds i16, i16* %55, i64 1
%57 = load i16, i16* %56
%58 = sext i16 %57 to i32
%59 = add nsw i32 %43, %58
%60 = load i16*, i16** %4
%61 = load i32, i32* %7
%62 = mul nsw i32 2, %61
%63 = load i32, i32* %9
%64 = mul nsw i32 %62, %63
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds i16, i16* %60, i64 %65
%67 = load i32, i32* %8
%68 = mul nsw i32 2, %67
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds i16, i16* %66, i64 %69
%71 = load i16, i16* %70
%72 = sext i16 %71 to i32
%73 = add nsw i32 %59, %72
%74 = load i16*, i16** %4
%75 = load i32, i32* %7
%76 = mul nsw i32 2, %75
%77 = load i32, i32* %9
%78 = mul nsw i32 %76, %77
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds i16, i16* %74, i64 %79
%81 = load i32, i32* %8
%82 = mul nsw i32 2, %81
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds i16, i16* %80, i64 %83
%85 = getelementptr inbounds i16, i16* %84, i64 1
%86 = load i16, i16* %85
%87 = sext i16 %86 to i32
%88 = add nsw i32 %73, %87
%89 = sdiv i32 %88, 4
%90 = trunc i32 %89 to i16
%91 = load i16*, i16** %4
%92 = load i32, i32* %7
%93 = load i32*, i32** %5
%94 = load i32, i32* %93
%95 = mul nsw i32 %92, %94
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds i16, i16* %91, i64 %96
%98 = load i32, i32* %8
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds i16, i16* %97, i64 %99
store i16 %90, i16* %100
br label %101
101: ; preds = %29
%102 = load i32, i32* %8
%103 = add nsw i32 %102, 1
store i32 %103, i32* %8
br label %24
104: ; preds = %24
br label %105
105: ; preds = %104
%106 = load i32, i32* %7
%107 = add nsw i32 %106, 1
store i32 %107, i32* %7
br label %18
108: ; preds = %18
ret void
}

LLVM IR (-O2):
define void @compress4(i16* nocapture %0, i32* nocapture %1, i32* nocapture %2) {
%4 = load i32, i32* %1
%5 = sdiv i32 %4, 2
store i32 %5, i32* %1
%6 = load i32, i32* %2
%7 = sdiv i32 %6, 2
store i32 %7, i32* %2
%8 = icmp sgt i32 %6, 1
br i1 %8, label %9, label %53
9: ; preds = %3
%10 = load i32, i32* %1
%11 = icmp sgt i32 %10, 0
%12 = sext i32 %4 to i64
%13 = sext i32 %10 to i64
%14 = zext i32 %7 to i64
%15 = zext i32 %10 to i64
br label %16
16: ; preds = %9, %50
%17 = phi i64 [ 0, %9 ], [ %51, %50 ]
br i1 %11, label %18, label %50
18: ; preds = %16
%19 = shl nuw nsw i64 %17, 1
%20 = or i64 %19, 1
%21 = mul nsw i64 %20, %12
%22 = getelementptr inbounds i16, i16* %0, i64 %21
%23 = mul nsw i64 %19, %12
%24 = getelementptr inbounds i16, i16* %0, i64 %23
%25 = mul nsw i64 %17, %13
%26 = getelementptr inbounds i16, i16* %0, i64 %25
br label %27
27: ; preds = %18, %27
%28 = phi i64 [ 0, %18 ], [ %48, %27 ]
%29 = shl nuw nsw i64 %28, 1
%30 = getelementptr inbounds i16, i16* %22, i64 %29
%31 = load i16, i16* %30
%32 = sext i16 %31 to i32
%33 = getelementptr inbounds i16, i16* %30, i64 1
%34 = load i16, i16* %33
%35 = sext i16 %34 to i32
%36 = add nsw i32 %35, %32
%37 = getelementptr inbounds i16, i16* %24, i64 %29
%38 = load i16, i16* %37
%39 = sext i16 %38 to i32
%40 = add nsw i32 %36, %39
%41 = getelementptr inbounds i16, i16* %37, i64 1
%42 = load i16, i16* %41
%43 = sext i16 %42 to i32
%44 = add nsw i32 %40, %43
%45 = sdiv i32 %44, 4
%46 = trunc i32 %45 to i16
%47 = getelementptr inbounds i16, i16* %26, i64 %28
store i16 %46, i16* %47
%48 = add nuw nsw i64 %28, 1
%49 = icmp eq i64 %48, %15
br i1 %49, label %50, label %27
50: ; preds = %27, %16
%51 = add nuw nsw i64 %17, 1
%52 = icmp eq i64 %51, %14
br i1 %52, label %53, label %16
53: ; preds = %50, %3
ret void
}