double _sum_row_sparse(double *X, int n, int i, int *J, int from, int to) {
  int j;
  double sum = 0.0;
  for (j = from; j < to; j++)
    sum += X[i * n + J[j]];
  return sum;
}


x86 assembly (clang, -O0):
<_sum_row_sparse>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9d,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    116a <_sum_row_sparse+0x6a>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0xc(%rbp),%ecx
mov    -0x18(%rbp),%rdx
movslq -0x24(%rbp),%rsi
add    (%rdx,%rsi,4),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1128 <_sum_row_sparse+0x28>
movsd  -0x30(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<_sum_row_sparse>:
xorpd  %xmm0,%xmm0
cmp    %r9d,%r8d
jge    118c <_sum_row_sparse+0x8c>
imul   %esi,%edx
movslq %r8d,%rax
movslq %r9d,%r10
sub    %r8d,%r9d
mov    %rax,%r8
not    %r8
add    %r10,%r8
movslq %edx,%rdx
xorpd  %xmm0,%xmm0
and    $0x3,%r9
je     1146 <_sum_row_sparse+0x46>
nop
movslq (%rcx,%rax,4),%rsi
add    %rdx,%rsi
addsd  (%rdi,%rsi,8),%xmm0
add    $0x1,%rax
add    $0xffffffffffffffff,%r9
jne    1130 <_sum_row_sparse+0x30>
cmp    $0x3,%r8
jb     118c <_sum_row_sparse+0x8c>
nopl   0x0(%rax)
movslq (%rcx,%rax,4),%rsi
add    %rdx,%rsi
addsd  (%rdi,%rsi,8),%xmm0
movslq 0x4(%rcx,%rax,4),%rsi
add    %rdx,%rsi
addsd  (%rdi,%rsi,8),%xmm0
movslq 0x8(%rcx,%rax,4),%rsi
add    %rdx,%rsi
addsd  (%rdi,%rsi,8),%xmm0
movslq 0xc(%rcx,%rax,4),%rsi
add    %rdx,%rsi
addsd  (%rdi,%rsi,8),%xmm0
add    $0x4,%rax
cmp    %rax,%r10
jne    1150 <_sum_row_sparse+0x50>
ret

LLVM IR (-O0):
define double @_sum_row_sparse(double* %0, i32 %1, i32 %2, i32* %3, i32 %4, i32 %5) {
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca double
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 i32 %5, i32* %12
store double 0.000000e+00, double* %14
%15 = load i32, i32* %11
store i32 %15, i32* %13
br label %16
16: ; preds = %36, %6
%17 = load i32, i32* %13
%18 = load i32, i32* %12
%19 = icmp slt i32 %17, %18
br i1 %19, label %20, label %39
20: ; preds = %16
%21 = load double*, double** %7
%22 = load i32, i32* %9
%23 = load i32, i32* %8
%24 = mul nsw i32 %22, %23
%25 = load i32*, i32** %10
%26 = load i32, i32* %13
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds i32, i32* %25, i64 %27
%29 = load i32, i32* %28
%30 = add nsw i32 %24, %29
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds double, double* %21, i64 %31
%33 = load double, double* %32
%34 = load double, double* %14
%35 = fadd double %34, %33
store double %35, double* %14
br label %36
36: ; preds = %20
%37 = load i32, i32* %13
%38 = add nsw i32 %37, 1
store i32 %38, i32* %13
br label %16
39: ; preds = %16
%40 = load double, double* %14
ret double %40
}

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