void formatMatrix(double *ma, int dim) {
  int i, j;
  for (i = 0; i < dim; i++) {
    for (j = i + 1; j < dim; j++) {
      ma[i * dim + j] = 0.0;
    }
  }
}


x86 assembly (clang, -O0):
<formatMatrix>:
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    116d <formatMatrix+0x6d>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    115a <formatMatrix+0x5a>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1127 <formatMatrix+0x27>
jmp    115f <formatMatrix+0x5f>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1112 <formatMatrix+0x12>
pop    %rbp
ret

x86 assembly (clang, -O2):
<formatMatrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rdi,(%rsp)
test   %esi,%esi
jle    122a <formatMatrix+0x11a>
mov    %esi,%r12d
cmp    $0x1,%esi
jne    1143 <formatMatrix+0x33>
xor    %ebx,%ebx
test   $0x1,%r12b
jne    11fa <formatMatrix+0xea>
jmp    122a <formatMatrix+0x11a>
lea    -0x3(%rsi),%r13d
lea    0x2(%rsi),%ebp
mov    %rsi,0x8(%rsp)
lea    (%rsi,%rsi,1),%r14d
add    $0x2,%r14d
mov    %r12d,%eax
and    $0xfffffffe,%eax
neg    %rax
mov    %rax,0x10(%rsp)
mov    $0x1,%ebx
mov    $0x1,%r15d
jmp    11a0 <formatMatrix+0x90>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0xfffffffe,%r13d
add    %r14d,%ebp
add    %r14d,%r15d
mov    0x10(%rsp),%rax
add    %rbx,%rax
add    $0x2,%rax
add    $0x2,%rbx
cmp    $0x1,%rax
je     11eb <formatMatrix+0xdb>
cmp    %r12,%rbx
jae    11c3 <formatMatrix+0xb3>
lea    0x1(%r13),%eax
lea    0x8(,%rax,8),%rdx
movslq %r15d,%rax
mov    (%rsp),%rcx
lea    (%rcx,%rax,8),%rdi
xor    %esi,%esi
call   1030 <memset@plt>
lea    0x1(%rbx),%rax
cmp    %r12,%rax
jae    1180 <formatMatrix+0x70>
mov    %r13d,%eax
lea    0x8(,%rax,8),%rdx
movslq %ebp,%rax
mov    (%rsp),%rcx
lea    (%rcx,%rax,8),%rdi
xor    %esi,%esi
call   1030 <memset@plt>
jmp    1180 <formatMatrix+0x70>
add    $0xffffffffffffffff,%rbx
mov    0x8(%rsp),%rsi
test   $0x1,%r12b
je     122a <formatMatrix+0x11a>
lea    0x1(%rbx),%rax
cmp    %r12,%rax
jae    122a <formatMatrix+0x11a>
lea    0x1(%rsi),%eax
add    $0xfffffffe,%esi
sub    %ebx,%esi
lea    0x8(,%rsi,8),%rdx
imul   %ebx,%eax
add    $0x1,%eax
cltq
mov    (%rsp),%rcx
lea    (%rcx,%rax,8),%rdi
xor    %esi,%esi
call   1030 <memset@plt>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @formatMatrix(double* %0, i32 %1) {
%3 = alloca double*
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
store double* %0, double** %3
store i32 %1, i32* %4
store i32 0, i32* %5
br label %7
7: ; preds = %31, %2
%8 = load i32, i32* %5
%9 = load i32, i32* %4
%10 = icmp slt i32 %8, %9
br i1 %10, label %11, label %34
11: ; preds = %7
%12 = load i32, i32* %5
%13 = add nsw i32 %12, 1
store i32 %13, i32* %6
br label %14
14: ; preds = %27, %11
%15 = load i32, i32* %6
%16 = load i32, i32* %4
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %30
18: ; preds = %14
%19 = load double*, double** %3
%20 = load i32, i32* %5
%21 = load i32, i32* %4
%22 = mul nsw i32 %20, %21
%23 = load i32, i32* %6
%24 = add nsw i32 %22, %23
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds double, double* %19, i64 %25
store double 0.000000e+00, double* %26
br label %27
27: ; preds = %18
%28 = load i32, i32* %6
%29 = add nsw i32 %28, 1
store i32 %29, i32* %6
br label %14
30: ; preds = %14
br label %31
31: ; preds = %30
%32 = load i32, i32* %5
%33 = add nsw i32 %32, 1
store i32 %33, i32* %5
br label %7
34: ; preds = %7
ret void
}

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