void make_zero_rows(double **q, const int dim) {
  for (int i = 0; i < dim; i++) {
    double sum = 0.0;
    for (int j = 0; j < dim; j++) {
      if (i != j) {
        sum += q[i][j];
      }
    }
    q[i][i] = -sum;
  }
}


x86 assembly (clang, -O0):
<make_zero_rows>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11b6 <make_zero_rows+0xb6>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1177 <make_zero_rows+0x77>
mov    -0x10(%rbp),%eax
cmp    -0x1c(%rbp),%eax
je     1164 <make_zero_rows+0x64>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
jmp    1169 <make_zero_rows+0x69>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    112d <make_zero_rows+0x2d>
movsd  -0x18(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1112 <make_zero_rows+0x12>
pop    %rbp
ret

x86 assembly (clang, -O2):
<make_zero_rows>:
test   %esi,%esi
jle    11ad <make_zero_rows+0xad>
mov    %esi,%r8d
mov    %r8d,%r9d
and    $0xfffffffe,%r9d
xor    %r10d,%r10d
movapd 0xee3(%rip),%xmm0
jmp    1137 <make_zero_rows+0x37>
nop
xorpd  %xmm0,%xmm1
mov    (%rdi,%r10,8),%rax
movlpd %xmm1,(%rax,%r10,8)
add    $0x1,%r10
cmp    %r8,%r10
je     11ad <make_zero_rows+0xad>
xorpd  %xmm1,%xmm1
cmp    $0x1,%esi
jne    1150 <make_zero_rows+0x50>
xor    %ecx,%ecx
test   $0x1,%r8b
je     1120 <make_zero_rows+0x20>
jmp    119a <make_zero_rows+0x9a>
nopw   0x0(%rax,%rax,1)
xor    %eax,%eax
jmp    1169 <make_zero_rows+0x69>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%rax
cmp    %r9,%rax
je     1190 <make_zero_rows+0x90>
mov    %rax,%rcx
cmp    %rax,%r10
je     117a <make_zero_rows+0x7a>
mov    (%rdi,%r10,8),%rax
addsd  (%rax,%rcx,8),%xmm1
lea    0x1(%rcx),%rax
cmp    %rax,%r10
je     1160 <make_zero_rows+0x60>
mov    (%rdi,%r10,8),%rdx
addsd  0x8(%rdx,%rcx,8),%xmm1
jmp    1160 <make_zero_rows+0x60>
nop
add    $0x2,%rcx
test   $0x1,%r8b
je     1120 <make_zero_rows+0x20>
cmp    %rcx,%r10
je     1120 <make_zero_rows+0x20>
mov    (%rdi,%r10,8),%rax
addsd  (%rax,%rcx,8),%xmm1
jmp    1120 <make_zero_rows+0x20>
ret

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

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