#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);
double compute_error(double solution[][500 + 2], double u[][500 + 2],
                     const int m) {
  double error = 0.0;
  int i, j;
  for (i = 1; i < m + 1; ++i) {
    for (j = 1; j < m + 1; ++j) {
      const double abs_diff = fabs(solution[i][j] - u[i][j]);
      if (error < abs_diff) {
        error = abs_diff;
      }
    }
  }
  return error;
}


x86 assembly (clang, -O0):
<compute_error>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    11d1 <compute_error+0xd1>
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    11be <compute_error+0xbe>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0xfb0,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
imul   $0xfb0,%rcx,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movaps 0xe7c(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x20(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11ab <compute_error+0xab>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    11b0 <compute_error+0xb0>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1136 <compute_error+0x36>
jmp    11c3 <compute_error+0xc3>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111e <compute_error+0x1e>
movsd  -0x20(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<compute_error>:
push   %r14
push   %rbx
test   %edx,%edx
jle    11f3 <compute_error+0xf3>
mov    %edx,%r8d
lea    0x1(%r8),%r9
mov    %r8d,%r11d
and    $0xfffffffe,%r11d
lea    0xfc0(%rsi),%rax
lea    0xfc0(%rdi),%rcx
xorpd  %xmm0,%xmm0
mov    $0x1,%r10d
movapd 0xec7(%rip),%xmm1
movapd 0xebf(%rip),%xmm2
jmp    116a <compute_error+0x6a>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r10
add    $0xfb0,%rax
add    $0xfb0,%rcx
cmp    %r9,%r10
je     11f7 <compute_error+0xf7>
cmp    $0x2,%r9d
jne    1180 <compute_error+0x80>
mov    $0x1,%edx
test   $0x1,%r8b
je     1150 <compute_error+0x50>
jmp    11c9 <compute_error+0xc9>
nopl   (%rax)
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  -0x8(%rcx,%rdx,8),%xmm3
subsd  -0x8(%rax,%rdx,8),%xmm3
andpd  %xmm2,%xmm3
maxsd  %xmm0,%xmm3
movsd  (%rcx,%rdx,8),%xmm0
subsd  (%rax,%rdx,8),%xmm0
andpd  %xmm2,%xmm0
maxsd  %xmm3,%xmm0
add    $0x2,%rdx
cmp    %rdx,%r11
jne    1190 <compute_error+0x90>
add    $0x1,%rdx
test   $0x1,%r8b
je     1150 <compute_error+0x50>
imul   $0xfb0,%r10,%r14
lea    (%rdi,%r14,1),%rbx
movsd  (%rbx,%rdx,8),%xmm3
add    %rsi,%r14
subsd  (%r14,%rdx,8),%xmm3
andpd  %xmm1,%xmm3
maxsd  %xmm0,%xmm3
movapd %xmm3,%xmm0
jmp    1150 <compute_error+0x50>
xorpd  %xmm0,%xmm0
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define double @compute_error([502 x double]* %0, [502 x double]* %1, i32 %2) {
%4 = alloca [502 x double]*
%5 = alloca [502 x double]*
%6 = alloca i32
%7 = alloca double
%8 = alloca i32
%9 = alloca i32
%10 = alloca double
store [502 x double]* %0, [502 x double]** %4
store [502 x double]* %1, [502 x double]** %5
store i32 %2, i32* %6
store double 0.000000e+00, double* %7
store i32 1, i32* %8
br label %11
11: ; preds = %51, %3
%12 = load i32, i32* %8
%13 = load i32, i32* %6
%14 = add nsw i32 %13, 1
%15 = icmp slt i32 %12, %14
br i1 %15, label %16, label %54
16: ; preds = %11
store i32 1, i32* %9
br label %17
17: ; preds = %47, %16
%18 = load i32, i32* %9
%19 = load i32, i32* %6
%20 = add nsw i32 %19, 1
%21 = icmp slt i32 %18, %20
br i1 %21, label %22, label %50
22: ; preds = %17
%23 = load [502 x double]*, [502 x double]** %4
%24 = load i32, i32* %8
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds [502 x double], [502 x double]* %23, i64 %25
%27 = load i32, i32* %9
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds [502 x double], [502 x double]* %26, i64 0, i64 %28
%30 = load double, double* %29
%31 = load [502 x double]*, [502 x double]** %5
%32 = load i32, i32* %8
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds [502 x double], [502 x double]* %31, i64 %33
%35 = load i32, i32* %9
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds [502 x double], [502 x double]* %34, i64 0, i64 %36
%38 = load double, double* %37
%39 = fsub double %30, %38
%40 = call double @llvm.fabs.f64(double %39)
store double %40, double* %10
%41 = load double, double* %7
%42 = load double, double* %10
%43 = fcmp olt double %41, %42
br i1 %43, label %44, label %46
44: ; preds = %22
%45 = load double, double* %10
store double %45, double* %7
br label %46
46: ; preds = %44, %22
br label %47
47: ; preds = %46
%48 = load i32, i32* %9
%49 = add nsw i32 %48, 1
store i32 %49, i32* %9
br label %17
50: ; preds = %17
br label %51
51: ; preds = %50
%52 = load i32, i32* %8
%53 = add nsw i32 %52, 1
store i32 %53, i32* %8
br label %11
54: ; preds = %11
%55 = load double, double* %7
ret double %55
}

LLVM IR (-O2):
define double @compute_error([502 x double]* nocapture %0, [502 x double]* nocapture %1, i32 %2) {
%4 = icmp slt i32 %2, 1
br i1 %4, label %57, label %5
5: ; preds = %3
%6 = add nuw i32 %2, 1
%7 = zext i32 %6 to i64
%8 = add nsw i64 %7, -1
%9 = and i64 %8, 1
%10 = icmp eq i32 %6, 2
%11 = and i64 %8, -2
%12 = icmp eq i64 %9, 0
br label %13
13: ; preds = %53, %5
%14 = phi i64 [ 1, %5 ], [ %55, %53 ]
%15 = phi double [ 0.000000e+00, %5 ], [ %54, %53 ]
br i1 %10, label %40, label %16
16: ; preds = %13, %16
%17 = phi i64 [ %37, %16 ], [ 1, %13 ]
%18 = phi double [ %36, %16 ], [ %15, %13 ]
%19 = phi i64 [ %38, %16 ], [ 0, %13 ]
%20 = getelementptr inbounds [502 x double], [502 x double]* %0, i64 %14, i64 %17
%21 = load double, double* %20
%22 = getelementptr inbounds [502 x double], [502 x double]* %1, i64 %14, i64 %17
%23 = load double, double* %22
%24 = fsub double %21, %23
%25 = tail call double @llvm.fabs.f64(double %24)
%26 = fcmp olt double %18, %25
%27 = select i1 %26, double %25, double %18
%28 = add nuw nsw i64 %17, 1
%29 = getelementptr inbounds [502 x double], [502 x double]* %0, i64 %14, i64 %28
%30 = load double, double* %29
%31 = getelementptr inbounds [502 x double], [502 x double]* %1, i64 %14, i64 %28
%32 = load double, double* %31
%33 = fsub double %30, %32
%34 = tail call double @llvm.fabs.f64(double %33)
%35 = fcmp olt double %27, %34
%36 = select i1 %35, double %34, double %27
%37 = add nuw nsw i64 %17, 2
%38 = add i64 %19, 2
%39 = icmp eq i64 %38, %11
br i1 %39, label %40, label %16
40: ; preds = %16, %13
%41 = phi double [ undef, %13 ], [ %36, %16 ]
%42 = phi i64 [ 1, %13 ], [ %37, %16 ]
%43 = phi double [ %15, %13 ], [ %36, %16 ]
br i1 %12, label %53, label %44
44: ; preds = %40
%45 = getelementptr inbounds [502 x double], [502 x double]* %0, i64 %14, i64 %42
%46 = load double, double* %45
%47 = getelementptr inbounds [502 x double], [502 x double]* %1, i64 %14, i64 %42
%48 = load double, double* %47
%49 = fsub double %46, %48
%50 = tail call double @llvm.fabs.f64(double %49)
%51 = fcmp olt double %43, %50
%52 = select i1 %51, double %50, double %43
br label %53
53: ; preds = %40, %44
%54 = phi double [ %41, %40 ], [ %52, %44 ]
%55 = add nuw nsw i64 %14, 1
%56 = icmp eq i64 %55, %7
br i1 %56, label %57, label %13
57: ; preds = %53, %3
%58 = phi double [ 0.000000e+00, %3 ], [ %54, %53 ]
ret double %58
}