void s_copy_pmat(int row, int col, int bs, float *A, int sda, float *B,
                 int sdb) {
  int i, ii, j, row2;
  for (ii = 0; ii < row; ii += bs) {
    row2 = row - ii;
    if (bs < row2)
      row2 = bs;
    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):
<s_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   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11e6 <s_copy_pmat+0xe6>
mov    -0x4(%rbp),%eax
sub    -0x30(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    114a <s_copy_pmat+0x4a>
mov    -0xc(%rbp),%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11d3 <s_copy_pmat+0xd3>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x38(%rbp),%eax
jge    11c0 <s_copy_pmat+0xc0>
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x34(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
mov    -0x34(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1164 <s_copy_pmat+0x64>
jmp    11c5 <s_copy_pmat+0xc5>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1151 <s_copy_pmat+0x51>
jmp    11d8 <s_copy_pmat+0xd8>
mov    -0xc(%rbp),%eax
add    -0x30(%rbp),%eax
mov    %eax,-0x30(%rbp)
jmp    1123 <s_copy_pmat+0x23>
pop    %rbp
ret

x86 assembly (clang, -O2):
<s_copy_pmat>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %esi,-0x7c(%rsp)
test   %edi,%edi
jle    13cb <s_copy_pmat+0x2cb>
movslq %edx,%r14
mov    %edi,%eax
movslq %r8d,%rdx
movslq 0x80(%rsp),%rsi
mov    -0x7c(%rsp),%r11d
mov    %rsi,-0x38(%rsp)
imul   %r14,%rsi
mov    %rdx,-0x30(%rsp)
mov    %rdx,%rdi
imul   %r14,%rdi
lea    0x30(%r9),%rdx
mov    %rdx,-0x58(%rsp)
mov    %rsi,-0x40(%rsp)
lea    0x0(,%rsi,4),%rdx
mov    %rdx,-0x20(%rsp)
lea    0x0(,%r14,4),%rdx
lea    0x30(%rcx),%rsi
mov    %rsi,-0x60(%rsp)
mov    %rdi,-0x48(%rsp)
lea    0x0(,%rdi,4),%rsi
mov    %rsi,-0x28(%rsp)
xor    %esi,%esi
mov    %rsi,-0x68(%rsp)
mov    %rcx,-0x70(%rsp)
mov    %r9,-0x78(%rsp)
xor    %edi,%edi
mov    %rax,-0x18(%rsp)
mov    %r14,0x8(%rsp)
jmp    11ea <s_copy_pmat+0xea>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    -0x50(%rsp),%rdi
add    %r14,%rdi
addq   $0x1,-0x68(%rsp)
mov    -0x20(%rsp),%rax
add    %rax,-0x58(%rsp)
mov    -0x28(%rsp),%rsi
add    %rsi,-0x60(%rsp)
add    %rax,-0x78(%rsp)
add    %rsi,-0x70(%rsp)
mov    -0x18(%rsp),%rax
cmp    %rax,%rdi
jge    13cb <s_copy_pmat+0x2cb>
mov    %rdi,-0x50(%rsp)
sub    %rdi,%rax
cmp    %r14,%rax
mov    %r14d,%r15d
cmovle %eax,%r15d
cmpl   $0x0,-0x7c(%rsp)
jle    11b0 <s_copy_pmat+0xb0>
mov    -0x40(%rsp),%rsi
mov    -0x68(%rsp),%rax
imul   %rax,%rsi
mov    -0x48(%rsp),%rdi
imul   %rax,%rdi
mov    -0x50(%rsp),%rax
mov    %rax,%rbp
imul   -0x30(%rsp),%rbp
mov    %rbp,(%rsp)
imul   -0x38(%rsp),%rax
mov    %rax,-0x8(%rsp)
movslq %r15d,%rbp
mov    %rsi,0x40(%rsp)
lea    (%rsi,%rbp,1),%rax
mov    %rax,0x28(%rsp)
mov    %rdi,0x30(%rsp)
lea    (%rdi,%rbp,1),%rax
mov    %rax,0x38(%rsp)
lea    -0x8(%rbp),%rax
mov    %rax,0x20(%rsp)
mov    %rax,%rsi
shr    $0x3,%rsi
add    $0x1,%rsi
mov    %rbp,%rax
and    $0xfffffffffffffff8,%rax
mov    %rax,0x10(%rsp)
mov    %rsi,%rax
mov    %rsi,0x18(%rsp)
and    $0xfffffffffffffffe,%rsi
mov    %rsi,-0x10(%rsp)
mov    -0x70(%rsp),%r13
mov    -0x78(%rsp),%r8
mov    -0x60(%rsp),%rax
mov    -0x58(%rsp),%rbx
xor    %r12d,%r12d
jmp    12b9 <s_copy_pmat+0x1b9>
add    $0x1,%r12
add    %rdx,%rbx
add    %rdx,%rax
add    %rdx,%r8
add    %rdx,%r13
cmp    %r11,%r12
je     11b0 <s_copy_pmat+0xb0>
test   %r15d,%r15d
jle    12a0 <s_copy_pmat+0x1a0>
cmp    $0x8,%r15d
jae    12d0 <s_copy_pmat+0x1d0>
xor    %esi,%esi
jmp    13b0 <s_copy_pmat+0x2b0>
nopl   0x0(%rax,%rax,1)
mov    %r12,%r10
imul   %r14,%r10
mov    0x40(%rsp),%rsi
add    %r10,%rsi
mov    0x38(%rsp),%rdi
add    %r10,%rdi
lea    (%r9,%rsi,4),%rsi
lea    (%rcx,%rdi,4),%rdi
cmp    %rdi,%rsi
jae    1318 <s_copy_pmat+0x218>
mov    0x28(%rsp),%rsi
add    %r10,%rsi
mov    0x30(%rsp),%rdi
add    %r10,%rdi
lea    (%r9,%rsi,4),%rsi
lea    (%rcx,%rdi,4),%rdi
cmp    %rsi,%rdi
jae    1318 <s_copy_pmat+0x218>
xor    %esi,%esi
jmp    13b0 <s_copy_pmat+0x2b0>
cmpq   $0x8,0x20(%rsp)
jae    1325 <s_copy_pmat+0x225>
xor    %r14d,%r14d
jmp    1368 <s_copy_pmat+0x268>
mov    -0x10(%rsp),%rsi
xor    %r14d,%r14d
nopl   (%rax)
movups -0x30(%rax,%r14,4),%xmm0
movups -0x20(%rax,%r14,4),%xmm1
movups %xmm0,-0x30(%rbx,%r14,4)
movups %xmm1,-0x20(%rbx,%r14,4)
movups -0x10(%rax,%r14,4),%xmm0
movups (%rax,%r14,4),%xmm1
movups %xmm0,-0x10(%rbx,%r14,4)
movups %xmm1,(%rbx,%r14,4)
add    $0x10,%r14
add    $0xfffffffffffffffe,%rsi
jne    1330 <s_copy_pmat+0x230>
testb  $0x1,0x18(%rsp)
je     1398 <s_copy_pmat+0x298>
add    %r10,%r14
mov    (%rsp),%rsi
add    %r14,%rsi
mov    %rcx,%rdi
movups (%rcx,%rsi,4),%xmm0
movups 0x10(%rcx,%rsi,4),%xmm1
add    -0x8(%rsp),%r14
mov    %r9,%rsi
movups %xmm0,(%r9,%r14,4)
movups %xmm1,0x10(%r9,%r14,4)
mov    0x10(%rsp),%rdi
mov    %rdi,%rsi
cmp    %rbp,%rdi
mov    0x8(%rsp),%r14
je     12a0 <s_copy_pmat+0x1a0>
xchg   %ax,%ax
movss  0x0(%r13,%rsi,4),%xmm0
movss  %xmm0,(%r8,%rsi,4)
add    $0x1,%rsi
cmp    %rbp,%rsi
jl     13b0 <s_copy_pmat+0x2b0>
jmp    12a0 <s_copy_pmat+0x1a0>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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