void matmult_vpur_notrans(double *C, const double *A, const double *B,
                          double alpha, int ldC, int ldA, int ldB, int M, int N,
                          int P, int S, int L, int R, int E, int vlen) {
  int j, k, vpS, vpL;
  int SLcnt, REcnt, REmod;
  const double *Bc, *Br, *Ac, *AvpS;
  double *Cc;
  register const double *Ar;
  register int i;
  register double *Cr;
  register double coeff, f0, f1, f2, f3;
  SLcnt = L - S;
  REcnt = E - R;
  REmod = REcnt % 4;
  REcnt -= REmod;
  if (vlen == 0) {
    vlen = 32;
  }
  vpS = 0;
  vpL = vlen < P ? vlen : P;
  while (vpS < P) {
    Cc = &C[S * ldC + R];
    Bc = &B[S * ldB + vpS];
    AvpS = &A[vpS * ldA + R];
    for (j = S; j < L; j++) {
      Ac = AvpS;
      Br = Bc;
      for (k = vpS; k < vpL; k++) {
        Cr = Cc;
        Ar = Ac;
        coeff = *Br;
        if (coeff != 0.0) {
          coeff *= alpha;
          for (i = 0; i < REmod; i++) {
            f0 = Ar[0] * coeff;
            Cr[0] += f0;
            Cr++;
            Ar++;
          }
          for (i = 0; i < REcnt; i += 4) {
            f0 = Ar[0] * coeff;
            Cr[0] += f0;
            f1 = Ar[1] * coeff;
            Cr[1] += f1;
            f2 = Ar[2] * coeff;
            Cr[2] += f2;
            f3 = Ar[3] * coeff;
            Cr[3] += f3;
            Cr += 4;
            Ar += 4;
          }
        }
        Br++;
        Ac += ldA;
      }
      Cc += ldC;
      Bc += ldB;
    }
    vpS = vpL;
    vpL += vlen;
    if (vpL > P) {
      vpL = P;
    }
  }
  return;
}


x86 assembly (clang, -O0):
<matmult_vpur_notrans>:
push   %rbp
mov    %rsp,%rbp
sub    $0x38,%rsp
mov    0x48(%rbp),%eax
mov    0x40(%rbp),%eax
mov    0x38(%rbp),%eax
mov    0x30(%rbp),%eax
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    %ecx,-0x24(%rbp)
mov    %r8d,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
mov    0x30(%rbp),%eax
sub    0x28(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    0x40(%rbp),%eax
sub    0x38(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %edx,-0x48(%rbp)
mov    -0x48(%rbp),%ecx
mov    -0x44(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x44(%rbp)
cmpl   $0x0,0x48(%rbp)
jne    1178 <matmult_vpur_notrans+0x78>
movl   $0x20,0x48(%rbp)
movl   $0x0,-0x38(%rbp)
mov    0x48(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    1199 <matmult_vpur_notrans+0x99>
mov    0x48(%rbp),%eax
mov    %eax,-0xb4(%rbp)
jmp    11a2 <matmult_vpur_notrans+0xa2>
mov    0x20(%rbp),%eax
mov    %eax,-0xb4(%rbp)
mov    -0xb4(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x38(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    149b <matmult_vpur_notrans+0x39b>
mov    -0x8(%rbp),%rax
mov    0x28(%rbp),%ecx
imul   -0x24(%rbp),%ecx
add    0x38(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x70(%rbp)
mov    -0x18(%rbp),%rax
mov    0x28(%rbp),%ecx
imul   -0x2c(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x28(%rbp),%ecx
add    0x38(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x68(%rbp)
mov    0x28(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    0x30(%rbp),%eax
jge    1475 <matmult_vpur_notrans+0x375>
mov    -0x68(%rbp),%rax
mov    %rax,-0x60(%rbp)
mov    -0x50(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x38(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    143d <matmult_vpur_notrans+0x33d>
mov    -0x70(%rbp),%rax
mov    %rax,-0x88(%rbp)
mov    -0x60(%rbp),%rax
mov    %rax,-0x78(%rbp)
mov    -0x58(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1282 <matmult_vpur_notrans+0x182>
jp     1282 <matmult_vpur_notrans+0x182>
jmp    140e <matmult_vpur_notrans+0x30e>
movsd  -0x20(%rbp),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x90(%rbp)
ff
movl   $0x0,-0x7c(%rbp)
mov    -0x7c(%rbp),%eax
cmp    -0x48(%rbp),%eax
jge    1305 <matmult_vpur_notrans+0x205>
mov    -0x78(%rbp),%rax
movsd  (%rax),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
mov    -0x88(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x88(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x88(%rbp)
mov    -0x78(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x78(%rbp)
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
jmp    129e <matmult_vpur_notrans+0x19e>
movl   $0x0,-0x7c(%rbp)
mov    -0x7c(%rbp),%eax
cmp    -0x44(%rbp),%eax
jge    1409 <matmult_vpur_notrans+0x309>
mov    -0x78(%rbp),%rax
movsd  (%rax),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
mov    -0x88(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x78(%rbp),%rax
movsd  0x8(%rax),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0xa0(%rbp),%xmm0
ff
mov    -0x88(%rbp),%rax
addsd  0x8(%rax),%xmm0
movsd  %xmm0,0x8(%rax)
mov    -0x78(%rbp),%rax
movsd  0x10(%rax),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
mov    -0x88(%rbp),%rax
addsd  0x10(%rax),%xmm0
movsd  %xmm0,0x10(%rax)
mov    -0x78(%rbp),%rax
movsd  0x18(%rax),%xmm0
mulsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
mov    -0x88(%rbp),%rax
addsd  0x18(%rax),%xmm0
movsd  %xmm0,0x18(%rax)
mov    -0x88(%rbp),%rax
add    $0x20,%rax
mov    %rax,-0x88(%rbp)
mov    -0x78(%rbp),%rax
add    $0x20,%rax
mov    %rax,-0x78(%rbp)
mov    -0x7c(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x7c(%rbp)
jmp    130c <matmult_vpur_notrans+0x20c>
jmp    140e <matmult_vpur_notrans+0x30e>
mov    -0x58(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x58(%rbp)
mov    -0x28(%rbp),%ecx
mov    -0x60(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x60(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1233 <matmult_vpur_notrans+0x133>
mov    -0x24(%rbp),%ecx
mov    -0x70(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x70(%rbp)
mov    -0x2c(%rbp),%ecx
mov    -0x50(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1211 <matmult_vpur_notrans+0x111>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    0x48(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    0x20(%rbp),%eax
jle    1496 <matmult_vpur_notrans+0x396>
mov    0x20(%rbp),%eax
mov    %eax,-0x3c(%rbp)
jmp    11ab <matmult_vpur_notrans+0xab>
add    $0x38,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<matmult_vpur_notrans>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r9d,%ebp
mov    %r8d,-0x7c(%rsp)
mov    %rdx,-0x20(%rsp)
mov    %rsi,-0x28(%rsp)
mov    0x90(%rsp),%r8d
mov    0x98(%rsp),%eax
mov    %r8d,%r9d
sub    0x88(%rsp),%r9d
lea    0x3(%r9),%ebx
test   %r9d,%r9d
cmovns %r9d,%ebx
test   %eax,%eax
mov    $0x20,%esi
cmovne %eax,%esi
mov    %rsi,-0x58(%rsp)
cmpl   $0x0,0x70(%rsp)
jle    1579 <matmult_vpur_notrans+0x479>
mov    0x78(%rsp),%edx
and    $0xfffffffc,%ebx
sub    %ebx,%r9d
mov    0x70(%rsp),%eax
mov    -0x58(%rsp),%rsi
cmp    %eax,%esi
mov    %eax,%r11d
cmovl  %esi,%r11d
mov    %edx,%eax
imul   %ecx,%eax
mov    0x88(%rsp),%esi
add    %esi,%eax
cltq
lea    (%rdi,%rax,8),%rax
mov    %rax,-0x30(%rsp)
mov    %edx,%eax
imul   %ebp,%eax
mov    %rax,-0x38(%rsp)
movslq -0x7c(%rsp),%r13
movslq %ecx,%rax
mov    %rax,-0x8(%rsp)
movslq %ebp,%rax
mov    %rax,-0x10(%rsp)
mov    %r9d,%eax
not    %eax
add    %r8d,%eax
sub    %esi,%eax
mov    %eax,%ecx
and    $0xfffffffc,%ecx
add    $0x4,%rcx
mov    %rcx,0x20(%rsp)
mov    %eax,-0x5c(%rsp)
shr    $0x2,%eax
add    $0x1,%eax
mov    %rax,0x10(%rsp)
and    $0xfffffffe,%eax
lea    0x0(,%rax,4),%rcx
mov    %rcx,(%rsp)
mov    %rax,-0x50(%rsp)
lea    0x0(,%rax,4),%eax
mov    %eax,-0x64(%rsp)
xor    %eax,%eax
mov    %rax,-0x78(%rsp)
xorpd  %xmm8,%xmm8
mov    %r13,0x8(%rsp)
jmp    1235 <matmult_vpur_notrans+0x135>
nopl   0x0(%rax)
mov    -0x58(%rsp),%rax
add    %r11d,%eax
mov    0x70(%rsp),%ecx
cmp    %ecx,%eax
cmovg  %ecx,%eax
mov    %r11d,%edx
mov    %rdx,-0x78(%rsp)
cmp    %ecx,%r11d
mov    %eax,%r11d
jge    1579 <matmult_vpur_notrans+0x479>
mov    0x78(%rsp),%eax
cmp    0x80(%rsp),%eax
jge    1210 <matmult_vpur_notrans+0x110>
mov    -0x78(%rsp),%rdx
mov    %edx,%eax
imul   -0x7c(%rsp),%eax
add    0x88(%rsp),%eax
cltq
mov    -0x28(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
mov    %rax,-0x18(%rsp)
mov    -0x38(%rsp),%rax
add    %edx,%eax
cltq
mov    -0x20(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
mov    %rax,-0x70(%rsp)
mov    0x78(%rsp),%eax
mov    -0x30(%rsp),%r10
jmp    12c0 <matmult_vpur_notrans+0x1c0>
nopw   0x0(%rax,%rax,1)
mov    -0x8(%rsp),%rax
lea    (%r10,%rax,8),%r10
mov    -0x10(%rsp),%rax
mov    -0x70(%rsp),%rcx
lea    (%rcx,%rax,8),%rcx
mov    %rcx,-0x70(%rsp)
mov    -0x60(%rsp),%eax
add    $0x1,%eax
cmp    0x80(%rsp),%eax
je     1210 <matmult_vpur_notrans+0x110>
mov    %eax,-0x60(%rsp)
cmp    %r11d,-0x78(%rsp)
jge    1290 <matmult_vpur_notrans+0x190>
lea    0x8(%r10),%rax
mov    %rax,-0x40(%rsp)
lea    0x10(%r10),%rax
mov    %rax,-0x48(%rsp)
lea    0x18(%r10),%rax
mov    %rax,0x18(%rsp)
mov    -0x78(%rsp),%rax
mov    %eax,%r14d
mov    -0x18(%rsp),%rdx
mov    -0x70(%rsp),%r12
jmp    1315 <matmult_vpur_notrans+0x215>
nopw   0x0(%rax,%rax,1)
add    $0x8,%r12
lea    (%rdx,%r13,8),%rdx
add    $0x1,%r14d
cmp    %r11d,%r14d
je     1290 <matmult_vpur_notrans+0x190>
movsd  (%r12),%xmm10
ucomisd %xmm8,%xmm10
jne    1324 <matmult_vpur_notrans+0x224>
jnp    1300 <matmult_vpur_notrans+0x200>
mulsd  %xmm0,%xmm10
mov    %rdx,%r15
mov    %r10,%r8
test   %r9d,%r9d
jle    1390 <matmult_vpur_notrans+0x290>
movsd  (%rdx),%xmm1
mulsd  %xmm10,%xmm1
addsd  (%r10),%xmm1
movsd  %xmm1,(%r10)
cmp    $0x1,%r9d
jne    135c <matmult_vpur_notrans+0x25c>
lea    0x8(%rdx),%r15
mov    -0x40(%rsp),%r8
test   %ebx,%ebx
jg     1398 <matmult_vpur_notrans+0x298>
jmp    1300 <matmult_vpur_notrans+0x200>
movsd  0x8(%rdx),%xmm1
mulsd  %xmm10,%xmm1
mov    -0x40(%rsp),%rax
addsd  (%rax),%xmm1
movsd  %xmm1,(%rax)
cmp    $0x2,%r9d
jne    154c <matmult_vpur_notrans+0x44c>
lea    0x10(%rdx),%r15
mov    -0x48(%rsp),%r8
cs nopw 0x0(%rax,%rax,1)
test   %ebx,%ebx
jle    1300 <matmult_vpur_notrans+0x200>
xor    %eax,%eax
cmpl   $0x4,-0x5c(%rsp)
jb     13c0 <matmult_vpur_notrans+0x2c0>
mov    0x20(%rsp),%rsi
lea    (%r15,%rsi,8),%rcx
cmp    %rcx,%r8
jae    1439 <matmult_vpur_notrans+0x339>
lea    (%r8,%rsi,8),%rcx
cmp    %rcx,%r15
jae    1439 <matmult_vpur_notrans+0x339>
mov    %r8,%rcx
mov    %r15,%rsi
mov    %eax,%ebp
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
movsd  (%rsi,%rax,8),%xmm1
mulsd  %xmm10,%xmm1
addsd  (%rcx,%rax,8),%xmm1
movsd  %xmm1,(%rcx,%rax,8)
movsd  0x8(%rsi,%rax,8),%xmm1
mulsd  %xmm10,%xmm1
addsd  0x8(%rcx,%rax,8),%xmm1
movsd  %xmm1,0x8(%rcx,%rax,8)
movsd  0x10(%rsi,%rax,8),%xmm1
mulsd  %xmm10,%xmm1
addsd  0x10(%rcx,%rax,8),%xmm1
movsd  %xmm1,0x10(%rcx,%rax,8)
movsd  0x18(%rsi,%rax,8),%xmm1
mulsd  %xmm10,%xmm1
addsd  0x18(%rcx,%rax,8),%xmm1
movsd  %xmm1,0x18(%rcx,%rax,8)
add    $0x4,%rax
lea    (%rax,%rbp,1),%edi
cmp    %ebx,%edi
jl     13d0 <matmult_vpur_notrans+0x2d0>
jmp    1300 <matmult_vpur_notrans+0x200>
mov    %r11,%rdi
mov    (%rsp),%rax
lea    (%r8,%rax,8),%rcx
lea    (%r15,%rax,8),%rsi
movapd %xmm10,%xmm9
unpcklpd %xmm10,%xmm9
mov    -0x50(%rsp),%r13
xor    %r11d,%r11d
nopw   0x0(%rax,%rax,1)
movupd (%r15,%r11,1),%xmm4
movupd 0x10(%r15,%r11,1),%xmm5
movupd 0x20(%r15,%r11,1),%xmm6
movupd 0x30(%r15,%r11,1),%xmm7
movapd %xmm4,%xmm1
unpcklpd %xmm6,%xmm1
unpckhpd %xmm6,%xmm4
movapd %xmm5,%xmm6
unpcklpd %xmm7,%xmm6
unpckhpd %xmm7,%xmm5
mulpd  %xmm9,%xmm1
movupd (%r8,%r11,1),%xmm2
movupd 0x10(%r8,%r11,1),%xmm3
movupd 0x20(%r8,%r11,1),%xmm12
movupd 0x30(%r8,%r11,1),%xmm11
movapd %xmm2,%xmm7
unpcklpd %xmm12,%xmm7
addpd  %xmm1,%xmm7
unpckhpd %xmm12,%xmm2
movapd %xmm3,%xmm1
unpcklpd %xmm11,%xmm1
unpckhpd %xmm11,%xmm3
mulpd  %xmm9,%xmm4
addpd  %xmm2,%xmm4
mulpd  %xmm9,%xmm6
addpd  %xmm1,%xmm6
mulpd  %xmm9,%xmm5
addpd  %xmm3,%xmm5
movapd %xmm6,%xmm1
unpcklpd %xmm5,%xmm1
movapd %xmm7,%xmm2
unpcklpd %xmm4,%xmm2
movhlps %xmm6,%xmm5
movhlps %xmm7,%xmm4
movups %xmm4,0x20(%r8,%r11,1)
movups %xmm5,0x30(%r8,%r11,1)
movupd %xmm2,(%r8,%r11,1)
movupd %xmm1,0x10(%r8,%r11,1)
add    $0x40,%r11
add    $0xfffffffffffffffe,%r13
jne    1460 <matmult_vpur_notrans+0x360>
mov    -0x64(%rsp),%eax
mov    -0x50(%rsp),%rbp
cmp    0x10(%rsp),%rbp
mov    %rdi,%r11
mov    0x8(%rsp),%r13
je     1300 <matmult_vpur_notrans+0x200>
jmp    13c6 <matmult_vpur_notrans+0x2c6>
movsd  0x10(%rdx),%xmm1
mulsd  %xmm10,%xmm1
mov    -0x48(%rsp),%rax
addsd  (%rax),%xmm1
movsd  %xmm1,(%rax)
lea    0x18(%rdx),%r15
mov    0x18(%rsp),%r8
test   %ebx,%ebx
jg     1398 <matmult_vpur_notrans+0x298>
jmp    1300 <matmult_vpur_notrans+0x200>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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