void copy_dmatrix(double **m1, long nrl1, long nrh1, long ncl1, long nch1,
                  double **m2, long nrl2, long nrh2, long ncl2, long nch2) {
  int i, j, ro, co;
  ro = nrl2 - nrl1;
  co = ncl2 - ncl1;
  for (i = nrl1; i <= nrh1; i++) {
    for (j = ncl1; j <= nch1; j++) {
      m2[i + ro][j + co] = m1[i][j];
    }
  }
}


x86 assembly (clang, -O0):
<copy_dmatrix>:
push   %rbp
mov    %rsp,%rbp
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
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    0x10(%rbp),%rax
sub    -0x10(%rbp),%rax
mov    %eax,-0x3c(%rbp)
mov    0x20(%rbp),%rax
sub    -0x20(%rbp),%rax
mov    %eax,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    %eax,-0x34(%rbp)
movslq -0x34(%rbp),%rax
cmp    -0x18(%rbp),%rax
jg     11c1 <copy_dmatrix+0xc1>
mov    -0x20(%rbp),%rax
mov    %eax,-0x38(%rbp)
movslq -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jg     11ae <copy_dmatrix+0xae>
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
mov    -0x34(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x38(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    115e <copy_dmatrix+0x5e>
jmp    11b3 <copy_dmatrix+0xb3>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1149 <copy_dmatrix+0x49>
pop    %rbp
ret

x86 assembly (clang, -O2):
<copy_dmatrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
movslq %esi,%rax
cmp    %rdx,%rax
jle    111d <copy_dmatrix+0x1d>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    0x48(%rsp),%rbp
mov    0x38(%rsp),%rbx
sub    %ecx,%ebp
sub    %esi,%ebx
movslq %ecx,%r15
movslq %ebp,%r13
movslq %ebx,%rbx
lea    (%r15,%r13,1),%r11
lea    0x1(%r8),%rcx
lea    0x1(%r15),%r10
cmp    %r10,%rcx
cmovg  %rcx,%r10
lea    (%r10,%r13,1),%rcx
mov    %rcx,-0x10(%rsp)
mov    %r10,%rcx
mov    %r10,-0x8(%rsp)
sub    %r15,%r10
lea    -0x4(%r10),%rcx
mov    %rcx,-0x20(%rsp)
mov    %rcx,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %r10,%rcx
and    $0xfffffffffffffffc,%rcx
mov    %rcx,-0x30(%rsp)
add    %r15,%rcx
mov    %rcx,-0x38(%rsp)
mov    %rsi,%rcx
mov    %rsi,-0x28(%rsp)
and    $0xfffffffffffffffe,%rsi
mov    %rsi,-0x40(%rsp)
lea    (%r15,%r13,1),%rcx
add    $0x6,%rcx
mov    %rcx,-0x48(%rsp)
lea    0x6(%r15),%rcx
mov    %rcx,-0x50(%rsp)
mov    %r11,-0x18(%rsp)
jmp    11cd <copy_dmatrix+0xcd>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%rax
cmp    %rdx,%rax
jg     1112 <copy_dmatrix+0x12>
cmp    %r8,%r15
jg     11c0 <copy_dmatrix+0xc0>
lea    (%rax,%rbx,1),%rcx
mov    (%rdi,%rax,8),%r12
mov    (%r9,%rcx,8),%r14
mov    %r15,%rcx
cmp    $0x4,%r10
jb     12c0 <copy_dmatrix+0x1c0>
lea    (%r14,%r11,8),%rcx
mov    -0x8(%rsp),%rsi
lea    (%r12,%rsi,8),%rbp
cmp    %rbp,%rcx
jae    1216 <copy_dmatrix+0x116>
mov    -0x10(%rsp),%rcx
lea    (%r14,%rcx,8),%rbp
lea    (%r12,%r15,8),%rsi
mov    %r15,%rcx
cmp    %rbp,%rsi
jb     12c0 <copy_dmatrix+0x1c0>
mov    %r13,%rsi
cmpq   $0x4,-0x20(%rsp)
jae    1226 <copy_dmatrix+0x126>
xor    %r11d,%r11d
jmp    127a <copy_dmatrix+0x17a>
mov    -0x48(%rsp),%rcx
lea    (%r14,%rcx,8),%rbp
mov    -0x50(%rsp),%rcx
lea    (%r12,%rcx,8),%r13
mov    -0x40(%rsp),%rcx
xor    %r11d,%r11d
movups -0x30(%r13,%r11,8),%xmm0
movups -0x20(%r13,%r11,8),%xmm1
movups %xmm0,-0x30(%rbp,%r11,8)
movups %xmm1,-0x20(%rbp,%r11,8)
movups -0x10(%r13,%r11,8),%xmm0
movups 0x0(%r13,%r11,8),%xmm1
movups %xmm0,-0x10(%rbp,%r11,8)
movups %xmm1,0x0(%rbp,%r11,8)
add    $0x8,%r11
add    $0xfffffffffffffffe,%rcx
jne    1240 <copy_dmatrix+0x140>
testb  $0x1,-0x28(%rsp)
mov    %rsi,%r13
je     12a0 <copy_dmatrix+0x1a0>
add    %r15,%r11
movups (%r12,%r11,8),%xmm0
movups 0x10(%r12,%r11,8),%xmm1
add    %r13,%r11
movups %xmm0,(%r14,%r11,8)
movups %xmm1,0x10(%r14,%r11,8)
mov    -0x38(%rsp),%rcx
cmp    -0x30(%rsp),%r10
mov    -0x18(%rsp),%r11
je     11c0 <copy_dmatrix+0xc0>
cs nopw 0x0(%rax,%rax,1)
nop
lea    (%r14,%r13,8),%rbp
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%r12,%rcx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rcx,8)
add    $0x1,%rcx
cmp    %r8,%rcx
jle    12d0 <copy_dmatrix+0x1d0>
jmp    11c0 <copy_dmatrix+0xc0>

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

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