void Array2D_double_copy(int M, int N, double **B, double **A) {
  int remainder = N & 3;
  int i = 0;
  int j = 0;
  for (i = 0; i < M; i++) {
    double *Bi = B[i];
    double *Ai = A[i];
    for (j = 0; j < remainder; j++)
      Bi[j] = Ai[j];
    for (j = remainder; j < N; j += 4) {
      Bi[j] = Ai[j];
      Bi[j + 1] = Ai[j + 1];
      Bi[j + 2] = Ai[j + 2];
      Bi[j + 3] = Ai[j + 3];
    }
  }
}


x86 assembly (clang, -O0):
<Array2D_double_copy>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    -0x8(%rbp),%eax
and    $0x3,%eax
mov    %eax,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1250 <Array2D_double_copy+0x150>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1197 <Array2D_double_copy+0x97>
mov    -0x38(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1163 <Array2D_double_copy+0x63>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    123d <Array2D_double_copy+0x13d>
mov    -0x38(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x24(%rbp)
jmp    119d <Array2D_double_copy+0x9d>
jmp    1242 <Array2D_double_copy+0x142>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1130 <Array2D_double_copy+0x30>
pop    %rbp
ret

x86 assembly (clang, -O2):
<Array2D_double_copy>:
push   %rbx
test   %edi,%edi
jle    11bc <Array2D_double_copy+0xbc>
mov    %esi,%r8d
and    $0x3,%r8d
movslq %esi,%r11
mov    %edi,%r9d
xor    %r10d,%r10d
jmp    112d <Array2D_double_copy+0x2d>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r10
cmp    %r9,%r10
je     11bc <Array2D_double_copy+0xbc>
mov    (%rdx,%r10,8),%rdi
mov    (%rcx,%r10,8),%rax
test   %r8d,%r8d
je     1170 <Array2D_double_copy+0x70>
movsd  (%rax),%xmm0
movsd  %xmm0,(%rdi)
cmp    $0x1,%r8d
je     1170 <Array2D_double_copy+0x70>
movsd  0x8(%rax),%xmm0
movsd  %xmm0,0x8(%rdi)
cmp    $0x2,%r8d
je     1170 <Array2D_double_copy+0x70>
movsd  0x10(%rax),%xmm0
movsd  %xmm0,0x10(%rdi)
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %r8,%rbx
cmp    $0x4,%esi
jl     1120 <Array2D_double_copy+0x20>
nopl   0x0(%rax,%rax,1)
movsd  (%rax,%rbx,8),%xmm0
movsd  %xmm0,(%rdi,%rbx,8)
movsd  0x8(%rax,%rbx,8),%xmm0
movsd  %xmm0,0x8(%rdi,%rbx,8)
movsd  0x10(%rax,%rbx,8),%xmm0
movsd  %xmm0,0x10(%rdi,%rbx,8)
movsd  0x18(%rax,%rbx,8),%xmm0
movsd  %xmm0,0x18(%rdi,%rbx,8)
add    $0x4,%rbx
cmp    %r11,%rbx
jl     1180 <Array2D_double_copy+0x80>
jmp    1120 <Array2D_double_copy+0x20>
pop    %rbx
ret

LLVM IR (-O0):
define void @Array2D_double_copy(i32 %0, i32 %1, double** %2, double** %3) {
%5 = alloca i32
%6 = alloca i32
%7 = alloca double**
%8 = alloca double**
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca double*
%13 = alloca double*
store i32 %0, i32* %5
store i32 %1, i32* %6
store double** %2, double*** %7
store double** %3, double*** %8
%14 = load i32, i32* %6
%15 = and i32 %14, 3
store i32 %15, i32* %9
store i32 0, i32* %10
store i32 0, i32* %11
store i32 0, i32* %10
br label %16
16: ; preds = %101, %4
%17 = load i32, i32* %10
%18 = load i32, i32* %5
%19 = icmp slt i32 %17, %18
br i1 %19, label %20, label %104
20: ; preds = %16
%21 = load double**, double*** %7
%22 = load i32, i32* %10
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double*, double** %21, i64 %23
%25 = load double*, double** %24
store double* %25, double** %12
%26 = load double**, double*** %8
%27 = load i32, i32* %10
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds double*, double** %26, i64 %28
%30 = load double*, double** %29
store double* %30, double** %13
store i32 0, i32* %11
br label %31
31: ; preds = %45, %20
%32 = load i32, i32* %11
%33 = load i32, i32* %9
%34 = icmp slt i32 %32, %33
br i1 %34, label %35, label %48
35: ; preds = %31
%36 = load double*, double** %13
%37 = load i32, i32* %11
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds double, double* %36, i64 %38
%40 = load double, double* %39
%41 = load double*, double** %12
%42 = load i32, i32* %11
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds double, double* %41, i64 %43
store double %40, double* %44
br label %45
45: ; preds = %35
%46 = load i32, i32* %11
%47 = add nsw i32 %46, 1
store i32 %47, i32* %11
br label %31
48: ; preds = %31
%49 = load i32, i32* %9
store i32 %49, i32* %11
br label %50
50: ; preds = %97, %48
%51 = load i32, i32* %11
%52 = load i32, i32* %6
%53 = icmp slt i32 %51, %52
br i1 %53, label %54, label %100
54: ; preds = %50
%55 = load double*, double** %13
%56 = load i32, i32* %11
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds double, double* %55, i64 %57
%59 = load double, double* %58
%60 = load double*, double** %12
%61 = load i32, i32* %11
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %60, i64 %62
store double %59, double* %63
%64 = load double*, double** %13
%65 = load i32, i32* %11
%66 = add nsw i32 %65, 1
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds double, double* %64, i64 %67
%69 = load double, double* %68
%70 = load double*, double** %12
%71 = load i32, i32* %11
%72 = add nsw i32 %71, 1
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %70, i64 %73
store double %69, double* %74
%75 = load double*, double** %13
%76 = load i32, i32* %11
%77 = add nsw i32 %76, 2
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds double, double* %75, i64 %78
%80 = load double, double* %79
%81 = load double*, double** %12
%82 = load i32, i32* %11
%83 = add nsw i32 %82, 2
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %81, i64 %84
store double %80, double* %85
%86 = load double*, double** %13
%87 = load i32, i32* %11
%88 = add nsw i32 %87, 3
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds double, double* %86, i64 %89
%91 = load double, double* %90
%92 = load double*, double** %12
%93 = load i32, i32* %11
%94 = add nsw i32 %93, 3
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds double, double* %92, i64 %95
store double %91, double* %96
br label %97
97: ; preds = %54
%98 = load i32, i32* %11
%99 = add nsw i32 %98, 4
store i32 %99, i32* %11
br label %50
100: ; preds = %50
br label %101
101: ; preds = %100
%102 = load i32, i32* %10
%103 = add nsw i32 %102, 1
store i32 %103, i32* %10
br label %16
104: ; preds = %16
ret void
}

LLVM IR (-O2):
define void @Array2D_double_copy(i32 %0, i32 %1, double** nocapture %2, double** nocapture %3) {
%5 = and i32 %1, 3
%6 = icmp sgt i32 %0, 0
br i1 %6, label %7, label %55
7: ; preds = %4
%8 = icmp eq i32 %5, 0
%9 = icmp sgt i32 %1, 3
%10 = and i32 %1, 3
%11 = zext i32 %10 to i64
%12 = sext i32 %1 to i64
%13 = zext i32 %0 to i64
%14 = icmp eq i32 %5, 1
%15 = icmp eq i32 %5, 2
br label %16
16: ; preds = %7, %52
%17 = phi i64 [ 0, %7 ], [ %53, %52 ]
%18 = getelementptr inbounds double*, double** %2, i64 %17
%19 = load double*, double** %18
%20 = getelementptr inbounds double*, double** %3, i64 %17
%21 = load double*, double** %20
br i1 %8, label %22, label %23
22: ; preds = %23, %25, %29, %16
br i1 %9, label %33, label %52
23: ; preds = %16
%24 = load double, double* %21
store double %24, double* %19
br i1 %14, label %22, label %25
25: ; preds = %23
%26 = getelementptr inbounds double, double* %21, i64 1
%27 = load double, double* %26
%28 = getelementptr inbounds double, double* %19, i64 1
store double %27, double* %28
br i1 %15, label %22, label %29
29: ; preds = %25
%30 = getelementptr inbounds double, double* %21, i64 2
%31 = load double, double* %30
%32 = getelementptr inbounds double, double* %19, i64 2
store double %31, double* %32
br label %22
33: ; preds = %22, %33
%34 = phi i64 [ %50, %33 ], [ %11, %22 ]
%35 = getelementptr inbounds double, double* %21, i64 %34
%36 = load double, double* %35
%37 = getelementptr inbounds double, double* %19, i64 %34
store double %36, double* %37
%38 = add nuw nsw i64 %34, 1
%39 = getelementptr inbounds double, double* %21, i64 %38
%40 = load double, double* %39
%41 = getelementptr inbounds double, double* %19, i64 %38
store double %40, double* %41
%42 = add nuw nsw i64 %34, 2
%43 = getelementptr inbounds double, double* %21, i64 %42
%44 = load double, double* %43
%45 = getelementptr inbounds double, double* %19, i64 %42
store double %44, double* %45
%46 = add nuw nsw i64 %34, 3
%47 = getelementptr inbounds double, double* %21, i64 %46
%48 = load double, double* %47
%49 = getelementptr inbounds double, double* %19, i64 %46
store double %48, double* %49
%50 = add nuw nsw i64 %34, 4
%51 = icmp slt i64 %50, %12
br i1 %51, label %33, label %52
52: ; preds = %33, %22
%53 = add nuw nsw i64 %17, 1
%54 = icmp eq i64 %53, %13
br i1 %54, label %55, label %16
55: ; preds = %52, %4
ret void
}