#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 */
scalar_t__ malloc(size_t);
double **matrix(int nrl, int nrh, int ncl, int nch) {
  int i;
  int nrow = nrh - nrl + 1;
  int ncol = nch - ncl + 1;
  double **pArray = (double **)malloc((size_t)(nrow * sizeof(double *)));
  double *pMatrix = (double *)malloc((size_t)(nrow * ncol * sizeof(double)));
  pArray[0] = pMatrix - ncl;
  for (i = 1; i < nrow; i++) {
    pArray[i] = pArray[i - 1] + ncol;
  }
  return pArray - nrl;
}


x86 assembly (clang, -O0):
<matrix>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    -0x8(%rbp),%eax
sub    -0x4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x10(%rbp),%eax
sub    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
movslq -0x18(%rbp),%rdi
shl    $0x3,%rdi
call   1030 <malloc@plt>
mov    %rax,-0x28(%rbp)
mov    -0x18(%rbp),%eax
imul   -0x1c(%rbp),%eax
movslq %eax,%rdi
shl    $0x3,%rdi
call   1030 <malloc@plt>
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rcx
movslq -0xc(%rbp),%rdx
xor    %eax,%eax
sub    %rdx,%rax
shl    $0x3,%rax
add    %rax,%rcx
mov    -0x28(%rbp),%rax
mov    %rcx,(%rax)
movl   $0x1,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11c8 <matrix+0xb8>
mov    -0x28(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rdx
movslq -0x1c(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdx
mov    -0x28(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1186 <matrix+0x76>
mov    -0x28(%rbp),%rax
movslq -0x4(%rbp),%rdx
xor    %ecx,%ecx
sub    %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<matrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %ecx,%r15d
mov    %edx,%r12d
mov    %esi,%ebp
mov    %edi,%r14d
sub    %edi,%ebp
sub    %edx,%r15d
add    $0x1,%r15d
movslq %ebp,%r13
add    $0x1,%ebp
lea    0x8(,%r13,8),%rdi
call   1030 <malloc@plt>
mov    %rax,%rbx
mov    %r15d,%eax
imul   %ebp,%eax
movslq %eax,%rdi
shl    $0x3,%rdi
call   1030 <malloc@plt>
movslq %r12d,%rcx
shl    $0x3,%rcx
sub    %rcx,%rax
mov    %rax,(%rbx)
test   %r13d,%r13d
jle    1220 <matrix+0x110>
movslq %r15d,%r9
lea    -0x1(%rbp),%rdi
add    $0xfffffffffffffffe,%rbp
mov    %edi,%edx
and    $0x7,%edx
mov    $0x1,%esi
cmp    $0x7,%rbp
jb     11f6 <matrix+0xe6>
and    $0xfffffffffffffff8,%rdi
mov    %r9,%r8
shl    $0x6,%r8
lea    0x0(,%r9,8),%rbp
xor    %esi,%esi
nopl   (%rax)
mov    %rax,%rcx
add    %rbp,%rax
mov    %rax,0x8(%rbx,%rsi,8)
add    %rbp,%rax
mov    %rax,0x10(%rbx,%rsi,8)
add    %rbp,%rax
mov    %rax,0x18(%rbx,%rsi,8)
add    %rbp,%rax
mov    %rax,0x20(%rbx,%rsi,8)
add    %rbp,%rax
mov    %rax,0x28(%rbx,%rsi,8)
add    %rbp,%rax
mov    %rax,0x30(%rbx,%rsi,8)
add    %rbp,%rax
mov    %rax,0x38(%rbx,%rsi,8)
add    %rbp,%rax
mov    %rax,0x40(%rbx,%rsi,8)
add    $0x8,%rsi
cmp    %rsi,%rdi
jne    11a0 <matrix+0x90>
add    %r8,%rcx
add    $0x1,%rsi
mov    %rcx,%rax
test   %rdx,%rdx
je     1220 <matrix+0x110>
lea    (%rax,%r9,8),%rax
shl    $0x3,%r9
lea    (%rbx,%rsi,8),%rcx
xor    %esi,%esi
nopl   0x0(%rax)
mov    %rax,(%rcx,%rsi,8)
add    $0x1,%rsi
add    %r9,%rax
cmp    %rsi,%rdx
jne    1210 <matrix+0x100>
movslq %r14d,%rax
shl    $0x3,%rax
sub    %rax,%rbx
mov    %rbx,%rax
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double** @matrix(i32 %0, i32 %1, i32 %2, i32 %3) {
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca double**
%13 = alloca double*
store i32 %0, i32* %5
store i32 %1, i32* %6
store i32 %2, i32* %7
store i32 %3, i32* %8
%14 = load i32, i32* %6
%15 = load i32, i32* %5
%16 = sub nsw i32 %14, %15
%17 = add nsw i32 %16, 1
store i32 %17, i32* %10
%18 = load i32, i32* %8
%19 = load i32, i32* %7
%20 = sub nsw i32 %18, %19
%21 = add nsw i32 %20, 1
store i32 %21, i32* %11
%22 = load i32, i32* %10
%23 = sext i32 %22 to i64
%24 = mul i64 %23, 8
%25 = call i64 @malloc(i64 %24)
%26 = inttoptr i64 %25 to double**
store double** %26, double*** %12
%27 = load i32, i32* %10
%28 = load i32, i32* %11
%29 = mul nsw i32 %27, %28
%30 = sext i32 %29 to i64
%31 = mul i64 %30, 8
%32 = call i64 @malloc(i64 %31)
%33 = inttoptr i64 %32 to double*
store double* %33, double** %13
%34 = load double*, double** %13
%35 = load i32, i32* %7
%36 = sext i32 %35 to i64
%37 = sub i64 0, %36
%38 = getelementptr inbounds double, double* %34, i64 %37
%39 = load double**, double*** %12
%40 = getelementptr inbounds double*, double** %39, i64 0
store double* %38, double** %40
store i32 1, i32* %9
br label %41
41: ; preds = %59, %4
%42 = load i32, i32* %9
%43 = load i32, i32* %10
%44 = icmp slt i32 %42, %43
br i1 %44, label %45, label %62
45: ; preds = %41
%46 = load double**, double*** %12
%47 = load i32, i32* %9
%48 = sub nsw i32 %47, 1
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double*, double** %46, i64 %49
%51 = load double*, double** %50
%52 = load i32, i32* %11
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds double, double* %51, i64 %53
%55 = load double**, double*** %12
%56 = load i32, i32* %9
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds double*, double** %55, i64 %57
store double* %54, double** %58
br label %59
59: ; preds = %45
%60 = load i32, i32* %9
%61 = add nsw i32 %60, 1
store i32 %61, i32* %9
br label %41
62: ; preds = %41
%63 = load double**, double*** %12
%64 = load i32, i32* %5
%65 = sext i32 %64 to i64
%66 = sub i64 0, %65
%67 = getelementptr inbounds double*, double** %63, i64 %66
ret double** %67
}

LLVM IR (-O2):
define double** @matrix(i32 %0, i32 %1, i32 %2, i32 %3) {
%5 = sub nsw i32 %1, %0
%6 = add nsw i32 %5, 1
%7 = sub nsw i32 %3, %2
%8 = add nsw i32 %7, 1
%9 = sext i32 %6 to i64
%10 = shl nsw i64 %9, 3
%11 = tail call i64 @malloc(i64 %10)
%12 = inttoptr i64 %11 to double**
%13 = mul nsw i32 %8, %6
%14 = sext i32 %13 to i64
%15 = shl nsw i64 %14, 3
%16 = tail call i64 @malloc(i64 %15)
%17 = inttoptr i64 %16 to double*
%18 = sext i32 %2 to i64
%19 = sub nsw i64 0, %18
%20 = getelementptr inbounds double, double* %17, i64 %19
store double* %20, double** %12
%21 = icmp slt i32 %5, 1
br i1 %21, label %76, label %22
22: ; preds = %4
%23 = sext i32 %8 to i64
%24 = add i32 %1, 1
%25 = sub i32 %24, %0
%26 = zext i32 %25 to i64
%27 = add nsw i64 %26, -1
%28 = add nsw i64 %26, -2
%29 = and i64 %27, 7
%30 = icmp ult i64 %28, 7
br i1 %30, label %63, label %31
31: ; preds = %22
%32 = and i64 %27, -8
br label %33
33: ; preds = %33, %31
%34 = phi double* [ %20, %31 ], [ %58, %33 ]
%35 = phi i64 [ 1, %31 ], [ %60, %33 ]
%36 = phi i64 [ 0, %31 ], [ %61, %33 ]
%37 = getelementptr inbounds double, double* %34, i64 %23
%38 = getelementptr inbounds double*, double** %12, i64 %35
store double* %37, double** %38
%39 = add nuw nsw i64 %35, 1
%40 = getelementptr inbounds double, double* %37, i64 %23
%41 = getelementptr inbounds double*, double** %12, i64 %39
store double* %40, double** %41
%42 = add nuw nsw i64 %35, 2
%43 = getelementptr inbounds double, double* %40, i64 %23
%44 = getelementptr inbounds double*, double** %12, i64 %42
store double* %43, double** %44
%45 = add nuw nsw i64 %35, 3
%46 = getelementptr inbounds double, double* %43, i64 %23
%47 = getelementptr inbounds double*, double** %12, i64 %45
store double* %46, double** %47
%48 = add nuw nsw i64 %35, 4
%49 = getelementptr inbounds double, double* %46, i64 %23
%50 = getelementptr inbounds double*, double** %12, i64 %48
store double* %49, double** %50
%51 = add nuw nsw i64 %35, 5
%52 = getelementptr inbounds double, double* %49, i64 %23
%53 = getelementptr inbounds double*, double** %12, i64 %51
store double* %52, double** %53
%54 = add nuw nsw i64 %35, 6
%55 = getelementptr inbounds double, double* %52, i64 %23
%56 = getelementptr inbounds double*, double** %12, i64 %54
store double* %55, double** %56
%57 = add nuw nsw i64 %35, 7
%58 = getelementptr inbounds double, double* %55, i64 %23
%59 = getelementptr inbounds double*, double** %12, i64 %57
store double* %58, double** %59
%60 = add nuw nsw i64 %35, 8
%61 = add i64 %36, 8
%62 = icmp eq i64 %61, %32
br i1 %62, label %63, label %33
63: ; preds = %33, %22
%64 = phi double* [ %20, %22 ], [ %58, %33 ]
%65 = phi i64 [ 1, %22 ], [ %60, %33 ]
%66 = icmp eq i64 %29, 0
br i1 %66, label %76, label %67
67: ; preds = %63, %67
%68 = phi double* [ %71, %67 ], [ %64, %63 ]
%69 = phi i64 [ %73, %67 ], [ %65, %63 ]
%70 = phi i64 [ %74, %67 ], [ 0, %63 ]
%71 = getelementptr inbounds double, double* %68, i64 %23
%72 = getelementptr inbounds double*, double** %12, i64 %69
store double* %71, double** %72
%73 = add nuw nsw i64 %69, 1
%74 = add i64 %70, 1
%75 = icmp eq i64 %74, %29
br i1 %75, label %76, label %67
76: ; preds = %63, %67, %4
%77 = sext i32 %0 to i64
%78 = sub nsw i64 0, %77
%79 = getelementptr inbounds double*, double** %12, i64 %78
ret double** %79
}