void setXYind2add(int nnz, double dd, double *y, double *x, int *s, int *invp) {
  int i;
  for (i = 0; i < nnz; ++i) {
    x[i] += dd * y[invp[s[i]]];
    y[invp[s[i]]] = 0.0;
  }
}


x86 assembly (clang, -O0):
<setXYind2add>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1199 <setXYind2add+0x99>
movsd  -0x10(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    -0x28(%rbp),%rdx
movslq -0x34(%rbp),%rsi
movslq (%rdx,%rsi,4),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    -0x28(%rbp),%rdx
movslq -0x34(%rbp),%rsi
movslq (%rdx,%rsi,4),%rdx
movslq (%rcx,%rdx,4),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1123 <setXYind2add+0x23>
pop    %rbp
ret

x86 assembly (clang, -O2):
<setXYind2add>:
test   %edi,%edi
jle    1197 <setXYind2add+0x97>
mov    %edi,%r9d
cmp    $0x1,%edi
jne    1114 <setXYind2add+0x14>
xor    %edi,%edi
jmp    1172 <setXYind2add+0x72>
mov    %r9d,%r10d
and    $0xfffffffe,%r10d
xor    %edi,%edi
nopl   (%rax)
movslq (%rcx,%rdi,4),%rax
movslq (%r8,%rax,4),%rax
movsd  (%rsi,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdx,%rdi,8),%xmm1
movsd  %xmm1,(%rdx,%rdi,8)
movq   $0x0,(%rsi,%rax,8)
movslq 0x4(%rcx,%rdi,4),%rax
movslq (%r8,%rax,4),%rax
movsd  (%rsi,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x8(%rdx,%rdi,8),%xmm1
movsd  %xmm1,0x8(%rdx,%rdi,8)
movq   $0x0,(%rsi,%rax,8)
add    $0x2,%rdi
cmp    %rdi,%r10
jne    1120 <setXYind2add+0x20>
test   $0x1,%r9b
je     1197 <setXYind2add+0x97>
movslq (%rcx,%rdi,4),%rax
movslq (%r8,%rax,4),%rax
mulsd  (%rsi,%rax,8),%xmm0
addsd  (%rdx,%rdi,8),%xmm0
movsd  %xmm0,(%rdx,%rdi,8)
movq   $0x0,(%rsi,%rax,8)
ret

LLVM IR (-O0):
define void @setXYind2add(i32 %0, double %1, double* %2, double* %3, i32* %4, i32* %5) {
%7 = alloca i32
%8 = alloca double
%9 = alloca double*
%10 = alloca double*
%11 = alloca i32*
%12 = alloca i32*
%13 = alloca i32
store i32 %0, i32* %7
store double %1, double* %8
store double* %2, double** %9
store double* %3, double** %10
store i32* %4, i32** %11
store i32* %5, i32** %12
store i32 0, i32* %13
br label %14
14: ; preds = %51, %6
%15 = load i32, i32* %13
%16 = load i32, i32* %7
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %54
18: ; preds = %14
%19 = load double, double* %8
%20 = load double*, double** %9
%21 = load i32*, i32** %12
%22 = load i32*, i32** %11
%23 = load i32, i32* %13
%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 i32, i32* %21, i64 %27
%29 = load i32, i32* %28
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %20, i64 %30
%32 = load double, double* %31
%33 = load double*, double** %10
%34 = load i32, i32* %13
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double, double* %33, i64 %35
%37 = load double, double* %36
%38 = call double @llvm.fmuladd.f64(double %19, double %32, double %37)
store double %38, double* %36
%39 = load double*, double** %9
%40 = load i32*, i32** %12
%41 = load i32*, i32** %11
%42 = load i32, i32* %13
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds i32, i32* %41, i64 %43
%45 = load i32, i32* %44
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds i32, i32* %40, i64 %46
%48 = load i32, i32* %47
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %39, i64 %49
store double 0.000000e+00, double* %50
br label %51
51: ; preds = %18
%52 = load i32, i32* %13
%53 = add nsw i32 %52, 1
store i32 %53, i32* %13
br label %14
54: ; preds = %14
ret void
}

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