void window_cell(int nwpv, double *array1, int gridID1, double *array2,
                 long gridsize2, int *cellidx) {
  long i;
  if (nwpv == 2) {
    for (i = 0; i < gridsize2; ++i) {
      array2[i * 2] = array1[cellidx[i] * 2];
      array2[i * 2 + 1] = array1[cellidx[i] * 2 + 1];
    }
  } else {
    for (i = 0; i < gridsize2; ++i)
      array2[i] = array1[cellidx[i]];
  }
}


x86 assembly (clang, -O0):
<window_cell>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
cmpl   $0x2,-0x4(%rbp)
jne    11aa <window_cell+0xaa>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    11a5 <window_cell+0xa5>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    -0x38(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
shl    $0x1,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    -0x38(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
shl    $0x1,%rcx
movsd  %xmm0,0x8(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    112c <window_cell+0x2c>
jmp    11f8 <window_cell+0xf8>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    11f3 <window_cell+0xf3>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    -0x38(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11b2 <window_cell+0xb2>
jmp    11f8 <window_cell+0xf8>
pop    %rbp
ret

x86 assembly (clang, -O2):
<window_cell>:
cmp    $0x2,%edi
jne    1157 <window_cell+0x57>
test   %r8,%r8
jle    11f7 <window_cell+0xf7>
add    $0x8,%rcx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movslq (%r9,%rax,4),%rdx
lea    (%rdx,%rdx,1),%rdi
shl    $0x4,%rdx
movsd  (%rsi,%rdx,1),%xmm0
movsd  %xmm0,-0x8(%rcx)
or     $0x1,%edi
movslq %edi,%rdx
movsd  (%rsi,%rdx,8),%xmm0
movsd  %xmm0,(%rcx)
add    $0x1,%rax
add    $0x10,%rcx
cmp    %rax,%r8
jne    1120 <window_cell+0x20>
jmp    11f7 <window_cell+0xf7>
test   %r8,%r8
jle    11f7 <window_cell+0xf7>
lea    -0x1(%r8),%rdx
mov    %r8d,%r10d
and    $0x3,%r10d
cmp    $0x3,%rdx
jae    1175 <window_cell+0x75>
xor    %edx,%edx
jmp    11c7 <window_cell+0xc7>
and    $0xfffffffffffffffc,%r8
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movslq (%r9,%rdx,4),%rdi
movsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,(%rcx,%rdx,8)
movslq 0x4(%r9,%rdx,4),%rdi
movsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,0x8(%rcx,%rdx,8)
movslq 0x8(%r9,%rdx,4),%rdi
movsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,0x10(%rcx,%rdx,8)
movslq 0xc(%r9,%rdx,4),%rdi
movsd  (%rsi,%rdi,8),%xmm0
movsd  %xmm0,0x18(%rcx,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%r8
jne    1180 <window_cell+0x80>
test   %r10,%r10
je     11f7 <window_cell+0xf7>
lea    (%rcx,%rdx,8),%rcx
lea    (%r9,%rdx,4),%rdx
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
movslq (%rdx,%rdi,4),%rax
movsd  (%rsi,%rax,8),%xmm0
movsd  %xmm0,(%rcx,%rdi,8)
add    $0x1,%rdi
cmp    %rdi,%r10
jne    11e0 <window_cell+0xe0>
ret

LLVM IR (-O0):
define void @window_cell(i32 %0, double* %1, i32 %2, double* %3, i64 %4, i32* %5) {
%7 = alloca i32
%8 = alloca double*
%9 = alloca i32
%10 = alloca double*
%11 = alloca i64
%12 = alloca i32*
%13 = alloca i64
store i32 %0, i32* %7
store double* %1, double** %8
store i32 %2, i32* %9
store double* %3, double** %10
store i64 %4, i64* %11
store i32* %5, i32** %12
%14 = load i32, i32* %7
%15 = icmp eq i32 %14, 2
br i1 %15, label %16, label %54
16: ; preds = %6
store i64 0, i64* %13
br label %17
17: ; preds = %50, %16
%18 = load i64, i64* %13
%19 = load i64, i64* %11
%20 = icmp slt i64 %18, %19
br i1 %20, label %21, label %53
21: ; preds = %17
%22 = load double*, double** %8
%23 = load i32*, i32** %12
%24 = load i64, i64* %13
%25 = getelementptr inbounds i32, i32* %23, i64 %24
%26 = load i32, i32* %25
%27 = mul nsw i32 %26, 2
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds double, double* %22, i64 %28
%30 = load double, double* %29
%31 = load double*, double** %10
%32 = load i64, i64* %13
%33 = mul nsw i64 %32, 2
%34 = getelementptr inbounds double, double* %31, i64 %33
store double %30, double* %34
%35 = load double*, double** %8
%36 = load i32*, i32** %12
%37 = load i64, i64* %13
%38 = getelementptr inbounds i32, i32* %36, i64 %37
%39 = load i32, i32* %38
%40 = mul nsw i32 %39, 2
%41 = add nsw i32 %40, 1
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds double, double* %35, i64 %42
%44 = load double, double* %43
%45 = load double*, double** %10
%46 = load i64, i64* %13
%47 = mul nsw i64 %46, 2
%48 = add nsw i64 %47, 1
%49 = getelementptr inbounds double, double* %45, i64 %48
store double %44, double* %49
br label %50
50: ; preds = %21
%51 = load i64, i64* %13
%52 = add nsw i64 %51, 1
store i64 %52, i64* %13
br label %17
53: ; preds = %17
br label %75
54: ; preds = %6
store i64 0, i64* %13
br label %55
55: ; preds = %71, %54
%56 = load i64, i64* %13
%57 = load i64, i64* %11
%58 = icmp slt i64 %56, %57
br i1 %58, label %59, label %74
59: ; preds = %55
%60 = load double*, double** %8
%61 = load i32*, i32** %12
%62 = load i64, i64* %13
%63 = getelementptr inbounds i32, i32* %61, i64 %62
%64 = load i32, i32* %63
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %60, i64 %65
%67 = load double, double* %66
%68 = load double*, double** %10
%69 = load i64, i64* %13
%70 = getelementptr inbounds double, double* %68, i64 %69
store double %67, double* %70
br label %71
71: ; preds = %59
%72 = load i64, i64* %13
%73 = add nsw i64 %72, 1
store i64 %73, i64* %13
br label %55
74: ; preds = %55
br label %75
75: ; preds = %74, %53
ret void
}

LLVM IR (-O2):
define void @window_cell(i32 %0, double* nocapture %1, i32 %2, double* nocapture writeonly %3, i64 %4, i32* nocapture %5) {
%7 = icmp eq i32 %0, 2
%8 = icmp sgt i64 %4, 0
br i1 %7, label %16, label %9
9: ; preds = %6
br i1 %8, label %10, label %83
10: ; preds = %9
%11 = add i64 %4, -1
%12 = and i64 %4, 3
%13 = icmp ult i64 %11, 3
br i1 %13, label %68, label %14
14: ; preds = %10
%15 = and i64 %4, -4
br label %35
16: ; preds = %6
br i1 %8, label %17, label %83
17: ; preds = %16, %17
%18 = phi i64 [ %33, %17 ], [ 0, %16 ]
%19 = getelementptr inbounds i32, i32* %5, i64 %18
%20 = load i32, i32* %19
%21 = shl nsw i32 %20, 1
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds double, double* %1, i64 %22
%24 = load double, double* %23
%25 = shl nuw nsw i64 %18, 1
%26 = getelementptr inbounds double, double* %3, i64 %25
store double %24, double* %26
%27 = or i32 %21, 1
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds double, double* %1, i64 %28
%30 = load double, double* %29
%31 = or i64 %25, 1
%32 = getelementptr inbounds double, double* %3, i64 %31
store double %30, double* %32
%33 = add nuw nsw i64 %18, 1
%34 = icmp eq i64 %33, %4
br i1 %34, label %83, label %17
35: ; preds = %35, %14
%36 = phi i64 [ 0, %14 ], [ %65, %35 ]
%37 = phi i64 [ 0, %14 ], [ %66, %35 ]
%38 = getelementptr inbounds i32, i32* %5, i64 %36
%39 = load i32, i32* %38
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds double, double* %1, i64 %40
%42 = load double, double* %41
%43 = getelementptr inbounds double, double* %3, i64 %36
store double %42, double* %43
%44 = or i64 %36, 1
%45 = getelementptr inbounds i32, i32* %5, i64 %44
%46 = load i32, i32* %45
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %1, i64 %47
%49 = load double, double* %48
%50 = getelementptr inbounds double, double* %3, i64 %44
store double %49, double* %50
%51 = or i64 %36, 2
%52 = getelementptr inbounds i32, i32* %5, i64 %51
%53 = load i32, i32* %52
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds double, double* %1, i64 %54
%56 = load double, double* %55
%57 = getelementptr inbounds double, double* %3, i64 %51
store double %56, double* %57
%58 = or i64 %36, 3
%59 = getelementptr inbounds i32, i32* %5, i64 %58
%60 = load i32, i32* %59
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double, double* %1, i64 %61
%63 = load double, double* %62
%64 = getelementptr inbounds double, double* %3, i64 %58
store double %63, double* %64
%65 = add nuw nsw i64 %36, 4
%66 = add i64 %37, 4
%67 = icmp eq i64 %66, %15
br i1 %67, label %68, label %35
68: ; preds = %35, %10
%69 = phi i64 [ 0, %10 ], [ %65, %35 ]
%70 = icmp eq i64 %12, 0
br i1 %70, label %83, label %71
71: ; preds = %68, %71
%72 = phi i64 [ %80, %71 ], [ %69, %68 ]
%73 = phi i64 [ %81, %71 ], [ 0, %68 ]
%74 = getelementptr inbounds i32, i32* %5, i64 %72
%75 = load i32, i32* %74
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %1, i64 %76
%78 = load double, double* %77
%79 = getelementptr inbounds double, double* %3, i64 %72
store double %78, double* %79
%80 = add nuw nsw i64 %72, 1
%81 = add i64 %73, 1
%82 = icmp eq i64 %81, %12
br i1 %82, label %83, label %71
83: ; preds = %68, %71, %17, %9, %16
ret void
}