void matsub(double **a, double **b, double **c, int n) {
  int i, j;
  for (j = 1; j <= n; j++)
    for (i = 1; i <= n; i++)
      c[i][j] = a[i][j] - b[i][j];
}


x86 assembly (clang, -O0):
<matsub>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     1199 <matsub+0x99>
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     1186 <matsub+0x86>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    112d <matsub+0x2d>
jmp    118b <matsub+0x8b>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111a <matsub+0x1a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<matsub>:
test   %ecx,%ecx
jle    11ca <matsub+0xca>
mov    %ecx,%r8d
lea    0x1(%r8),%r9
mov    %r8d,%r10d
and    $0xfffffffe,%r10d
neg    %r10
mov    $0x1,%r11d
jmp    113d <matsub+0x3d>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%r11
cmp    %r9,%r11
je     11ca <matsub+0xca>
mov    $0x1,%eax
cmp    $0x2,%r9d
je     11a1 <matsub+0xa1>
mov    $0x1,%eax
nopl   (%rax)
mov    (%rdi,%rax,8),%rcx
movsd  (%rcx,%r11,8),%xmm0
mov    (%rsi,%rax,8),%rcx
subsd  (%rcx,%r11,8),%xmm0
mov    (%rdx,%rax,8),%rcx
movsd  %xmm0,(%rcx,%r11,8)
mov    0x8(%rdi,%rax,8),%rcx
movsd  (%rcx,%r11,8),%xmm0
mov    0x8(%rsi,%rax,8),%rcx
subsd  (%rcx,%r11,8),%xmm0
mov    0x8(%rdx,%rax,8),%rcx
movsd  %xmm0,(%rcx,%r11,8)
lea    (%r10,%rax,1),%rcx
add    $0x2,%rcx
add    $0x2,%rax
cmp    $0x1,%rcx
jne    1150 <matsub+0x50>
test   $0x1,%r8b
je     1130 <matsub+0x30>
mov    (%rdi,%rax,8),%rcx
movsd  (%rcx,%r11,8),%xmm0
mov    (%rsi,%rax,8),%rcx
subsd  (%rcx,%r11,8),%xmm0
mov    (%rdx,%rax,8),%rax
movsd  %xmm0,(%rax,%r11,8)
jmp    1130 <matsub+0x30>
ret

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

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