void pp2rr_sorted(const int *ii, const float *fdata, const int n,
                  const int stride, float *rr) {
  for (int j = 0; j < n; ++j) {
    const int i = ii[j];
    const float *r = fdata + j * stride;
    for (int c = 0; c < 3; ++c)
      rr[3 * i + c] = r[c];
  }
}


x86 assembly (clang, -O0):
<pp2rr_sorted>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11a2 <pp2rr_sorted+0xa2>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
cmpl   $0x3,-0x34(%rbp)
jge    118f <pp2rr_sorted+0x8f>
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
imul   $0x3,-0x28(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1157 <pp2rr_sorted+0x57>
jmp    1194 <pp2rr_sorted+0x94>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111d <pp2rr_sorted+0x1d>
pop    %rbp
ret

x86 assembly (clang, -O2):
<pp2rr_sorted>:
test   %edx,%edx
jle    1156 <pp2rr_sorted+0x56>
movslq %ecx,%rax
mov    %edx,%r9d
add    $0x8,%rsi
shl    $0x2,%rax
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movslq (%rdi,%rdx,4),%rcx
movss  -0x8(%rsi),%xmm0
lea    (%rcx,%rcx,2),%rcx
movss  %xmm0,(%r8,%rcx,4)
movss  -0x4(%rsi),%xmm0
movss  %xmm0,0x4(%r8,%rcx,4)
movss  (%rsi),%xmm0
movss  %xmm0,0x8(%r8,%rcx,4)
add    $0x1,%rdx
add    %rax,%rsi
cmp    %rdx,%r9
jne    1120 <pp2rr_sorted+0x20>
ret

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

LLVM IR (-O2):
define void @pp2rr_sorted(i32* nocapture %0, float* nocapture %1, i32 %2, i32 %3, float* nocapture writeonly %4) {
%6 = icmp sgt i32 %2, 0
br i1 %6, label %7, label %10
7: ; preds = %5
%8 = sext i32 %3 to i64
%9 = zext i32 %2 to i64
br label %11
10: ; preds = %11, %5
ret void
11: ; preds = %7, %11
%12 = phi i64 [ 0, %7 ], [ %29, %11 ]
%13 = getelementptr inbounds i32, i32* %0, i64 %12
%14 = load i32, i32* %13
%15 = mul nsw i64 %12, %8
%16 = getelementptr inbounds float, float* %1, i64 %15
%17 = mul nsw i32 %14, 3
%18 = sext i32 %17 to i64
%19 = load float, float* %16
%20 = getelementptr inbounds float, float* %4, i64 %18
store float %19, float* %20
%21 = getelementptr inbounds float, float* %16, i64 1
%22 = load float, float* %21
%23 = add nsw i64 %18, 1
%24 = getelementptr inbounds float, float* %4, i64 %23
store float %22, float* %24
%25 = getelementptr inbounds float, float* %16, i64 2
%26 = load float, float* %25
%27 = add nsw i64 %18, 2
%28 = getelementptr inbounds float, float* %4, i64 %27
store float %26, float* %28
%29 = add nuw nsw i64 %12, 1
%30 = icmp eq i64 %29, %9
br i1 %30, label %10, label %11
}