void d_copy_pmat_l(int row, int bs_dummy, double *A, int sda, double *B,
                   int sdb) {
  const int bs = 4;
  int i, ii, j, row2, row0;
  ii = 0;
  for (; ii < row - 3; ii += bs) {
    j = 0;
    for (; j < ii; 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];
    }
    for (; j < ii + bs; j++) {
      row0 = j - ii;
      if (row0 < 0)
        row0 = 0;
      for (i = row0; i < bs; i++) {
        B[i + j * bs + ii * sdb] = A[i + j * bs + ii * sda];
      }
    }
  }
  if (ii < row) {
    row2 = row - ii;
    if (bs < row2)
      row2 = bs;
    for (j = 0; j < ii + row2; j++) {
      row0 = j - ii;
      if (row0 < 0)
        row0 = 0;
      for (i = row0; i < row2; i++) {
        B[i + j * bs + ii * sdb] = A[i + j * bs + ii * sda];
      }
    }
  }
}


x86 assembly (clang, -O0):
<d_copy_pmat_l>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
movl   $0x4,-0x28(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x3,%ecx
cmp    %ecx,%eax
jge    12fd <d_copy_pmat_l+0x1fd>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x30(%rbp),%eax
jge    1249 <d_copy_pmat_l+0x149>
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x0,%ecx
mov    -0x30(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x0,%ecx
mov    -0x30(%rbp),%edx
imul   -0x24(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x1,%ecx
mov    -0x30(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x1,%ecx
mov    -0x30(%rbp),%edx
imul   -0x24(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x2,%ecx
mov    -0x30(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x2,%ecx
mov    -0x30(%rbp),%edx
imul   -0x24(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
mov    -0x30(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
mov    -0x30(%rbp),%edx
imul   -0x24(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    113f <d_copy_pmat_l+0x3f>
jmp    124e <d_copy_pmat_l+0x14e>
mov    -0x34(%rbp),%eax
mov    -0x30(%rbp),%ecx
add    $0x4,%ecx
cmp    %ecx,%eax
jge    12ea <d_copy_pmat_l+0x1ea>
mov    -0x34(%rbp),%eax
sub    -0x30(%rbp),%eax
mov    %eax,-0x3c(%rbp)
cmpl   $0x0,-0x3c(%rbp)
jge    1279 <d_copy_pmat_l+0x179>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x4,-0x2c(%rbp)
jge    12d7 <d_copy_pmat_l+0x1d7>
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x34(%rbp),%edx
shl    $0x2,%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x34(%rbp),%edx
shl    $0x2,%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
imul   -0x24(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    127f <d_copy_pmat_l+0x17f>
jmp    12dc <d_copy_pmat_l+0x1dc>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    124e <d_copy_pmat_l+0x14e>
jmp    12ef <d_copy_pmat_l+0x1ef>
mov    -0x30(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x30(%rbp)
jmp    1127 <d_copy_pmat_l+0x27>
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    13d1 <d_copy_pmat_l+0x2d1>
mov    -0x4(%rbp),%eax
sub    -0x30(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    $0x4,%eax
cmp    -0x38(%rbp),%eax
jge    1327 <d_copy_pmat_l+0x227>
movl   $0x4,-0x38(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x30(%rbp),%ecx
add    -0x38(%rbp),%ecx
cmp    %ecx,%eax
jge    13cc <d_copy_pmat_l+0x2cc>
mov    -0x34(%rbp),%eax
sub    -0x30(%rbp),%eax
mov    %eax,-0x3c(%rbp)
cmpl   $0x0,-0x3c(%rbp)
jge    1359 <d_copy_pmat_l+0x259>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    13b9 <d_copy_pmat_l+0x2b9>
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x34(%rbp),%edx
shl    $0x2,%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x34(%rbp),%edx
shl    $0x2,%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
imul   -0x24(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    135f <d_copy_pmat_l+0x25f>
jmp    13be <d_copy_pmat_l+0x2be>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    132e <d_copy_pmat_l+0x22e>
jmp    13d1 <d_copy_pmat_l+0x2d1>
pop    %rbp
ret

x86 assembly (clang, -O2):
<d_copy_pmat_l>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
xor    %r12d,%r12d
cmp    $0x4,%edi
jl     13b8 <d_copy_pmat_l+0x2b8>
mov    %rdi,-0x58(%rsp)
lea    -0x3(%rdi),%eax
mov    %rax,-0x10(%rsp)
mov    %ecx,-0x60(%rsp)
movslq %ecx,%rax
mov    %r9d,-0x5c(%rsp)
movslq %r9d,%rcx
lea    0x0(,%rcx,4),%rsi
mov    %rsi,-0x28(%rsp)
mov    %rcx,-0x40(%rsp)
shl    $0x5,%rcx
mov    %rcx,-0x18(%rsp)
add    $0x80,%rcx
mov    %rcx,-0x48(%rsp)
lea    0x0(,%rax,4),%rcx
mov    %rcx,-0x50(%rsp)
mov    %rax,-0x38(%rsp)
shl    $0x5,%rax
mov    %rax,-0x20(%rsp)
add    $0x80,%rax
mov    %rax,-0x30(%rsp)
xor    %r15d,%r15d
xor    %edi,%edi
mov    %rdx,%r14
mov    %r8,%r9
xor    %ebx,%ebx
xor    %r12d,%r12d
jmp    11c9 <d_copy_pmat_l+0xc9>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rcx,%rbx
add    $0x1,%rbx
add    -0x18(%rsp),%r9
add    -0x20(%rsp),%r14
mov    -0x8(%rsp),%rdi
sub    $0xffffffffffffff80,%rdi
add    $0xfffffffffffffffc,%r15
cmp    -0x10(%rsp),%r12
jae    13aa <d_copy_pmat_l+0x2aa>
mov    %r12,%rax
mov    $0x0,%esi
test   %r12,%r12
je     12ba <d_copy_pmat_l+0x1ba>
mov    -0x30(%rsp),%rcx
imul   %rbx,%rcx
add    %rdx,%rcx
mov    -0x28(%rsp),%rsi
imul   %rbx,%rsi
lea    (%r8,%rsi,8),%rsi
cmp    %rsi,%rcx
jbe    1270 <d_copy_pmat_l+0x170>
mov    -0x50(%rsp),%rcx
imul   %rbx,%rcx
lea    (%rdx,%rcx,8),%rcx
mov    -0x48(%rsp),%rsi
imul   %rbx,%rsi
add    %r8,%rsi
cmp    %rcx,%rsi
jbe    1270 <d_copy_pmat_l+0x170>
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movsd  (%r14,%rcx,1),%xmm0
movsd  %xmm0,(%r9,%rcx,1)
movsd  0x8(%r14,%rcx,1),%xmm0
movsd  %xmm0,0x8(%r9,%rcx,1)
movsd  0x10(%r14,%rcx,1),%xmm0
movsd  %xmm0,0x10(%r9,%rcx,1)
movsd  0x18(%r14,%rcx,1),%xmm0
movsd  %xmm0,0x18(%r9,%rcx,1)
add    $0x20,%rcx
cmp    %rcx,%rdi
jne    1220 <d_copy_pmat_l+0x120>
jmp    12b8 <d_copy_pmat_l+0x1b8>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rax,%rcx
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nop
movups (%r14,%rbp,1),%xmm0
movups 0x10(%r14,%rbp,1),%xmm1
movups 0x20(%r14,%rbp,1),%xmm2
movups 0x30(%r14,%rbp,1),%xmm3
movups %xmm2,0x20(%r9,%rbp,1)
movups %xmm3,0x30(%r9,%rbp,1)
movups %xmm0,(%r9,%rbp,1)
movups %xmm1,0x10(%r9,%rbp,1)
add    $0x40,%rbp
add    $0xfffffffffffffffe,%rcx
jne    1280 <d_copy_pmat_l+0x180>
mov    %eax,%esi
mov    %rbx,%rcx
mov    %rdi,-0x8(%rsp)
lea    0x4(%rax),%r12
mov    %esi,%r11d
cmp    %r11,%r12
jbe    11a0 <d_copy_pmat_l+0xa0>
mov    %rax,%r10
imul   -0x38(%rsp),%r10
imul   -0x40(%rsp),%rax
lea    0x0(,%r11,4),%r13d
jmp    1301 <d_copy_pmat_l+0x201>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r11
add    $0x4,%r13d
cmp    %r12,%r11
jae    11a0 <d_copy_pmat_l+0xa0>
mov    %r15,%rbp
add    %r11,%rbp
mov    $0x0,%esi
cmovns %rbp,%rsi
test   $0xfffffffc,%esi
jne    12f0 <d_copy_pmat_l+0x1f0>
test   %ebp,%ebp
mov    $0x0,%esi
cmovg  %ebp,%esi
movslq %r13d,%rbx
add    %rsi,%rbx
mov    %r15,%rdi
lea    (%rbx,%r10,1),%r15
movsd  (%rdx,%r15,8),%xmm0
mov    %rdi,%r15
lea    (%rbx,%rax,1),%rdi
movsd  %xmm0,(%r8,%rdi,8)
cmp    $0x2,%esi
ja     12f0 <d_copy_pmat_l+0x1f0>
lea    (%r10,%rbx,1),%rdi
add    $0x1,%rdi
movsd  (%rdx,%rdi,8),%xmm0
lea    (%rax,%rbx,1),%rdi
add    $0x1,%rdi
movsd  %xmm0,(%r8,%rdi,8)
cmp    $0x2,%esi
je     12f0 <d_copy_pmat_l+0x1f0>
lea    (%r10,%rbx,1),%rsi
add    $0x2,%rsi
movsd  (%rdx,%rsi,8),%xmm0
lea    (%rax,%rbx,1),%rsi
add    $0x2,%rsi
movsd  %xmm0,(%r8,%rsi,8)
test   %ebp,%ebp
jg     12f0 <d_copy_pmat_l+0x1f0>
lea    (%r10,%rbx,1),%rsi
add    $0x3,%rsi
movsd  (%rdx,%rsi,8),%xmm0
lea    (%rax,%rbx,1),%rsi
add    $0x3,%rsi
movsd  %xmm0,(%r8,%rsi,8)
jmp    12f0 <d_copy_pmat_l+0x1f0>
mov    -0x5c(%rsp),%r9d
mov    -0x60(%rsp),%ecx
mov    -0x58(%rsp),%rdi
cmp    %edi,%r12d
jge    14a7 <d_copy_pmat_l+0x3a7>
sub    %r12d,%edi
cmp    $0x4,%edi
mov    $0x4,%r13d
cmovl  %edi,%r13d
lea    (%r12,%r13,1),%eax
test   %eax,%eax
jle    14a7 <d_copy_pmat_l+0x3a7>
imul   %r12d,%ecx
imul   %r12d,%r9d
movslq %ecx,%r10
movslq %r9d,%r11
movslq %r13d,%r14
movslq %r12d,%rcx
mov    %eax,%r15d
neg    %rcx
shl    $0x2,%r15
xor    %ebx,%ebx
jmp    1411 <d_copy_pmat_l+0x311>
nop
add    $0x1,%rcx
add    $0x4,%rbx
cmp    %rbx,%r15
je     14a7 <d_copy_pmat_l+0x3a7>
test   %rcx,%rcx
mov    $0x0,%eax
cmovg  %ecx,%eax
cmp    %eax,%r13d
jle    1400 <d_copy_pmat_l+0x300>
test   %ecx,%ecx
mov    $0x0,%edi
cmovg  %ecx,%edi
movslq %ebx,%rbp
lea    (%rdi,%rbp,1),%rsi
lea    (%rsi,%r10,1),%rax
movsd  (%rdx,%rax,8),%xmm0
add    %r11,%rsi
movsd  %xmm0,(%r8,%rsi,8)
lea    0x1(%rdi),%rsi
cmp    %r14,%rsi
jge    1400 <d_copy_pmat_l+0x300>
add    %rbp,%rsi
lea    (%rsi,%r10,1),%rax
movsd  (%rdx,%rax,8),%xmm0
add    %r11,%rsi
movsd  %xmm0,(%r8,%rsi,8)
lea    0x2(%rdi),%rsi
cmp    %r14,%rsi
jge    1400 <d_copy_pmat_l+0x300>
add    %rbp,%rsi
lea    (%rsi,%r10,1),%rax
movsd  (%rdx,%rax,8),%xmm0
add    %r11,%rsi
movsd  %xmm0,(%r8,%rsi,8)
add    $0x3,%rdi
cmp    %r14,%rdi
jge    1400 <d_copy_pmat_l+0x300>
add    %rbp,%rdi
lea    (%rdi,%r10,1),%rax
movsd  (%rdx,%rax,8),%xmm0
add    %r11,%rdi
movsd  %xmm0,(%r8,%rdi,8)
jmp    1400 <d_copy_pmat_l+0x300>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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