void d_copy_pmat(int row, int col, int bs_dummy, double *A, int sda, double *B,
                 int sdb) {
  const int bs = 4;
  int i, ii, j, row2;
  ii = 0;
  for (; ii < row - 3; ii += bs) {
    for (j = 0; j < col; j++) {
      B[0 + j * bs + ii * sdb] = A[0 + j * bs + ii * sda];
      B[1 + j * bs + ii * sdb] = A[1 + j * bs + ii * sda];
      B[2 + j * bs + ii * sdb] = A[2 + j * bs + ii * sda];
      B[3 + j * bs + ii * sdb] = A[3 + j * bs + ii * sda];
    }
  }
  if (ii < row) {
    row2 = row - ii;
    for (j = 0; j < col; j++) {
      for (i = 0; i < row2; i++) {
        B[i + j * bs + ii * sdb] = A[i + j * bs + ii * sda];
      }
    }
  }
}


x86 assembly (clang, -O0):
<d_copy_pmat>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%eax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9,-0x28(%rbp)
movl   $0x4,-0x2c(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x3,%ecx
cmp    %ecx,%eax
jge    125f <d_copy_pmat+0x15f>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    124c <d_copy_pmat+0x14c>
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
shl    $0x2,%ecx
add    $0x0,%ecx
mov    -0x34(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
shl    $0x2,%ecx
add    $0x0,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
shl    $0x2,%ecx
add    $0x1,%ecx
mov    -0x34(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
shl    $0x2,%ecx
add    $0x1,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
shl    $0x2,%ecx
add    $0x2,%ecx
mov    -0x34(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
shl    $0x2,%ecx
add    $0x2,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
mov    -0x34(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1142 <d_copy_pmat+0x42>
jmp    1251 <d_copy_pmat+0x151>
mov    -0x34(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x34(%rbp)
jmp    112a <d_copy_pmat+0x2a>
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1300 <d_copy_pmat+0x200>
mov    -0x4(%rbp),%eax
sub    -0x34(%rbp),%eax
mov    %eax,-0x3c(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    12fb <d_copy_pmat+0x1fb>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    12e8 <d_copy_pmat+0x1e8>
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
mov    -0x38(%rbp),%edx
shl    $0x2,%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%ecx
mov    -0x38(%rbp),%edx
shl    $0x2,%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    128e <d_copy_pmat+0x18e>
jmp    12ed <d_copy_pmat+0x1ed>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    127b <d_copy_pmat+0x17b>
jmp    1300 <d_copy_pmat+0x200>
pop    %rbp
ret

x86 assembly (clang, -O2):
<d_copy_pmat>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    0x38(%rsp),%ebp
xor    %eax,%eax
cmp    $0x4,%edi
mov    %r9,-0x68(%rsp)
mov    %rcx,-0x70(%rsp)
jl     138f <d_copy_pmat+0x28f>
mov    %rdi,-0x10(%rsp)
lea    -0x3(%rdi),%edx
mov    %r8d,-0x14(%rsp)
movslq %r8d,%r11
movslq %ebp,%r14
mov    %esi,%r12d
lea    0x0(,%r14,4),%rax
mov    %rax,-0x20(%rsp)
lea    0x0(,%r12,4),%r10
lea    0x0(,%r11,4),%rax
mov    %rax,-0x28(%rsp)
mov    %r12d,%eax
and    $0xfffffffe,%eax
mov    %rax,-0x38(%rsp)
add    $0xfffffffffffffffe,%rax
mov    %rax,-0x40(%rsp)
mov    %rax,%rdi
shr    %rdi
add    $0x1,%rdi
mov    %rdi,%rax
mov    %rdi,-0x48(%rsp)
and    $0xfffffffffffffffe,%rdi
mov    %rdi,-0x50(%rsp)
mov    %r14,-0x8(%rsp)
shl    $0x5,%r14
mov    %r11,-0x58(%rsp)
shl    $0x5,%r11
xor    %r8d,%r8d
mov    %rcx,%rdi
mov    %r9,%r15
xor    %eax,%eax
mov    %r12,-0x30(%rsp)
mov    %rdx,-0x60(%rsp)
jmp    11d7 <d_copy_pmat+0xd7>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x4,%rax
add    $0x1,%r8
add    %r14,%r15
add    %r11,%rdi
cmp    %rdx,%rax
jae    1381 <d_copy_pmat+0x281>
test   %esi,%esi
jle    11c0 <d_copy_pmat+0xc0>
cmp    $0x1,%esi
jne    11f0 <d_copy_pmat+0xf0>
xor    %ebp,%ebp
jmp    1330 <d_copy_pmat+0x230>
nopw   0x0(%rax,%rax,1)
mov    -0x20(%rsp),%rbp
imul   %r8,%rbp
mov    -0x28(%rsp),%rbx
imul   %r8,%rbx
lea    (%r10,%rbx,1),%rdx
mov    %r10,%r13
lea    (%r9,%rbp,8),%r10
lea    (%rcx,%rdx,8),%rdx
cmp    %rdx,%r10
jae    1235 <d_copy_pmat+0x135>
mov    %r13,%r10
add    %r13,%rbp
lea    (%r9,%rbp,8),%rdx
lea    (%rcx,%rbx,8),%rbp
cmp    %rdx,%rbp
jae    1235 <d_copy_pmat+0x135>
xor    %ebp,%ebp
mov    -0x60(%rsp),%rdx
jmp    1330 <d_copy_pmat+0x230>
cmpq   $0x0,-0x40(%rsp)
je     12b6 <d_copy_pmat+0x1b6>
mov    $0x20,%ebp
mov    -0x50(%rsp),%r12
xor    %r10d,%r10d
nopw   0x0(%rax,%rax,1)
movups -0x40(%rdi,%rbp,2),%xmm0
movups -0x30(%rdi,%rbp,2),%xmm1
movups -0x20(%rdi,%rbp,2),%xmm2
movups -0x10(%rdi,%rbp,2),%xmm3
movups %xmm3,-0x10(%r15,%rbp,2)
movups %xmm0,-0x40(%r15,%rbp,2)
movups %xmm1,-0x30(%r15,%rbp,2)
movups %xmm2,-0x20(%r15,%rbp,2)
movups (%rdi,%rbp,2),%xmm0
movups 0x10(%rdi,%rbp,2),%xmm1
movups 0x20(%rdi,%rbp,2),%xmm2
movups 0x30(%rdi,%rbp,2),%xmm3
movups %xmm2,0x20(%r15,%rbp,2)
movups %xmm3,0x30(%r15,%rbp,2)
movups %xmm0,(%r15,%rbp,2)
movups %xmm1,0x10(%r15,%rbp,2)
add    $0x4,%r10
add    $0x40,%rbp
add    $0xfffffffffffffffe,%r12
jne    1250 <d_copy_pmat+0x150>
jmp    12b9 <d_copy_pmat+0x1b9>
xor    %r10d,%r10d
testb  $0x1,-0x48(%rsp)
mov    -0x68(%rsp),%r9
mov    -0x70(%rsp),%rcx
je     130e <d_copy_pmat+0x20e>
mov    %rax,%rdx
imul   -0x58(%rsp),%rdx
mov    %rax,%rbx
imul   -0x8(%rsp),%rbx
lea    (%rdx,%r10,4),%rdx
movups (%rcx,%rdx,8),%xmm0
movups 0x10(%rcx,%rdx,8),%xmm1
movups 0x20(%rcx,%rdx,8),%xmm2
movups 0x30(%rcx,%rdx,8),%xmm3
lea    (%rbx,%r10,4),%rdx
movups %xmm2,0x20(%r9,%rdx,8)
movups %xmm3,0x30(%r9,%rdx,8)
movups %xmm0,(%r9,%rdx,8)
movups %xmm1,0x10(%r9,%rdx,8)
mov    -0x38(%rsp),%rdx
mov    %rdx,%rbp
mov    -0x30(%rsp),%r12
cmp    %r12,%rdx
mov    %r13,%r10
mov    -0x60(%rsp),%rdx
je     11c0 <d_copy_pmat+0xc0>
nopl   0x0(%rax)
mov    %r12,%rbx
sub    %rbp,%rbx
shl    $0x5,%rbp
or     $0x10,%rbp
xchg   %ax,%ax
movsd  -0x10(%rdi,%rbp,1),%xmm0
movsd  %xmm0,-0x10(%r15,%rbp,1)
movsd  -0x8(%rdi,%rbp,1),%xmm0
movsd  %xmm0,-0x8(%r15,%rbp,1)
movsd  (%rdi,%rbp,1),%xmm0
movsd  %xmm0,(%r15,%rbp,1)
movsd  0x8(%rdi,%rbp,1),%xmm0
movsd  %xmm0,0x8(%r15,%rbp,1)
add    $0x20,%rbp
add    $0xffffffffffffffff,%rbx
jne    1340 <d_copy_pmat+0x240>
jmp    11c0 <d_copy_pmat+0xc0>
mov    -0x10(%rsp),%rdi
mov    -0x14(%rsp),%r8d
mov    0x38(%rsp),%ebp
cmp    %edi,%eax
jge    15dc <d_copy_pmat+0x4dc>
test   %esi,%esi
jle    15dc <d_copy_pmat+0x4dc>
sub    %eax,%edi
imul   %eax,%r8d
imul   %eax,%ebp
movslq %r8d,%r12
movslq %ebp,%r15
mov    %esi,%r13d
mov    %edi,%ebp
mov    -0x68(%rsp),%rax
lea    (%rax,%rbp,8),%rcx
mov    %rcx,-0x28(%rsp)
mov    -0x70(%rsp),%rcx
lea    (%rcx,%rbp,8),%r10
mov    %ebp,%edx
and    $0xfffffffc,%edx
mov    %rdx,-0x30(%rsp)
add    $0xfffffffffffffffc,%rdx
mov    %rdx,-0x38(%rsp)
mov    %rdx,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %rsi,%rdx
mov    %rsi,-0x40(%rsp)
and    $0xfffffffffffffffe,%rsi
mov    %rsi,-0x48(%rsp)
mov    %rbp,%r14
neg    %r14
lea    0x30(%rax),%rdx
mov    %rdx,-0x50(%rsp)
lea    0x30(%rcx),%rdx
mov    %rdx,-0x58(%rsp)
lea    0x8(%rax,%r15,8),%rax
mov    %rax,-0x20(%rsp)
lea    (%rcx,%r12,8),%rax
add    $0x8,%rax
mov    %rax,-0x60(%rsp)
xor    %r11d,%r11d
xor    %r8d,%r8d
jmp    1441 <d_copy_pmat+0x341>
add    $0x1,%r8
add    $0x4,%r11d
cmp    %r13,%r8
je     15dc <d_copy_pmat+0x4dc>
test   %edi,%edi
jle    1430 <d_copy_pmat+0x330>
movslq %r11d,%r9
lea    0x0(,%r8,4),%eax
movslq %eax,%rcx
cmp    $0x4,%edi
jb     1496 <d_copy_pmat+0x396>
lea    0x0(,%r8,4),%eax
cltq
lea    (%r15,%rax,1),%rdx
add    %r12,%rax
mov    -0x68(%rsp),%rsi
lea    (%rsi,%rdx,8),%rbx
lea    (%r10,%rax,8),%rsi
cmp    %rsi,%rbx
jae    1514 <d_copy_pmat+0x414>
mov    -0x28(%rsp),%rsi
lea    (%rsi,%rdx,8),%rdx
mov    -0x70(%rsp),%rsi
lea    (%rsi,%rax,8),%rax
cmp    %rdx,%rax
jae    1514 <d_copy_pmat+0x414>
xor    %eax,%eax
mov    %rax,%rbx
test   $0x1,%bpl
je     14c6 <d_copy_pmat+0x3c6>
add    %rax,%rcx
lea    (%rcx,%r12,1),%rdx
mov    -0x70(%rsp),%rsi
movsd  (%rsi,%rdx,8),%xmm0
add    %r15,%rcx
mov    -0x68(%rsp),%rdx
movsd  %xmm0,(%rdx,%rcx,8)
mov    %rax,%rbx
or     $0x1,%rbx
not    %rax
cmp    %r14,%rax
je     1430 <d_copy_pmat+0x330>
mov    -0x20(%rsp),%rax
lea    (%rax,%r9,8),%rax
mov    -0x60(%rsp),%rcx
lea    (%rcx,%r9,8),%rcx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  -0x8(%rcx,%rbx,8),%xmm0
movsd  %xmm0,-0x8(%rax,%rbx,8)
movsd  (%rcx,%rbx,8),%xmm0
movsd  %xmm0,(%rax,%rbx,8)
add    $0x2,%rbx
cmp    %rbx,%rbp
jne    14f0 <d_copy_pmat+0x3f0>
jmp    1430 <d_copy_pmat+0x330>
mov    %r10,%rsi
cmpq   $0x0,-0x38(%rsp)
je     15d0 <d_copy_pmat+0x4d0>
lea    (%r15,%r9,1),%rax
lea    (%r12,%r9,1),%rdx
mov    -0x50(%rsp),%rbx
lea    (%rbx,%rax,8),%rax
mov    -0x58(%rsp),%rbx
lea    (%rbx,%rdx,8),%rbx
mov    -0x48(%rsp),%rdx
xor    %r10d,%r10d
cs nopw 0x0(%rax,%rax,1)
nop
movups -0x30(%rbx,%r10,8),%xmm0
movups -0x20(%rbx,%r10,8),%xmm1
movups %xmm0,-0x30(%rax,%r10,8)
movups %xmm1,-0x20(%rax,%r10,8)
movups -0x10(%rbx,%r10,8),%xmm0
movups (%rbx,%r10,8),%xmm1
movups %xmm0,-0x10(%rax,%r10,8)
movups %xmm1,(%rax,%r10,8)
add    $0x8,%r10
add    $0xfffffffffffffffe,%rdx
jne    1550 <d_copy_pmat+0x450>
testb  $0x1,-0x40(%rsp)
je     15b7 <d_copy_pmat+0x4b7>
add    %rcx,%r10
lea    (%r10,%r12,1),%rax
mov    -0x70(%rsp),%rdx
movups (%rdx,%rax,8),%xmm0
movups 0x10(%rdx,%rax,8),%xmm1
add    %r15,%r10
mov    -0x68(%rsp),%rax
movups %xmm0,(%rax,%r10,8)
movups %xmm1,0x10(%rax,%r10,8)
mov    -0x30(%rsp),%rdx
mov    %rdx,%rax
cmp    %rbp,%rdx
mov    %rsi,%r10
je     1430 <d_copy_pmat+0x330>
jmp    1498 <d_copy_pmat+0x398>
xor    %r10d,%r10d
testb  $0x1,-0x40(%rsp)
jne    158f <d_copy_pmat+0x48f>
jmp    15b7 <d_copy_pmat+0x4b7>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @d_copy_pmat(i32 %0, i32 %1, i32 %2, double* %3, i32 %4, double* %5, i32 %6) {
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca double*
%12 = alloca i32
%13 = alloca double*
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
store i32 %0, i32* %8
store i32 %1, i32* %9
store i32 %2, i32* %10
store double* %3, double** %11
store i32 %4, i32* %12
store double* %5, double** %13
store i32 %6, i32* %14
store i32 4, i32* %15
store i32 0, i32* %17
br label %20
20: ; preds = %119, %7
%21 = load i32, i32* %17
%22 = load i32, i32* %8
%23 = sub nsw i32 %22, 3
%24 = icmp slt i32 %21, %23
br i1 %24, label %25, label %122
25: ; preds = %20
store i32 0, i32* %18
br label %26
26: ; preds = %115, %25
%27 = load i32, i32* %18
%28 = load i32, i32* %9
%29 = icmp slt i32 %27, %28
br i1 %29, label %30, label %118
30: ; preds = %26
%31 = load double*, double** %11
%32 = load i32, i32* %18
%33 = mul nsw i32 %32, 4
%34 = add nsw i32 0, %33
%35 = load i32, i32* %17
%36 = load i32, i32* %12
%37 = mul nsw i32 %35, %36
%38 = add nsw i32 %34, %37
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds double, double* %31, i64 %39
%41 = load double, double* %40
%42 = load double*, double** %13
%43 = load i32, i32* %18
%44 = mul nsw i32 %43, 4
%45 = add nsw i32 0, %44
%46 = load i32, i32* %17
%47 = load i32, i32* %14
%48 = mul nsw i32 %46, %47
%49 = add nsw i32 %45, %48
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double, double* %42, i64 %50
store double %41, double* %51
%52 = load double*, double** %11
%53 = load i32, i32* %18
%54 = mul nsw i32 %53, 4
%55 = add nsw i32 1, %54
%56 = load i32, i32* %17
%57 = load i32, i32* %12
%58 = mul nsw i32 %56, %57
%59 = add nsw i32 %55, %58
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %52, i64 %60
%62 = load double, double* %61
%63 = load double*, double** %13
%64 = load i32, i32* %18
%65 = mul nsw i32 %64, 4
%66 = add nsw i32 1, %65
%67 = load i32, i32* %17
%68 = load i32, i32* %14
%69 = mul nsw i32 %67, %68
%70 = add nsw i32 %66, %69
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds double, double* %63, i64 %71
store double %62, double* %72
%73 = load double*, double** %11
%74 = load i32, i32* %18
%75 = mul nsw i32 %74, 4
%76 = add nsw i32 2, %75
%77 = load i32, i32* %17
%78 = load i32, i32* %12
%79 = mul nsw i32 %77, %78
%80 = add nsw i32 %76, %79
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %73, i64 %81
%83 = load double, double* %82
%84 = load double*, double** %13
%85 = load i32, i32* %18
%86 = mul nsw i32 %85, 4
%87 = add nsw i32 2, %86
%88 = load i32, i32* %17
%89 = load i32, i32* %14
%90 = mul nsw i32 %88, %89
%91 = add nsw i32 %87, %90
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds double, double* %84, i64 %92
store double %83, double* %93
%94 = load double*, double** %11
%95 = load i32, i32* %18
%96 = mul nsw i32 %95, 4
%97 = add nsw i32 3, %96
%98 = load i32, i32* %17
%99 = load i32, i32* %12
%100 = mul nsw i32 %98, %99
%101 = add nsw i32 %97, %100
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds double, double* %94, i64 %102
%104 = load double, double* %103
%105 = load double*, double** %13
%106 = load i32, i32* %18
%107 = mul nsw i32 %106, 4
%108 = add nsw i32 3, %107
%109 = load i32, i32* %17
%110 = load i32, i32* %14
%111 = mul nsw i32 %109, %110
%112 = add nsw i32 %108, %111
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds double, double* %105, i64 %113
store double %104, double* %114
br label %115
115: ; preds = %30
%116 = load i32, i32* %18
%117 = add nsw i32 %116, 1
store i32 %117, i32* %18
br label %26
118: ; preds = %26
br label %119
119: ; preds = %118
%120 = load i32, i32* %17
%121 = add nsw i32 %120, 4
store i32 %121, i32* %17
br label %20
122: ; preds = %20
%123 = load i32, i32* %17
%124 = load i32, i32* %8
%125 = icmp slt i32 %123, %124
br i1 %125, label %126, label %171
126: ; preds = %122
%127 = load i32, i32* %8
%128 = load i32, i32* %17
%129 = sub nsw i32 %127, %128
store i32 %129, i32* %19
store i32 0, i32* %18
br label %130
130: ; preds = %167, %126
%131 = load i32, i32* %18
%132 = load i32, i32* %9
%133 = icmp slt i32 %131, %132
br i1 %133, label %134, label %170
134: ; preds = %130
store i32 0, i32* %16
br label %135
135: ; preds = %163, %134
%136 = load i32, i32* %16
%137 = load i32, i32* %19
%138 = icmp slt i32 %136, %137
br i1 %138, label %139, label %166
139: ; preds = %135
%140 = load double*, double** %11
%141 = load i32, i32* %16
%142 = load i32, i32* %18
%143 = mul nsw i32 %142, 4
%144 = add nsw i32 %141, %143
%145 = load i32, i32* %17
%146 = load i32, i32* %12
%147 = mul nsw i32 %145, %146
%148 = add nsw i32 %144, %147
%149 = sext i32 %148 to i64
%150 = getelementptr inbounds double, double* %140, i64 %149
%151 = load double, double* %150
%152 = load double*, double** %13
%153 = load i32, i32* %16
%154 = load i32, i32* %18
%155 = mul nsw i32 %154, 4
%156 = add nsw i32 %153, %155
%157 = load i32, i32* %17
%158 = load i32, i32* %14
%159 = mul nsw i32 %157, %158
%160 = add nsw i32 %156, %159
%161 = sext i32 %160 to i64
%162 = getelementptr inbounds double, double* %152, i64 %161
store double %151, double* %162
br label %163
163: ; preds = %139
%164 = load i32, i32* %16
%165 = add nsw i32 %164, 1
store i32 %165, i32* %16
br label %135
166: ; preds = %135
br label %167
167: ; preds = %166
%168 = load i32, i32* %18
%169 = add nsw i32 %168, 1
store i32 %169, i32* %18
br label %130
170: ; preds = %130
br label %171
171: ; preds = %170, %122
ret void
}

LLVM IR (-O2):
define void @d_copy_pmat(i32 %0, i32 %1, i32 %2, double* nocapture %3, i32 %4, double* nocapture writeonly %5, i32 %6) {
%8 = icmp sgt i32 %0, 3
br i1 %8, label %9, label %126
9: ; preds = %7
%10 = add nsw i32 %0, -3
%11 = icmp sgt i32 %1, 0
%12 = sext i32 %4 to i64
%13 = sext i32 %6 to i64
%14 = zext i32 %10 to i64
%15 = zext i32 %1 to i64
%16 = shl nsw i64 %13, 2
%17 = shl nuw nsw i64 %15, 2
%18 = shl nsw i64 %12, 2
%19 = and i64 %15, 4294967294
%20 = add nsw i64 %19, -2
%21 = lshr exact i64 %20, 1
%22 = add nuw i64 %21, 1
%23 = icmp eq i32 %1, 1
%24 = and i64 %15, 4294967294
%25 = and i64 %22, 1
%26 = icmp eq i64 %20, 0
%27 = and i64 %22, -2
%28 = getelementptr inbounds double, double* %5, i64 -3
%29 = getelementptr inbounds double, double* %5, i64 -3
%30 = icmp eq i64 %25, 0
%31 = getelementptr inbounds double, double* %5, i64 -3
%32 = icmp eq i64 %24, %15
br label %33
33: ; preds = %9, %120
%34 = phi i64 [ 0, %9 ], [ %123, %120 ]
%35 = phi i64 [ 0, %9 ], [ %121, %120 ]
%36 = mul i64 %16, %34
%37 = getelementptr double, double* %5, i64 %36
%38 = add i64 %17, %36
%39 = getelementptr double, double* %5, i64 %38
%40 = mul i64 %18, %34
%41 = getelementptr double, double* %3, i64 %40
%42 = add i64 %17, %40
%43 = getelementptr double, double* %3, i64 %42
br i1 %11, label %44, label %120
44: ; preds = %33
%45 = mul nsw i64 %35, %12
%46 = mul nsw i64 %35, %13
br i1 %23, label %90, label %47
47: ; preds = %44
%48 = icmp ult double* %37, %43
%49 = icmp ult double* %41, %39
%50 = and i1 %48, %49
br i1 %50, label %90, label %51
51: ; preds = %47
br i1 %26, label %77, label %52
52: ; preds = %51, %52
%53 = phi i64 [ %74, %52 ], [ 0, %51 ]
%54 = phi i64 [ %75, %52 ], [ 0, %51 ]
%55 = shl nsw i64 %53, 2
%56 = add nsw i64 %55, %45
%57 = getelementptr inbounds double, double* %3, i64 %56
%58 = bitcast double* %57 to <8 x double>*
%59 = load <8 x double>, <8 x double>* %58
%60 = or i64 %55, 3
%61 = add nsw i64 %60, %46
%62 = getelementptr inbounds double, double* %28, i64 %61
%63 = bitcast double* %62 to <8 x double>*
store <8 x double> %59, <8 x double>* %63
%64 = shl i64 %53, 2
%65 = or i64 %64, 8
%66 = add nsw i64 %65, %45
%67 = getelementptr inbounds double, double* %3, i64 %66
%68 = bitcast double* %67 to <8 x double>*
%69 = load <8 x double>, <8 x double>* %68
%70 = or i64 %64, 11
%71 = add nsw i64 %70, %46
%72 = getelementptr inbounds double, double* %29, i64 %71
%73 = bitcast double* %72 to <8 x double>*
store <8 x double> %69, <8 x double>* %73
%74 = add nuw i64 %53, 4
%75 = add nuw i64 %54, 2
%76 = icmp eq i64 %75, %27
br i1 %76, label %77, label %52
77: ; preds = %52, %51
%78 = phi i64 [ 0, %51 ], [ %74, %52 ]
br i1 %30, label %89, label %79
79: ; preds = %77
%80 = shl nsw i64 %78, 2
%81 = add nsw i64 %80, %45
%82 = getelementptr inbounds double, double* %3, i64 %81
%83 = bitcast double* %82 to <8 x double>*
%84 = load <8 x double>, <8 x double>* %83
%85 = or i64 %80, 3
%86 = add nsw i64 %85, %46
%87 = getelementptr inbounds double, double* %31, i64 %86
%88 = bitcast double* %87 to <8 x double>*
store <8 x double> %84, <8 x double>* %88
br label %89
89: ; preds = %77, %79
br i1 %32, label %120, label %90
90: ; preds = %47, %44, %89
%91 = phi i64 [ 0, %47 ], [ 0, %44 ], [ %24, %89 ]
br label %92
92: ; preds = %90, %92
%93 = phi i64 [ %118, %92 ], [ %91, %90 ]
%94 = shl nsw i64 %93, 2
%95 = add nsw i64 %94, %45
%96 = getelementptr inbounds double, double* %3, i64 %95
%97 = load double, double* %96
%98 = add nsw i64 %94, %46
%99 = getelementptr inbounds double, double* %5, i64 %98
store double %97, double* %99
%100 = or i64 %94, 1
%101 = add nsw i64 %100, %45
%102 = getelementptr inbounds double, double* %3, i64 %101
%103 = load double, double* %102
%104 = add nsw i64 %100, %46
%105 = getelementptr inbounds double, double* %5, i64 %104
store double %103, double* %105
%106 = or i64 %94, 2
%107 = add nsw i64 %106, %45
%108 = getelementptr inbounds double, double* %3, i64 %107
%109 = load double, double* %108
%110 = add nsw i64 %106, %46
%111 = getelementptr inbounds double, double* %5, i64 %110
store double %109, double* %111
%112 = or i64 %94, 3
%113 = add nsw i64 %112, %45
%114 = getelementptr inbounds double, double* %3, i64 %113
%115 = load double, double* %114
%116 = add nsw i64 %112, %46
%117 = getelementptr inbounds double, double* %5, i64 %116
store double %115, double* %117
%118 = add nuw nsw i64 %93, 1
%119 = icmp eq i64 %118, %15
br i1 %119, label %120, label %92
120: ; preds = %92, %89, %33
%121 = add nuw nsw i64 %35, 4
%122 = icmp ult i64 %121, %14
%123 = add i64 %34, 1
br i1 %122, label %33, label %124
124: ; preds = %120
%125 = trunc i64 %121 to i32
br label %126
126: ; preds = %124, %7
%127 = phi i32 [ 0, %7 ], [ %125, %124 ]
%128 = icmp slt i32 %127, %0
br i1 %128, label %129, label %260
129: ; preds = %126
%130 = sub i32 %0, %127
%131 = icmp sgt i32 %1, 0
br i1 %131, label %132, label %260
132: ; preds = %129
%133 = icmp sgt i32 %130, 0
%134 = mul i32 %127, %4
%135 = mul i32 %127, %6
%136 = sext i32 %134 to i64
%137 = sext i32 %135 to i64
%138 = zext i32 %1 to i64
%139 = zext i32 %130 to i64
%140 = getelementptr double, double* %5, i64 %139
%141 = getelementptr double, double* %3, i64 %139
%142 = and i64 %139, 4294967292
%143 = add nsw i64 %142, -4
%144 = lshr exact i64 %143, 2
%145 = add nuw nsw i64 %144, 1
%146 = icmp ult i32 %130, 4
%147 = and i64 %139, 4294967292
%148 = and i64 %145, 1
%149 = icmp eq i64 %143, 0
%150 = and i64 %145, 9223372036854775806
%151 = icmp eq i64 %148, 0
%152 = icmp eq i64 %147, %139
%153 = and i64 %139, 1
%154 = icmp eq i64 %153, 0
%155 = sub nsw i64 0, %139
br label %156
156: ; preds = %132, %257
%157 = phi i64 [ 0, %132 ], [ %258, %257 ]
%158 = trunc i64 %157 to i32
%159 = shl i32 %158, 2
%160 = sext i32 %159 to i64
%161 = add nsw i64 %137, %160
%162 = getelementptr double, double* %5, i64 %161
%163 = getelementptr double, double* %140, i64 %161
%164 = add nsw i64 %136, %160
%165 = getelementptr double, double* %3, i64 %164
%166 = getelementptr double, double* %141, i64 %164
br i1 %133, label %167, label %257
167: ; preds = %156
%168 = trunc i64 %157 to i32
%169 = shl nsw i32 %168, 2
%170 = sext i32 %169 to i64
br i1 %146, label %226, label %171
171: ; preds = %167
%172 = icmp ult double* %162, %166
%173 = icmp ult double* %165, %163
%174 = and i1 %172, %173
br i1 %174, label %226, label %175
175: ; preds = %171
br i1 %149, label %209, label %176
176: ; preds = %175, %176
%177 = phi i64 [ %206, %176 ], [ 0, %175 ]
%178 = phi i64 [ %207, %176 ], [ 0, %175 ]
%179 = add nuw nsw i64 %177, %170
%180 = add nsw i64 %179, %136
%181 = getelementptr inbounds double, double* %3, i64 %180
%182 = bitcast double* %181 to <2 x double>*
%183 = load <2 x double>, <2 x double>* %182
%184 = getelementptr inbounds double, double* %181, i64 2
%185 = bitcast double* %184 to <2 x double>*
%186 = load <2 x double>, <2 x double>* %185
%187 = add nsw i64 %179, %137
%188 = getelementptr inbounds double, double* %5, i64 %187
%189 = bitcast double* %188 to <2 x double>*
store <2 x double> %183, <2 x double>* %189
%190 = getelementptr inbounds double, double* %188, i64 2
%191 = bitcast double* %190 to <2 x double>*
store <2 x double> %186, <2 x double>* %191
%192 = or i64 %177, 4
%193 = add nuw nsw i64 %192, %170
%194 = add nsw i64 %193, %136
%195 = getelementptr inbounds double, double* %3, i64 %194
%196 = bitcast double* %195 to <2 x double>*
%197 = load <2 x double>, <2 x double>* %196
%198 = getelementptr inbounds double, double* %195, i64 2
%199 = bitcast double* %198 to <2 x double>*
%200 = load <2 x double>, <2 x double>* %199
%201 = add nsw i64 %193, %137
%202 = getelementptr inbounds double, double* %5, i64 %201
%203 = bitcast double* %202 to <2 x double>*
store <2 x double> %197, <2 x double>* %203
%204 = getelementptr inbounds double, double* %202, i64 2
%205 = bitcast double* %204 to <2 x double>*
store <2 x double> %200, <2 x double>* %205
%206 = add nuw i64 %177, 8
%207 = add i64 %178, 2
%208 = icmp eq i64 %207, %150
br i1 %208, label %209, label %176
209: ; preds = %176, %175
%210 = phi i64 [ 0, %175 ], [ %206, %176 ]
br i1 %151, label %225, label %211
211: ; preds = %209
%212 = add nuw nsw i64 %210, %170
%213 = add nsw i64 %212, %136
%214 = getelementptr inbounds double, double* %3, i64 %213
%215 = bitcast double* %214 to <2 x double>*
%216 = load <2 x double>, <2 x double>* %215
%217 = getelementptr inbounds double, double* %214, i64 2
%218 = bitcast double* %217 to <2 x double>*
%219 = load <2 x double>, <2 x double>* %218
%220 = add nsw i64 %212, %137
%221 = getelementptr inbounds double, double* %5, i64 %220
%222 = bitcast double* %221 to <2 x double>*
store <2 x double> %216, <2 x double>* %222
%223 = getelementptr inbounds double, double* %221, i64 2
%224 = bitcast double* %223 to <2 x double>*
store <2 x double> %219, <2 x double>* %224
br label %225
225: ; preds = %209, %211
br i1 %152, label %257, label %226
226: ; preds = %171, %167, %225
%227 = phi i64 [ 0, %171 ], [ 0, %167 ], [ %147, %225 ]
%228 = xor i64 %227, -1
br i1 %154, label %237, label %229
229: ; preds = %226
%230 = add nuw nsw i64 %227, %170
%231 = add nsw i64 %230, %136
%232 = getelementptr inbounds double, double* %3, i64 %231
%233 = load double, double* %232
%234 = add nsw i64 %230, %137
%235 = getelementptr inbounds double, double* %5, i64 %234
store double %233, double* %235
%236 = or i64 %227, 1
br label %237
237: ; preds = %229, %226
%238 = phi i64 [ %227, %226 ], [ %236, %229 ]
%239 = icmp eq i64 %228, %155
br i1 %239, label %257, label %240
240: ; preds = %237, %240
%241 = phi i64 [ %255, %240 ], [ %238, %237 ]
%242 = add nuw nsw i64 %241, %170
%243 = add nsw i64 %242, %136
%244 = getelementptr inbounds double, double* %3, i64 %243
%245 = load double, double* %244
%246 = add nsw i64 %242, %137
%247 = getelementptr inbounds double, double* %5, i64 %246
store double %245, double* %247
%248 = add nuw nsw i64 %241, 1
%249 = add nuw nsw i64 %248, %170
%250 = add nsw i64 %249, %136
%251 = getelementptr inbounds double, double* %3, i64 %250
%252 = load double, double* %251
%253 = add nsw i64 %249, %137
%254 = getelementptr inbounds double, double* %5, i64 %253
store double %252, double* %254
%255 = add nuw nsw i64 %241, 2
%256 = icmp eq i64 %255, %139
br i1 %256, label %257, label %240
257: ; preds = %237, %240, %225, %156
%258 = add nuw nsw i64 %157, 1
%259 = icmp eq i64 %258, %138
br i1 %259, label %260, label %156
260: ; preds = %257, %129, %126
ret void
}