void getDInvValues(double *D, double *values, unsigned long M, unsigned long N,
                   unsigned long long nz) {
  for (int i = 0; i < M; i++) {
    for (int j = 0; j < N; j++) {
      if (i == j) {
        D[i * N + j] = 1.0 / values[i * N + j];
      } else {
        D[i * N + j] = 0.0;
      }
    }
  }
}


x86 assembly (clang, -O0):
<getDInvValues>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
movslq -0x2c(%rbp),%rax
cmp    -0x18(%rbp),%rax
jae    11cf <getDInvValues+0xcf>
movl   $0x0,-0x30(%rbp)
movslq -0x30(%rbp),%rax
cmp    -0x20(%rbp),%rax
jae    11bc <getDInvValues+0xbc>
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jne    118d <getDInvValues+0x8d>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   -0x20(%rbp),%rcx
movslq -0x30(%rbp),%rdx
add    %rdx,%rcx
movsd  0xe96(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   -0x20(%rbp),%rcx
movslq -0x30(%rbp),%rdx
add    %rdx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    11a9 <getDInvValues+0xa9>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   -0x20(%rbp),%rcx
movslq -0x30(%rbp),%rdx
add    %rdx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    11ae <getDInvValues+0xae>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1134 <getDInvValues+0x34>
jmp    11c1 <getDInvValues+0xc1>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    111f <getDInvValues+0x1f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<getDInvValues>:
test   %rdx,%rdx
je     1168 <getDInvValues+0x68>
lea    0x0(,%rcx,8),%r8
xor    %r10d,%r10d
movsd  0xee8(%rip),%xmm0
mov    %rdi,%r9
jmp    112c <getDInvValues+0x2c>
nopl   (%rax)
add    $0x1,%r10
add    %r8,%r9
cmp    %rdx,%r10
je     1168 <getDInvValues+0x68>
test   %rcx,%rcx
je     1120 <getDInvValues+0x20>
mov    %r10,%r11
imul   %rcx,%r11
add    %r10,%r11
xor    %eax,%eax
jmp    1151 <getDInvValues+0x51>
nop
movq   $0x0,(%r9,%rax,8)
add    $0x1,%rax
cmp    %rax,%rcx
je     1120 <getDInvValues+0x20>
cmp    %rax,%r10
jne    1140 <getDInvValues+0x40>
movapd %xmm0,%xmm1
divsd  (%rsi,%r11,8),%xmm1
movsd  %xmm1,(%rdi,%r11,8)
jmp    1148 <getDInvValues+0x48>
ret

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

LLVM IR (-O2):
define void @getDInvValues(double* nocapture writeonly %0, double* nocapture %1, i64 %2, i64 %3, i64 %4) {
%6 = icmp eq i64 %2, 0
br i1 %6, label %16, label %7
7: ; preds = %5
%8 = icmp eq i64 %3, 0
br label %9
9: ; preds = %7, %17
%10 = phi i64 [ 0, %7 ], [ %18, %17 ]
br i1 %8, label %17, label %11
11: ; preds = %9
%12 = mul i64 %10, %3
%13 = add i64 %10, %12
%14 = getelementptr inbounds double, double* %1, i64 %13
%15 = getelementptr inbounds double, double* %0, i64 %13
br label %20
16: ; preds = %17, %5
ret void
17: ; preds = %29, %9
%18 = add nuw nsw i64 %10, 1
%19 = icmp eq i64 %18, %2
br i1 %19, label %16, label %9
20: ; preds = %11, %29
%21 = phi i64 [ 0, %11 ], [ %30, %29 ]
%22 = icmp eq i64 %10, %21
br i1 %22, label %23, label %26
23: ; preds = %20
%24 = load double, double* %14
%25 = fdiv double 1.000000e+00, %24
store double %25, double* %15
br label %29
26: ; preds = %20
%27 = add i64 %21, %12
%28 = getelementptr inbounds double, double* %0, i64 %27
store double 0.000000e+00, double* %28
br label %29
29: ; preds = %23, %26
%30 = add nuw nsw i64 %21, 1
%31 = icmp eq i64 %30, %3
br i1 %31, label %17, label %20
}