#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double fabs(double);
int compare_matrix(double *a, double *b, int n) {
  int i, j;
  double eps = 1e-6;
  for (i = 0; i < n; ++i) {
    for (j = 0; j < n; ++j) {
      if (fabs(a[i * n + j] - b[i * n + j]) > eps) {
        return 0;
      }
    }
  }
  return 1;
}


x86 assembly (clang, -O0):
<compare_matrix>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
movsd  0xee9(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11b6 <compare_matrix+0xb6>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11a3 <compare_matrix+0xa3>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movaps 0xe9b(%rip),%xmm1
pand   %xmm1,%xmm0
ucomisd -0x30(%rbp),%xmm0
jbe    1190 <compare_matrix+0x90>
movl   $0x0,-0x4(%rbp)
jmp    11bd <compare_matrix+0xbd>
jmp    1195 <compare_matrix+0x95>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1136 <compare_matrix+0x36>
jmp    11a8 <compare_matrix+0xa8>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1123 <compare_matrix+0x23>
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<compare_matrix>:
mov    $0x1,%eax
test   %edx,%edx
jle    115c <compare_matrix+0x5c>
mov    %edx,%ecx
lea    0x0(,%rcx,8),%r8
xor    %r9d,%r9d
movapd 0xee2(%rip),%xmm0
movsd  0xeea(%rip),%xmm1
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%rdx,8),%xmm2
subsd  (%rsi,%rdx,8),%xmm2
andpd  %xmm0,%xmm2
ucomisd %xmm1,%xmm2
ja     115d <compare_matrix+0x5d>
add    $0x1,%rdx
cmp    %rdx,%rcx
jne    1130 <compare_matrix+0x30>
add    $0x1,%r9
add    %r8,%rdi
add    %r8,%rsi
cmp    %rcx,%r9
jne    1126 <compare_matrix+0x26>
ret
xor    %eax,%eax
ret

LLVM IR (-O0):
define i32 @compare_matrix(double* %0, double* %1, i32 %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca double*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca double
store double* %0, double** %5
store double* %1, double** %6
store i32 %2, i32* %7
store double 0x3EB0C6F7A0B5ED8D, double* %10
store i32 0, i32* %8
br label %11
11: ; preds = %49, %3
%12 = load i32, i32* %8
%13 = load i32, i32* %7
%14 = icmp slt i32 %12, %13
br i1 %14, label %15, label %52
15: ; preds = %11
store i32 0, i32* %9
br label %16
16: ; preds = %45, %15
%17 = load i32, i32* %9
%18 = load i32, i32* %7
%19 = icmp slt i32 %17, %18
br i1 %19, label %20, label %48
20: ; preds = %16
%21 = load double*, double** %5
%22 = load i32, i32* %8
%23 = load i32, i32* %7
%24 = mul nsw i32 %22, %23
%25 = load i32, i32* %9
%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** %6
%31 = load i32, i32* %8
%32 = load i32, i32* %7
%33 = mul nsw i32 %31, %32
%34 = load i32, i32* %9
%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 %29, %38
%40 = call double @llvm.fabs.f64(double %39)
%41 = load double, double* %10
%42 = fcmp ogt double %40, %41
br i1 %42, label %43, label %44
43: ; preds = %20
store i32 0, i32* %4
br label %53
44: ; preds = %20
br label %45
45: ; preds = %44
%46 = load i32, i32* %9
%47 = add nsw i32 %46, 1
store i32 %47, i32* %9
br label %16
48: ; preds = %16
br label %49
49: ; preds = %48
%50 = load i32, i32* %8
%51 = add nsw i32 %50, 1
store i32 %51, i32* %8
br label %11
52: ; preds = %11
store i32 1, i32* %4
br label %53
53: ; preds = %52, %43
%54 = load i32, i32* %4
ret i32 %54
}

LLVM IR (-O2):
define i32 @compare_matrix(double* nocapture %0, double* nocapture %1, i32 %2) {
%4 = icmp sgt i32 %2, 0
br i1 %4, label %5, label %27
5: ; preds = %3
%6 = zext i32 %2 to i64
%7 = zext i32 %2 to i64
br label %8
8: ; preds = %24, %5
%9 = phi i64 [ 0, %5 ], [ %25, %24 ]
%10 = mul nsw i64 %9, %6
br label %13
11: ; preds = %13
%12 = icmp eq i64 %23, %7
br i1 %12, label %24, label %13
13: ; preds = %8, %11
%14 = phi i64 [ 0, %8 ], [ %23, %11 ]
%15 = add nuw nsw i64 %14, %10
%16 = getelementptr inbounds double, double* %0, i64 %15
%17 = load double, double* %16
%18 = getelementptr inbounds double, double* %1, i64 %15
%19 = load double, double* %18
%20 = fsub double %17, %19
%21 = tail call double @llvm.fabs.f64(double %20)
%22 = fcmp ogt double %21, 0x3EB0C6F7A0B5ED8D
%23 = add nuw nsw i64 %14, 1
br i1 %22, label %27, label %11
24: ; preds = %11
%25 = add nuw nsw i64 %9, 1
%26 = icmp eq i64 %25, %7
br i1 %26, label %27, label %8
27: ; preds = %24, %13, %3
%28 = phi i32 [ 1, %3 ], [ 0, %13 ], [ 1, %24 ]
ret i32 %28
}