void _normalize_all_sparse(double *X, int *I, int *J, int n, int n_idx) {
  int k;
  double sum = 0.0;
  for (k = 0; k < n_idx; k++)
    sum += X[I[k] * n + J[k]];
  for (k = 0; k < n_idx; k++)
    X[I[k] * n + J[k]] /= sum;
}


x86 assembly (clang, -O0):
<_normalize_all_sparse>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1170 <_normalize_all_sparse+0x70>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x24(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
imul   -0x1c(%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    1126 <_normalize_all_sparse+0x26>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11c5 <_normalize_all_sparse+0xc5>
movsd  -0x30(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x24(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
imul   -0x1c(%rbp),%ecx
mov    -0x18(%rbp),%rdx
movslq -0x24(%rbp),%rsi
add    (%rdx,%rsi,4),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1177 <_normalize_all_sparse+0x77>
pop    %rbp
ret

x86 assembly (clang, -O2):
<_normalize_all_sparse>:
push   %r14
push   %rbx
test   %r8d,%r8d
jle    1208 <_normalize_all_sparse+0x108>
mov    %r8d,%r9d
mov    %r9,%r10
movslq %ecx,%r14
add    $0xffffffffffffffff,%r10
je     120c <_normalize_all_sparse+0x10c>
mov    %r9d,%r11d
and    $0xfffffffe,%r11d
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
nopl   0x0(%rax)
movslq (%rsi,%rcx,4),%rbx
imul   %r14,%rbx
movslq (%rdx,%rcx,4),%rax
add    %rbx,%rax
addsd  (%rdi,%rax,8),%xmm0
movslq 0x4(%rsi,%rcx,4),%rax
imul   %r14,%rax
movslq 0x4(%rdx,%rcx,4),%rbx
add    %rax,%rbx
addsd  (%rdi,%rbx,8),%xmm0
add    $0x2,%rcx
cmp    %rcx,%r11
jne    1130 <_normalize_all_sparse+0x30>
test   $0x1,%r9b
je     117d <_normalize_all_sparse+0x7d>
movslq (%rsi,%rcx,4),%r11
imul   %r14,%r11
movslq (%rdx,%rcx,4),%rcx
add    %r11,%rcx
addsd  (%rdi,%rcx,8),%xmm0
test   %r8d,%r8d
jle    1208 <_normalize_all_sparse+0x108>
test   %r10,%r10
je     1221 <_normalize_all_sparse+0x121>
mov    %r9d,%r8d
and    $0xfffffffe,%r8d
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movslq (%rsi,%rcx,4),%rax
imul   %r14,%rax
movslq (%rdx,%rcx,4),%rbx
add    %rax,%rbx
movsd  (%rdi,%rbx,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi,%rbx,8)
movslq 0x4(%rsi,%rcx,4),%rax
imul   %r14,%rax
movslq 0x4(%rdx,%rcx,4),%rbx
add    %rax,%rbx
movsd  (%rdi,%rbx,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi,%rbx,8)
add    $0x2,%rcx
cmp    %rcx,%r8
jne    11a0 <_normalize_all_sparse+0xa0>
test   $0x1,%r9b
je     1208 <_normalize_all_sparse+0x108>
movslq (%rsi,%rcx,4),%rax
imul   %r14,%rax
movslq (%rdx,%rcx,4),%rcx
add    %rax,%rcx
movsd  (%rdi,%rcx,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi,%rcx,8)
pop    %rbx
pop    %r14
ret
xor    %ecx,%ecx
xorpd  %xmm0,%xmm0
test   $0x1,%r9b
jne    1169 <_normalize_all_sparse+0x69>
jmp    117d <_normalize_all_sparse+0x7d>
xor    %ecx,%ecx
test   $0x1,%r9b
jne    11eb <_normalize_all_sparse+0xeb>
jmp    1208 <_normalize_all_sparse+0x108>

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

LLVM IR (-O2):
define void @_normalize_all_sparse(double* nocapture %0, i32* nocapture %1, i32* nocapture %2, i32 %3, i32 %4) {
%6 = icmp sgt i32 %4, 0
br i1 %6, label %7, label %106
7: ; preds = %5
%8 = zext i32 %4 to i64
%9 = add nsw i64 %8, -1
%10 = and i64 %8, 1
%11 = icmp eq i64 %9, 0
br i1 %11, label %14, label %12
12: ; preds = %7
%13 = and i64 %8, 4294967294
br label %37
14: ; preds = %37, %7
%15 = phi double [ undef, %7 ], [ %61, %37 ]
%16 = phi i64 [ 0, %7 ], [ %62, %37 ]
%17 = phi double [ 0.000000e+00, %7 ], [ %61, %37 ]
%18 = icmp eq i64 %10, 0
br i1 %18, label %30, label %19
19: ; preds = %14
%20 = getelementptr inbounds i32, i32* %1, i64 %16
%21 = load i32, i32* %20
%22 = mul nsw i32 %21, %3
%23 = getelementptr inbounds i32, i32* %2, i64 %16
%24 = load i32, i32* %23
%25 = add nsw i32 %22, %24
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds double, double* %0, i64 %26
%28 = load double, double* %27
%29 = fadd double %17, %28
br label %30
30: ; preds = %14, %19
%31 = phi double [ %15, %14 ], [ %29, %19 ]
br i1 %6, label %32, label %106
32: ; preds = %30
%33 = and i64 %8, 1
%34 = icmp eq i64 %9, 0
br i1 %34, label %92, label %35
35: ; preds = %32
%36 = and i64 %8, 4294967294
br label %65
37: ; preds = %37, %12
%38 = phi i64 [ 0, %12 ], [ %62, %37 ]
%39 = phi double [ 0.000000e+00, %12 ], [ %61, %37 ]
%40 = phi i64 [ 0, %12 ], [ %63, %37 ]
%41 = getelementptr inbounds i32, i32* %1, i64 %38
%42 = load i32, i32* %41
%43 = mul nsw i32 %42, %3
%44 = getelementptr inbounds i32, i32* %2, i64 %38
%45 = load i32, i32* %44
%46 = add nsw i32 %43, %45
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %0, i64 %47
%49 = load double, double* %48
%50 = fadd double %39, %49
%51 = or i64 %38, 1
%52 = getelementptr inbounds i32, i32* %1, i64 %51
%53 = load i32, i32* %52
%54 = mul nsw i32 %53, %3
%55 = getelementptr inbounds i32, i32* %2, i64 %51
%56 = load i32, i32* %55
%57 = add nsw i32 %54, %56
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %0, i64 %58
%60 = load double, double* %59
%61 = fadd double %50, %60
%62 = add nuw nsw i64 %38, 2
%63 = add i64 %40, 2
%64 = icmp eq i64 %63, %13
br i1 %64, label %14, label %37
65: ; preds = %65, %35
%66 = phi i64 [ 0, %35 ], [ %89, %65 ]
%67 = phi i64 [ 0, %35 ], [ %90, %65 ]
%68 = getelementptr inbounds i32, i32* %1, i64 %66
%69 = load i32, i32* %68
%70 = mul nsw i32 %69, %3
%71 = getelementptr inbounds i32, i32* %2, i64 %66
%72 = load i32, i32* %71
%73 = add nsw i32 %70, %72
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds double, double* %0, i64 %74
%76 = load double, double* %75
%77 = fdiv double %76, %31
store double %77, double* %75
%78 = or i64 %66, 1
%79 = getelementptr inbounds i32, i32* %1, i64 %78
%80 = load i32, i32* %79
%81 = mul nsw i32 %80, %3
%82 = getelementptr inbounds i32, i32* %2, i64 %78
%83 = load i32, i32* %82
%84 = add nsw i32 %81, %83
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %0, i64 %85
%87 = load double, double* %86
%88 = fdiv double %87, %31
store double %88, double* %86
%89 = add nuw nsw i64 %66, 2
%90 = add i64 %67, 2
%91 = icmp eq i64 %90, %36
br i1 %91, label %92, label %65
92: ; preds = %65, %32
%93 = phi i64 [ 0, %32 ], [ %89, %65 ]
%94 = icmp eq i64 %33, 0
br i1 %94, label %106, label %95
95: ; preds = %92
%96 = getelementptr inbounds i32, i32* %1, i64 %93
%97 = load i32, i32* %96
%98 = mul nsw i32 %97, %3
%99 = getelementptr inbounds i32, i32* %2, i64 %93
%100 = load i32, i32* %99
%101 = add nsw i32 %98, %100
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds double, double* %0, i64 %102
%104 = load double, double* %103
%105 = fdiv double %104, %31
store double %105, double* %103
br label %106
106: ; preds = %95, %92, %5, %30
ret void
}