void cholmult_sparse_neardiag(double **a, int N, int *Nvec, double *p,
                              double *b, double *x) {
  int i, k;
  double sum;
  for (i = N - 1; i >= 0; i--) {
    for (sum = b[i] * p[i], k = i - 1; k >= 0; k--) {
      if (Nvec[k] < i)
        break;
      sum += a[i][k] * b[k];
    }
    x[i] = sum;
  }
}


x86 assembly (clang, -O0):
<cholmult_sparse_neardiag>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x34(%rbp)
jl     11db <cholmult_sparse_neardiag+0xdb>
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x34(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x0,-0x38(%rbp)
jl     11bb <cholmult_sparse_neardiag+0xbb>
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    -0x34(%rbp),%eax
jge    1179 <cholmult_sparse_neardiag+0x79>
jmp    11bb <cholmult_sparse_neardiag+0xbb>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x38(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x38(%rbp)
jmp    1156 <cholmult_sparse_neardiag+0x56>
movsd  -0x40(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    1124 <cholmult_sparse_neardiag+0x24>
pop    %rbp
ret

x86 assembly (clang, -O2):
<cholmult_sparse_neardiag>:
push   %r14
push   %rbx
test   %esi,%esi
jle    1168 <cholmult_sparse_neardiag+0x68>
mov    %esi,%r10d
add    $0xffffffff,%esi
jmp    1122 <cholmult_sparse_neardiag+0x22>
nop
movsd  %xmm0,(%r9,%r14,8)
add    $0xffffffff,%esi
cmp    $0x1,%r10
mov    %r11,%r10
jbe    1168 <cholmult_sparse_neardiag+0x68>
lea    -0x1(%r10),%r11
mov    %r11d,%r14d
movsd  (%r8,%r14,8),%xmm0
mulsd  (%rcx,%r14,8),%xmm0
cmp    $0x1,%r10
jbe    1162 <cholmult_sparse_neardiag+0x62>
mov    %esi,%ebx
nopl   (%rax)
add    $0xffffffff,%ebx
cmp    %r11d,(%rdx,%rbx,4)
jl     1110 <cholmult_sparse_neardiag+0x10>
mov    (%rdi,%r14,8),%rax
movsd  (%rax,%rbx,8),%xmm1
mulsd  (%r8,%rbx,8),%xmm1
addsd  %xmm1,%xmm0
test   %ebx,%ebx
jg     1140 <cholmult_sparse_neardiag+0x40>
jmp    1110 <cholmult_sparse_neardiag+0x10>
movsd  %xmm0,(%r9,%r14,8)
pop    %rbx
pop    %r14
ret

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

LLVM IR (-O2):
define void @cholmult_sparse_neardiag(double** nocapture %0, i32 %1, i32* nocapture %2, double* nocapture %3, double* nocapture %4, double* nocapture writeonly %5) {
%7 = icmp sgt i32 %1, 0
br i1 %7, label %8, label %46
8: ; preds = %6
%9 = zext i32 %1 to i64
br label %10
10: ; preds = %8, %43
%11 = phi i64 [ %9, %8 ], [ %12, %43 ]
%12 = add nsw i64 %11, -1
%13 = trunc i64 %12 to i32
%14 = and i64 %12, 4294967295
%15 = getelementptr inbounds double, double* %4, i64 %14
%16 = load double, double* %15
%17 = getelementptr inbounds double, double* %3, i64 %14
%18 = load double, double* %17
%19 = fmul double %16, %18
%20 = icmp ugt i64 %11, 1
br i1 %20, label %23, label %21
21: ; preds = %10
%22 = getelementptr inbounds double, double* %5, i64 %14
store double %19, double* %22
br label %46
23: ; preds = %10
%24 = trunc i64 %11 to i32
%25 = add nsw i32 %24, -2
%26 = getelementptr inbounds double*, double** %0, i64 %14
br label %27
27: ; preds = %23, %34
%28 = phi double [ %19, %23 ], [ %40, %34 ]
%29 = phi i32 [ %25, %23 ], [ %41, %34 ]
%30 = zext i32 %29 to i64
%31 = getelementptr inbounds i32, i32* %2, i64 %30
%32 = load i32, i32* %31
%33 = icmp slt i32 %32, %13
br i1 %33, label %43, label %34
34: ; preds = %27
%35 = load double*, double** %26
%36 = getelementptr inbounds double, double* %35, i64 %30
%37 = load double, double* %36
%38 = getelementptr inbounds double, double* %4, i64 %30
%39 = load double, double* %38
%40 = tail call double @llvm.fmuladd.f64(double %37, double %39, double %28)
%41 = add nsw i32 %29, -1
%42 = icmp sgt i32 %29, 0
br i1 %42, label %27, label %43
43: ; preds = %34, %27
%44 = phi double [ %40, %34 ], [ %28, %27 ]
%45 = getelementptr inbounds double, double* %5, i64 %14
store double %44, double* %45
br i1 %20, label %10, label %46
46: ; preds = %43, %21, %6
ret void
}