void extract_smatrix(short **src, long nrl, long nrh, long ncl, long nch,
                     short **dst) {
  int i, j, ncol = nch - ncl + 1;
  for (i = nrl; i <= nrh; i++) {
    for (j = ncl; j <= nch; j++) {
      dst[i - nrl][j - ncl] = src[i][j];
    }
  }
}


x86 assembly (clang, -O0):
<extract_smatrix>:
push   %rbp
mov    %rsp,%rbp
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    -0x28(%rbp),%rax
sub    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %eax,-0x3c(%rbp)
mov    -0x10(%rbp),%rax
mov    %eax,-0x34(%rbp)
movslq -0x34(%rbp),%rax
cmp    -0x18(%rbp),%rax
jg     11a6 <extract_smatrix+0xa6>
mov    -0x20(%rbp),%rax
mov    %eax,-0x38(%rbp)
movslq -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jg     1193 <extract_smatrix+0x93>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,2),%dx
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
sub    -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
sub    -0x20(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1147 <extract_smatrix+0x47>
jmp    1198 <extract_smatrix+0x98>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1132 <extract_smatrix+0x32>
pop    %rbp
ret

x86 assembly (clang, -O2):
<extract_smatrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rcx,-0x58(%rsp)
movslq %esi,%rax
cmp    %rdx,%rax
jle    1122 <extract_smatrix+0x22>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %rdi,%r14
mov    -0x58(%rsp),%rdi
movslq %edi,%r10
mov    %r10,%r13
sub    %rdi,%r13
lea    0x1(%r8),%rcx
lea    0x1(%r10),%r15
cmp    %r15,%rcx
cmovg  %rcx,%r15
mov    %r15,%rcx
sub    %rdi,%rcx
mov    %rcx,-0x10(%rsp)
mov    %r15,%rcx
mov    %r15,-0x8(%rsp)
sub    %r10,%r15
lea    -0x10(%r15),%rcx
mov    %rcx,-0x20(%rsp)
mov    %rcx,%rbp
shr    $0x4,%rbp
add    $0x1,%rbp
mov    %r15,%rcx
and    $0xfffffffffffffff0,%rcx
mov    %rcx,-0x30(%rsp)
add    %r10,%rcx
mov    %rcx,-0x38(%rsp)
mov    %rbp,%rcx
mov    %rbp,-0x28(%rsp)
and    $0xfffffffffffffffe,%rbp
mov    %rbp,-0x40(%rsp)
lea    0x18(%r10),%rbp
mov    %rbp,%rcx
mov    %rbp,-0x48(%rsp)
sub    %rdi,%rbp
mov    %rbp,-0x50(%rsp)
neg    %rdi
mov    %r14,-0x18(%rsp)
jmp    11bd <extract_smatrix+0xbd>
nop
add    $0x1,%rax
cmp    %rdx,%rax
jg     1117 <extract_smatrix+0x17>
cmp    %r8,%r10
jg     11b0 <extract_smatrix+0xb0>
mov    %rax,%rcx
sub    %rsi,%rcx
mov    (%r14,%rax,8),%r12
mov    (%r9,%rcx,8),%rbp
mov    %r10,%rcx
cmp    $0x10,%r15
jb     12d0 <extract_smatrix+0x1d0>
lea    0x0(,%r13,2),%rcx
add    %rbp,%rcx
mov    -0x8(%rsp),%rbx
lea    (%r12,%rbx,2),%rbx
cmp    %rbx,%rcx
jae    1216 <extract_smatrix+0x116>
mov    -0x10(%rsp),%rcx
lea    0x0(,%rcx,2),%r11
add    %rbp,%r11
lea    (%r12,%r10,2),%rbx
mov    %r10,%rcx
cmp    %r11,%rbx
jb     12d0 <extract_smatrix+0x1d0>
mov    %r13,%rbx
cmpq   $0x10,-0x20(%rsp)
jae    1226 <extract_smatrix+0x126>
xor    %r14d,%r14d
jmp    1289 <extract_smatrix+0x189>
mov    -0x50(%rsp),%rcx
lea    0x0(,%rcx,2),%r13
add    %rbp,%r13
mov    -0x48(%rsp),%rcx
lea    (%r12,%rcx,2),%r11
mov    -0x40(%rsp),%rcx
xor    %r14d,%r14d
nopw   0x0(%rax,%rax,1)
movups -0x30(%r11,%r14,2),%xmm0
movups -0x20(%r11,%r14,2),%xmm1
movups %xmm0,-0x30(%r13,%r14,2)
movups %xmm1,-0x20(%r13,%r14,2)
movups -0x10(%r11,%r14,2),%xmm0
movups (%r11,%r14,2),%xmm1
movups %xmm0,-0x10(%r13,%r14,2)
movups %xmm1,0x0(%r13,%r14,2)
add    $0x20,%r14
add    $0xfffffffffffffffe,%rcx
jne    1250 <extract_smatrix+0x150>
testb  $0x1,-0x28(%rsp)
mov    %rbx,%r13
je     12b2 <extract_smatrix+0x1b2>
add    %r10,%r14
movups (%r12,%r14,2),%xmm0
movups 0x10(%r12,%r14,2),%xmm1
sub    -0x58(%rsp),%r14
movups %xmm0,0x0(%rbp,%r14,2)
movups %xmm1,0x10(%rbp,%r14,2)
mov    -0x38(%rsp),%rcx
cmp    -0x30(%rsp),%r15
mov    -0x18(%rsp),%r14
je     11b0 <extract_smatrix+0xb0>
nopw   0x0(%rax,%rax,1)
lea    0x0(%rbp,%rdi,2),%rbp
cs nopw 0x0(%rax,%rax,1)
nop
movzwl (%r12,%rcx,2),%ebx
mov    %bx,0x0(%rbp,%rcx,2)
add    $0x1,%rcx
cmp    %r8,%rcx
jle    12e0 <extract_smatrix+0x1e0>
jmp    11b0 <extract_smatrix+0xb0>

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

LLVM IR (-O2):
define void @extract_smatrix(i16** nocapture %0, i64 %1, i64 %2, i64 %3, i64 %4, i16** nocapture %5) {
%7 = shl i64 %1, 32
%8 = ashr exact i64 %7, 32
%9 = icmp sgt i64 %8, %2
br i1 %9, label %111, label %10
10: ; preds = %6
%11 = shl i64 %3, 32
%12 = ashr exact i64 %11, 32
%13 = icmp sgt i64 %12, %4
%14 = sub i64 %12, %3
%15 = add nsw i64 %4, 1
%16 = add nsw i64 %12, 1
%17 = call i64 @llvm.smax.i64(i64 %15, i64 %16)
%18 = sub i64 %17, %3
%19 = add nsw i64 %4, 1
%20 = add nsw i64 %12, 1
%21 = call i64 @llvm.smax.i64(i64 %19, i64 %20)
%22 = sub i64 %21, %12
%23 = add i64 %22, -16
%24 = lshr i64 %23, 4
%25 = add nuw nsw i64 %24, 1
%26 = icmp ult i64 %22, 16
%27 = and i64 %22, -16
%28 = add i64 %12, %27
%29 = and i64 %25, 1
%30 = icmp ult i64 %23, 16
%31 = and i64 %25, 2305843009213693950
%32 = icmp eq i64 %29, 0
%33 = icmp eq i64 %22, %27
br label %34
34: ; preds = %10, %108
%35 = phi i64 [ %8, %10 ], [ %109, %108 ]
br i1 %13, label %108, label %36
36: ; preds = %34
%37 = getelementptr inbounds i16*, i16** %0, i64 %35
%38 = load i16*, i16** %37
%39 = sub nsw i64 %35, %1
%40 = getelementptr inbounds i16*, i16** %5, i64 %39
%41 = load i16*, i16** %40
br i1 %26, label %98, label %42
42: ; preds = %36
%43 = getelementptr i16, i16* %41, i64 %14
%44 = getelementptr i16, i16* %41, i64 %18
%45 = getelementptr i16, i16* %38, i64 %12
%46 = getelementptr i16, i16* %38, i64 %17
%47 = icmp ult i16* %43, %46
%48 = icmp ult i16* %45, %44
%49 = and i1 %47, %48
br i1 %49, label %98, label %50
50: ; preds = %42
br i1 %30, label %82, label %51
51: ; preds = %50, %51
%52 = phi i64 [ %79, %51 ], [ 0, %50 ]
%53 = phi i64 [ %80, %51 ], [ 0, %50 ]
%54 = add i64 %12, %52
%55 = getelementptr inbounds i16, i16* %38, i64 %54
%56 = bitcast i16* %55 to <8 x i16>*
%57 = load <8 x i16>, <8 x i16>* %56
%58 = getelementptr inbounds i16, i16* %55, i64 8
%59 = bitcast i16* %58 to <8 x i16>*
%60 = load <8 x i16>, <8 x i16>* %59
%61 = sub nsw i64 %54, %3
%62 = getelementptr inbounds i16, i16* %41, i64 %61
%63 = bitcast i16* %62 to <8 x i16>*
store <8 x i16> %57, <8 x i16>* %63
%64 = getelementptr inbounds i16, i16* %62, i64 8
%65 = bitcast i16* %64 to <8 x i16>*
store <8 x i16> %60, <8 x i16>* %65
%66 = or i64 %52, 16
%67 = add i64 %12, %66
%68 = getelementptr inbounds i16, i16* %38, i64 %67
%69 = bitcast i16* %68 to <8 x i16>*
%70 = load <8 x i16>, <8 x i16>* %69
%71 = getelementptr inbounds i16, i16* %68, i64 8
%72 = bitcast i16* %71 to <8 x i16>*
%73 = load <8 x i16>, <8 x i16>* %72
%74 = sub nsw i64 %67, %3
%75 = getelementptr inbounds i16, i16* %41, i64 %74
%76 = bitcast i16* %75 to <8 x i16>*
store <8 x i16> %70, <8 x i16>* %76
%77 = getelementptr inbounds i16, i16* %75, i64 8
%78 = bitcast i16* %77 to <8 x i16>*
store <8 x i16> %73, <8 x i16>* %78
%79 = add nuw i64 %52, 32
%80 = add i64 %53, 2
%81 = icmp eq i64 %80, %31
br i1 %81, label %82, label %51
82: ; preds = %51, %50
%83 = phi i64 [ 0, %50 ], [ %79, %51 ]
br i1 %32, label %97, label %84
84: ; preds = %82
%85 = add i64 %12, %83
%86 = getelementptr inbounds i16, i16* %38, i64 %85
%87 = bitcast i16* %86 to <8 x i16>*
%88 = load <8 x i16>, <8 x i16>* %87
%89 = getelementptr inbounds i16, i16* %86, i64 8
%90 = bitcast i16* %89 to <8 x i16>*
%91 = load <8 x i16>, <8 x i16>* %90
%92 = sub nsw i64 %85, %3
%93 = getelementptr inbounds i16, i16* %41, i64 %92
%94 = bitcast i16* %93 to <8 x i16>*
store <8 x i16> %88, <8 x i16>* %94
%95 = getelementptr inbounds i16, i16* %93, i64 8
%96 = bitcast i16* %95 to <8 x i16>*
store <8 x i16> %91, <8 x i16>* %96
br label %97
97: ; preds = %82, %84
br i1 %33, label %108, label %98
98: ; preds = %42, %36, %97
%99 = phi i64 [ %12, %42 ], [ %12, %36 ], [ %28, %97 ]
br label %100
100: ; preds = %98, %100
%101 = phi i64 [ %106, %100 ], [ %99, %98 ]
%102 = getelementptr inbounds i16, i16* %38, i64 %101
%103 = load i16, i16* %102
%104 = sub nsw i64 %101, %3
%105 = getelementptr inbounds i16, i16* %41, i64 %104
store i16 %103, i16* %105
%106 = add i64 %101, 1
%107 = icmp sgt i64 %106, %4
br i1 %107, label %108, label %100
108: ; preds = %100, %97, %34
%109 = add i64 %35, 1
%110 = icmp sgt i64 %109, %2
br i1 %110, label %111, label %34
111: ; preds = %108, %6
ret void
}