#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 errorcheck3d(int n1, int n2, int n3, double scale, double ***a) {
  int j1, j2, j3, seed = 0;
  double err = 0, e;
  for (j1 = 0; j1 <= n1 - 1; j1++) {
    for (j2 = 0; j2 <= n2 - 1; j2++) {
      for (j3 = 0; j3 <= n3 - 1; j3++) {
        e = ((*(&seed) = (*(&seed) * 7141 + 54773) % 259200) * (1.0 / 259200)) -
            a[j1][j2][j3] * scale;
        err = ((err) > (fabs(e)) ? (err) : (fabs(e)));
      }
    }
  }
  return err;
}


x86 assembly (clang, -O0):
<errorcheck3d>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1257 <errorcheck3d+0x157>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1244 <errorcheck3d+0x144>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1231 <errorcheck3d+0x131>
imul   $0x1be5,-0x30(%rbp),%eax
add    $0xd5f5,%eax
mov    $0x3f480,%ecx
cltd
idiv   %ecx
mov    %edx,-0x30(%rbp)
cvtsi2sd %edx,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x18(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
movsd  0xe47(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm1
movaps 0xe19(%rip),%xmm2
pand   %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jbe    1204 <errorcheck3d+0x104>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
jmp    1219 <errorcheck3d+0x119>
movsd  -0x40(%rbp),%xmm0
movaps 0xdf0(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    115c <errorcheck3d+0x5c>
jmp    1236 <errorcheck3d+0x136>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1144 <errorcheck3d+0x44>
jmp    1249 <errorcheck3d+0x149>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    112c <errorcheck3d+0x2c>
movsd  -0x38(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<errorcheck3d>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
test   %edi,%edi
jle    11e8 <errorcheck3d+0xe8>
mov    %edi,%r8d
mov    %esi,%r9d
mov    %edx,%r14d
xor    %r10d,%r10d
xorpd  %xmm1,%xmm1
movsd  0xed8(%rip),%xmm2
movapd 0xee0(%rip),%xmm3
xor    %edi,%edi
jmp    114d <errorcheck3d+0x4d>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%r10
cmp    %r8,%r10
je     11ec <errorcheck3d+0xec>
test   %esi,%esi
jle    1140 <errorcheck3d+0x40>
xor    %r11d,%r11d
jmp    1169 <errorcheck3d+0x69>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r11
cmp    %r9,%r11
je     1140 <errorcheck3d+0x40>
test   %edx,%edx
jle    1160 <errorcheck3d+0x60>
mov    (%rcx,%r10,8),%rax
mov    (%rax,%r11,8),%r15
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
imul   $0x1be5,%edi,%r12d
mov    %r12d,%edi
add    $0xd5f5,%edi
movslq %edi,%rdi
imul   $0xffffffff81742e05,%rdi,%rbx
shr    $0x20,%rbx
add    %r12d,%ebx
add    $0xd5f5,%ebx
mov    %ebx,%ebp
shr    $0x1f,%ebp
sar    $0x11,%ebx
add    %ebp,%ebx
imul   $0x3f480,%ebx,%ebx
sub    %ebx,%edi
xorps  %xmm4,%xmm4
cvtsi2sd %edi,%xmm4
movsd  (%r15,%rax,8),%xmm5
mulsd  %xmm0,%xmm5
mulsd  %xmm2,%xmm4
subsd  %xmm5,%xmm4
andpd  %xmm3,%xmm4
maxsd  %xmm4,%xmm1
add    $0x1,%rax
cmp    %rax,%r14
jne    1180 <errorcheck3d+0x80>
jmp    1160 <errorcheck3d+0x60>
xorpd  %xmm1,%xmm1
movapd %xmm1,%xmm0
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @errorcheck3d(i32 %0, i32 %1, i32 %2, double %3, double*** %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca double
%10 = alloca double***
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca double
%16 = alloca double
store i32 %0, i32* %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store double %3, double* %9
store double*** %4, double**** %10
store i32 0, i32* %14
store double 0.000000e+00, double* %15
store i32 0, i32* %11
br label %17
17: ; preds = %76, %5
%18 = load i32, i32* %11
%19 = load i32, i32* %6
%20 = sub nsw i32 %19, 1
%21 = icmp sle i32 %18, %20
br i1 %21, label %22, label %79
22: ; preds = %17
store i32 0, i32* %12
br label %23
23: ; preds = %72, %22
%24 = load i32, i32* %12
%25 = load i32, i32* %7
%26 = sub nsw i32 %25, 1
%27 = icmp sle i32 %24, %26
br i1 %27, label %28, label %75
28: ; preds = %23
store i32 0, i32* %13
br label %29
29: ; preds = %68, %28
%30 = load i32, i32* %13
%31 = load i32, i32* %8
%32 = sub nsw i32 %31, 1
%33 = icmp sle i32 %30, %32
br i1 %33, label %34, label %71
34: ; preds = %29
%35 = load i32, i32* %14
%36 = mul nsw i32 %35, 7141
%37 = add nsw i32 %36, 54773
%38 = srem i32 %37, 259200
store i32 %38, i32* %14
%39 = sitofp i32 %38 to double
%40 = load double***, double**** %10
%41 = load i32, i32* %11
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds double**, double*** %40, i64 %42
%44 = load double**, double*** %43
%45 = load i32, i32* %12
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds double*, double** %44, i64 %46
%48 = load double*, double** %47
%49 = load i32, i32* %13
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double, double* %48, i64 %50
%52 = load double, double* %51
%53 = load double, double* %9
%54 = fmul double %52, %53
%55 = fneg double %54
%56 = call double @llvm.fmuladd.f64(double %39, double 0x3ED02E85C0898B71, double %55)
store double %56, double* %16
%57 = load double, double* %15
%58 = load double, double* %16
%59 = call double @llvm.fabs.f64(double %58)
%60 = fcmp ogt double %57, %59
br i1 %60, label %61, label %63
61: ; preds = %34
%62 = load double, double* %15
br label %66
63: ; preds = %34
%64 = load double, double* %16
%65 = call double @llvm.fabs.f64(double %64)
br label %66
66: ; preds = %63, %61
%67 = phi double [ %62, %61 ], [ %65, %63 ]
store double %67, double* %15
br label %68
68: ; preds = %66
%69 = load i32, i32* %13
%70 = add nsw i32 %69, 1
store i32 %70, i32* %13
br label %29
71: ; preds = %29
br label %72
72: ; preds = %71
%73 = load i32, i32* %12
%74 = add nsw i32 %73, 1
store i32 %74, i32* %12
br label %23
75: ; preds = %23
br label %76
76: ; preds = %75
%77 = load i32, i32* %11
%78 = add nsw i32 %77, 1
store i32 %78, i32* %11
br label %17
79: ; preds = %17
%80 = load double, double* %15
ret double %80
}

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