double diff_mats(int m, int n, double *a, int a_row_stride, int a_col_stride,
                 double *b, int b_row_stride, int b_col_stride) {
  int i, j;
  double diff = 0.0;
  for (i = 0; i < m; i++) {
    for (j = 0; j < n; j++) {
      diff += abs(a[(i)*a_row_stride + (j)*a_col_stride] -
                  b[(i)*b_row_stride + (j)*b_col_stride]);
    }
  }
  return diff;
}


x86 assembly (clang, -O0):
<diff_mats>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11d1 <diff_mats+0xc1>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11be <diff_mats+0xae>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
mov    -0x28(%rbp),%edx
imul   -0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   0x10(%rbp),%ecx
mov    -0x28(%rbp),%edx
imul   0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <abs@plt>
cvtsi2sd %eax,%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1155 <diff_mats+0x45>
jmp    11c3 <diff_mats+0xb3>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1142 <diff_mats+0x32>
movsd  -0x30(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<diff_mats>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edi,%edi
jle    127c <diff_mats+0x17c>
movslq %r8d,%r8
movslq 0x40(%rsp),%r13
movslq %ecx,%rax
movslq 0x38(%rsp),%rcx
mov    %edi,%edi
mov    %rdi,-0x10(%rsp)
mov    %esi,%ebx
mov    %rbx,-0x28(%rsp)
and    $0xfffffffe,%ebx
mov    %rax,-0x40(%rsp)
lea    0x0(,%rax,8),%rax
mov    %rax,-0x18(%rsp)
mov    %r8,%rdi
shl    $0x4,%rdi
mov    %rcx,-0x48(%rsp)
lea    0x0(,%rcx,8),%rax
mov    %rax,-0x20(%rsp)
mov    %r13,%rax
shl    $0x4,%rax
xorpd  %xmm0,%xmm0
xor    %r10d,%r10d
mov    %r9,-0x30(%rsp)
mov    %rdx,-0x38(%rsp)
mov    %rdx,%r14
mov    %esi,-0x4c(%rsp)
jmp    1199 <diff_mats+0x99>
nop
add    $0x1,%r10
add    -0x18(%rsp),%r14
add    -0x20(%rsp),%r9
cmp    -0x10(%rsp),%r10
je     1280 <diff_mats+0x180>
test   %esi,%esi
jle    1180 <diff_mats+0x80>
cmp    $0x1,%esi
mov    %r9,-0x8(%rsp)
jne    11b0 <diff_mats+0xb0>
xor    %r11d,%r11d
jmp    1215 <diff_mats+0x115>
nopl   0x0(%rax)
mov    %r14,%rdx
xor    %r11d,%r11d
cs nopw 0x0(%rax,%rax,1)
movsd  (%rdx),%xmm1
subsd  (%r9),%xmm1
cvttsd2si %xmm1,%r12d
mov    %r12d,%r15d
neg    %r15d
cmovs  %r12d,%r15d
movsd  (%rdx,%r8,8),%xmm1
subsd  (%r9,%r13,8),%xmm1
cvttsd2si %xmm1,%ebp
xorps  %xmm1,%xmm1
cvtsi2sd %r15d,%xmm1
mov    %ebp,%ecx
neg    %ecx
cmovs  %ebp,%ecx
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
addsd  %xmm1,%xmm0
add    $0x2,%r11
add    %rdi,%rdx
add    %rax,%r9
cmp    %r11,%rbx
jne    11c0 <diff_mats+0xc0>
testb  $0x1,-0x28(%rsp)
mov    -0x4c(%rsp),%esi
mov    -0x8(%rsp),%r9
je     1180 <diff_mats+0x80>
mov    %r10,%rcx
imul   -0x40(%rsp),%rcx
mov    %r10,%rdx
imul   -0x48(%rsp),%rdx
mov    %r11,%rbp
imul   %r8,%rbp
add    %rcx,%rbp
mov    -0x38(%rsp),%rcx
movsd  (%rcx,%rbp,8),%xmm1
imul   %r13,%r11
add    %rdx,%r11
mov    -0x30(%rsp),%rcx
subsd  (%rcx,%r11,8),%xmm1
cvttsd2si %xmm1,%ecx
mov    %ecx,%edx
neg    %edx
cmovs  %ecx,%edx
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
addsd  %xmm1,%xmm0
jmp    1180 <diff_mats+0x80>
xorpd  %xmm0,%xmm0
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @diff_mats(i32 %0, i32 %1, double* %2, i32 %3, i32 %4, double* %5, i32 %6, i32 %7) {
%9 = alloca i32
%10 = alloca i32
%11 = alloca double*
%12 = alloca i32
%13 = alloca i32
%14 = alloca double*
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca double
store i32 %0, i32* %9
store i32 %1, i32* %10
store double* %2, double** %11
store i32 %3, i32* %12
store i32 %4, i32* %13
store double* %5, double** %14
store i32 %6, i32* %15
store i32 %7, i32* %16
store double 0.000000e+00, double* %19
store i32 0, i32* %17
br label %20
20: ; preds = %62, %8
%21 = load i32, i32* %17
%22 = load i32, i32* %9
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %65
24: ; preds = %20
store i32 0, i32* %18
br label %25
25: ; preds = %58, %24
%26 = load i32, i32* %18
%27 = load i32, i32* %10
%28 = icmp slt i32 %26, %27
br i1 %28, label %29, label %61
29: ; preds = %25
%30 = load double*, double** %11
%31 = load i32, i32* %17
%32 = load i32, i32* %12
%33 = mul nsw i32 %31, %32
%34 = load i32, i32* %18
%35 = load i32, i32* %13
%36 = mul nsw i32 %34, %35
%37 = add nsw i32 %33, %36
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds double, double* %30, i64 %38
%40 = load double, double* %39
%41 = load double*, double** %14
%42 = load i32, i32* %17
%43 = load i32, i32* %15
%44 = mul nsw i32 %42, %43
%45 = load i32, i32* %18
%46 = load i32, i32* %16
%47 = mul nsw i32 %45, %46
%48 = add nsw i32 %44, %47
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %41, i64 %49
%51 = load double, double* %50
%52 = fsub double %40, %51
%53 = fptosi double %52 to i32
%54 = call i32 @abs(i32 %53)
%55 = sitofp i32 %54 to double
%56 = load double, double* %19
%57 = fadd double %56, %55
store double %57, double* %19
br label %58
58: ; preds = %29
%59 = load i32, i32* %18
%60 = add nsw i32 %59, 1
store i32 %60, i32* %18
br label %25
61: ; preds = %25
br label %62
62: ; preds = %61
%63 = load i32, i32* %17
%64 = add nsw i32 %63, 1
store i32 %64, i32* %17
br label %20
65: ; preds = %20
%66 = load double, double* %19
ret double %66
}

LLVM IR (-O2):
define double @diff_mats(i32 %0, i32 %1, double* nocapture %2, i32 %3, i32 %4, double* nocapture %5, i32 %6, i32 %7) {
%9 = icmp sgt i32 %0, 0
br i1 %9, label %10, label %84
10: ; preds = %8
%11 = icmp sgt i32 %1, 0
%12 = sext i32 %4 to i64
%13 = sext i32 %7 to i64
%14 = sext i32 %3 to i64
%15 = sext i32 %6 to i64
%16 = zext i32 %0 to i64
%17 = zext i32 %1 to i64
%18 = and i64 %17, 1
%19 = icmp eq i32 %1, 1
%20 = and i64 %17, 4294967294
%21 = icmp eq i64 %18, 0
br label %22
22: ; preds = %10, %80
%23 = phi i64 [ 0, %10 ], [ %82, %80 ]
%24 = phi double [ 0.000000e+00, %10 ], [ %81, %80 ]
br i1 %11, label %25, label %80
25: ; preds = %22
%26 = mul nsw i64 %23, %14
%27 = mul nsw i64 %23, %15
br i1 %19, label %62, label %28
28: ; preds = %25, %28
%29 = phi i64 [ %59, %28 ], [ 0, %25 ]
%30 = phi double [ %58, %28 ], [ %24, %25 ]
%31 = phi i64 [ %60, %28 ], [ 0, %25 ]
%32 = mul nsw i64 %29, %12
%33 = add nsw i64 %32, %26
%34 = getelementptr inbounds double, double* %2, i64 %33
%35 = load double, double* %34
%36 = mul nsw i64 %29, %13
%37 = add nsw i64 %36, %27
%38 = getelementptr inbounds double, double* %5, i64 %37
%39 = load double, double* %38
%40 = fsub double %35, %39
%41 = fptosi double %40 to i32
%42 = tail call i32 @llvm.abs.i32(i32 %41, i1 true)
%43 = sitofp i32 %42 to double
%44 = fadd double %30, %43
%45 = or i64 %29, 1
%46 = mul nsw i64 %45, %12
%47 = add nsw i64 %46, %26
%48 = getelementptr inbounds double, double* %2, i64 %47
%49 = load double, double* %48
%50 = mul nsw i64 %45, %13
%51 = add nsw i64 %50, %27
%52 = getelementptr inbounds double, double* %5, i64 %51
%53 = load double, double* %52
%54 = fsub double %49, %53
%55 = fptosi double %54 to i32
%56 = tail call i32 @llvm.abs.i32(i32 %55, i1 true)
%57 = sitofp i32 %56 to double
%58 = fadd double %44, %57
%59 = add nuw nsw i64 %29, 2
%60 = add i64 %31, 2
%61 = icmp eq i64 %60, %20
br i1 %61, label %62, label %28
62: ; preds = %28, %25
%63 = phi double [ undef, %25 ], [ %58, %28 ]
%64 = phi i64 [ 0, %25 ], [ %59, %28 ]
%65 = phi double [ %24, %25 ], [ %58, %28 ]
br i1 %21, label %80, label %66
66: ; preds = %62
%67 = mul nsw i64 %64, %12
%68 = add nsw i64 %67, %26
%69 = getelementptr inbounds double, double* %2, i64 %68
%70 = load double, double* %69
%71 = mul nsw i64 %64, %13
%72 = add nsw i64 %71, %27
%73 = getelementptr inbounds double, double* %5, i64 %72
%74 = load double, double* %73
%75 = fsub double %70, %74
%76 = fptosi double %75 to i32
%77 = tail call i32 @llvm.abs.i32(i32 %76, i1 true)
%78 = sitofp i32 %77 to double
%79 = fadd double %65, %78
br label %80
80: ; preds = %66, %62, %22
%81 = phi double [ %24, %22 ], [ %63, %62 ], [ %79, %66 ]
%82 = add nuw nsw i64 %23, 1
%83 = icmp eq i64 %82, %16
br i1 %83, label %84, label %22
84: ; preds = %80, %8
%85 = phi double [ 0.000000e+00, %8 ], [ %81, %80 ]
ret double %85
}