#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 */
int /*<<< orphan*/ printf(char *, ...);
void compare(int n, double wref[][n], double w[][n]) {
  double maxdiff, this_diff;
  int numdiffs;
  int i, j;
  numdiffs = 0;
  maxdiff = 0;
  for (i = 0; i < n; i++)
    for (j = 0; j < n; j++) {
      this_diff = wref[i][j] - w[i][j];
      if (this_diff < 0)
        this_diff = -1.0 * this_diff;
      if (this_diff > 0.0000001) {
        numdiffs++;
        if (this_diff > maxdiff)
          maxdiff = this_diff;
      }
    }
  if (numdiffs > 0)
    printf("%d Diffs found over threshold %f; Max Diff = %f\n", numdiffs,
           0.0000001, maxdiff);
  else
    printf("Passed Correctness Check\n");
}


x86 assembly (clang, -O0):
<compare>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x4(%rbp),%eax
mov    %rax,-0x48(%rbp)
mov    -0x4(%rbp),%eax
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x2c(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1230 <compare+0x120>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    121d <compare+0x10d>
mov    -0x40(%rbp),%rdx
mov    -0x48(%rbp),%rsi
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
imul   %rsi,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    11cb <compare+0xbb>
movsd  0xe47(%rip),%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  0xe28(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    120a <compare+0xfa>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
movsd  -0x28(%rbp),%xmm0
ucomisd -0x20(%rbp),%xmm0
jbe    1205 <compare+0xf5>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    120a <compare+0xfa>
jmp    120f <compare+0xff>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    115a <compare+0x4a>
jmp    1222 <compare+0x112>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1147 <compare+0x37>
cmpl   $0x0,-0x2c(%rbp)
jle    125d <compare+0x14d>
mov    -0x2c(%rbp),%esi
movsd  -0x20(%rbp),%xmm1
lea    0xdc7(%rip),%rdi
movsd  0xdaf(%rip),%xmm0
mov    $0x2,%al
call   1030 <printf@plt>
jmp    126b <compare+0x15b>
lea    0xddd(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<compare>:
push   %r15
push   %r14
push   %r12
push   %rbx
test   %edi,%edi
jle    128b <compare+0x16b>
mov    %edi,%r8d
mov    %r8d,%r9d
and    $0xfffffffe,%r9d
lea    0x8(%rsi),%rcx
lea    0x0(,%r8,8),%r10
lea    0x8(%rdx),%rax
xorpd  %xmm0,%xmm0
xor    %r14d,%r14d
movapd 0xea8(%rip),%xmm2
movsd  0xeb0(%rip),%xmm3
xor    %r11d,%r11d
xorpd  %xmm1,%xmm1
jmp    1183 <compare+0x63>
nopl   0x0(%rax)
add    $0x1,%r14
add    %r10,%rcx
add    %r10,%rax
cmp    %r8,%r14
je     1266 <compare+0x146>
xor    %ebx,%ebx
cmp    $0x1,%edi
jne    11e9 <compare+0xc9>
test   $0x1,%r8b
je     1170 <compare+0x50>
mov    %r14,%r15
imul   %r8,%r15
lea    (%rsi,%r15,8),%r12
lea    (%rdx,%r15,8),%r15
movsd  (%r12,%rbx,8),%xmm5
subsd  (%r15,%rbx,8),%xmm5
movapd %xmm5,%xmm6
xorpd  %xmm2,%xmm6
movapd %xmm5,%xmm4
cmpltsd %xmm0,%xmm4
andpd  %xmm4,%xmm6
andnpd %xmm5,%xmm4
orpd   %xmm6,%xmm4
ucomisd %xmm3,%xmm4
jbe    1170 <compare+0x50>
add    $0x1,%r11d
ucomisd %xmm1,%xmm4
jbe    1170 <compare+0x50>
movapd %xmm4,%xmm1
jmp    1170 <compare+0x50>
xchg   %ax,%ax
add    $0x2,%rbx
cmp    %rbx,%r9
je     118a <compare+0x6a>
movsd  -0x8(%rcx,%rbx,8),%xmm5
subsd  -0x8(%rax,%rbx,8),%xmm5
movapd %xmm5,%xmm6
xorpd  %xmm2,%xmm6
movapd %xmm5,%xmm4
cmpltsd %xmm0,%xmm4
andpd  %xmm4,%xmm6
andnpd %xmm5,%xmm4
orpd   %xmm6,%xmm4
ucomisd %xmm3,%xmm4
jbe    1226 <compare+0x106>
add    $0x1,%r11d
ucomisd %xmm1,%xmm4
jbe    1226 <compare+0x106>
movapd %xmm4,%xmm1
movsd  (%rcx,%rbx,8),%xmm5
subsd  (%rax,%rbx,8),%xmm5
movapd %xmm5,%xmm6
xorpd  %xmm2,%xmm6
movapd %xmm5,%xmm4
cmpltsd %xmm0,%xmm4
andpd  %xmm4,%xmm6
andnpd %xmm5,%xmm4
orpd   %xmm6,%xmm4
ucomisd %xmm3,%xmm4
jbe    11e0 <compare+0xc0>
add    $0x1,%r11d
ucomisd %xmm1,%xmm4
jbe    11e0 <compare+0xc0>
movapd %xmm4,%xmm1
jmp    11e0 <compare+0xc0>
test   %r11d,%r11d
jle    128b <compare+0x16b>
lea    0xda6(%rip),%rdi
movsd  0xd96(%rip),%xmm0
mov    %r11d,%esi
mov    $0x2,%al
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
jmp    1040 <printf@plt>
lea    0xdb7(%rip),%rdi
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
jmp    1030 <puts@plt>

LLVM IR (-O0):
@.str = [49 x i8] c"%d Diffs found over threshold %f; Max Diff = %f\0A\00"
@.str.1 = [26 x i8] c"Passed Correctness Check\0A\00"
define void @compare(i32 %0, double* %1, double* %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca double*
%7 = alloca double
%8 = alloca double
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
store i32 %0, i32* %4
store double* %1, double** %5
store double* %2, double** %6
%12 = load i32, i32* %4
%13 = zext i32 %12 to i64
%14 = load i32, i32* %4
%15 = zext i32 %14 to i64
store i32 0, i32* %9
store double 0.000000e+00, double* %7
store i32 0, i32* %10
br label %16
16: ; preds = %67, %3
%17 = load i32, i32* %10
%18 = load i32, i32* %4
%19 = icmp slt i32 %17, %18
br i1 %19, label %20, label %70
20: ; preds = %16
store i32 0, i32* %11
br label %21
21: ; preds = %63, %20
%22 = load i32, i32* %11
%23 = load i32, i32* %4
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %66
25: ; preds = %21
%26 = load double*, double** %5
%27 = load i32, i32* %10
%28 = sext i32 %27 to i64
%29 = mul nsw i64 %28, %13
%30 = getelementptr inbounds double, double* %26, i64 %29
%31 = load i32, i32* %11
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %30, i64 %32
%34 = load double, double* %33
%35 = load double*, double** %6
%36 = load i32, i32* %10
%37 = sext i32 %36 to i64
%38 = mul nsw i64 %37, %15
%39 = getelementptr inbounds double, double* %35, i64 %38
%40 = load i32, i32* %11
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds double, double* %39, i64 %41
%43 = load double, double* %42
%44 = fsub double %34, %43
store double %44, double* %8
%45 = load double, double* %8
%46 = fcmp olt double %45, 0.000000e+00
br i1 %46, label %47, label %50
47: ; preds = %25
%48 = load double, double* %8
%49 = fmul double -1.000000e+00, %48
store double %49, double* %8
br label %50
50: ; preds = %47, %25
%51 = load double, double* %8
%52 = fcmp ogt double %51, 0x3E7AD7F29ABCAF48
br i1 %52, label %53, label %62
53: ; preds = %50
%54 = load i32, i32* %9
%55 = add nsw i32 %54, 1
store i32 %55, i32* %9
%56 = load double, double* %8
%57 = load double, double* %7
%58 = fcmp ogt double %56, %57
br i1 %58, label %59, label %61
59: ; preds = %53
%60 = load double, double* %8
store double %60, double* %7
br label %61
61: ; preds = %59, %53
br label %62
62: ; preds = %61, %50
br label %63
63: ; preds = %62
%64 = load i32, i32* %11
%65 = add nsw i32 %64, 1
store i32 %65, i32* %11
br label %21
66: ; preds = %21
br label %67
67: ; preds = %66
%68 = load i32, i32* %10
%69 = add nsw i32 %68, 1
store i32 %69, i32* %10
br label %16
70: ; preds = %16
%71 = load i32, i32* %9
%72 = icmp sgt i32 %71, 0
br i1 %72, label %73, label %77
73: ; preds = %70
%74 = load i32, i32* %9
%75 = load double, double* %7
%76 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([49 x i8], [49 x i8]* @.str, i64 0, i64 0), i32 %74, double 0x3E7AD7F29ABCAF48, double %75)
br label %79
77: ; preds = %70
%78 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([26 x i8], [26 x i8]* @.str.1, i64 0, i64 0))
br label %79
79: ; preds = %77, %73
ret void
}

LLVM IR (-O2):
@.str = [49 x i8] c"%d Diffs found over threshold %f; Max Diff = %f\0A\00"
@str = [25 x i8] c"Passed Correctness Check\00"
define void @compare(i32 %0, double* nocapture %1, double* nocapture %2) {
%4 = zext i32 %0 to i64
%5 = icmp sgt i32 %0, 0
br i1 %5, label %6, label %88
6: ; preds = %3
%7 = and i64 %4, 1
%8 = icmp eq i32 %0, 1
%9 = and i64 %4, 4294967294
%10 = icmp eq i64 %7, 0
br label %11
11: ; preds = %6, %79
%12 = phi i64 [ %82, %79 ], [ 0, %6 ]
%13 = phi i32 [ %81, %79 ], [ 0, %6 ]
%14 = phi double [ %80, %79 ], [ 0.000000e+00, %6 ]
%15 = mul nuw nsw i64 %12, %4
%16 = getelementptr inbounds double, double* %1, i64 %15
%17 = getelementptr inbounds double, double* %2, i64 %15
br i1 %8, label %59, label %18
18: ; preds = %11, %53
%19 = phi i64 [ %56, %53 ], [ 0, %11 ]
%20 = phi i32 [ %55, %53 ], [ %13, %11 ]
%21 = phi double [ %54, %53 ], [ %14, %11 ]
%22 = phi i64 [ %57, %53 ], [ 0, %11 ]
%23 = getelementptr inbounds double, double* %16, i64 %19
%24 = load double, double* %23
%25 = getelementptr inbounds double, double* %17, i64 %19
%26 = load double, double* %25
%27 = fsub double %24, %26
%28 = fcmp olt double %27, 0.000000e+00
%29 = fneg double %27
%30 = select i1 %28, double %29, double %27
%31 = fcmp ogt double %30, 0x3E7AD7F29ABCAF48
br i1 %31, label %32, label %36
32: ; preds = %18
%33 = add nsw i32 %20, 1
%34 = fcmp ogt double %30, %21
br i1 %34, label %35, label %36
35: ; preds = %32
br label %36
36: ; preds = %18, %35, %32
%37 = phi double [ %30, %35 ], [ %21, %32 ], [ %21, %18 ]
%38 = phi i32 [ %33, %35 ], [ %33, %32 ], [ %20, %18 ]
%39 = or i64 %19, 1
%40 = getelementptr inbounds double, double* %16, i64 %39
%41 = load double, double* %40
%42 = getelementptr inbounds double, double* %17, i64 %39
%43 = load double, double* %42
%44 = fsub double %41, %43
%45 = fcmp olt double %44, 0.000000e+00
%46 = fneg double %44
%47 = select i1 %45, double %46, double %44
%48 = fcmp ogt double %47, 0x3E7AD7F29ABCAF48
br i1 %48, label %49, label %53
49: ; preds = %36
%50 = add nsw i32 %38, 1
%51 = fcmp ogt double %47, %37
br i1 %51, label %52, label %53
52: ; preds = %49
br label %53
53: ; preds = %52, %49, %36
%54 = phi double [ %47, %52 ], [ %37, %49 ], [ %37, %36 ]
%55 = phi i32 [ %50, %52 ], [ %50, %49 ], [ %38, %36 ]
%56 = add nuw nsw i64 %19, 2
%57 = add i64 %22, 2
%58 = icmp eq i64 %57, %9
br i1 %58, label %59, label %18
59: ; preds = %53, %11
%60 = phi double [ undef, %11 ], [ %54, %53 ]
%61 = phi i32 [ undef, %11 ], [ %55, %53 ]
%62 = phi i64 [ 0, %11 ], [ %56, %53 ]
%63 = phi i32 [ %13, %11 ], [ %55, %53 ]
%64 = phi double [ %14, %11 ], [ %54, %53 ]
br i1 %10, label %79, label %65
65: ; preds = %59
%66 = getelementptr inbounds double, double* %16, i64 %62
%67 = load double, double* %66
%68 = getelementptr inbounds double, double* %17, i64 %62
%69 = load double, double* %68
%70 = fsub double %67, %69
%71 = fcmp olt double %70, 0.000000e+00
%72 = fneg double %70
%73 = select i1 %71, double %72, double %70
%74 = fcmp ogt double %73, 0x3E7AD7F29ABCAF48
br i1 %74, label %75, label %79
75: ; preds = %65
%76 = add nsw i32 %63, 1
%77 = fcmp ogt double %73, %64
br i1 %77, label %78, label %79
78: ; preds = %75
br label %79
79: ; preds = %65, %75, %78, %59
%80 = phi double [ %60, %59 ], [ %73, %78 ], [ %64, %75 ], [ %64, %65 ]
%81 = phi i32 [ %61, %59 ], [ %76, %78 ], [ %76, %75 ], [ %63, %65 ]
%82 = add nuw nsw i64 %12, 1
%83 = icmp eq i64 %82, %4
br i1 %83, label %84, label %11
84: ; preds = %79
%85 = icmp sgt i32 %81, 0
br i1 %85, label %86, label %88
86: ; preds = %84
%87 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([49 x i8], [49 x i8]* @.str, i64 0, i64 0), i32 %81, double 0x3E7AD7F29ABCAF48, double %80)
br label %90
88: ; preds = %3, %84
%89 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([25 x i8], [25 x i8]* @str, i64 0, i64 0))
br label %90
90: ; preds = %88, %86
ret void
}