double trc_dmat(int n, double **a) {
  double _sum = 0;
  while (n--)
    _sum += a[n][n];
  return _sum;
}


x86 assembly (clang, -O0):
<trc_dmat>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x4(%rbp),%eax
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x4(%rbp)
cmp    $0x0,%eax
je     114b <trc_dmat+0x4b>
mov    -0x10(%rbp),%rax
movslq -0x4(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
jmp    1113 <trc_dmat+0x13>
movsd  -0x18(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<trc_dmat>:
test   %edi,%edi
je     1175 <trc_dmat+0x75>
movslq %edi,%rax
lea    -0x1(%rdi),%r8d
xorpd  %xmm0,%xmm0
and    $0x3,%edi
je     1137 <trc_dmat+0x37>
mov    %rax,%rdx
nopw   0x0(%rax,%rax,1)
lea    -0x1(%rdx),%rax
mov    -0x8(%rsi,%rdx,8),%rcx
addsd  -0x8(%rcx,%rdx,8),%xmm0
mov    %rax,%rdx
add    $0xffffffff,%edi
jne    1120 <trc_dmat+0x20>
cmp    $0x3,%r8d
jb     1174 <trc_dmat+0x74>
nopl   (%rax)
mov    -0x8(%rsi,%rax,8),%rcx
addsd  -0x8(%rcx,%rax,8),%xmm0
mov    -0x10(%rsi,%rax,8),%rcx
addsd  -0x10(%rcx,%rax,8),%xmm0
mov    -0x20(%rsi,%rax,8),%rcx
mov    -0x18(%rsi,%rax,8),%rdx
addsd  -0x18(%rdx,%rax,8),%xmm0
addsd  -0x20(%rcx,%rax,8),%xmm0
add    $0xfffffffffffffffc,%rax
test   %eax,%eax
jne    1140 <trc_dmat+0x40>
ret
xorps  %xmm0,%xmm0
ret

LLVM IR (-O0):
define double @trc_dmat(i32 %0, double** %1) {
%3 = alloca i32
%4 = alloca double**
%5 = alloca double
store i32 %0, i32* %3
store double** %1, double*** %4
store double 0.000000e+00, double* %5
br label %6
6: ; preds = %10, %2
%7 = load i32, i32* %3
%8 = add nsw i32 %7, -1
store i32 %8, i32* %3
%9 = icmp ne i32 %7, 0
br i1 %9, label %10, label %22
10: ; preds = %6
%11 = load double**, double*** %4
%12 = load i32, i32* %3
%13 = sext i32 %12 to i64
%14 = getelementptr inbounds double*, double** %11, i64 %13
%15 = load double*, double** %14
%16 = load i32, i32* %3
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds double, double* %15, i64 %17
%19 = load double, double* %18
%20 = load double, double* %5
%21 = fadd double %20, %19
store double %21, double* %5
br label %6
22: ; preds = %6
%23 = load double, double* %5
ret double %23
}

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