void wkeep_1D_center(double *sigIn, int sigInLength, double *sigOut,
                     int sigOutLength) {
  int count = 0;
  for (count = 0; count < sigOutLength; count++)
    sigOut[count] = sigIn[(sigInLength - sigOutLength) / 2 + count];
  return;
}


x86 assembly (clang, -O0):
<wkeep_1D_center>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    116e <wkeep_1D_center+0x6e>
mov    -0x8(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0xc(%rbp),%eax
sub    -0x1c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x28(%rbp),%rax
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1120 <wkeep_1D_center+0x20>
pop    %rbp
ret

x86 assembly (clang, -O2):
<wkeep_1D_center>:
test   %ecx,%ecx
jle    12a6 <wkeep_1D_center+0x1a6>
sub    %ecx,%esi
mov    %esi,%eax
shr    $0x1f,%eax
add    %esi,%eax
sar    %eax
movslq %eax,%r8
mov    %ecx,%r10d
cmp    $0x4,%ecx
jb     1140 <wkeep_1D_center+0x40>
lea    (%r8,%r10,1),%rax
lea    (%rdi,%rax,8),%rax
cmp    %rdx,%rax
jbe    11cc <wkeep_1D_center+0xcc>
lea    (%rdx,%r10,8),%rax
lea    (%rdi,%r8,8),%rcx
cmp    %rax,%rcx
jae    11cc <wkeep_1D_center+0xcc>
xor    %ecx,%ecx
mov    %rcx,%r9
not    %r9
add    %r10,%r9
mov    %r10,%rsi
and    $0x3,%rsi
je     1174 <wkeep_1D_center+0x74>
lea    (%rdi,%r8,8),%rax
nopl   0x0(%rax,%rax,1)
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rdx,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    1160 <wkeep_1D_center+0x60>
cmp    $0x3,%r9
jb     12a6 <wkeep_1D_center+0x1a6>
lea    (%rdi,%r8,8),%rax
add    $0x18,%rax
cs nopw 0x0(%rax,%rax,1)
movsd  -0x18(%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rdx,%rcx,8)
movsd  -0x10(%rax,%rcx,8),%xmm0
movsd  %xmm0,0x8(%rdx,%rcx,8)
movsd  -0x8(%rax,%rcx,8),%xmm0
movsd  %xmm0,0x10(%rdx,%rcx,8)
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,0x18(%rdx,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%r10
jne    1190 <wkeep_1D_center+0x90>
jmp    12a6 <wkeep_1D_center+0x1a6>
mov    %r10d,%ecx
and    $0xfffffffc,%ecx
lea    -0x4(%rcx),%rax
mov    %rax,%r11
shr    $0x2,%r11
add    $0x1,%r11
mov    %r11d,%r9d
and    $0x3,%r9d
cmp    $0xc,%rax
jae    11f2 <wkeep_1D_center+0xf2>
xor    %eax,%eax
jmp    1258 <wkeep_1D_center+0x158>
and    $0xfffffffffffffffc,%r11
lea    (%rdi,%r8,8),%rsi
add    $0x70,%rsi
xor    %eax,%eax
movups -0x70(%rsi,%rax,8),%xmm0
movups -0x60(%rsi,%rax,8),%xmm1
movups %xmm0,(%rdx,%rax,8)
movups %xmm1,0x10(%rdx,%rax,8)
movups -0x50(%rsi,%rax,8),%xmm0
movups -0x40(%rsi,%rax,8),%xmm1
movups %xmm0,0x20(%rdx,%rax,8)
movups %xmm1,0x30(%rdx,%rax,8)
movups -0x30(%rsi,%rax,8),%xmm0
movups -0x20(%rsi,%rax,8),%xmm1
movups %xmm0,0x40(%rdx,%rax,8)
movups %xmm1,0x50(%rdx,%rax,8)
movups -0x10(%rsi,%rax,8),%xmm0
movups (%rsi,%rax,8),%xmm1
movups %xmm0,0x60(%rdx,%rax,8)
movups %xmm1,0x70(%rdx,%rax,8)
add    $0x10,%rax
add    $0xfffffffffffffffc,%r11
jne    1200 <wkeep_1D_center+0x100>
test   %r9,%r9
je     129d <wkeep_1D_center+0x19d>
lea    (%rdx,%rax,8),%r11
add    $0x10,%r11
add    %r8,%rax
lea    (%rdi,%rax,8),%rax
add    $0x10,%rax
shl    $0x5,%r9
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
movups -0x10(%rax,%rsi,1),%xmm0
movups (%rax,%rsi,1),%xmm1
movups %xmm0,-0x10(%r11,%rsi,1)
movups %xmm1,(%r11,%rsi,1)
add    $0x20,%rsi
cmp    %rsi,%r9
jne    1280 <wkeep_1D_center+0x180>
cmp    %r10,%rcx
jne    1142 <wkeep_1D_center+0x42>
ret

LLVM IR (-O0):
define void @wkeep_1D_center(double* %0, i32 %1, double* %2, i32 %3) {
%5 = alloca double*
%6 = alloca i32
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
store double* %0, double** %5
store i32 %1, i32* %6
store double* %2, double** %7
store i32 %3, i32* %8
store i32 0, i32* %9
store i32 0, i32* %9
br label %10
10: ; preds = %29, %4
%11 = load i32, i32* %9
%12 = load i32, i32* %8
%13 = icmp slt i32 %11, %12
br i1 %13, label %14, label %32
14: ; preds = %10
%15 = load double*, double** %5
%16 = load i32, i32* %6
%17 = load i32, i32* %8
%18 = sub nsw i32 %16, %17
%19 = sdiv i32 %18, 2
%20 = load i32, i32* %9
%21 = add nsw i32 %19, %20
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds double, double* %15, i64 %22
%24 = load double, double* %23
%25 = load double*, double** %7
%26 = load i32, i32* %9
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
store double %24, double* %28
br label %29
29: ; preds = %14
%30 = load i32, i32* %9
%31 = add nsw i32 %30, 1
store i32 %31, i32* %9
br label %10
32: ; preds = %10
ret void
}

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