float **submatrix(a, oldrl, oldrh, oldcl, oldch, newrl, newcl)
float **a;
long newcl, newrl, oldch, oldcl, oldrh, oldrl;
{
  long i, j, nrow = oldrh - oldrl + 1, ncol = oldcl - newcl;
  float **m;
  m = (float **)malloc((unsigned int)((nrow + 1) * sizeof(float *)));
  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),%rax
add    $0x1,%rax
shl    $0x3,%rax
mov    %eax,%eax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x58(%rbp)
cmpq   $0x0,-0x58(%rbp)
jne    1192 <submatrix+0x72>
lea    0xe75(%rip),%rdi
mov    $0x0,%al
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     1214 <submatrix+0xf4>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
mov    -0x50(%rbp),%rax
shl    $0x2,%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    11c6 <submatrix+0xa6>
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,%r15
mov    %rcx,%r14
mov    %rdx,%rbx
mov    %rsi,%r12
mov    %rdi,%r13
mov    %ebx,%eax
sub    %r12d,%eax
lea    0x10(,%rax,8),%edi
call   1040 <malloc@plt>
mov    %rax,%rbp
test   %rax,%rax
jne    1161 <submatrix+0x41>
lea    0xea6(%rip),%rdi
xor    %eax,%eax
call   1030 <nrerror@plt>
lea    0x0(,%r15,8),%r9
mov    %rbp,%rax
sub    %r9,%rax
add    $0x8,%rax
mov    %rbx,%rdx
sub    %r12,%rdx
jl     12b4 <submatrix+0x194>
sub    0x40(%rsp),%r14
add    $0x1,%rdx
cmp    $0x4,%rdx
jb     126c <submatrix+0x14c>
mov    %rdx,%r8
and    $0xfffffffffffffffc,%r8
lea    -0x4(%r8),%rcx
mov    %rcx,%r10
shr    $0x2,%r10
add    $0x1,%r10
movq   %r14,%xmm0
test   %rcx,%rcx
je     12c3 <submatrix+0x1a3>
mov    %r10,%rsi
and    $0xfffffffffffffffe,%rsi
lea    0x30(,%r12,8),%rcx
add    %r13,%rcx
xor    %edi,%edi
pshufd $0x44,%xmm0,%xmm1
psllq  $0x2,%xmm1
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movdqu -0x30(%rcx,%rdi,8),%xmm2
movdqu -0x20(%rcx,%rdi,8),%xmm3
paddq  %xmm1,%xmm2
paddq  %xmm1,%xmm3
movdqu %xmm2,0x8(%rbp,%rdi,8)
movdqu %xmm3,0x18(%rbp,%rdi,8)
movdqu -0x10(%rcx,%rdi,8),%xmm2
movdqu (%rcx,%rdi,8),%xmm3
paddq  %xmm1,%xmm2
paddq  %xmm1,%xmm3
movdqu %xmm2,0x28(%rbp,%rdi,8)
movdqu %xmm3,0x38(%rbp,%rdi,8)
add    $0x8,%rdi
add    $0xfffffffffffffffe,%rsi
jne    11e0 <submatrix+0xc0>
test   $0x1,%r10b
je     1261 <submatrix+0x141>
lea    (%rdi,%r12,1),%rcx
add    %r15,%rdi
movdqu 0x0(%r13,%rcx,8),%xmm1
movdqu 0x10(%r13,%rcx,8),%xmm2
pshufd $0x44,%xmm0,%xmm0
psllq  $0x2,%xmm0
paddq  %xmm0,%xmm1
paddq  %xmm2,%xmm0
movdqu %xmm1,(%rax,%rdi,8)
movdqu %xmm0,0x10(%rax,%rdi,8)
cmp    %r8,%rdx
je     12b4 <submatrix+0x194>
add    %r8,%r12
add    %r8,%r15
sub    %r12,%rbx
add    $0x1,%rbx
lea    0x0(,%r12,8),%rcx
add    %r13,%rcx
shl    $0x3,%r15
sub    %r9,%r15
lea    (%r15,%rbp,1),%rdx
add    $0x8,%rdx
xor    %esi,%esi
shl    $0x2,%r14
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    (%rcx,%rsi,8),%rdi
add    %r14,%rdi
mov    %rdi,(%rdx,%rsi,8)
add    $0x1,%rsi
cmp    %rsi,%rbx
jne    12a0 <submatrix+0x180>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %edi,%edi
test   $0x1,%r10b
jne    122f <submatrix+0x10f>
jmp    1261 <submatrix+0x141>

LLVM IR (-O0):
@.str = [34 x i8] c"allocation failure in submatrix()\00"
define float** @submatrix(float** %0, i64 %1, i64 %2, i64 %3, i64 %4, i64 %5, i64 %6) {
%8 = alloca float**
%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 float**
store float** %0, float*** %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 = trunc i64 %29 to i32
%31 = zext i32 %30 to i64
%32 = call i8* @malloc(i64 %31)
%33 = bitcast i8* %32 to float**
store float** %33, float*** %19
%34 = load float**, float*** %19
%35 = icmp ne float** %34, null
br i1 %35, label %38, label %36
36: ; preds = %7
%37 = call i32 (i8*, ...) bitcast (i32 (...)* @nrerror to i32 (i8*, ...)*)(i8* getelementptr inbounds ([34 x i8], [34 x i8]* @.str, i64 0, i64 0))
br label %38
38: ; preds = %36, %7
%39 = load float**, float*** %19
%40 = getelementptr inbounds float*, float** %39, i64 1
store float** %40, float*** %19
%41 = load i64, i64* %13
%42 = load float**, float*** %19
%43 = sub i64 0, %41
%44 = getelementptr inbounds float*, float** %42, i64 %43
store float** %44, float*** %19
%45 = load i64, i64* %9
store i64 %45, i64* %15
%46 = load i64, i64* %13
store i64 %46, i64* %16
br label %47
47: ; preds = %61, %38
%48 = load i64, i64* %15
%49 = load i64, i64* %10
%50 = icmp sle i64 %48, %49
br i1 %50, label %51, label %66
51: ; preds = %47
%52 = load float**, float*** %8
%53 = load i64, i64* %15
%54 = getelementptr inbounds float*, float** %52, i64 %53
%55 = load float*, float** %54
%56 = load i64, i64* %18
%57 = getelementptr inbounds float, float* %55, i64 %56
%58 = load float**, float*** %19
%59 = load i64, i64* %16
%60 = getelementptr inbounds float*, float** %58, i64 %59
store float* %57, float** %60
br label %61
61: ; preds = %51
%62 = load i64, i64* %15
%63 = add nsw i64 %62, 1
store i64 %63, i64* %15
%64 = load i64, i64* %16
%65 = add nsw i64 %64, 1
store i64 %65, i64* %16
br label %47
66: ; preds = %47
%67 = load float**, float*** %19
ret float** %67
}

LLVM IR (-O2):
@.str = [34 x i8] c"allocation failure in submatrix()\00"
define float** @submatrix(float** 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 = and i64 %11, 4294967288
%13 = tail call i8* @malloc(i64 %12)
%14 = bitcast i8* %13 to float**
%15 = icmp eq i8* %13, null
br i1 %15, label %16, label %18
16: ; preds = %7
%17 = tail call i32 (i8*, ...) bitcast (i32 (...)* @nrerror to i32 (i8*, ...)*)(i8* getelementptr inbounds ([34 x i8], [34 x i8]* @.str, i64 0, i64 0))
br label %18
18: ; preds = %16, %7
%19 = getelementptr inbounds float*, float** %14, i64 1
%20 = sub i64 0, %5
%21 = getelementptr inbounds float*, float** %19, i64 %20
%22 = icmp slt i64 %2, %1
br i1 %22, label %106, label %23
23: ; preds = %18
%24 = add i64 %2, 1
%25 = sub i64 %24, %1
%26 = icmp ult i64 %25, 4
br i1 %26, label %93, label %27
27: ; preds = %23
%28 = and i64 %25, -4
%29 = add i64 %28, %1
%30 = add i64 %28, %5
%31 = add i64 %28, -4
%32 = lshr exact i64 %31, 2
%33 = add nuw nsw i64 %32, 1
%34 = and i64 %33, 1
%35 = icmp eq i64 %31, 0
br i1 %35, label %73, label %36
36: ; preds = %27
%37 = and i64 %33, 9223372036854775806
br label %38
38: ; preds = %38, %36
%39 = phi i64 [ 0, %36 ], [ %70, %38 ]
%40 = phi i64 [ 0, %36 ], [ %71, %38 ]
%41 = add i64 %39, %1
%42 = add i64 %39, %5
%43 = getelementptr inbounds float*, float** %0, i64 %41
%44 = bitcast float** %43 to <2 x float*>*
%45 = load <2 x float*>, <2 x float*>* %44
%46 = getelementptr inbounds float*, float** %43, i64 2
%47 = bitcast float** %46 to <2 x float*>*
%48 = load <2 x float*>, <2 x float*>* %47
%49 = getelementptr inbounds float, <2 x float*> %45, i64 %9
%50 = getelementptr inbounds float, <2 x float*> %48, i64 %9
%51 = getelementptr inbounds float*, float** %21, i64 %42
%52 = bitcast float** %51 to <2 x float*>*
store <2 x float*> %49, <2 x float*>* %52
%53 = getelementptr inbounds float*, float** %51, i64 2
%54 = bitcast float** %53 to <2 x float*>*
store <2 x float*> %50, <2 x float*>* %54
%55 = or i64 %39, 4
%56 = add i64 %55, %1
%57 = add i64 %55, %5
%58 = getelementptr inbounds float*, float** %0, i64 %56
%59 = bitcast float** %58 to <2 x float*>*
%60 = load <2 x float*>, <2 x float*>* %59
%61 = getelementptr inbounds float*, float** %58, i64 2
%62 = bitcast float** %61 to <2 x float*>*
%63 = load <2 x float*>, <2 x float*>* %62
%64 = getelementptr inbounds float, <2 x float*> %60, i64 %9
%65 = getelementptr inbounds float, <2 x float*> %63, i64 %9
%66 = getelementptr inbounds float*, float** %21, i64 %57
%67 = bitcast float** %66 to <2 x float*>*
store <2 x float*> %64, <2 x float*>* %67
%68 = getelementptr inbounds float*, float** %66, i64 2
%69 = bitcast float** %68 to <2 x float*>*
store <2 x float*> %65, <2 x float*>* %69
%70 = add nuw i64 %39, 8
%71 = add i64 %40, 2
%72 = icmp eq i64 %71, %37
br i1 %72, label %73, label %38
73: ; preds = %38, %27
%74 = phi i64 [ 0, %27 ], [ %70, %38 ]
%75 = icmp eq i64 %34, 0
br i1 %75, label %91, label %76
76: ; preds = %73
%77 = add i64 %74, %1
%78 = add i64 %74, %5
%79 = getelementptr inbounds float*, float** %0, i64 %77
%80 = bitcast float** %79 to <2 x float*>*
%81 = load <2 x float*>, <2 x float*>* %80
%82 = getelementptr inbounds float*, float** %79, i64 2
%83 = bitcast float** %82 to <2 x float*>*
%84 = load <2 x float*>, <2 x float*>* %83
%85 = getelementptr inbounds float, <2 x float*> %81, i64 %9
%86 = getelementptr inbounds float, <2 x float*> %84, i64 %9
%87 = getelementptr inbounds float*, float** %21, i64 %78
%88 = bitcast float** %87 to <2 x float*>*
store <2 x float*> %85, <2 x float*>* %88
%89 = getelementptr inbounds float*, float** %87, i64 2
%90 = bitcast float** %89 to <2 x float*>*
store <2 x float*> %86, <2 x float*>* %90
br label %91
91: ; preds = %73, %76
%92 = icmp eq i64 %25, %28
br i1 %92, label %106, label %93
93: ; preds = %23, %91
%94 = phi i64 [ %1, %23 ], [ %29, %91 ]
%95 = phi i64 [ %5, %23 ], [ %30, %91 ]
br label %96
96: ; preds = %93, %96
%97 = phi i64 [ %103, %96 ], [ %94, %93 ]
%98 = phi i64 [ %104, %96 ], [ %95, %93 ]
%99 = getelementptr inbounds float*, float** %0, i64 %97
%100 = load float*, float** %99
%101 = getelementptr inbounds float, float* %100, i64 %9
%102 = getelementptr inbounds float*, float** %21, i64 %98
store float* %101, float** %102
%103 = add i64 %97, 1
%104 = add nsw i64 %98, 1
%105 = icmp eq i64 %97, %2
br i1 %105, label %106, label %96
106: ; preds = %96, %91, %18
ret float** %21
}