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


x86 assembly (clang, -O0):
<window_cell>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x18(%rbp),%rax
jge    115d <window_cell+0x5d>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    -0x28(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    111c <window_cell+0x1c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<window_cell>:
test   %rdx,%rdx
jle    11a7 <window_cell+0xa7>
lea    -0x1(%rdx),%rax
mov    %edx,%r8d
and    $0x3,%r8d
cmp    $0x3,%rax
jae    111f <window_cell+0x1f>
xor    %r9d,%r9d
jmp    117b <window_cell+0x7b>
and    $0xfffffffffffffffc,%rdx
xor    %r9d,%r9d
cs nopw 0x0(%rax,%rax,1)
movslq (%rcx,%r9,4),%rax
movsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,(%rsi,%r9,8)
movslq 0x4(%rcx,%r9,4),%rax
movsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,0x8(%rsi,%r9,8)
movslq 0x8(%rcx,%r9,4),%rax
movsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,0x10(%rsi,%r9,8)
movslq 0xc(%rcx,%r9,4),%rax
movsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,0x18(%rsi,%r9,8)
add    $0x4,%r9
cmp    %r9,%rdx
jne    1130 <window_cell+0x30>
test   %r8,%r8
je     11a7 <window_cell+0xa7>
lea    (%rsi,%r9,8),%rdx
lea    (%rcx,%r9,4),%rax
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movslq (%rax,%rcx,4),%rsi
movsd  (%rdi,%rsi,8),%xmm0
movsd  %xmm0,(%rdx,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%r8
jne    1190 <window_cell+0x90>
ret

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

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