#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, float wref[][n], float w[][n]) {
  float 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    $0x40,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x4(%rbp),%eax
mov    %rax,-0x40(%rbp)
mov    -0x4(%rbp),%eax
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    123e <compare+0x12e>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    122b <compare+0x11b>
mov    -0x38(%rbp),%rdx
mov    -0x40(%rbp),%rsi
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
imul   %rsi,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
imul   %rdx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x20(%rbp)
xorps  %xmm0,%xmm0
ucomiss -0x20(%rbp),%xmm0
jbe    11d6 <compare+0xc6>
movss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe3f(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe19(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1218 <compare+0x108>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
movss  -0x20(%rbp),%xmm0
ucomiss -0x1c(%rbp),%xmm0
jbe    1213 <compare+0x103>
movss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x1c(%rbp)
jmp    1218 <compare+0x108>
jmp    121d <compare+0x10d>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    115a <compare+0x4a>
jmp    1230 <compare+0x120>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1147 <compare+0x37>
cmpl   $0x0,-0x24(%rbp)
jle    126f <compare+0x15f>
mov    -0x24(%rbp),%esi
movss  -0x1c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
lea    0xdb5(%rip),%rdi
movsd  0xd9d(%rip),%xmm0
mov    $0x2,%al
call   1030 <printf@plt>
jmp    127d <compare+0x16d>
lea    0xdcb(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<compare>:
test   %edi,%edi
jle    11d7 <compare+0xb7>
mov    %edi,%ecx
lea    0x0(,%rcx,4),%r8
xorps  %xmm0,%xmm0
xor    %r9d,%r9d
movaps 0xec1(%rip),%xmm1
movsd  0xec9(%rip),%xmm2
xor    %eax,%eax
xorps  %xmm3,%xmm3
jmp    115f <compare+0x3f>
xchg   %ax,%ax
add    $0x1,%r9
add    %r8,%rdx
add    %r8,%rsi
cmp    %rcx,%r9
je     11b4 <compare+0x94>
xor    %edi,%edi
jmp    1179 <compare+0x59>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rdi
cmp    %rdi,%rcx
je     1150 <compare+0x30>
movss  (%rsi,%rdi,4),%xmm5
subss  (%rdx,%rdi,4),%xmm5
movaps %xmm5,%xmm6
xorps  %xmm1,%xmm6
movaps %xmm5,%xmm4
cmpltss %xmm0,%xmm4
andps  %xmm4,%xmm6
andnps %xmm5,%xmm4
orps   %xmm6,%xmm4
xorps  %xmm5,%xmm5
cvtss2sd %xmm4,%xmm5
ucomisd %xmm2,%xmm5
jbe    1170 <compare+0x50>
add    $0x1,%eax
ucomiss %xmm3,%xmm4
jbe    1170 <compare+0x50>
movaps %xmm4,%xmm3
jmp    1170 <compare+0x50>
test   %eax,%eax
jle    11d7 <compare+0xb7>
xorps  %xmm1,%xmm1
cvtss2sd %xmm3,%xmm1
lea    0xe52(%rip),%rdi
movsd  0xe42(%rip),%xmm0
mov    %eax,%esi
mov    $0x2,%al
jmp    1040 <printf@plt>
lea    0xe6b(%rip),%rdi
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, float* %1, float* %2) {
%4 = alloca i32
%5 = alloca float*
%6 = alloca float*
%7 = alloca float
%8 = alloca float
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
store i32 %0, i32* %4
store float* %1, float** %5
store float* %2, float** %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 float 0.000000e+00, float* %7
store i32 0, i32* %10
br label %16
16: ; preds = %70, %3
%17 = load i32, i32* %10
%18 = load i32, i32* %4
%19 = icmp slt i32 %17, %18
br i1 %19, label %20, label %73
20: ; preds = %16
store i32 0, i32* %11
br label %21
21: ; preds = %66, %20
%22 = load i32, i32* %11
%23 = load i32, i32* %4
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %69
25: ; preds = %21
%26 = load float*, float** %5
%27 = load i32, i32* %10
%28 = sext i32 %27 to i64
%29 = mul nsw i64 %28, %13
%30 = getelementptr inbounds float, float* %26, i64 %29
%31 = load i32, i32* %11
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds float, float* %30, i64 %32
%34 = load float, float* %33
%35 = load float*, float** %6
%36 = load i32, i32* %10
%37 = sext i32 %36 to i64
%38 = mul nsw i64 %37, %15
%39 = getelementptr inbounds float, float* %35, i64 %38
%40 = load i32, i32* %11
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds float, float* %39, i64 %41
%43 = load float, float* %42
%44 = fsub float %34, %43
store float %44, float* %8
%45 = load float, float* %8
%46 = fcmp olt float %45, 0.000000e+00
br i1 %46, label %47, label %52
47: ; preds = %25
%48 = load float, float* %8
%49 = fpext float %48 to double
%50 = fmul double -1.000000e+00, %49
%51 = fptrunc double %50 to float
store float %51, float* %8
br label %52
52: ; preds = %47, %25
%53 = load float, float* %8
%54 = fpext float %53 to double
%55 = fcmp ogt double %54, 0x3E7AD7F29ABCAF48
br i1 %55, label %56, label %65
56: ; preds = %52
%57 = load i32, i32* %9
%58 = add nsw i32 %57, 1
store i32 %58, i32* %9
%59 = load float, float* %8
%60 = load float, float* %7
%61 = fcmp ogt float %59, %60
br i1 %61, label %62, label %64
62: ; preds = %56
%63 = load float, float* %8
store float %63, float* %7
br label %64
64: ; preds = %62, %56
br label %65
65: ; preds = %64, %52
br label %66
66: ; preds = %65
%67 = load i32, i32* %11
%68 = add nsw i32 %67, 1
store i32 %68, i32* %11
br label %21
69: ; preds = %21
br label %70
70: ; preds = %69
%71 = load i32, i32* %10
%72 = add nsw i32 %71, 1
store i32 %72, i32* %10
br label %16
73: ; preds = %16
%74 = load i32, i32* %9
%75 = icmp sgt i32 %74, 0
br i1 %75, label %76, label %81
76: ; preds = %73
%77 = load i32, i32* %9
%78 = load float, float* %7
%79 = fpext float %78 to double
%80 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([49 x i8], [49 x i8]* @.str, i64 0, i64 0), i32 %77, double 0x3E7AD7F29ABCAF48, double %79)
br label %83
81: ; preds = %73
%82 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([26 x i8], [26 x i8]* @.str.1, i64 0, i64 0))
br label %83
83: ; preds = %81, %76
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, float* nocapture %1, float* nocapture %2) {
%4 = zext i32 %0 to i64
%5 = icmp sgt i32 %0, 0
br i1 %5, label %6, label %44
6: ; preds = %3, %36
%7 = phi i64 [ %37, %36 ], [ 0, %3 ]
%8 = phi i32 [ %33, %36 ], [ 0, %3 ]
%9 = phi float [ %32, %36 ], [ 0.000000e+00, %3 ]
%10 = mul nuw nsw i64 %7, %4
%11 = getelementptr inbounds float, float* %1, i64 %10
%12 = getelementptr inbounds float, float* %2, i64 %10
br label %13
13: ; preds = %6, %31
%14 = phi i64 [ 0, %6 ], [ %34, %31 ]
%15 = phi i32 [ %8, %6 ], [ %33, %31 ]
%16 = phi float [ %9, %6 ], [ %32, %31 ]
%17 = getelementptr inbounds float, float* %11, i64 %14
%18 = load float, float* %17
%19 = getelementptr inbounds float, float* %12, i64 %14
%20 = load float, float* %19
%21 = fsub float %18, %20
%22 = fcmp olt float %21, 0.000000e+00
%23 = fneg float %21
%24 = select i1 %22, float %23, float %21
%25 = fpext float %24 to double
%26 = fcmp ogt double %25, 0x3E7AD7F29ABCAF48
br i1 %26, label %27, label %31
27: ; preds = %13
%28 = add nsw i32 %15, 1
%29 = fcmp ogt float %24, %16
br i1 %29, label %30, label %31
30: ; preds = %27
br label %31
31: ; preds = %13, %30, %27
%32 = phi float [ %24, %30 ], [ %16, %27 ], [ %16, %13 ]
%33 = phi i32 [ %28, %30 ], [ %28, %27 ], [ %15, %13 ]
%34 = add nuw nsw i64 %14, 1
%35 = icmp eq i64 %34, %4
br i1 %35, label %36, label %13
36: ; preds = %31
%37 = add nuw nsw i64 %7, 1
%38 = icmp eq i64 %37, %4
br i1 %38, label %39, label %6
39: ; preds = %36
%40 = icmp sgt i32 %33, 0
br i1 %40, label %41, label %44
41: ; preds = %39
%42 = fpext float %32 to double
%43 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([49 x i8], [49 x i8]* @.str, i64 0, i64 0), i32 %33, double 0x3E7AD7F29ABCAF48, double %42)
br label %46
44: ; preds = %3, %39
%45 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([25 x i8], [25 x i8]* @str, i64 0, i64 0))
br label %46
46: ; preds = %44, %41
ret void
}