void dyaddown_2D_keep_odd_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 + col * matrixInRow * 2];
  }
  return;
}


x86 assembly (clang, -O0):
<dyaddown_2D_keep_odd_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    1194 <dyaddown_2D_keep_odd_col+0x94>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1181 <dyaddown_2D_keep_odd_col+0x81>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    -0x28(%rbp),%edx
imul   -0xc(%rbp),%edx
shl    $0x1,%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_odd_col+0x34>
jmp    1186 <dyaddown_2D_keep_odd_col+0x86>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1121 <dyaddown_2D_keep_odd_col+0x21>
pop    %rbp
ret

x86 assembly (clang, -O2):
<dyaddown_2D_keep_odd_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    136e <dyaddown_2D_keep_odd_col+0x26e>
movslq %esi,%rdx
mov    %r9d,%r9d
mov    %edx,%r11d
lea    (%rdx,%rdx,1),%r15d
mov    -0x58(%rsp),%rax
lea    (%rax,%r11,8),%rcx
mov    %rcx,-0x10(%rsp)
mov    %r11d,%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    %r11d,%r8d
and    $0x3,%r8d
lea    0x30(%rax),%rcx
mov    %rcx,-0x40(%rsp)
mov    -0x50(%rsp),%rcx
lea    0x30(%rcx),%rbp
mov    %rdx,-0x8(%rsp)
lea    0x0(,%rdx,8),%r13
add    $0x18,%rax
mov    %rax,-0x18(%rsp)
xor    %r14d,%r14d
xor    %r10d,%r10d
mov    %r8,-0x48(%rsp)
jmp    11c6 <dyaddown_2D_keep_odd_col+0xc6>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%r10
add    %r15d,%r14d
add    %r13,%rbp
add    %r13,%rcx
cmp    %r9,%r10
je     136e <dyaddown_2D_keep_odd_col+0x26e>
test   %esi,%esi
jle    11b0 <dyaddown_2D_keep_odd_col+0xb0>
movslq %r14d,%r12
cmp    $0x4,%esi
jae    11e0 <dyaddown_2D_keep_odd_col+0xe0>
xor    %ebx,%ebx
jmp    12d0 <dyaddown_2D_keep_odd_col+0x1d0>
nopl   0x0(%rax)
mov    %r10,%r8
imul   -0x8(%rsp),%r8
mov    %r15d,%eax
imul   %r10d,%eax
mov    -0x50(%rsp),%rbx
lea    (%rbx,%r8,8),%rdx
cltq
mov    -0x10(%rsp),%rdi
lea    (%rdi,%rax,8),%rdi
cmp    %rdi,%rdx
jae    122b <dyaddown_2D_keep_odd_col+0x12b>
lea    (%r8,%r11,1),%rdx
lea    (%rbx,%rdx,8),%rdx
mov    -0x58(%rsp),%rdi
lea    (%rdi,%rax,8),%rax
cmp    %rdx,%rax
jae    122b <dyaddown_2D_keep_odd_col+0x12b>
xor    %ebx,%ebx
mov    -0x48(%rsp),%r8
jmp    12d0 <dyaddown_2D_keep_odd_col+0x1d0>
mov    %r9,%rdi
cmpq   $0x0,-0x28(%rsp)
je     135c <dyaddown_2D_keep_odd_col+0x25c>
mov    -0x40(%rsp),%rax
lea    (%rax,%r12,8),%rdx
mov    -0x38(%rsp),%r9
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
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_odd_col+0x150>
testb  $0x1,-0x30(%rsp)
je     12b0 <dyaddown_2D_keep_odd_col+0x1b0>
lea    (%r8,%r8,1),%eax
cltq
add    %rbx,%rax
mov    -0x58(%rsp),%rdx
movups (%rdx,%rax,8),%xmm0
movups 0x10(%rdx,%rax,8),%xmm1
add    %r8,%rbx
mov    -0x50(%rsp),%rax
movups %xmm0,(%rax,%rbx,8)
movups %xmm1,0x10(%rax,%rbx,8)
mov    -0x20(%rsp),%rax
mov    %rax,%rbx
cmp    %r11,%rax
mov    %rdi,%r9
mov    -0x48(%rsp),%r8
je     11b0 <dyaddown_2D_keep_odd_col+0xb0>
nopl   0x0(%rax)
mov    %rbx,%rdx
not    %rdx
add    %r11,%rdx
test   %r8,%r8
je     1304 <dyaddown_2D_keep_odd_col+0x204>
mov    -0x58(%rsp),%rax
lea    (%rax,%r12,8),%rdi
mov    %r8,%rax
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rbx,8),%xmm0
movsd  %xmm0,(%rcx,%rbx,8)
add    $0x1,%rbx
add    $0xffffffffffffffff,%rax
jne    12f0 <dyaddown_2D_keep_odd_col+0x1f0>
cmp    $0x3,%rdx
jb     11b0 <dyaddown_2D_keep_odd_col+0xb0>
mov    -0x18(%rsp),%rax
lea    (%rax,%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,%r11
jne    1320 <dyaddown_2D_keep_odd_col+0x220>
jmp    11b0 <dyaddown_2D_keep_odd_col+0xb0>
xor    %ebx,%ebx
testb  $0x1,-0x30(%rsp)
jne    1288 <dyaddown_2D_keep_odd_col+0x188>
jmp    12b0 <dyaddown_2D_keep_odd_col+0x1b0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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