#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 max_abs_value(double **M, int nx, int ny) {
  int i, j;
  double tmp;
  double max = -1.0e19;
  for (j = 1; j <= ny; j++) {
    for (i = 1; i <= nx; i++) {
      tmp = fabs(M[i][j]);
      if (tmp > max)
        max = tmp;
    }
  }
  return max;
}


x86 assembly (clang, -O0):
<max_abs_value>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movsd  0xeea(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x1,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     11a6 <max_abs_value+0xa6>
movl   $0x1,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1193 <max_abs_value+0x93>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x18(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xeb3(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    1180 <max_abs_value+0x80>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    1185 <max_abs_value+0x85>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1135 <max_abs_value+0x35>
jmp    1198 <max_abs_value+0x98>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1122 <max_abs_value+0x22>
movsd  -0x28(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<max_abs_value>:
push   %rbx
test   %edx,%edx
jle    11b7 <max_abs_value+0xb7>
lea    0x1(%rsi),%r8d
add    $0x1,%edx
lea    -0x1(%r8),%r9
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
neg    %r10
movsd  0xeda(%rip),%xmm0
mov    $0x1,%r11d
movapd 0xedc(%rip),%xmm1
movapd 0xed4(%rip),%xmm2
jmp    1149 <max_abs_value+0x49>
xchg   %ax,%ax
add    $0x1,%r11
cmp    %rdx,%r11
je     11b5 <max_abs_value+0xb5>
test   %esi,%esi
jle    1140 <max_abs_value+0x40>
mov    $0x1,%ecx
cmp    $0x2,%r8d
je     1197 <max_abs_value+0x97>
mov    $0x1,%ecx
nopl   (%rax)
mov    (%rdi,%rcx,8),%rax
mov    0x8(%rdi,%rcx,8),%rbx
movsd  (%rax,%r11,8),%xmm3
andpd  %xmm2,%xmm3
maxsd  %xmm0,%xmm3
movsd  (%rbx,%r11,8),%xmm0
andpd  %xmm2,%xmm0
maxsd  %xmm3,%xmm0
lea    (%r10,%rcx,1),%rax
add    $0x2,%rax
add    $0x2,%rcx
cmp    $0x1,%rax
jne    1160 <max_abs_value+0x60>
test   $0x1,%r9b
je     1140 <max_abs_value+0x40>
mov    (%rdi,%rcx,8),%rax
movsd  (%rax,%r11,8),%xmm3
andpd  %xmm1,%xmm3
maxsd  %xmm0,%xmm3
movapd %xmm3,%xmm0
jmp    1140 <max_abs_value+0x40>
pop    %rbx
ret
movsd  0xe41(%rip),%xmm0
pop    %rbx
ret

LLVM IR (-O0):
define double @max_abs_value(double** %0, i32 %1, i32 %2) {
%4 = alloca double**
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca double
%10 = alloca double
store double** %0, double*** %4
store i32 %1, i32* %5
store i32 %2, i32* %6
store double -1.000000e+19, double* %10
store i32 1, i32* %8
br label %11
11: ; preds = %41, %3
%12 = load i32, i32* %8
%13 = load i32, i32* %6
%14 = icmp sle i32 %12, %13
br i1 %14, label %15, label %44
15: ; preds = %11
store i32 1, i32* %7
br label %16
16: ; preds = %37, %15
%17 = load i32, i32* %7
%18 = load i32, i32* %5
%19 = icmp sle i32 %17, %18
br i1 %19, label %20, label %40
20: ; preds = %16
%21 = load double**, double*** %4
%22 = load i32, i32* %7
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double*, double** %21, i64 %23
%25 = load double*, double** %24
%26 = load i32, i32* %8
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
%29 = load double, double* %28
%30 = call double @llvm.fabs.f64(double %29)
store double %30, double* %9
%31 = load double, double* %9
%32 = load double, double* %10
%33 = fcmp ogt double %31, %32
br i1 %33, label %34, label %36
34: ; preds = %20
%35 = load double, double* %9
store double %35, double* %10
br label %36
36: ; preds = %34, %20
br label %37
37: ; preds = %36
%38 = load i32, i32* %7
%39 = add nsw i32 %38, 1
store i32 %39, i32* %7
br label %16
40: ; preds = %16
br label %41
41: ; preds = %40
%42 = load i32, i32* %8
%43 = add nsw i32 %42, 1
store i32 %43, i32* %8
br label %11
44: ; preds = %11
%45 = load double, double* %10
ret double %45
}

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