void dyaddown_2D_keep_even_col(double *matrixIn, int matrixInRow,
                               int matrixInCol, double *matrixOut,
                               int matrixOutRow, int matrixOutCol) {
  int row, col;
  for (col = 0; col < matrixOutCol; col++) {
    for (row = 0; row < matrixInRow; row++)
      matrixOut[row + col * matrixInRow] =
          matrixIn[row + (2 * col + 1) * matrixInRow];
  }
  return;
}


x86 assembly (clang, -O0):
<dyaddown_2D_keep_even_col>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9d,-0x20(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1197 <dyaddown_2D_keep_even_col+0x97>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1184 <dyaddown_2D_keep_even_col+0x84>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x28(%rbp),%edx
shl    $0x1,%edx
add    $0x1,%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x28(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1134 <dyaddown_2D_keep_even_col+0x34>
jmp    1189 <dyaddown_2D_keep_even_col+0x89>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1121 <dyaddown_2D_keep_even_col+0x21>
pop    %rbp
ret

x86 assembly (clang, -O2):
<dyaddown_2D_keep_even_col>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rcx,-0x50(%rsp)
mov    %rdi,-0x58(%rsp)
test   %r9d,%r9d
jle    137e <dyaddown_2D_keep_even_col+0x27e>
mov    %esi,%r11d
movslq %esi,%rsi
mov    %r9d,%r9d
mov    %esi,%eax
lea    (%rsi,%rsi,1),%r15d
mov    -0x58(%rsp),%rdx
lea    (%rdx,%rax,8),%rcx
mov    %rcx,-0x10(%rsp)
mov    %eax,%ecx
and    $0xfffffffc,%ecx
mov    %rcx,-0x20(%rsp)
add    $0xfffffffffffffffc,%rcx
mov    %rcx,-0x28(%rsp)
mov    %rcx,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
mov    %rdi,%rcx
mov    %rdi,-0x30(%rsp)
and    $0xfffffffffffffffe,%rdi
mov    %rdi,-0x38(%rsp)
mov    %eax,%r8d
and    $0x3,%r8d
lea    0x30(%rdx),%rcx
mov    %rcx,-0x40(%rsp)
mov    -0x50(%rsp),%rcx
lea    0x30(%rcx),%rbp
mov    %rsi,-0x8(%rsp)
lea    0x0(,%rsi,8),%r13
add    $0x18,%rdx
mov    %rdx,-0x18(%rsp)
xor    %r14d,%r14d
mov    %r11d,%r10d
mov    %r8,-0x48(%rsp)
jmp    11c6 <dyaddown_2D_keep_even_col+0xc6>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r14
add    %r15d,%r10d
add    %r13,%rbp
add    %r13,%rcx
cmp    %r9,%r14
je     137e <dyaddown_2D_keep_even_col+0x27e>
test   %r11d,%r11d
jle    11b0 <dyaddown_2D_keep_even_col+0xb0>
movslq %r10d,%r12
cmp    $0x4,%r11d
jae    11e0 <dyaddown_2D_keep_even_col+0xe0>
xor    %ebx,%ebx
jmp    12e0 <dyaddown_2D_keep_even_col+0x1e0>
nopl   0x0(%rax,%rax,1)
mov    %r14,%r8
imul   -0x8(%rsp),%r8
mov    %r15d,%edx
imul   %r14d,%edx
add    %r11d,%edx
mov    -0x50(%rsp),%rbx
lea    (%rbx,%r8,8),%rsi
movslq %edx,%rdx
mov    -0x10(%rsp),%rdi
lea    (%rdi,%rdx,8),%rdi
cmp    %rdi,%rsi
jae    122f <dyaddown_2D_keep_even_col+0x12f>
lea    (%r8,%rax,1),%rsi
lea    (%rbx,%rsi,8),%rsi
mov    -0x58(%rsp),%rdi
lea    (%rdi,%rdx,8),%rdx
cmp    %rsi,%rdx
jae    122f <dyaddown_2D_keep_even_col+0x12f>
xor    %ebx,%ebx
mov    -0x48(%rsp),%r8
jmp    12e0 <dyaddown_2D_keep_even_col+0x1e0>
mov    %r9,%rdi
cmpq   $0x0,-0x28(%rsp)
je     136c <dyaddown_2D_keep_even_col+0x26c>
mov    -0x40(%rsp),%rdx
lea    (%rdx,%r12,8),%rdx
mov    -0x38(%rsp),%r9
xor    %ebx,%ebx
xchg   %ax,%ax
movups -0x30(%rdx,%rbx,8),%xmm0
movups -0x20(%rdx,%rbx,8),%xmm1
movups %xmm0,-0x30(%rbp,%rbx,8)
movups %xmm1,-0x20(%rbp,%rbx,8)
movups -0x10(%rdx,%rbx,8),%xmm0
movups (%rdx,%rbx,8),%xmm1
movups %xmm0,-0x10(%rbp,%rbx,8)
movups %xmm1,0x0(%rbp,%rbx,8)
add    $0x8,%rbx
add    $0xfffffffffffffffe,%r9
jne    1250 <dyaddown_2D_keep_even_col+0x150>
testb  $0x1,-0x30(%rsp)
je     12ba <dyaddown_2D_keep_even_col+0x1ba>
lea    (%r14,%r14,1),%edx
add    $0x1,%edx
imul   %r11d,%edx
movslq %edx,%rdx
add    %rbx,%rdx
mov    -0x58(%rsp),%rsi
movups (%rsi,%rdx,8),%xmm0
movups 0x10(%rsi,%rdx,8),%xmm1
add    %rbx,%r8
mov    -0x50(%rsp),%rdx
movups %xmm0,(%rdx,%r8,8)
movups %xmm1,0x10(%rdx,%r8,8)
mov    -0x20(%rsp),%rdx
mov    %rdx,%rbx
cmp    %rax,%rdx
mov    %rdi,%r9
mov    -0x48(%rsp),%r8
je     11b0 <dyaddown_2D_keep_even_col+0xb0>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %rbx,%rdx
not    %rdx
add    %rax,%rdx
test   %r8,%r8
je     1314 <dyaddown_2D_keep_even_col+0x214>
mov    -0x58(%rsp),%rsi
lea    (%rsi,%r12,8),%rdi
mov    %r8,%rsi
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rbx,8),%xmm0
movsd  %xmm0,(%rcx,%rbx,8)
add    $0x1,%rbx
add    $0xffffffffffffffff,%rsi
jne    1300 <dyaddown_2D_keep_even_col+0x200>
cmp    $0x3,%rdx
jb     11b0 <dyaddown_2D_keep_even_col+0xb0>
mov    -0x18(%rsp),%rdx
lea    (%rdx,%r12,8),%rdx
nopw   0x0(%rax,%rax,1)
movsd  -0x18(%rdx,%rbx,8),%xmm0
movsd  %xmm0,(%rcx,%rbx,8)
movsd  -0x10(%rdx,%rbx,8),%xmm0
movsd  %xmm0,0x8(%rcx,%rbx,8)
movsd  -0x8(%rdx,%rbx,8),%xmm0
movsd  %xmm0,0x10(%rcx,%rbx,8)
movsd  (%rdx,%rbx,8),%xmm0
movsd  %xmm0,0x18(%rcx,%rbx,8)
add    $0x4,%rbx
cmp    %rbx,%rax
jne    1330 <dyaddown_2D_keep_even_col+0x230>
jmp    11b0 <dyaddown_2D_keep_even_col+0xb0>
xor    %ebx,%ebx
testb  $0x1,-0x30(%rsp)
jne    1288 <dyaddown_2D_keep_even_col+0x188>
jmp    12ba <dyaddown_2D_keep_even_col+0x1ba>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @dyaddown_2D_keep_even_col(double* %0, i32 %1, i32 %2, double* %3, i32 %4, i32 %5) {
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca double*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
store double* %0, double** %7
store i32 %1, i32* %8
store i32 %2, i32* %9
store double* %3, double** %10
store i32 %4, i32* %11
store i32 %5, i32* %12
store i32 0, i32* %14
br label %15
15: ; preds = %48, %6
%16 = load i32, i32* %14
%17 = load i32, i32* %12
%18 = icmp slt i32 %16, %17
br i1 %18, label %19, label %51
19: ; preds = %15
store i32 0, i32* %13
br label %20
20: ; preds = %44, %19
%21 = load i32, i32* %13
%22 = load i32, i32* %8
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %47
24: ; preds = %20
%25 = load double*, double** %7
%26 = load i32, i32* %13
%27 = load i32, i32* %14
%28 = mul nsw i32 2, %27
%29 = add nsw i32 %28, 1
%30 = load i32, i32* %8
%31 = mul nsw i32 %29, %30
%32 = add nsw i32 %26, %31
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds double, double* %25, i64 %33
%35 = load double, double* %34
%36 = load double*, double** %10
%37 = load i32, i32* %13
%38 = load i32, i32* %14
%39 = load i32, i32* %8
%40 = mul nsw i32 %38, %39
%41 = add nsw i32 %37, %40
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds double, double* %36, i64 %42
store double %35, double* %43
br label %44
44: ; preds = %24
%45 = load i32, i32* %13
%46 = add nsw i32 %45, 1
store i32 %46, i32* %13
br label %20
47: ; preds = %20
br label %48
48: ; preds = %47
%49 = load i32, i32* %14
%50 = add nsw i32 %49, 1
store i32 %50, i32* %14
br label %15
51: ; preds = %15
ret void
}

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