void matNd_identity(double m[], int n) {
  for (int row = 0; row < n; row++)
    for (int col = 0; col < n; col++)
      if (row == col)
        m[row + col * n] = 1.0;
      else
        m[row + col * n] = 0.0;
}


x86 assembly (clang, -O0):
<matNd_identity>:
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    11a3 <matNd_identity+0xa3>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1190 <matNd_identity+0x90>
mov    -0x10(%rbp),%eax
cmp    -0x14(%rbp),%eax
jne    1162 <matNd_identity+0x62>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
mov    -0x14(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  0xea8(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    117d <matNd_identity+0x7d>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
mov    -0x14(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1182 <matNd_identity+0x82>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1125 <matNd_identity+0x25>
jmp    1195 <matNd_identity+0x95>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1112 <matNd_identity+0x12>
pop    %rbp
ret

x86 assembly (clang, -O2):
<matNd_identity>:
push   %r15
push   %r14
push   %rbx
test   %esi,%esi
jle    11db <matNd_identity+0xdb>
mov    %esi,%r8d
mov    %r8d,%r11d
and    $0xfffffffe,%r11d
mov    %r8,%r14
shl    $0x4,%r14
xor    %eax,%eax
movabs $0x3ff0000000000000,%r9
mov    %rdi,%r10
jmp    1145 <matNd_identity+0x45>
nop
mov    %r9,(%rdi,%r15,8)
add    $0x1,%rax
add    $0x8,%r10
cmp    %r8,%rax
je     11db <matNd_identity+0xdb>
mov    %rax,%r15
imul   %r8,%r15
add    %rax,%r15
cmp    $0x1,%esi
jne    1160 <matNd_identity+0x60>
xor    %ecx,%ecx
test   $0x1,%r8b
je     1134 <matNd_identity+0x34>
jmp    11be <matNd_identity+0xbe>
xchg   %ax,%ax
mov    %r10,%rdx
xor    %ebx,%ebx
jmp    1180 <matNd_identity+0x80>
nopw   0x0(%rax,%rax,1)
mov    %r9,(%rdi,%r15,8)
add    %r14,%rdx
add    $0x1,%rbx
cmp    %r11,%rbx
je     11b0 <matNd_identity+0xb0>
mov    %rbx,%rcx
cmp    %rbx,%rax
jne    1190 <matNd_identity+0x90>
mov    %r9,(%rdi,%r15,8)
jmp    1197 <matNd_identity+0x97>
xchg   %ax,%ax
movq   $0x0,(%rdx)
lea    0x1(%rcx),%rbx
cmp    %rbx,%rax
je     1170 <matNd_identity+0x70>
movq   $0x0,(%rdx,%r8,8)
jmp    1174 <matNd_identity+0x74>
nopw   0x0(%rax,%rax,1)
add    $0x2,%rcx
test   $0x1,%r8b
je     1134 <matNd_identity+0x34>
cmp    %rcx,%rax
je     1130 <matNd_identity+0x30>
imul   %r8,%rcx
add    %rax,%rcx
movq   $0x0,(%rdi,%rcx,8)
jmp    1134 <matNd_identity+0x34>
pop    %rbx
pop    %r14
pop    %r15
ret

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

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