double dw_ddoti(const int nz, const double *x, const int *ind_x,
                const double *y) {
  double ret = 0.0;
  int i;
  for (i = 0; i < nz; i++) {
    ret += x[i] * y[ind_x[i]];
  }
  return ret;
}


x86 assembly (clang, -O0):
<dw_ddoti>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1170 <dw_ddoti+0x70>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x2c(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1122 <dw_ddoti+0x22>
movsd  -0x28(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<dw_ddoti>:
test   %edi,%edi
jle    1120 <dw_ddoti+0x20>
mov    %edi,%r9d
lea    -0x1(%r9),%rdi
mov    %r9d,%r8d
and    $0x3,%r8d
cmp    $0x3,%rdi
jae    1124 <dw_ddoti+0x24>
xorpd  %xmm0,%xmm0
xor    %edi,%edi
jmp    1187 <dw_ddoti+0x87>
xorps  %xmm0,%xmm0
ret
and    $0xfffffffc,%r9d
xorpd  %xmm0,%xmm0
xor    %edi,%edi
xchg   %ax,%ax
movslq (%rdx,%rdi,4),%rax
movsd  (%rsi,%rdi,8),%xmm1
movsd  0x8(%rsi,%rdi,8),%xmm2
mulsd  (%rcx,%rax,8),%xmm1
addsd  %xmm0,%xmm1
movslq 0x4(%rdx,%rdi,4),%rax
mulsd  (%rcx,%rax,8),%xmm2
addsd  %xmm1,%xmm2
movsd  0x10(%rsi,%rdi,8),%xmm1
movslq 0x8(%rdx,%rdi,4),%rax
mulsd  (%rcx,%rax,8),%xmm1
movsd  0x18(%rsi,%rdi,8),%xmm0
movslq 0xc(%rdx,%rdi,4),%rax
mulsd  (%rcx,%rax,8),%xmm0
addsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
add    $0x4,%rdi
cmp    %rdi,%r9
jne    1130 <dw_ddoti+0x30>
test   %r8,%r8
je     11bb <dw_ddoti+0xbb>
lea    (%rdx,%rdi,4),%rax
lea    (%rsi,%rdi,8),%rdx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
movsd  (%rdx,%rsi,8),%xmm1
movslq (%rax,%rsi,4),%rdi
mulsd  (%rcx,%rdi,8),%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rsi
cmp    %rsi,%r8
jne    11a0 <dw_ddoti+0xa0>
ret

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

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