float **convert_matrix(a, nrl, nrh, ncl, nch)
float *a;
long nch, ncl, nrh, nrl;
{
  long i, j, nrow = nrh - nrl + 1, ncol = nch - ncl + 1;
  float **m;
  m = (float **)malloc((unsigned int)((nrow + 1) * sizeof(float *)));
  if (!m)
    nrerror("allocation failure in convert_matrix()");
  m += 1;
  m -= nrl;
  m[nrl] = a - ncl;
  for (i = 1, j = nrl + 1; i < nrow; i++, j++)
    m[j] = m[j - 1] + ncol;
  return m;
}


x86 assembly (clang, -O0):
<convert_matrix>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x18(%rbp),%rax
sub    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%rax
sub    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
shl    $0x3,%rax
mov    %eax,%eax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x50(%rbp)
cmpq   $0x0,-0x50(%rbp)
jne    118e <convert_matrix+0x6e>
lea    0xe79(%rip),%rdi
mov    $0x0,%al
call   1030 <nrerror@plt>
mov    -0x50(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x50(%rbp)
mov    -0x10(%rbp),%rdx
mov    -0x50(%rbp),%rax
xor    %ecx,%ecx
sub    %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x8(%rbp),%rdx
xor    %eax,%eax
sub    -0x20(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdx
mov    -0x50(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
movq   $0x1,-0x30(%rbp)
mov    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x40(%rbp),%rax
jge    1235 <convert_matrix+0x115>
mov    -0x50(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    $0x1,%rcx
mov    (%rax,%rcx,8),%rdx
mov    -0x48(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdx
mov    -0x50(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11e3 <convert_matrix+0xc3>
mov    -0x50(%rbp),%rax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<convert_matrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %r8,%r12
mov    %rcx,%rbp
mov    %rdx,%r15
mov    %rsi,%r14
mov    %rdi,%rbx
mov    %rcx,0x8(%rsp)
neg    %rbp
mov    %rdx,(%rsp)
sub    %rsi,%r15
lea    0x10(,%r15,8),%edi
call   1040 <malloc@plt>
mov    %rax,%r13
test   %rax,%rax
jne    116f <convert_matrix+0x4f>
lea    0xe98(%rip),%rdi
xor    %eax,%eax
call   1030 <nrerror@plt>
lea    0x0(,%r14,8),%rdi
mov    %r13,%rax
sub    %rdi,%rax
lea    (%rbx,%rbp,4),%rcx
mov    %rcx,0x8(%rax,%r14,8)
add    $0x8,%rax
test   %r15,%r15
jle    12b0 <convert_matrix+0x190>
mov    %r14,%rdx
not    %rdx
add    (%rsp),%rdx
mov    %r15d,%r10d
and    $0x3,%r10d
cmp    $0x3,%rdx
jb     125f <convert_matrix+0x13f>
mov    %rdi,(%rsp)
and    $0xfffffffffffffffc,%r15
mov    %r12,%r11
shl    $0x4,%r11
add    $0x10,%r11
mov    0x8(%rsp),%r8
mov    %r8,%rdx
shl    $0x4,%rdx
mov    %r11,%r9
sub    %rdx,%r9
lea    0x0(,%r8,4),%rsi
lea    (%rsi,%rsi,4),%rcx
sub    %rcx,%r11
lea    0x0(,%r12,4),%rdi
lea    (%rdi,%rdi,2),%rcx
sub    %rdx,%rcx
add    $0xc,%rcx
lea    0x0(,%r12,8),%rbp
mov    %rsi,0x10(%rsp)
lea    (%rsi,%rsi,2),%rdx
sub    %rdx,%rbp
add    $0x8,%rbp
lea    0x0(,%r8,8),%rdx
sub    %rdx,%rdi
add    $0x4,%rdi
xor    %edx,%edx
xchg   %ax,%ax
lea    (%rbx,%rdi,1),%r8
mov    %r8,0x10(%r13,%rdx,8)
lea    (%rbx,%rbp,1),%rsi
mov    %rsi,0x18(%r13,%rdx,8)
lea    (%rbx,%rcx,1),%rsi
mov    %rsi,0x20(%r13,%rdx,8)
lea    (%rbx,%r11,1),%rsi
mov    %rsi,0x28(%r13,%rdx,8)
add    $0x4,%rdx
add    %r9,%rbx
cmp    %rdx,%r15
jne    1220 <convert_matrix+0x100>
sub    0x10(%rsp),%rbx
add    %rdx,%r14
mov    %rbx,%rcx
mov    (%rsp),%rdi
test   %r10,%r10
je     12b0 <convert_matrix+0x190>
mov    %r12,%rdx
mov    0x8(%rsp),%rsi
sub    %rsi,%rdx
add    $0x1,%rdx
lea    (%rcx,%rdx,4),%rcx
shl    $0x2,%r12
shl    $0x2,%rsi
sub    %rsi,%r12
add    $0x4,%r12
shl    $0x3,%r14
sub    %rdi,%r14
lea    (%r14,%r13,1),%rdx
add    $0x10,%rdx
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
mov    %rcx,(%rdx,%rsi,8)
add    $0x1,%rsi
add    %r12,%rcx
cmp    %rsi,%r10
jne    12a0 <convert_matrix+0x180>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@.str = [39 x i8] c"allocation failure in convert_matrix()\00"
define float** @convert_matrix(float* %0, i64 %1, i64 %2, i64 %3, i64 %4) {
%6 = sub i64 0, %3
%7 = sub nsw i64 %2, %1
%8 = sub i64 1, %3
%9 = add i64 %8, %4
%10 = shl i64 %7, 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 = %5
%17 = tail call i32 (i8*, ...) bitcast (i32 (...)* @nrerror to i32 (i8*, ...)*)(i8* getelementptr inbounds ([39 x i8], [39 x i8]* @.str, i64 0, i64 0))
br label %18
18: ; preds = %16, %5
%19 = getelementptr inbounds float*, float** %14, i64 1
%20 = sub i64 0, %1
%21 = getelementptr inbounds float*, float** %19, i64 %20
%22 = getelementptr inbounds float, float* %0, i64 %6
%23 = getelementptr inbounds float*, float** %21, i64 %1
store float* %22, float** %23
%24 = icmp slt i64 %7, 1
br i1 %24, label %63, label %25
25: ; preds = %18
%26 = xor i64 %1, -1
%27 = add i64 %26, %2
%28 = and i64 %7, 3
%29 = icmp ult i64 %27, 3
br i1 %29, label %50, label %30
30: ; preds = %25
%31 = and i64 %7, -4
br label %32
32: ; preds = %32, %30
%33 = phi float* [ %22, %30 ], [ %46, %32 ]
%34 = phi i64 [ %1, %30 ], [ %45, %32 ]
%35 = phi i64 [ 0, %30 ], [ %48, %32 ]
%36 = add nsw i64 %34, 1
%37 = getelementptr inbounds float, float* %33, i64 %9
%38 = getelementptr inbounds float*, float** %21, i64 %36
store float* %37, float** %38
%39 = add nsw i64 %34, 2
%40 = getelementptr inbounds float, float* %37, i64 %9
%41 = getelementptr inbounds float*, float** %21, i64 %39
store float* %40, float** %41
%42 = add nsw i64 %34, 3
%43 = getelementptr inbounds float, float* %40, i64 %9
%44 = getelementptr inbounds float*, float** %21, i64 %42
store float* %43, float** %44
%45 = add nsw i64 %34, 4
%46 = getelementptr inbounds float, float* %43, i64 %9
%47 = getelementptr inbounds float*, float** %21, i64 %45
store float* %46, float** %47
%48 = add i64 %35, 4
%49 = icmp eq i64 %48, %31
br i1 %49, label %50, label %32
50: ; preds = %32, %25
%51 = phi float* [ %22, %25 ], [ %46, %32 ]
%52 = phi i64 [ %1, %25 ], [ %45, %32 ]
%53 = icmp eq i64 %28, 0
br i1 %53, label %63, label %54
54: ; preds = %50, %54
%55 = phi float* [ %59, %54 ], [ %51, %50 ]
%56 = phi i64 [ %58, %54 ], [ %52, %50 ]
%57 = phi i64 [ %61, %54 ], [ 0, %50 ]
%58 = add nsw i64 %56, 1
%59 = getelementptr inbounds float, float* %55, i64 %9
%60 = getelementptr inbounds float*, float** %21, i64 %58
store float* %59, float** %60
%61 = add i64 %57, 1
%62 = icmp eq i64 %61, %28
br i1 %62, label %63, label %54
63: ; preds = %50, %54, %18
ret float** %21
}