#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);
int /*<<< orphan*/ nrerror(char *);
double **submatrix(double **a, long oldrl, long oldrh, long oldcl, long oldch,
                   long newrl, long newcl) {
  long i, j, nrow = oldrh - oldrl + 1, ncol = oldcl - newcl;
  double **m;
  m = (double **)malloc((size_t)((nrow + 1) * sizeof(double *)));
  if (!m)
    nrerror("allocation failure in submatrix()");
  m += 1;
  m -= newrl;
  for (i = oldrl, j = newrl; i <= oldrh; i++, j++)
    m[j] = a[i] + ncol;
  return m;
}


x86 assembly (clang, -O0):
<submatrix>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x18(%rbp),%rax
sub    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
mov    -0x20(%rbp),%rax
sub    0x10(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x48(%rbp),%rdi
add    $0x1,%rdi
shl    $0x3,%rdi
call   1040 <malloc@plt>
mov    %rax,-0x58(%rbp)
cmpq   $0x0,-0x58(%rbp)
jne    118c <submatrix+0x6c>
lea    0xe79(%rip),%rdi
call   1030 <nrerror@plt>
mov    -0x58(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x58(%rbp)
mov    -0x30(%rbp),%rdx
mov    -0x58(%rbp),%rax
xor    %ecx,%ecx
sub    %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x30(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x18(%rbp),%rax
jg     120e <submatrix+0xee>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
mov    -0x50(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdx
mov    -0x58(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    11c0 <submatrix+0xa0>
mov    -0x58(%rbp),%rax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<submatrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %r9,%rbx
mov    %rcx,%r15
mov    %rdx,%r14
mov    %rsi,%r12
mov    %rdi,%rbp
mov    %rdx,%rax
sub    %rsi,%rax
lea    0x10(,%rax,8),%rdi
call   1040 <malloc@plt>
mov    %rax,%r13
test   %rax,%rax
jne    1161 <submatrix+0x41>
lea    0xea4(%rip),%rdi
call   1030 <nrerror@plt>
lea    0x8(%r13),%rcx
lea    0x0(,%rbx,8),%rdx
mov    %rcx,%rax
sub    %rdx,%rax
mov    %r14,%rdx
sub    %r12,%rdx
jl     1388 <submatrix+0x268>
sub    0x40(%rsp),%r15
lea    0x1(%rdx),%r9
cmp    $0x4,%r9
jb     11c1 <submatrix+0xa1>
lea    0x8(,%r14,8),%rsi
add    %rbp,%rsi
cmp    %rsi,%rcx
jae    12b0 <submatrix+0x190>
lea    0x10(,%rdx,8),%rcx
add    %r13,%rcx
lea    0x0(,%r12,8),%rdx
add    %rbp,%rdx
cmp    %rcx,%rdx
jae    12b0 <submatrix+0x190>
mov    %rbx,%r10
mov    %r14d,%esi
sub    %r12d,%esi
add    $0x1,%esi
mov    %r14,%r8
sub    %r12,%r8
and    $0x3,%rsi
je     1236 <submatrix+0x116>
mov    %rbp,(%rsp)
lea    0x0(,%r12,8),%r9
add    %rbp,%r9
mov    %r10,%rcx
mov    %rbx,%r11
sub    %rbx,%rcx
lea    0x0(,%rcx,8),%rdx
add    %r13,%rdx
neg    %rsi
mov    $0x8,%ebp
xor    %ecx,%ecx
lea    0x0(,%r15,8),%rbx
xchg   %ax,%ax
mov    -0x8(%r9,%rbp,1),%rdi
add    %rbx,%rdi
mov    %rdi,(%rdx,%rbp,1)
add    $0xffffffffffffffff,%rcx
add    $0x8,%rbp
cmp    %rcx,%rsi
jne    1210 <submatrix+0xf0>
sub    %rcx,%r12
sub    %rcx,%r10
mov    %r11,%rbx
mov    (%rsp),%rbp
cmp    $0x3,%r8
jb     1388 <submatrix+0x268>
sub    %r12,%r14
lea    0x18(,%r12,8),%rdx
add    %rbp,%rdx
sub    %rbx,%r10
lea    0x20(,%r10,8),%rcx
add    %r13,%rcx
mov    $0xffffffffffffffff,%rsi
shl    $0x3,%r15
nopw   0x0(%rax,%rax,1)
mov    -0x10(%rdx,%rsi,8),%rdi
add    %r15,%rdi
mov    %rdi,-0x10(%rcx,%rsi,8)
mov    -0x8(%rdx,%rsi,8),%rdi
add    %r15,%rdi
mov    %rdi,-0x8(%rcx,%rsi,8)
mov    (%rdx,%rsi,8),%rdi
add    %r15,%rdi
mov    %rdi,(%rcx,%rsi,8)
mov    0x8(%rdx,%rsi,8),%rdi
add    %r15,%rdi
mov    %rdi,0x8(%rcx,%rsi,8)
add    $0x4,%rsi
cmp    %rsi,%r14
jne    1270 <submatrix+0x150>
jmp    1388 <submatrix+0x268>
mov    %r9,%r10
and    $0xfffffffffffffffc,%r10
lea    -0x4(%r10),%rcx
mov    %rcx,%r8
shr    $0x2,%r8
add    $0x1,%r8
movq   %r15,%xmm0
test   %rcx,%rcx
je     13a2 <submatrix+0x282>
mov    %r8,%rdi
and    $0xfffffffffffffffe,%rdi
lea    0x30(,%r12,8),%rdx
add    %rbp,%rdx
xor    %esi,%esi
pshufd $0x44,%xmm0,%xmm1
psllq  $0x3,%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movdqu -0x30(%rdx,%rsi,8),%xmm2
movdqu -0x20(%rdx,%rsi,8),%xmm3
paddq  %xmm1,%xmm2
paddq  %xmm1,%xmm3
movdqu %xmm2,0x8(%r13,%rsi,8)
movdqu %xmm3,0x18(%r13,%rsi,8)
movdqu -0x10(%rdx,%rsi,8),%xmm2
movdqu (%rdx,%rsi,8),%xmm3
paddq  %xmm1,%xmm2
paddq  %xmm1,%xmm3
movdqu %xmm2,0x28(%r13,%rsi,8)
movdqu %xmm3,0x38(%r13,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1300 <submatrix+0x1e0>
test   $0x1,%r8b
je     1383 <submatrix+0x263>
lea    (%rsi,%r12,1),%rcx
add    %rbx,%rsi
movdqu 0x0(%rbp,%rcx,8),%xmm1
movdqu 0x10(%rbp,%rcx,8),%xmm2
pshufd $0x44,%xmm0,%xmm0
psllq  $0x3,%xmm0
paddq  %xmm0,%xmm1
paddq  %xmm2,%xmm0
movdqu %xmm1,(%rax,%rsi,8)
movdqu %xmm0,0x10(%rax,%rsi,8)
cmp    %r10,%r9
jne    1397 <submatrix+0x277>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
add    %r10,%r12
add    %rbx,%r10
jmp    11c4 <submatrix+0xa4>
xor    %esi,%esi
test   $0x1,%r8b
jne    1353 <submatrix+0x233>
jmp    1383 <submatrix+0x263>

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

LLVM IR (-O2):
@.str = [34 x i8] c"allocation failure in submatrix()\00"
define double** @submatrix(double** nocapture %0, i64 %1, i64 %2, i64 %3, i64 %4, i64 %5, i64 %6) {
%8 = sub nsw i64 %2, %1
%9 = sub nsw i64 %3, %6
%10 = shl i64 %8, 3
%11 = add i64 %10, 16
%12 = tail call i64 @malloc(i64 %11)
%13 = inttoptr i64 %12 to double**
%14 = icmp eq i64 %12, 0
br i1 %14, label %15, label %17
15: ; preds = %7
%16 = tail call i32 @nrerror(i8* getelementptr inbounds ([34 x i8], [34 x i8]* @.str, i64 0, i64 0))
br label %17
17: ; preds = %15, %7
%18 = getelementptr double*, double** %13, i64 1
%19 = sub i64 0, %5
%20 = getelementptr inbounds double*, double** %18, i64 %19
%21 = icmp slt i64 %2, %1
br i1 %21, label %154, label %22
22: ; preds = %17
%23 = add i64 %2, 1
%24 = sub i64 %23, %1
%25 = icmp ult i64 %24, 4
br i1 %25, label %102, label %26
26: ; preds = %22
%27 = add i64 %2, 2
%28 = sub i64 %27, %1
%29 = getelementptr double*, double** %13, i64 %28
%30 = getelementptr double*, double** %0, i64 %1
%31 = add i64 %2, 1
%32 = getelementptr double*, double** %0, i64 %31
%33 = icmp ult double** %18, %32
%34 = icmp ult double** %30, %29
%35 = and i1 %33, %34
br i1 %35, label %102, label %36
36: ; preds = %26
%37 = and i64 %24, -4
%38 = add i64 %37, %1
%39 = add i64 %37, %5
%40 = add i64 %37, -4
%41 = lshr exact i64 %40, 2
%42 = add nuw nsw i64 %41, 1
%43 = and i64 %42, 1
%44 = icmp eq i64 %40, 0
br i1 %44, label %82, label %45
45: ; preds = %36
%46 = and i64 %42, 9223372036854775806
br label %47
47: ; preds = %47, %45
%48 = phi i64 [ 0, %45 ], [ %79, %47 ]
%49 = phi i64 [ 0, %45 ], [ %80, %47 ]
%50 = add i64 %48, %1
%51 = add i64 %48, %5
%52 = getelementptr inbounds double*, double** %0, i64 %50
%53 = bitcast double** %52 to <2 x double*>*
%54 = load <2 x double*>, <2 x double*>* %53
%55 = getelementptr inbounds double*, double** %52, i64 2
%56 = bitcast double** %55 to <2 x double*>*
%57 = load <2 x double*>, <2 x double*>* %56
%58 = getelementptr inbounds double, <2 x double*> %54, i64 %9
%59 = getelementptr inbounds double, <2 x double*> %57, i64 %9
%60 = getelementptr inbounds double*, double** %20, i64 %51
%61 = bitcast double** %60 to <2 x double*>*
store <2 x double*> %58, <2 x double*>* %61
%62 = getelementptr inbounds double*, double** %60, i64 2
%63 = bitcast double** %62 to <2 x double*>*
store <2 x double*> %59, <2 x double*>* %63
%64 = or i64 %48, 4
%65 = add i64 %64, %1
%66 = add i64 %64, %5
%67 = getelementptr inbounds double*, double** %0, i64 %65
%68 = bitcast double** %67 to <2 x double*>*
%69 = load <2 x double*>, <2 x double*>* %68
%70 = getelementptr inbounds double*, double** %67, i64 2
%71 = bitcast double** %70 to <2 x double*>*
%72 = load <2 x double*>, <2 x double*>* %71
%73 = getelementptr inbounds double, <2 x double*> %69, i64 %9
%74 = getelementptr inbounds double, <2 x double*> %72, i64 %9
%75 = getelementptr inbounds double*, double** %20, i64 %66
%76 = bitcast double** %75 to <2 x double*>*
store <2 x double*> %73, <2 x double*>* %76
%77 = getelementptr inbounds double*, double** %75, i64 2
%78 = bitcast double** %77 to <2 x double*>*
store <2 x double*> %74, <2 x double*>* %78
%79 = add nuw i64 %48, 8
%80 = add i64 %49, 2
%81 = icmp eq i64 %80, %46
br i1 %81, label %82, label %47
82: ; preds = %47, %36
%83 = phi i64 [ 0, %36 ], [ %79, %47 ]
%84 = icmp eq i64 %43, 0
br i1 %84, label %100, label %85
85: ; preds = %82
%86 = add i64 %83, %1
%87 = add i64 %83, %5
%88 = getelementptr inbounds double*, double** %0, i64 %86
%89 = bitcast double** %88 to <2 x double*>*
%90 = load <2 x double*>, <2 x double*>* %89
%91 = getelementptr inbounds double*, double** %88, i64 2
%92 = bitcast double** %91 to <2 x double*>*
%93 = load <2 x double*>, <2 x double*>* %92
%94 = getelementptr inbounds double, <2 x double*> %90, i64 %9
%95 = getelementptr inbounds double, <2 x double*> %93, i64 %9
%96 = getelementptr inbounds double*, double** %20, i64 %87
%97 = bitcast double** %96 to <2 x double*>*
store <2 x double*> %94, <2 x double*>* %97
%98 = getelementptr inbounds double*, double** %96, i64 2
%99 = bitcast double** %98 to <2 x double*>*
store <2 x double*> %95, <2 x double*>* %99
br label %100
100: ; preds = %82, %85
%101 = icmp eq i64 %24, %37
br i1 %101, label %154, label %102
102: ; preds = %26, %22, %100
%103 = phi i64 [ %1, %26 ], [ %1, %22 ], [ %38, %100 ]
%104 = phi i64 [ %5, %26 ], [ %5, %22 ], [ %39, %100 ]
%105 = add i64 %2, 1
%106 = sub i64 %105, %103
%107 = sub i64 %2, %103
%108 = and i64 %106, 3
%109 = icmp eq i64 %108, 0
br i1 %109, label %122, label %110
110: ; preds = %102, %110
%111 = phi i64 [ %118, %110 ], [ %103, %102 ]
%112 = phi i64 [ %119, %110 ], [ %104, %102 ]
%113 = phi i64 [ %120, %110 ], [ 0, %102 ]
%114 = getelementptr inbounds double*, double** %0, i64 %111
%115 = load double*, double** %114
%116 = getelementptr inbounds double, double* %115, i64 %9
%117 = getelementptr inbounds double*, double** %20, i64 %112
store double* %116, double** %117
%118 = add i64 %111, 1
%119 = add nsw i64 %112, 1
%120 = add i64 %113, 1
%121 = icmp eq i64 %120, %108
br i1 %121, label %122, label %110
122: ; preds = %110, %102
%123 = phi i64 [ %103, %102 ], [ %118, %110 ]
%124 = phi i64 [ %104, %102 ], [ %119, %110 ]
%125 = icmp ult i64 %107, 3
br i1 %125, label %154, label %126
126: ; preds = %122, %126
%127 = phi i64 [ %151, %126 ], [ %123, %122 ]
%128 = phi i64 [ %152, %126 ], [ %124, %122 ]
%129 = getelementptr inbounds double*, double** %0, i64 %127
%130 = load double*, double** %129
%131 = getelementptr inbounds double, double* %130, i64 %9
%132 = getelementptr inbounds double*, double** %20, i64 %128
store double* %131, double** %132
%133 = add i64 %127, 1
%134 = add nsw i64 %128, 1
%135 = getelementptr inbounds double*, double** %0, i64 %133
%136 = load double*, double** %135
%137 = getelementptr inbounds double, double* %136, i64 %9
%138 = getelementptr inbounds double*, double** %20, i64 %134
store double* %137, double** %138
%139 = add i64 %127, 2
%140 = add nsw i64 %128, 2
%141 = getelementptr inbounds double*, double** %0, i64 %139
%142 = load double*, double** %141
%143 = getelementptr inbounds double, double* %142, i64 %9
%144 = getelementptr inbounds double*, double** %20, i64 %140
store double* %143, double** %144
%145 = add i64 %127, 3
%146 = add nsw i64 %128, 3
%147 = getelementptr inbounds double*, double** %0, i64 %145
%148 = load double*, double** %147
%149 = getelementptr inbounds double, double* %148, i64 %9
%150 = getelementptr inbounds double*, double** %20, i64 %146
store double* %149, double** %150
%151 = add i64 %127, 4
%152 = add nsw i64 %128, 4
%153 = icmp eq i64 %145, %2
br i1 %153, label %154, label %126
154: ; preds = %122, %126, %100, %17
ret double** %20
}