void cvt_d2s_mat2pmat(int row, int col, int offset, int bs, double *A, int lda,
                      float *pA, int sda) {
  int i, ii, j, row0, row1, row2;
  row0 = (bs - offset % bs) % bs;
  if (row0 > row)
    row0 = row;
  row1 = row - row0;
  if (row0 > 0) {
    for (j = 0; j < col; j++) {
      for (i = 0; i < row0; i++) {
        pA[i + j * bs] = (float)A[i + j * lda];
      }
    }
    A += row0;
    pA += row0 + bs * (sda - 1);
  }
  for (ii = 0; ii < row1; ii += bs) {
    row2 = row1 - ii;
    if (bs < row2)
      row2 = bs;
    for (j = 0; j < col; j++) {
      for (i = 0; i < row2; i++) {
        pA[i + j * bs + ii * sda] = (float)A[i + ii + j * lda];
      }
    }
  }
}


x86 assembly (clang, -O0):
<cvt_d2s_mat2pmat>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8,-0x18(%rbp)
mov    %r9d,-0x1c(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0xc(%rbp),%eax
cltd
idivl  -0x10(%rbp)
mov    -0x38(%rbp),%eax
sub    %edx,%eax
cltd
idivl  -0x10(%rbp)
mov    %edx,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jle    114a <cvt_d2s_mat2pmat+0x4a>
mov    -0x4(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x4(%rbp),%eax
sub    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jle    1210 <cvt_d2s_mat2pmat+0x110>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11d8 <cvt_d2s_mat2pmat+0xd8>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    11c5 <cvt_d2s_mat2pmat+0xc5>
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    -0x28(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    -0x28(%rbp),%edx
imul   -0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1177 <cvt_d2s_mat2pmat+0x77>
jmp    11ca <cvt_d2s_mat2pmat+0xca>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1164 <cvt_d2s_mat2pmat+0x64>
mov    -0x2c(%rbp),%ecx
mov    -0x18(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x2c(%rbp),%ecx
mov    -0x10(%rbp),%eax
mov    0x18(%rbp),%edx
sub    $0x1,%edx
imul   %edx,%eax
add    %eax,%ecx
mov    0x10(%rbp),%rax
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,0x10(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x30(%rbp),%eax
jge    12d8 <cvt_d2s_mat2pmat+0x1d8>
mov    -0x30(%rbp),%eax
sub    -0x24(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    123e <cvt_d2s_mat2pmat+0x13e>
mov    -0x10(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    12c5 <cvt_d2s_mat2pmat+0x1c5>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    12b2 <cvt_d2s_mat2pmat+0x1b2>
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x24(%rbp),%ecx
mov    -0x28(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    -0x28(%rbp),%edx
imul   -0x10(%rbp),%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
imul   0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1258 <cvt_d2s_mat2pmat+0x158>
jmp    12b7 <cvt_d2s_mat2pmat+0x1b7>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1245 <cvt_d2s_mat2pmat+0x145>
jmp    12ca <cvt_d2s_mat2pmat+0x1ca>
mov    -0x10(%rbp),%eax
add    -0x24(%rbp),%eax
mov    %eax,-0x24(%rbp)
jmp    1217 <cvt_d2s_mat2pmat+0x117>
pop    %rbp
ret

x86 assembly (clang, -O2):
<cvt_d2s_mat2pmat>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x20,%rsp
mov    %edx,%eax
cltd
idiv   %ecx
mov    %ecx,%eax
sub    %edx,%eax
cltd
idiv   %ecx
cmp    %edi,%edx
cmovg  %edi,%edx
mov    0x58(%rsp),%rax
sub    %edx,%edi
movslq %ecx,%rbp
mov    %rbp,-0x48(%rsp)
test   %edx,%edx
mov    %r9d,-0x7c(%rsp)
mov    %ecx,-0x58(%rsp)
mov    %esi,-0x54(%rsp)
mov    %edi,-0x68(%rsp)
jle    1309 <cvt_d2s_mat2pmat+0x209>
test   %esi,%esi
jle    12d9 <cvt_d2s_mat2pmat+0x1d9>
movslq -0x7c(%rsp),%rcx
mov    %esi,%r12d
mov    %edx,%ebp
mov    %ebp,%edi
and    $0xfffffffc,%edi
mov    %rdi,-0x50(%rsp)
add    $0xfffffffffffffffc,%rdi
mov    %rdi,-0x38(%rsp)
mov    %rdi,%rbx
shr    $0x2,%rbx
add    $0x1,%rbx
mov    %rbx,%rdi
mov    %rbx,-0x70(%rsp)
and    $0xfffffffffffffffe,%rbx
mov    %rbx,-0x78(%rsp)
lea    0x30(%r8),%r13
mov    %rcx,-0x40(%rsp)
lea    0x0(,%rcx,8),%r9
lea    0x18(%rax),%rbx
mov    -0x48(%rsp),%rcx
lea    0x0(,%rcx,4),%r10
xor    %r14d,%r14d
mov    %r8,%r11
mov    %rax,%r15
jmp    11d9 <cvt_d2s_mat2pmat+0xd9>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r14
add    %r9,%r13
add    %r10,%rbx
add    %r10,%r15
add    %r9,%r11
cmp    %r12,%r14
je     12d9 <cvt_d2s_mat2pmat+0x1d9>
cmp    $0x4,%edx
jae    11f0 <cvt_d2s_mat2pmat+0xf0>
xor    %ecx,%ecx
jmp    12b0 <cvt_d2s_mat2pmat+0x1b0>
cs nopw 0x0(%rax,%rax,1)
nop
cmpq   $0x0,-0x38(%rsp)
je     12ce <cvt_d2s_mat2pmat+0x1ce>
mov    -0x78(%rsp),%rcx
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd -0x30(%r13,%rdi,8),%xmm0
movupd -0x20(%r13,%rdi,8),%xmm1
cvtpd2ps %xmm0,%xmm0
cvtpd2ps %xmm1,%xmm1
unpcklpd %xmm1,%xmm0
movupd %xmm0,-0x18(%rbx,%rdi,4)
movupd -0x10(%r13,%rdi,8),%xmm0
movupd 0x0(%r13,%rdi,8),%xmm1
cvtpd2ps %xmm0,%xmm0
cvtpd2ps %xmm1,%xmm1
unpcklpd %xmm1,%xmm0
movupd %xmm0,-0x8(%rbx,%rdi,4)
add    $0x8,%rdi
add    $0xfffffffffffffffe,%rcx
jne    1210 <cvt_d2s_mat2pmat+0x110>
testb  $0x1,-0x70(%rsp)
je     129e <cvt_d2s_mat2pmat+0x19e>
mov    %r14,%rcx
imul   -0x40(%rsp),%rcx
mov    %r12,%rsi
mov    %r14,%r12
add    %rdi,%rcx
movupd (%r8,%rcx,8),%xmm0
movupd 0x10(%r8,%rcx,8),%xmm1
cvtpd2ps %xmm0,%xmm0
cvtpd2ps %xmm1,%xmm1
imul   -0x48(%rsp),%r12
unpcklpd %xmm1,%xmm0
add    %rdi,%r12
movupd %xmm0,(%rax,%r12,4)
mov    %rsi,%r12
mov    -0x50(%rsp),%rdi
mov    %rdi,%rcx
cmp    %rbp,%rdi
je     11c0 <cvt_d2s_mat2pmat+0xc0>
nop
movsd  (%r11,%rcx,8),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%r15,%rcx,4)
add    $0x1,%rcx
cmp    %rcx,%rbp
jne    12b0 <cvt_d2s_mat2pmat+0x1b0>
jmp    11c0 <cvt_d2s_mat2pmat+0xc0>
xor    %edi,%edi
testb  $0x1,-0x70(%rsp)
jne    1261 <cvt_d2s_mat2pmat+0x161>
jmp    129e <cvt_d2s_mat2pmat+0x19e>
movslq %edx,%rcx
lea    (%r8,%rcx,8),%r8
mov    0x60(%rsp),%ecx
lea    -0x1(%rcx),%edi
mov    -0x58(%rsp),%ecx
imul   %ecx,%edi
add    %edx,%edi
movslq %edi,%rdx
lea    (%rax,%rdx,4),%rax
mov    -0x7c(%rsp),%r10d
mov    -0x68(%rsp),%edi
test   %edi,%edi
jg     1314 <cvt_d2s_mat2pmat+0x214>
jmp    1612 <cvt_d2s_mat2pmat+0x512>
mov    %r9d,%r10d
test   %edi,%edi
jle    1612 <cvt_d2s_mat2pmat+0x512>
mov    0x60(%rsp),%edx
imul   %ecx,%edx
mov    %rdx,-0x20(%rsp)
mov    %ecx,%r12d
mov    %r10d,%edx
mov    %rdx,-0x38(%rsp)
xor    %edx,%edx
movl   $0x4,-0x60(%rsp)
xor    %ebp,%ebp
mov    %rbp,-0x28(%rsp)
mov    %edi,%r9d
movl   $0x4,-0x64(%rsp)
xor    %ebp,%ebp
jmp    1387 <cvt_d2s_mat2pmat+0x287>
nopl   0x0(%rax,%rax,1)
mov    -0x70(%rsp),%rdx
add    -0x48(%rsp),%rdx
mov    -0x5c(%rsp),%ebp
add    $0x1,%ebp
add    %ecx,-0x64(%rsp)
mov    -0x20(%rsp),%rdi
add    %edi,-0x60(%rsp)
mov    -0x30(%rsp),%r9d
sub    %ecx,%r9d
add    %rdi,-0x28(%rsp)
mov    -0x68(%rsp),%edi
cmp    %edx,%edi
jle    1612 <cvt_d2s_mat2pmat+0x512>
cmp    %r9d,%ecx
mov    %r9d,-0x30(%rsp)
cmovl  %ecx,%r9d
mov    %ebp,-0x5c(%rsp)
imul   %ecx,%ebp
mov    %edi,%ebx
sub    %ebp,%ebx
cmp    %ebx,%ecx
cmovl  %ecx,%ebx
mov    %rdx,-0x70(%rsp)
sub    %edx,%edi
cmp    %ecx,%edi
cmovg  %ecx,%edi
mov    %edi,-0x50(%rsp)
cmpl   $0x0,-0x54(%rsp)
jle    1350 <cvt_d2s_mat2pmat+0x250>
and    $0xfffffffc,%r9d
movslq %r9d,%rdx
add    $0xfffffffffffffffc,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
and    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x18(%rsp)
movslq %ebx,%rdx
and    $0xfffffffc,%ebx
movslq %ebx,%rsi
mov    %rdx,%rbx
add    $0xfffffffffffffffc,%rsi
mov    %rsi,(%rsp)
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %rsi,-0x8(%rsp)
add    $0xffffffffffffffff,%rdx
mov    %rdx,-0x40(%rsp)
mov    -0x20(%rsp),%rdx
imul   -0x5c(%rsp),%edx
mov    %rdx,0x8(%rsp)
mov    -0x70(%rsp),%r15
mov    %r15d,%esi
mov    0x60(%rsp),%edx
imul   %edx,%esi
mov    %esi,-0x2c(%rsp)
movslq -0x50(%rsp),%rbp
mov    %rbx,%rdx
and    $0xfffffffffffffffc,%rdx
mov    %rdx,-0x10(%rsp)
xor    %r11d,%r11d
mov    -0x28(%rsp),%r13
mov    -0x60(%rsp),%r14d
mov    -0x64(%rsp),%r9d
mov    %rbx,-0x78(%rsp)
jmp    146d <cvt_d2s_mat2pmat+0x36d>
nopl   0x0(%rax)
add    $0x1,%r11d
add    %r10d,%r9d
add    %ecx,%r14d
add    %r12,%r13
add    -0x38(%rsp),%r15
cmp    -0x54(%rsp),%r11d
je     1350 <cvt_d2s_mat2pmat+0x250>
cmpl   $0x0,-0x50(%rsp)
jle    1450 <cvt_d2s_mat2pmat+0x350>
cmp    $0x4,%ebx
jae    14b0 <cvt_d2s_mat2pmat+0x3b0>
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
lea    (%r15,%rdi,1),%edx
movslq %edx,%rdx
movsd  (%r8,%rdx,8),%xmm0
cvtsd2ss %xmm0,%xmm0
lea    (%rdi,%r13,1),%edx
movslq %edx,%rdx
movss  %xmm0,(%rax,%rdx,4)
add    $0x1,%rdi
cmp    %rbp,%rdi
jl     1480 <cvt_d2s_mat2pmat+0x380>
jmp    1450 <cvt_d2s_mat2pmat+0x350>
nopl   0x0(%rax,%rax,1)
mov    %r11d,%esi
imul   %r10d,%esi
add    -0x70(%rsp),%esi
mov    -0x40(%rsp),%rdi
lea    (%rsi,%rdi,1),%edx
cmp    %esi,%edx
jl     15f8 <cvt_d2s_mat2pmat+0x4f8>
mov    %rdi,%rdx
shr    $0x20,%rdx
jne    15ff <cvt_d2s_mat2pmat+0x4ff>
mov    %rsi,0x18(%rsp)
mov    %r11d,%esi
imul   %ecx,%esi
mov    0x8(%rsp),%rbx
mov    %rsi,0x10(%rsp)
add    %esi,%ebx
add    %ebx,%edi
cmp    %ebx,%edi
jl     1606 <cvt_d2s_mat2pmat+0x506>
mov    $0x0,%edi
test   %rdx,%rdx
mov    -0x78(%rsp),%rbx
jne    1480 <cvt_d2s_mat2pmat+0x380>
mov    %r12,%rbx
cmpq   $0x0,(%rsp)
je     1592 <cvt_d2s_mat2pmat+0x492>
mov    -0x18(%rsp),%r12
mov    %r14d,%ecx
mov    %r9d,%r10d
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
lea    -0x4(%r10),%edi
movslq %edi,%rdi
movupd (%r8,%rdi,8),%xmm0
movupd 0x10(%r8,%rdi,8),%xmm1
cvtpd2ps %xmm0,%xmm0
cvtpd2ps %xmm1,%xmm1
unpcklpd %xmm1,%xmm0
lea    -0x4(%rcx),%edi
movslq %edi,%rdi
movupd %xmm0,(%rax,%rdi,4)
movslq %r10d,%r10
movupd (%r8,%r10,8),%xmm0
movupd 0x10(%r8,%r10,8),%xmm1
cvtpd2ps %xmm0,%xmm0
cvtpd2ps %xmm1,%xmm1
unpcklpd %xmm1,%xmm0
movslq %ecx,%rcx
movupd %xmm0,(%rax,%rcx,4)
add    $0x8,%rdx
add    $0x8,%r10d
add    $0x8,%ecx
add    $0xfffffffffffffffe,%r12
jne    1530 <cvt_d2s_mat2pmat+0x430>
jmp    1594 <cvt_d2s_mat2pmat+0x494>
xor    %edx,%edx
testb  $0x1,-0x8(%rsp)
mov    -0x7c(%rsp),%r10d
mov    -0x58(%rsp),%ecx
mov    %rbx,%r12
je     15dd <cvt_d2s_mat2pmat+0x4dd>
mov    0x10(%rsp),%rsi
add    -0x2c(%rsp),%esi
mov    0x18(%rsp),%rdi
add    %edx,%edi
movslq %edi,%rdi
movupd (%r8,%rdi,8),%xmm0
movupd 0x10(%r8,%rdi,8),%xmm1
cvtpd2ps %xmm0,%xmm0
cvtpd2ps %xmm1,%xmm1
unpcklpd %xmm1,%xmm0
add    %edx,%esi
movslq %esi,%rdx
movupd %xmm0,(%rax,%rdx,4)
mov    -0x10(%rsp),%rdx
mov    %rdx,%rdi
mov    -0x78(%rsp),%rbx
cmp    %rbx,%rdx
jne    1480 <cvt_d2s_mat2pmat+0x380>
jmp    1450 <cvt_d2s_mat2pmat+0x350>
xor    %edi,%edi
jmp    1480 <cvt_d2s_mat2pmat+0x380>
xor    %edi,%edi
jmp    1480 <cvt_d2s_mat2pmat+0x380>
xor    %edi,%edi
mov    -0x78(%rsp),%rbx
jmp    1480 <cvt_d2s_mat2pmat+0x380>
add    $0x20,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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