float spfvec_inner_product(int *idx1, float *val1, int nz1, int *idx2,
                           float *val2, int nz2) {
  double s = 0;
  long i1 = 0, i2 = 0;
  while (1) {
    if (i1 == nz1)
      break;
    if (i2 == nz2)
      break;
    if (idx1[i1] == idx2[i2]) {
      s += val1[i1] * val2[i2];
      i1++;
      i2++;
    } else {
      if (idx1[i1] < idx2[i2]) {
        i1++;
      } else {
        i2++;
      }
    }
  }
  return (float)s;
}


x86 assembly (clang, -O0):
<spfvec_inner_product>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movq   $0x0,-0x40(%rbp)
movq   $0x0,-0x48(%rbp)
mov    -0x40(%rbp),%rax
movslq -0x14(%rbp),%rcx
cmp    %rcx,%rax
jne    1149 <spfvec_inner_product+0x49>
jmp    1203 <spfvec_inner_product+0x103>
mov    -0x48(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmp    %rcx,%rax
jne    115f <spfvec_inner_product+0x5f>
jmp    1203 <spfvec_inner_product+0x103>
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x20(%rbp),%rcx
mov    -0x48(%rbp),%rdx
cmp    (%rcx,%rdx,4),%eax
jne    11c0 <spfvec_inner_product+0xc0>
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x48(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    11fe <spfvec_inner_product+0xfe>
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x20(%rbp),%rcx
mov    -0x48(%rbp),%rdx
cmp    (%rcx,%rdx,4),%eax
jge    11ed <spfvec_inner_product+0xed>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    11f9 <spfvec_inner_product+0xf9>
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    11fe <spfvec_inner_product+0xfe>
jmp    1133 <spfvec_inner_product+0x33>
movsd  -0x38(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<spfvec_inner_product>:
xorpd  %xmm0,%xmm0
test   %edx,%edx
je     116e <spfvec_inner_product+0x6e>
test   %r9d,%r9d
je     116e <spfvec_inner_product+0x6e>
movslq %edx,%r10
movslq %r9d,%r9
xorpd  %xmm0,%xmm0
xor    %r11d,%r11d
xor    %eax,%eax
nopl   0x0(%rax)
mov    (%rcx,%r11,4),%edx
cmp    %edx,(%rdi,%rax,4)
jne    1150 <spfvec_inner_product+0x50>
movss  (%rsi,%rax,4),%xmm1
mulss  (%r8,%r11,4),%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rax
jmp    1160 <spfvec_inner_product+0x60>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
jge    1160 <spfvec_inner_product+0x60>
add    $0x1,%rax
cmp    %r10,%rax
jne    1169 <spfvec_inner_product+0x69>
jmp    116e <spfvec_inner_product+0x6e>
nopl   (%rax)
add    $0x1,%r11
cmp    %r10,%rax
je     116e <spfvec_inner_product+0x6e>
cmp    %r9,%r11
jne    1120 <spfvec_inner_product+0x20>
cvtsd2ss %xmm0,%xmm0
ret

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

LLVM IR (-O2):
define float @spfvec_inner_product(i32* nocapture %0, float* nocapture %1, i32 %2, i32* nocapture %3, float* nocapture %4, i32 %5) {
%7 = sext i32 %2 to i64
%8 = sext i32 %5 to i64
%9 = icmp eq i32 %2, 0
%10 = icmp eq i32 %5, 0
%11 = or i1 %9, %10
br i1 %11, label %44, label %12
12: ; preds = %6, %37
%13 = phi i64 [ %40, %37 ], [ 0, %6 ]
%14 = phi i64 [ %39, %37 ], [ 0, %6 ]
%15 = phi double [ %38, %37 ], [ 0.000000e+00, %6 ]
%16 = getelementptr inbounds i32, i32* %0, i64 %14
%17 = load i32, i32* %16
%18 = getelementptr inbounds i32, i32* %3, i64 %13
%19 = load i32, i32* %18
%20 = icmp eq i32 %17, %19
br i1 %20, label %21, label %31
21: ; preds = %12
%22 = getelementptr inbounds float, float* %1, i64 %14
%23 = load float, float* %22
%24 = getelementptr inbounds float, float* %4, i64 %13
%25 = load float, float* %24
%26 = fmul float %23, %25
%27 = fpext float %26 to double
%28 = fadd double %15, %27
%29 = add nsw i64 %14, 1
%30 = add nsw i64 %13, 1
br label %37
31: ; preds = %12
%32 = icmp slt i32 %17, %19
br i1 %32, label %33, label %35
33: ; preds = %31
%34 = add nsw i64 %14, 1
br label %37
35: ; preds = %31
%36 = add nsw i64 %13, 1
br label %37
37: ; preds = %33, %35, %21
%38 = phi double [ %28, %21 ], [ %15, %33 ], [ %15, %35 ]
%39 = phi i64 [ %29, %21 ], [ %34, %33 ], [ %14, %35 ]
%40 = phi i64 [ %30, %21 ], [ %13, %33 ], [ %36, %35 ]
%41 = icmp eq i64 %39, %7
%42 = icmp eq i64 %40, %8
%43 = select i1 %41, i1 true, i1 %42
br i1 %43, label %44, label %12
44: ; preds = %37, %6
%45 = phi double [ 0.000000e+00, %6 ], [ %38, %37 ]
%46 = fptrunc double %45 to float
ret float %46
}