int cd_mat_sub(double *dst, const double *src, int m, int n) {
  int i, j;
  for (i = 0; i < m; i++)
    for (j = 0; j < n; j++)
      dst[i * n + j] -= src[i * n + j];
  return 0;
}


x86 assembly (clang, -O0):
<cd_mat_sub>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    118e <cd_mat_sub+0x8e>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    117b <cd_mat_sub+0x7b>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    112c <cd_mat_sub+0x2c>
jmp    1180 <cd_mat_sub+0x80>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1119 <cd_mat_sub+0x19>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<cd_mat_sub>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edx,%edx
jle    1367 <cd_mat_sub+0x267>
movslq %ecx,%r15
mov    %edx,%r12d
mov    %r15d,%r8d
mov    %r8d,%eax
and    $0xfffffffc,%eax
mov    %rax,-0x10(%rsp)
add    $0xfffffffffffffffc,%rax
mov    %rax,-0x18(%rsp)
mov    %rax,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    %rdx,%rax
mov    %rdx,-0x20(%rsp)
and    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x28(%rsp)
mov    %r8,%rax
neg    %rax
mov    %rax,-0x8(%rsp)
lea    0x30(%rdi),%rbx
lea    0x0(,%r15,8),%r13
lea    0x30(%rsi),%rdx
lea    0x8(%rsi),%rbp
lea    0x8(%rdi),%r14
xor    %r10d,%r10d
mov    %r12,-0x30(%rsp)
jmp    1199 <cd_mat_sub+0x99>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r10
add    %r13,%rbx
add    %r13,%rdx
add    %r13,%rbp
add    %r13,%r14
cmp    %r12,%r10
je     1367 <cd_mat_sub+0x267>
test   %ecx,%ecx
jle    1180 <cd_mat_sub+0x80>
mov    %r10,%r11
imul   %r15,%r11
cmp    $0x4,%ecx
jae    11c0 <cd_mat_sub+0xc0>
xor    %r9d,%r9d
jmp    12e0 <cd_mat_sub+0x1e0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r15,%r12
lea    (%r11,%r8,1),%r15
lea    (%rdi,%r11,8),%r9
lea    (%rsi,%r15,8),%rax
cmp    %rax,%r9
jae    11f1 <cd_mat_sub+0xf1>
lea    (%rdi,%r15,8),%r9
lea    (%rsi,%r11,8),%rax
cmp    %r9,%rax
jae    11f1 <cd_mat_sub+0xf1>
xor    %r9d,%r9d
mov    %r12,%r15
mov    -0x30(%rsp),%r12
jmp    12e0 <cd_mat_sub+0x1e0>
mov    %rsi,%rax
cmpq   $0x0,-0x18(%rsp)
je     1354 <cd_mat_sub+0x254>
mov    -0x28(%rsp),%r15
xor    %r9d,%r9d
nopl   0x0(%rax,%rax,1)
movupd -0x30(%rdx,%r9,8),%xmm0
movupd -0x20(%rdx,%r9,8),%xmm1
movupd -0x30(%rbx,%r9,8),%xmm2
subpd  %xmm0,%xmm2
movupd -0x20(%rbx,%r9,8),%xmm0
subpd  %xmm1,%xmm0
movupd -0x10(%rbx,%r9,8),%xmm1
movupd (%rbx,%r9,8),%xmm3
movupd %xmm2,-0x30(%rbx,%r9,8)
movupd %xmm0,-0x20(%rbx,%r9,8)
movupd -0x10(%rdx,%r9,8),%xmm0
subpd  %xmm0,%xmm1
movupd (%rdx,%r9,8),%xmm0
subpd  %xmm0,%xmm3
movupd %xmm1,-0x10(%rbx,%r9,8)
movupd %xmm3,(%rbx,%r9,8)
add    $0x8,%r9
add    $0xfffffffffffffffe,%r15
jne    1210 <cd_mat_sub+0x110>
testb  $0x1,-0x20(%rsp)
je     12b7 <cd_mat_sub+0x1b7>
add    %r11,%r9
mov    %rax,%rsi
movupd (%rax,%r9,8),%xmm0
movupd 0x10(%rax,%r9,8),%xmm1
movupd (%rdi,%r9,8),%xmm2
subpd  %xmm0,%xmm2
movupd 0x10(%rdi,%r9,8),%xmm0
subpd  %xmm1,%xmm0
movupd %xmm2,(%rdi,%r9,8)
movupd %xmm0,0x10(%rdi,%r9,8)
mov    -0x10(%rsp),%rsi
mov    %rsi,%r9
cmp    %r8,%rsi
mov    %r12,%r15
mov    -0x30(%rsp),%r12
mov    %rax,%rsi
je     1180 <cd_mat_sub+0x80>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r9,%rax
test   $0x1,%r8b
je     1305 <cd_mat_sub+0x205>
add    %r9,%r11
movsd  (%rdi,%r11,8),%xmm0
subsd  (%rsi,%r11,8),%xmm0
movsd  %xmm0,(%rdi,%r11,8)
mov    %r9,%rax
or     $0x1,%rax
not    %r9
cmp    -0x8(%rsp),%r9
je     1180 <cd_mat_sub+0x80>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  -0x8(%r14,%rax,8),%xmm0
subsd  -0x8(%rbp,%rax,8),%xmm0
movsd  (%r14,%rax,8),%xmm1
movsd  %xmm0,-0x8(%r14,%rax,8)
subsd  0x0(%rbp,%rax,8),%xmm1
movsd  %xmm1,(%r14,%rax,8)
add    $0x2,%rax
cmp    %rax,%r8
jne    1320 <cd_mat_sub+0x220>
jmp    1180 <cd_mat_sub+0x80>
xor    %r9d,%r9d
testb  $0x1,-0x20(%rsp)
jne    1282 <cd_mat_sub+0x182>
jmp    12b7 <cd_mat_sub+0x1b7>
xor    %eax,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @cd_mat_sub(double* nocapture %0, double* nocapture %1, i32 %2, i32 %3) {
%5 = icmp sgt i32 %2, 0
br i1 %5, label %6, label %136
6: ; preds = %4
%7 = icmp sgt i32 %3, 0
%8 = sext i32 %3 to i64
%9 = zext i32 %2 to i64
%10 = zext i32 %3 to i64
%11 = and i64 %10, 4294967292
%12 = add nsw i64 %11, -4
%13 = lshr exact i64 %12, 2
%14 = add nuw nsw i64 %13, 1
%15 = icmp ult i32 %3, 4
%16 = and i64 %10, 4294967292
%17 = and i64 %14, 1
%18 = icmp eq i64 %12, 0
%19 = and i64 %14, 9223372036854775806
%20 = icmp eq i64 %17, 0
%21 = icmp eq i64 %16, %10
%22 = and i64 %10, 1
%23 = icmp eq i64 %22, 0
%24 = sub nsw i64 0, %10
br label %25
25: ; preds = %6, %133
%26 = phi i64 [ 0, %6 ], [ %134, %133 ]
%27 = mul i64 %26, %8
%28 = getelementptr double, double* %0, i64 %27
%29 = add i64 %27, %10
%30 = getelementptr double, double* %0, i64 %29
%31 = getelementptr double, double* %1, i64 %27
%32 = getelementptr double, double* %1, i64 %29
br i1 %7, label %33, label %133
33: ; preds = %25
%34 = mul nsw i64 %26, %8
br i1 %15, label %102, label %35
35: ; preds = %33
%36 = icmp ult double* %28, %32
%37 = icmp ult double* %31, %30
%38 = and i1 %36, %37
br i1 %38, label %102, label %39
39: ; preds = %35
br i1 %18, label %81, label %40
40: ; preds = %39, %40
%41 = phi i64 [ %78, %40 ], [ 0, %39 ]
%42 = phi i64 [ %79, %40 ], [ 0, %39 ]
%43 = add nsw i64 %41, %34
%44 = getelementptr inbounds double, double* %1, i64 %43
%45 = bitcast double* %44 to <2 x double>*
%46 = load <2 x double>, <2 x double>* %45
%47 = getelementptr inbounds double, double* %44, i64 2
%48 = bitcast double* %47 to <2 x double>*
%49 = load <2 x double>, <2 x double>* %48
%50 = getelementptr inbounds double, double* %0, i64 %43
%51 = bitcast double* %50 to <2 x double>*
%52 = load <2 x double>, <2 x double>* %51
%53 = getelementptr inbounds double, double* %50, i64 2
%54 = bitcast double* %53 to <2 x double>*
%55 = load <2 x double>, <2 x double>* %54
%56 = fsub <2 x double> %52, %46
%57 = fsub <2 x double> %55, %49
%58 = bitcast double* %50 to <2 x double>*
store <2 x double> %56, <2 x double>* %58
%59 = bitcast double* %53 to <2 x double>*
store <2 x double> %57, <2 x double>* %59
%60 = or i64 %41, 4
%61 = add nsw i64 %60, %34
%62 = getelementptr inbounds double, double* %1, i64 %61
%63 = bitcast double* %62 to <2 x double>*
%64 = load <2 x double>, <2 x double>* %63
%65 = getelementptr inbounds double, double* %62, i64 2
%66 = bitcast double* %65 to <2 x double>*
%67 = load <2 x double>, <2 x double>* %66
%68 = getelementptr inbounds double, double* %0, i64 %61
%69 = bitcast double* %68 to <2 x double>*
%70 = load <2 x double>, <2 x double>* %69
%71 = getelementptr inbounds double, double* %68, i64 2
%72 = bitcast double* %71 to <2 x double>*
%73 = load <2 x double>, <2 x double>* %72
%74 = fsub <2 x double> %70, %64
%75 = fsub <2 x double> %73, %67
%76 = bitcast double* %68 to <2 x double>*
store <2 x double> %74, <2 x double>* %76
%77 = bitcast double* %71 to <2 x double>*
store <2 x double> %75, <2 x double>* %77
%78 = add nuw i64 %41, 8
%79 = add i64 %42, 2
%80 = icmp eq i64 %79, %19
br i1 %80, label %81, label %40
81: ; preds = %40, %39
%82 = phi i64 [ 0, %39 ], [ %78, %40 ]
br i1 %20, label %101, label %83
83: ; preds = %81
%84 = add nsw i64 %82, %34
%85 = getelementptr inbounds double, double* %1, i64 %84
%86 = bitcast double* %85 to <2 x double>*
%87 = load <2 x double>, <2 x double>* %86
%88 = getelementptr inbounds double, double* %85, i64 2
%89 = bitcast double* %88 to <2 x double>*
%90 = load <2 x double>, <2 x double>* %89
%91 = getelementptr inbounds double, double* %0, i64 %84
%92 = bitcast double* %91 to <2 x double>*
%93 = load <2 x double>, <2 x double>* %92
%94 = getelementptr inbounds double, double* %91, i64 2
%95 = bitcast double* %94 to <2 x double>*
%96 = load <2 x double>, <2 x double>* %95
%97 = fsub <2 x double> %93, %87
%98 = fsub <2 x double> %96, %90
%99 = bitcast double* %91 to <2 x double>*
store <2 x double> %97, <2 x double>* %99
%100 = bitcast double* %94 to <2 x double>*
store <2 x double> %98, <2 x double>* %100
br label %101
101: ; preds = %81, %83
br i1 %21, label %133, label %102
102: ; preds = %35, %33, %101
%103 = phi i64 [ 0, %35 ], [ 0, %33 ], [ %16, %101 ]
%104 = xor i64 %103, -1
br i1 %23, label %113, label %105
105: ; preds = %102
%106 = add nsw i64 %103, %34
%107 = getelementptr inbounds double, double* %1, i64 %106
%108 = load double, double* %107
%109 = getelementptr inbounds double, double* %0, i64 %106
%110 = load double, double* %109
%111 = fsub double %110, %108
store double %111, double* %109
%112 = or i64 %103, 1
br label %113
113: ; preds = %105, %102
%114 = phi i64 [ %103, %102 ], [ %112, %105 ]
%115 = icmp eq i64 %104, %24
br i1 %115, label %133, label %116
116: ; preds = %113, %116
%117 = phi i64 [ %131, %116 ], [ %114, %113 ]
%118 = add nsw i64 %117, %34
%119 = getelementptr inbounds double, double* %1, i64 %118
%120 = load double, double* %119
%121 = getelementptr inbounds double, double* %0, i64 %118
%122 = load double, double* %121
%123 = fsub double %122, %120
store double %123, double* %121
%124 = add nuw nsw i64 %117, 1
%125 = add nsw i64 %124, %34
%126 = getelementptr inbounds double, double* %1, i64 %125
%127 = load double, double* %126
%128 = getelementptr inbounds double, double* %0, i64 %125
%129 = load double, double* %128
%130 = fsub double %129, %127
store double %130, double* %128
%131 = add nuw nsw i64 %117, 2
%132 = icmp eq i64 %131, %10
br i1 %132, label %133, label %116
133: ; preds = %113, %116, %101, %25
%134 = add nuw nsw i64 %26, 1
%135 = icmp eq i64 %134, %9
br i1 %135, label %136, label %25
136: ; preds = %133, %4
ret i32 0
}