void extract_dmatrix(double **src, long nrl, long nrh, long ncl, long nch,
                     double **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_dmatrix>:
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     11a8 <extract_dmatrix+0xa8>
mov    -0x20(%rbp),%rax
mov    %eax,-0x38(%rbp)
movslq -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jg     1195 <extract_dmatrix+0x95>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
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
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1147 <extract_dmatrix+0x47>
jmp    119a <extract_dmatrix+0x9a>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1132 <extract_dmatrix+0x32>
pop    %rbp
ret

x86 assembly (clang, -O2):
<extract_dmatrix>:
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_dmatrix+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    -0x4(%r15),%rcx
mov    %rcx,-0x20(%rsp)
mov    %rcx,%rbp
shr    $0x2,%rbp
add    $0x1,%rbp
mov    %r15,%rcx
and    $0xfffffffffffffffc,%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    0x6(%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_dmatrix+0xbd>
nop
add    $0x1,%rax
cmp    %rdx,%rax
jg     1117 <extract_dmatrix+0x17>
cmp    %r8,%r10
jg     11b0 <extract_dmatrix+0xb0>
mov    %rax,%rcx
sub    %rsi,%rcx
mov    (%r14,%rax,8),%r12
mov    (%r9,%rcx,8),%rbp
mov    %r10,%rcx
cmp    $0x4,%r15
jb     12d0 <extract_dmatrix+0x1d0>
lea    0x0(,%r13,8),%rcx
add    %rbp,%rcx
mov    -0x8(%rsp),%rbx
lea    (%r12,%rbx,8),%rbx
cmp    %rbx,%rcx
jae    1216 <extract_dmatrix+0x116>
mov    -0x10(%rsp),%rcx
lea    0x0(,%rcx,8),%r11
add    %rbp,%r11
lea    (%r12,%r10,8),%rbx
mov    %r10,%rcx
cmp    %r11,%rbx
jb     12d0 <extract_dmatrix+0x1d0>
mov    %r13,%rbx
cmpq   $0x4,-0x20(%rsp)
jae    1226 <extract_dmatrix+0x126>
xor    %r14d,%r14d
jmp    1289 <extract_dmatrix+0x189>
mov    -0x50(%rsp),%rcx
lea    0x0(,%rcx,8),%r13
add    %rbp,%r13
mov    -0x48(%rsp),%rcx
lea    (%r12,%rcx,8),%r11
mov    -0x40(%rsp),%rcx
xor    %r14d,%r14d
nopw   0x0(%rax,%rax,1)
movups -0x30(%r11,%r14,8),%xmm0
movups -0x20(%r11,%r14,8),%xmm1
movups %xmm0,-0x30(%r13,%r14,8)
movups %xmm1,-0x20(%r13,%r14,8)
movups -0x10(%r11,%r14,8),%xmm0
movups (%r11,%r14,8),%xmm1
movups %xmm0,-0x10(%r13,%r14,8)
movups %xmm1,0x0(%r13,%r14,8)
add    $0x8,%r14
add    $0xfffffffffffffffe,%rcx
jne    1250 <extract_dmatrix+0x150>
testb  $0x1,-0x28(%rsp)
mov    %rbx,%r13
je     12b2 <extract_dmatrix+0x1b2>
add    %r10,%r14
movups (%r12,%r14,8),%xmm0
movups 0x10(%r12,%r14,8),%xmm1
sub    -0x58(%rsp),%r14
movups %xmm0,0x0(%rbp,%r14,8)
movups %xmm1,0x10(%rbp,%r14,8)
mov    -0x38(%rsp),%rcx
cmp    -0x30(%rsp),%r15
mov    -0x18(%rsp),%r14
je     11b0 <extract_dmatrix+0xb0>
nopw   0x0(%rax,%rax,1)
lea    0x0(%rbp,%rdi,8),%rbp
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%r12,%rcx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rcx,8)
add    $0x1,%rcx
cmp    %r8,%rcx
jle    12e0 <extract_dmatrix+0x1e0>
jmp    11b0 <extract_dmatrix+0xb0>

LLVM IR (-O0):
define void @extract_dmatrix(double** %0, i64 %1, i64 %2, i64 %3, i64 %4, double** %5) {
%7 = alloca double**
%8 = alloca i64
%9 = alloca i64
%10 = alloca i64
%11 = alloca i64
%12 = alloca double**
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
store double** %0, double*** %7
store i64 %1, i64* %8
store i64 %2, i64* %9
store i64 %3, i64* %10
store i64 %4, i64* %11
store double** %5, double*** %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 double**, double*** %7
%38 = load i32, i32* %13
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds double*, double** %37, i64 %39
%41 = load double*, double** %40
%42 = load i32, i32* %14
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds double, double* %41, i64 %43
%45 = load double, double* %44
%46 = load double**, double*** %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 double*, double** %46, i64 %50
%52 = load double*, double** %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 double, double* %52, i64 %56
store double %45, double* %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_dmatrix(double** nocapture %0, i64 %1, i64 %2, i64 %3, i64 %4, double** 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, -4
%24 = lshr i64 %23, 2
%25 = add nuw nsw i64 %24, 1
%26 = icmp ult i64 %22, 4
%27 = and i64 %22, -4
%28 = add i64 %12, %27
%29 = and i64 %25, 1
%30 = icmp ult i64 %23, 4
%31 = and i64 %25, 9223372036854775806
%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 double*, double** %0, i64 %35
%38 = load double*, double** %37
%39 = sub nsw i64 %35, %1
%40 = getelementptr inbounds double*, double** %5, i64 %39
%41 = load double*, double** %40
br i1 %26, label %98, label %42
42: ; preds = %36
%43 = getelementptr double, double* %41, i64 %14
%44 = getelementptr double, double* %41, i64 %18
%45 = getelementptr double, double* %38, i64 %12
%46 = getelementptr double, double* %38, i64 %17
%47 = icmp ult double* %43, %46
%48 = icmp ult double* %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 double, double* %38, i64 %54
%56 = bitcast double* %55 to <2 x double>*
%57 = load <2 x double>, <2 x double>* %56
%58 = getelementptr inbounds double, double* %55, i64 2
%59 = bitcast double* %58 to <2 x double>*
%60 = load <2 x double>, <2 x double>* %59
%61 = sub nsw i64 %54, %3
%62 = getelementptr inbounds double, double* %41, i64 %61
%63 = bitcast double* %62 to <2 x double>*
store <2 x double> %57, <2 x double>* %63
%64 = getelementptr inbounds double, double* %62, i64 2
%65 = bitcast double* %64 to <2 x double>*
store <2 x double> %60, <2 x double>* %65
%66 = or i64 %52, 4
%67 = add i64 %12, %66
%68 = getelementptr inbounds double, double* %38, i64 %67
%69 = bitcast double* %68 to <2 x double>*
%70 = load <2 x double>, <2 x double>* %69
%71 = getelementptr inbounds double, double* %68, i64 2
%72 = bitcast double* %71 to <2 x double>*
%73 = load <2 x double>, <2 x double>* %72
%74 = sub nsw i64 %67, %3
%75 = getelementptr inbounds double, double* %41, i64 %74
%76 = bitcast double* %75 to <2 x double>*
store <2 x double> %70, <2 x double>* %76
%77 = getelementptr inbounds double, double* %75, i64 2
%78 = bitcast double* %77 to <2 x double>*
store <2 x double> %73, <2 x double>* %78
%79 = add nuw i64 %52, 8
%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 double, double* %38, i64 %85
%87 = bitcast double* %86 to <2 x double>*
%88 = load <2 x double>, <2 x double>* %87
%89 = getelementptr inbounds double, double* %86, i64 2
%90 = bitcast double* %89 to <2 x double>*
%91 = load <2 x double>, <2 x double>* %90
%92 = sub nsw i64 %85, %3
%93 = getelementptr inbounds double, double* %41, i64 %92
%94 = bitcast double* %93 to <2 x double>*
store <2 x double> %88, <2 x double>* %94
%95 = getelementptr inbounds double, double* %93, i64 2
%96 = bitcast double* %95 to <2 x double>*
store <2 x double> %91, <2 x double>* %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 double, double* %38, i64 %101
%103 = load double, double* %102
%104 = sub nsw i64 %101, %3
%105 = getelementptr inbounds double, double* %41, i64 %104
store double %103, double* %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
}