void remove2columns(float *pt, int Nnx, int Nnyz) {
  int ix, iyz;
  float *tmp1, *tmp2;
  for (iyz = 0; iyz < Nnyz; iyz++) {
    tmp1 = pt + Nnx * iyz;
    tmp2 = pt + (Nnx + 2) * iyz;
    for (ix = 0; ix < Nnx; ix++)
      *tmp1++ = *tmp2++;
  }
}


x86 assembly (clang, -O0):
<remove2columns>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    11b0 <remove2columns+0xb0>
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
add    $0x2,%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    119d <remove2columns+0x9d>
mov    -0x28(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x28(%rbp)
movss  (%rax),%xmm0
mov    -0x20(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x20(%rbp)
movss  %xmm0,(%rax)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    115d <remove2columns+0x5d>
jmp    11a2 <remove2columns+0xa2>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1115 <remove2columns+0x15>
pop    %rbp
ret

x86 assembly (clang, -O2):
<remove2columns>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edx,%edx
jle    13cb <remove2columns+0x2cb>
lea    0x2(%rsi),%eax
movslq %esi,%rcx
cltq
mov    %edx,%edx
mov    %rdx,-0x8(%rsp)
lea    -0x1(%rcx),%edx
mov    %rdx,-0x30(%rsp)
add    $0x1,%rdx
mov    %rdx,-0x58(%rsp)
and    $0xfffffffffffffff8,%rdx
mov    %rdx,-0x60(%rsp)
add    $0xfffffffffffffff8,%rdx
mov    %rdx,-0x40(%rsp)
shr    $0x3,%rdx
add    $0x1,%rdx
mov    %edx,%ebp
and    $0x3,%ebp
and    $0xfffffffffffffffc,%rdx
mov    %rdx,-0x50(%rsp)
lea    0x70(%rdi),%r10
mov    %rax,-0x28(%rsp)
lea    0x0(,%rax,4),%rax
mov    %rax,-0x10(%rsp)
mov    %rcx,-0x20(%rsp)
lea    0x0(,%rcx,4),%rax
mov    %rax,-0x18(%rsp)
mov    %rbp,-0x48(%rsp)
shl    $0x5,%rbp
xor    %r11d,%r11d
mov    %rdi,%r15
mov    %rdi,%r14
mov    %r10,%rbx
mov    %rdi,-0x38(%rsp)
jmp    11c5 <remove2columns+0xc5>
xchg   %ax,%ax
add    $0x1,%r11
mov    -0x10(%rsp),%rax
add    %rax,%rbx
mov    -0x18(%rsp),%rcx
add    %rcx,%r10
add    %rcx,%r14
add    %rax,%r15
cmp    -0x8(%rsp),%r11
je     13cb <remove2columns+0x2cb>
test   %esi,%esi
jle    11a0 <remove2columns+0xa0>
mov    %r11,%rax
imul   -0x20(%rsp),%rax
mov    %r11,%rcx
imul   -0x28(%rsp),%rcx
lea    (%rdi,%rax,4),%r8
lea    (%rdi,%rcx,4),%r9
xor    %r13d,%r13d
cmpl   $0x7,-0x30(%rsp)
jb     12f0 <remove2columns+0x1f0>
mov    -0x58(%rsp),%rdx
add    %rdx,%rcx
lea    (%rdi,%rcx,4),%rcx
cmp    %rcx,%r8
jae    1212 <remove2columns+0x112>
add    %rdx,%rax
lea    (%rdi,%rax,4),%rax
cmp    %rax,%r9
jb     12f0 <remove2columns+0x1f0>
cmpq   $0x18,-0x40(%rsp)
jae    121e <remove2columns+0x11e>
xor    %edx,%edx
jmp    1290 <remove2columns+0x190>
mov    -0x50(%rsp),%rax
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nop
movups -0x70(%rbx,%rdx,4),%xmm0
movups -0x60(%rbx,%rdx,4),%xmm1
movups %xmm0,-0x70(%r10,%rdx,4)
movups %xmm1,-0x60(%r10,%rdx,4)
movups -0x50(%rbx,%rdx,4),%xmm0
movups -0x40(%rbx,%rdx,4),%xmm1
movups %xmm0,-0x50(%r10,%rdx,4)
movups %xmm1,-0x40(%r10,%rdx,4)
movups -0x30(%rbx,%rdx,4),%xmm0
movups -0x20(%rbx,%rdx,4),%xmm1
movups %xmm0,-0x30(%r10,%rdx,4)
movups %xmm1,-0x20(%r10,%rdx,4)
movups -0x10(%rbx,%rdx,4),%xmm0
movups (%rbx,%rdx,4),%xmm1
movups %xmm0,-0x10(%r10,%rdx,4)
movups %xmm1,(%r10,%rdx,4)
add    $0x20,%rdx
add    $0xfffffffffffffffc,%rax
jne    1230 <remove2columns+0x130>
cmpq   $0x0,-0x48(%rsp)
je     12cb <remove2columns+0x1cb>
lea    (%r14,%rdx,4),%rax
lea    (%r15,%rdx,4),%rdx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movups (%rdx,%rcx,1),%xmm0
movups 0x10(%rdx,%rcx,1),%xmm1
movups %xmm0,(%rax,%rcx,1)
movups %xmm1,0x10(%rax,%rcx,1)
add    $0x20,%rcx
cmp    %rcx,%rbp
jne    12b0 <remove2columns+0x1b0>
mov    -0x60(%rsp),%rax
cmp    %rax,-0x58(%rsp)
je     11a0 <remove2columns+0xa0>
mov    -0x60(%rsp),%rax
lea    (%r9,%rax,4),%r9
lea    (%r8,%rax,4),%r8
mov    %eax,%r13d
nopl   0x0(%rax,%rax,1)
mov    %esi,%edx
sub    %r13d,%edx
mov    %r13d,%r12d
not    %r12d
add    %esi,%r12d
and    $0x7,%edx
je     132c <remove2columns+0x22c>
neg    %edx
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
movss  (%r9),%xmm0
add    $0x4,%r9
movss  %xmm0,(%r8)
add    $0x4,%r8
add    $0xffffffff,%eax
cmp    %eax,%edx
jne    1310 <remove2columns+0x210>
sub    %eax,%r13d
cmp    $0x7,%r12d
mov    -0x38(%rsp),%rdi
jb     11a0 <remove2columns+0xa0>
mov    %esi,%edx
sub    %r13d,%edx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%r9,%rax,4),%xmm0
movss  %xmm0,(%r8,%rax,4)
movss  0x4(%r9,%rax,4),%xmm0
movss  %xmm0,0x4(%r8,%rax,4)
movss  0x8(%r9,%rax,4),%xmm0
movss  %xmm0,0x8(%r8,%rax,4)
movss  0xc(%r9,%rax,4),%xmm0
movss  %xmm0,0xc(%r8,%rax,4)
movss  0x10(%r9,%rax,4),%xmm0
movss  %xmm0,0x10(%r8,%rax,4)
movss  0x14(%r9,%rax,4),%xmm0
movss  %xmm0,0x14(%r8,%rax,4)
movss  0x18(%r9,%rax,4),%xmm0
movss  %xmm0,0x18(%r8,%rax,4)
movss  0x1c(%r9,%rax,4),%xmm0
movss  %xmm0,0x1c(%r8,%rax,4)
add    $0x8,%rax
cmp    %eax,%edx
jne    1350 <remove2columns+0x250>
jmp    11a0 <remove2columns+0xa0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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