#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
void DoDiagonalOfQ(double *mat, const int n) {
  int i, j;
  if (NULL == mat)
    return;
  for (i = 0; i < n; i++) {
    mat[i * n + i] = 0.;
    double rsum = 0.;
    for (j = 0; j < n; j++)
      rsum += mat[i * n + j];
    mat[i * n + i] = -rsum;
  }
}


x86 assembly (clang, -O0):
<DoDiagonalOfQ>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
xor    %eax,%eax
cmp    -0x8(%rbp),%rax
jne    111c <DoDiagonalOfQ+0x1c>
jmp    11d1 <DoDiagonalOfQ+0xd1>
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11d1 <DoDiagonalOfQ+0xd1>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x10(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1191 <DoDiagonalOfQ+0x91>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1157 <DoDiagonalOfQ+0x57>
movsd  -0x20(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x10(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1123 <DoDiagonalOfQ+0x23>
pop    %rbp
ret

x86 assembly (clang, -O2):
<DoDiagonalOfQ>:
push   %r15
push   %r14
push   %r12
push   %rbx
test   %rdi,%rdi
je     11d0 <DoDiagonalOfQ+0xd0>
test   %esi,%esi
jle    11d0 <DoDiagonalOfQ+0xd0>
mov    %esi,%r9d
lea    -0x1(%r9),%r8
mov    %r9d,%edx
and    $0x3,%edx
mov    %r9d,%r12d
and    $0xfffffffc,%r12d
lea    0x18(%rdi),%rcx
lea    0x0(,%r9,8),%r10
xor    %r11d,%r11d
movapd 0xebd(%rip),%xmm0
mov    %rdi,%r15
jmp    1169 <DoDiagonalOfQ+0x69>
nopl   0x0(%rax,%rax,1)
xorpd  %xmm0,%xmm1
movlpd %xmm1,(%rdi,%r14,8)
add    $0x1,%r11
add    %r10,%rcx
add    %r10,%r15
cmp    %r9,%r11
je     11d0 <DoDiagonalOfQ+0xd0>
mov    %r11d,%eax
imul   %esi,%eax
movslq %eax,%r14
add    %r11,%r14
movq   $0x0,(%rdi,%r14,8)
xorpd  %xmm1,%xmm1
xor    %ebx,%ebx
cmp    $0x3,%r8
jb     11b0 <DoDiagonalOfQ+0xb0>
nopl   0x0(%rax)
addsd  -0x18(%rcx,%rbx,8),%xmm1
addsd  -0x10(%rcx,%rbx,8),%xmm1
addsd  -0x8(%rcx,%rbx,8),%xmm1
addsd  (%rcx,%rbx,8),%xmm1
add    $0x4,%rbx
cmp    %rbx,%r12
jne    1190 <DoDiagonalOfQ+0x90>
test   %rdx,%rdx
je     1150 <DoDiagonalOfQ+0x50>
lea    (%r15,%rbx,8),%rbx
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
addsd  (%rbx,%rax,8),%xmm1
add    $0x1,%rax
cmp    %rax,%rdx
jne    11c0 <DoDiagonalOfQ+0xc0>
jmp    1150 <DoDiagonalOfQ+0x50>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

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

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