void sp_restrict_sorted(double *y, unsigned int Rn, const unsigned int *Ri,
                        unsigned int xn, const unsigned int *xi,
                        const double *x) {
  const unsigned int *xe = xi + xn;
  double *ye = y + Rn;
  unsigned int iy;
  if (y == ye)
    return;
  iy = *Ri;
  for (; xi != xe; ++xi, ++x) {
    unsigned int ix = *xi;
    while (iy < ix) {
      *y++ = 0;
      if (y == ye)
        return;
      iy = *(++Ri);
    }
    if (iy == ix) {
      *y++ = *x;
      if (y == ye)
        return;
      iy = *(++Ri);
    }
  }
  while (y != ye)
    *y++ = 0;
}


x86 assembly (clang, -O0):
<sp_restrict_sorted>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x1c(%rbp),%ecx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x8(%rbp),%rax
cmp    -0x40(%rbp),%rax
jne    1151 <sp_restrict_sorted+0x51>
jmp    125f <sp_restrict_sorted+0x15f>
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x44(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x38(%rbp),%rax
je     1231 <sp_restrict_sorted+0x131>
mov    -0x28(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x48(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x48(%rbp),%eax
jae    11c0 <sp_restrict_sorted+0xc0>
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x8(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
cmp    -0x40(%rbp),%rax
jne    11a6 <sp_restrict_sorted+0xa6>
jmp    125f <sp_restrict_sorted+0x15f>
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
mov    0x4(%rax),%eax
mov    %eax,-0x44(%rbp)
jmp    1171 <sp_restrict_sorted+0x71>
mov    -0x44(%rbp),%eax
cmp    -0x48(%rbp),%eax
jne    120f <sp_restrict_sorted+0x10f>
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x8(%rbp)
movsd  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
cmp    -0x40(%rbp),%rax
jne    11fa <sp_restrict_sorted+0xfa>
jmp    125f <sp_restrict_sorted+0x15f>
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
mov    0x4(%rax),%eax
mov    %eax,-0x44(%rbp)
jmp    1214 <sp_restrict_sorted+0x114>
mov    -0x28(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x28(%rbp)
mov    -0x30(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x30(%rbp)
jmp    115a <sp_restrict_sorted+0x5a>
jmp    1236 <sp_restrict_sorted+0x136>
mov    -0x8(%rbp),%rax
cmp    -0x40(%rbp),%rax
je     125f <sp_restrict_sorted+0x15f>
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x8(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
jmp    1236 <sp_restrict_sorted+0x136>
pop    %rbp
ret

x86 assembly (clang, -O2):
<sp_restrict_sorted>:
push   %rax
test   %esi,%esi
je     119d <sp_restrict_sorted+0x8d>
mov    %rdx,%rax
mov    %esi,%edx
lea    (%rdi,%rdx,8),%rdx
test   %ecx,%ecx
je     118a <sp_restrict_sorted+0x7a>
mov    %ecx,%ecx
lea    (%r8,%rcx,4),%r10
mov    (%rax),%esi
jmp    113d <sp_restrict_sorted+0x2d>
add    $0x4,%r8
add    $0x8,%r9
cmp    %r10,%r8
je     118a <sp_restrict_sorted+0x7a>
mov    (%r8),%ecx
cmp    %ecx,%esi
jae    116b <sp_restrict_sorted+0x5b>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movq   $0x0,(%rdi)
add    $0x8,%rdi
cmp    %rdx,%rdi
je     119d <sp_restrict_sorted+0x8d>
mov    0x4(%rax),%esi
add    $0x4,%rax
cmp    %ecx,%esi
jb     1150 <sp_restrict_sorted+0x40>
cmp    %ecx,%esi
jne    1130 <sp_restrict_sorted+0x20>
movsd  (%r9),%xmm0
movsd  %xmm0,(%rdi)
add    $0x8,%rdi
cmp    %rdx,%rdi
je     119d <sp_restrict_sorted+0x8d>
mov    0x4(%rax),%esi
add    $0x4,%rax
jmp    1130 <sp_restrict_sorted+0x20>
cmp    %rdx,%rdi
je     119d <sp_restrict_sorted+0x8d>
sub    %rdi,%rdx
and    $0xfffffffffffffff8,%rdx
xor    %esi,%esi
call   1030 <memset@plt>
pop    %rax
ret

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

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