// File: /workspace/github_cpp_makefile_star2_O2/canzarlab[P]Trajan/./Eigen/src/Core/products/GeneralMatrixVector.h
// Function: Eigen::internal::general_matrix_vector_product<long, double, Eigen::internal::const_blas_data_mapper<double, long, 0>, 0, false, double, Eigen::internal::const_blas_data_mapper<double, long, 1>, false, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<double, long, 0> const&, Eigen::internal::const_blas_data_mapper<double, long, 1> const&, double*, long, double)
EIGEN_DONT_INLINE void general_matrix_vector_product<Index,LhsScalar,LhsMapper,ColMajor,ConjugateLhs,RhsScalar,RhsMapper,ConjugateRhs,Version>::run(
  Index rows, Index cols,
  const LhsMapper& lhs,
  const RhsMapper& rhs,
        ResScalar* res, Index resIncr,
  RhsScalar alpha)
{
  EIGEN_UNUSED_VARIABLE(resIncr);
  eigen_internal_assert(resIncr==1);
  #ifdef _EIGEN_ACCUMULATE_PACKETS
  #error _EIGEN_ACCUMULATE_PACKETS has already been defined
  #endif
  #define _EIGEN_ACCUMULATE_PACKETS(Alignment0,Alignment13,Alignment2) \
    pstore(&res[j], \
      padd(pload<ResPacket>(&res[j]), \
        padd( \
      padd(pcj.pmul(lhs0.template load<LhsPacket, Alignment0>(j),    ptmp0), \
      pcj.pmul(lhs1.template load<LhsPacket, Alignment13>(j),   ptmp1)),   \
      padd(pcj.pmul(lhs2.template load<LhsPacket, Alignment2>(j),    ptmp2), \
      pcj.pmul(lhs3.template load<LhsPacket, Alignment13>(j),   ptmp3)) )))

  typedef typename LhsMapper::VectorMapper LhsScalars;

  conj_helper<LhsScalar,RhsScalar,ConjugateLhs,ConjugateRhs> cj;
  conj_helper<LhsPacket,RhsPacket,ConjugateLhs,ConjugateRhs> pcj;
  if(ConjugateRhs)
    alpha = numext::conj(alpha);

  enum { AllAligned = 0, EvenAligned, FirstAligned, NoneAligned };
  const Index columnsAtOnce = 4;
  const Index peels = 2;
  const Index LhsPacketAlignedMask = LhsPacketSize-1;
  const Index ResPacketAlignedMask = ResPacketSize-1;
//  const Index PeelAlignedMask = ResPacketSize*peels-1;
  const Index size = rows;

  const Index lhsStride = lhs.stride();

  // How many coeffs of the result do we have to skip to be aligned.
  // Here we assume data are at least aligned on the base scalar type.
  Index alignedStart = internal::first_default_aligned(res,size);
  Index alignedSize = ResPacketSize>1 ? alignedStart + ((size-alignedStart) & ~ResPacketAlignedMask) : 0;
  const Index peeledSize = alignedSize - RhsPacketSize*peels - RhsPacketSize + 1;

  const Index alignmentStep = LhsPacketSize>1 ? (LhsPacketSize - lhsStride % LhsPacketSize) & LhsPacketAlignedMask : 0;
  Index alignmentPattern = alignmentStep==0 ? AllAligned
                       : alignmentStep==(LhsPacketSize/2) ? EvenAligned
                       : FirstAligned;

  // we cannot assume the first element is aligned because of sub-matrices
  const Index lhsAlignmentOffset = lhs.firstAligned(size);

  // find how many columns do we have to skip to be aligned with the result (if possible)
  Index skipColumns = 0;
  // if the data cannot be aligned (TODO add some compile time tests when possible, e.g. for floats)
  if( (lhsAlignmentOffset < 0) || (lhsAlignmentOffset == size) || (UIntPtr(res)%sizeof(ResScalar)) )
  {
    alignedSize = 0;
    alignedStart = 0;
    alignmentPattern = NoneAligned;
  }
  else if(LhsPacketSize > 4)
  {
    // TODO: extend the code to support aligned loads whenever possible when LhsPacketSize > 4.
    // Currently, it seems to be better to perform unaligned loads anyway
    alignmentPattern = NoneAligned;
  }
  else if (LhsPacketSize>1)
  {
  //    eigen_internal_assert(size_t(firstLhs+lhsAlignmentOffset)%sizeof(LhsPacket)==0 || size<LhsPacketSize);

    while (skipColumns<LhsPacketSize &&
          alignedStart != ((lhsAlignmentOffset + alignmentStep*skipColumns)%LhsPacketSize))
      ++skipColumns;
    if (skipColumns==LhsPacketSize)
    {
      // nothing can be aligned, no need to skip any column
      alignmentPattern = NoneAligned;
      skipColumns = 0;
    }
    else
    {
      skipColumns = (std::min)(skipColumns,cols);
      // note that the skiped columns are processed later.
    }

    /*    eigen_internal_assert(  (alignmentPattern==NoneAligned)
                      || (skipColumns + columnsAtOnce >= cols)
                      || LhsPacketSize > size
                      || (size_t(firstLhs+alignedStart+lhsStride*skipColumns)%sizeof(LhsPacket))==0);*/
  }
  else if(Vectorizable)
  {
    alignedStart = 0;
    alignedSize = size;
    alignmentPattern = AllAligned;
  }

  const Index offset1 = (alignmentPattern==FirstAligned && alignmentStep==1)?3:1;
  const Index offset3 = (alignmentPattern==FirstAligned && alignmentStep==1)?1:3;

  Index columnBound = ((cols-skipColumns)/columnsAtOnce)*columnsAtOnce + skipColumns;
  for (Index i=skipColumns; i<columnBound; i+=columnsAtOnce)
  {
    RhsPacket ptmp0 = pset1<RhsPacket>(alpha*rhs(i, 0)),
              ptmp1 = pset1<RhsPacket>(alpha*rhs(i+offset1, 0)),
              ptmp2 = pset1<RhsPacket>(alpha*rhs(i+2, 0)),
              ptmp3 = pset1<RhsPacket>(alpha*rhs(i+offset3, 0));

    // this helps a lot generating better binary code
    const LhsScalars lhs0 = lhs.getVectorMapper(0, i+0),   lhs1 = lhs.getVectorMapper(0, i+offset1),
                     lhs2 = lhs.getVectorMapper(0, i+2),   lhs3 = lhs.getVectorMapper(0, i+offset3);

    if (Vectorizable)
    {
      /* explicit vectorization */
      // process initial unaligned coeffs
      for (Index j=0; j<alignedStart; ++j)
      {
        res[j] = cj.pmadd(lhs0(j), pfirst(ptmp0), res[j]);
        res[j] = cj.pmadd(lhs1(j), pfirst(ptmp1), res[j]);
        res[j] = cj.pmadd(lhs2(j), pfirst(ptmp2), res[j]);
        res[j] = cj.pmadd(lhs3(j), pfirst(ptmp3), res[j]);
      }

      if (alignedSize>alignedStart)
      {
        switch(alignmentPattern)
        {
          case AllAligned:
            for (Index j = alignedStart; j<alignedSize; j+=ResPacketSize)
              _EIGEN_ACCUMULATE_PACKETS(Aligned,Aligned,Aligned);
            break;
          case EvenAligned:
            for (Index j = alignedStart; j<alignedSize; j+=ResPacketSize)
              _EIGEN_ACCUMULATE_PACKETS(Aligned,Unaligned,Aligned);
            break;
          case FirstAligned:
          {
            Index j = alignedStart;
            if(peels>1)
            {
              LhsPacket A00, A01, A02, A03, A10, A11, A12, A13;
              ResPacket T0, T1;

              A01 = lhs1.template load<LhsPacket, Aligned>(alignedStart-1);
              A02 = lhs2.template load<LhsPacket, Aligned>(alignedStart-2);
              A03 = lhs3.template load<LhsPacket, Aligned>(alignedStart-3);

              for (; j<peeledSize; j+=peels*ResPacketSize)
              {
                A11 = lhs1.template load<LhsPacket, Aligned>(j-1+LhsPacketSize);  palign<1>(A01,A11);
                A12 = lhs2.template load<LhsPacket, Aligned>(j-2+LhsPacketSize);  palign<2>(A02,A12);
                A13 = lhs3.template load<LhsPacket, Aligned>(j-3+LhsPacketSize);  palign<3>(A03,A13);

                A00 = lhs0.template load<LhsPacket, Aligned>(j);
                A10 = lhs0.template load<LhsPacket, Aligned>(j+LhsPacketSize);
                T0  = pcj.pmadd(A00, ptmp0, pload<ResPacket>(&res[j]));
                T1  = pcj.pmadd(A10, ptmp0, pload<ResPacket>(&res[j+ResPacketSize]));

                T0  = pcj.pmadd(A01, ptmp1, T0);
                A01 = lhs1.template load<LhsPacket, Aligned>(j-1+2*LhsPacketSize);  palign<1>(A11,A01);
                T0  = pcj.pmadd(A02, ptmp2, T0);
                A02 = lhs2.template load<LhsPacket, Aligned>(j-2+2*LhsPacketSize);  palign<2>(A12,A02);
                T0  = pcj.pmadd(A03, ptmp3, T0);
                pstore(&res[j],T0);
                A03 = lhs3.template load<LhsPacket, Aligned>(j-3+2*LhsPacketSize);  palign<3>(A13,A03);
                T1  = pcj.pmadd(A11, ptmp1, T1);
                T1  = pcj.pmadd(A12, ptmp2, T1);
                T1  = pcj.pmadd(A13, ptmp3, T1);
                pstore(&res[j+ResPacketSize],T1);
              }
            }
            for (; j<alignedSize; j+=ResPacketSize)
              _EIGEN_ACCUMULATE_PACKETS(Aligned,Unaligned,Unaligned);
            break;
          }
          default:
            for (Index j = alignedStart; j<alignedSize; j+=ResPacketSize)
              _EIGEN_ACCUMULATE_PACKETS(Unaligned,Unaligned,Unaligned);
            break;
        }
      }
    } // end explicit vectorization

    /* process remaining coeffs (or all if there is no explicit vectorization) */
    for (Index j=alignedSize; j<size; ++j)
    {
      res[j] = cj.pmadd(lhs0(j), pfirst(ptmp0), res[j]);
      res[j] = cj.pmadd(lhs1(j), pfirst(ptmp1), res[j]);
      res[j] = cj.pmadd(lhs2(j), pfirst(ptmp2), res[j]);
      res[j] = cj.pmadd(lhs3(j), pfirst(ptmp3), res[j]);
    }
  }

  // process remaining first and last columns (at most columnsAtOnce-1)
  Index end = cols;
  Index start = columnBound;
  do
  {
    for (Index k=start; k<end; ++k)
    {
      RhsPacket ptmp0 = pset1<RhsPacket>(alpha*rhs(k, 0));
      const LhsScalars lhs0 = lhs.getVectorMapper(0, k);

      if (Vectorizable)
      {
        /* explicit vectorization */
        // process first unaligned result's coeffs
        for (Index j=0; j<alignedStart; ++j)
          res[j] += cj.pmul(lhs0(j), pfirst(ptmp0));
        // process aligned result's coeffs
        if (lhs0.template aligned<LhsPacket>(alignedStart))
          for (Index i = alignedStart;i<alignedSize;i+=ResPacketSize)
            pstore(&res[i], pcj.pmadd(lhs0.template load<LhsPacket, Aligned>(i), ptmp0, pload<ResPacket>(&res[i])));
        else
          for (Index i = alignedStart;i<alignedSize;i+=ResPacketSize)
            pstore(&res[i], pcj.pmadd(lhs0.template load<LhsPacket, Unaligned>(i), ptmp0, pload<ResPacket>(&res[i])));
      }

      // process remaining scalars (or all if no explicit vectorization)
      for (Index i=alignedSize; i<size; ++i)
        res[i] += cj.pmul(lhs0(i), pfirst(ptmp0));
    }
    if (skipColumns)
    {
      start = 0;
      end = skipColumns;
      skipColumns = 0;
    }
    else
      break;
  } while(Vectorizable);
  #undef _EIGEN_ACCUMULATE_PACKETS
}

// Assembly:
movq %rdx, %r9
movq %rdi, %rax
testb $0x7, %r8b
jne 0x21e44
movl %r8d, %eax
shrl $0x3, %eax
andl $0x1, %eax
cmpq %rdi, %rax
cmovgeq %rdi, %rax
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
movq %rdi, %rdx
subq %rax, %rdx
movq (%r9), %rbx
testb $0x7, %bl
je 0x21e68
movabsq $-0x1, %rbx
jmp 0x21e75
shrl $0x3, %ebx
andl $0x1, %ebx
cmpq %rdi, %rbx
cmovgeq %rdi, %rbx
andq $-0x2, %rdx
xorl %r14d, %r14d
movabsq $0x3, %r10
movq %r10, -0x68(%rsp)
testb $0x7, %r8b
jne 0x21ead
testq %rbx, %rbx
js 0x21ead
cmpq %rdi, %rbx
je 0x21ead
movq 0x8(%r9), %r11
leaq (%rdx,%rax), %r10
cmpq %rbx, %rax
jne 0x21eb5
xorl %r14d, %r14d
jmp 0x21eca
xorl %r10d, %r10d
xorl %r15d, %r15d
jmp 0x21ee2
addl %r11d, %ebx
andl $0x1, %ebx
cmpq %rbx, %rax
jne 0x21edc
movabsq $0x1, %r14
andl $0x1, %r11d
cmpq %rsi, %r14
cmovgeq %rsi, %r14
movq %r11, -0x68(%rsp)
jmp 0x21edf
xorl %r14d, %r14d
movq %rax, %r15
addq %rdx, %rax
addq $-0x5, %rax
movq %rax, -0x8(%rsp)
movq %rsi, %rax
subq %r14, %rax
movabsq $0x4, %r11
cqto
idivq %r11
leaq (%r14,%rax,4), %rdx
xorl %eax, %eax
testq %r15, %r15
cmovgq %r15, %rax
leaq (,%r14,8), %rbp
leaq 0x8(,%r14,8), %r11
movq %r11, -0x50(%rsp)
leaq 0x10(,%r14,8), %r11
movq %r11, -0x58(%rsp)
leaq 0x18(,%r14,8), %r11
movq %r11, -0x60(%rsp)
movq %r14, %r13
movq %rsi, -0x38(%rsp)
movq %r14, -0x40(%rsp)
cmpq %rdx, %r13
jge 0x223e2
movq (%rcx), %rsi
movq 0x8(%rcx), %r11
movq %r11, %rbx
imulq %r13, %rbx
movsd (%rsi,%rbx,8), %xmm1
mulsd %xmm0, %xmm1
movapd %xmm1, %xmm5
leaq 0x1(%r13), %r14
movq %r11, %rbx
imulq %r14, %rbx
movsd (%rsi,%rbx,8), %xmm2
mulsd %xmm0, %xmm2
movapd %xmm2, %xmm6
leaq 0x2(%r13), %r12
movq %r11, %rbx
imulq %r12, %rbx
movsd (%rsi,%rbx,8), %xmm3
mulsd %xmm0, %xmm3
movapd %xmm3, %xmm7
movq %r13, -0x18(%rsp)
addq $0x3, %r13
imulq %r13, %r11
movsd (%rsi,%r11,8), %xmm4
mulsd %xmm0, %xmm4
movapd %xmm4, %xmm8
movq (%r9), %r11
movq 0x8(%r9), %rbx
imulq %rbx, %r14
leaq (%r11,%r14,8), %rsi
movq %rsi, -0x70(%rsp)
imulq %rbx, %r12
leaq (%r11,%r12,8), %rsi
movq %rsi, -0x78(%rsp)
imulq %rbx, %r13
leaq (%r11,%r13,8), %rsi
movq %rsi, -0x80(%rsp)
movq %rbx, %rsi
movq %rbp, -0x10(%rsp)
imulq %rbp, %rsi
movq %rsi, -0x20(%rsp)
leaq (%r11,%rsi), %r13
movq %rbx, %rsi
imulq -0x50(%rsp), %rsi
movq %rsi, -0x28(%rsp)
leaq (%r11,%rsi), %rbp
movq %rbx, %rsi
imulq -0x58(%rsp), %rsi
movq %rsi, -0x30(%rsp)
leaq (%r11,%rsi), %r14
imulq -0x60(%rsp), %rbx
leaq (%r11,%rbx), %r12
xorl %esi, %esi
cmpq %rsi, %rax
je 0x22094
movsd (%r13,%rsi,8), %xmm9
mulsd %xmm1, %xmm9
addsd (%r8,%rsi,8), %xmm9
movsd %xmm9, (%r8,%rsi,8)
movsd (%rbp,%rsi,8), %xmm10
mulsd %xmm2, %xmm10
addsd %xmm9, %xmm10
movsd %xmm10, (%r8,%rsi,8)
movsd (%r14,%rsi,8), %xmm9
mulsd %xmm3, %xmm9
addsd %xmm10, %xmm9
movsd %xmm9, (%r8,%rsi,8)
movsd (%r12,%rsi,8), %xmm10
mulsd %xmm4, %xmm10
addsd %xmm9, %xmm10
movsd %xmm10, (%r8,%rsi,8)
incq %rsi
jmp 0x2202f
cmpq %r15, %r10
jle 0x22299
movq %r11, -0x48(%rsp)
unpcklpd %xmm5, %xmm5            # xmm5 = xmm5[0,0]
unpcklpd %xmm6, %xmm6            # xmm6 = xmm6[0,0]
unpcklpd %xmm7, %xmm7            # xmm7 = xmm7[0,0]
unpcklpd %xmm8, %xmm8            # xmm8 = xmm8[0,0]
movq -0x68(%rsp), %rsi
leaq 0x319e1(%rip), %r11     # 0x53aa0
movslq (%r11,%rsi,4), %rsi
addq %r11, %rsi
jmpq *%rsi
movq %r15, %rsi
cmpq %r10, %rsi
jge 0x22299
movapd (%r13,%rsi,8), %xmm9
mulpd %xmm5, %xmm9
movapd (%rbp,%rsi,8), %xmm10
mulpd %xmm6, %xmm10
addpd %xmm9, %xmm10
movapd (%r14,%rsi,8), %xmm9
mulpd %xmm7, %xmm9
movapd (%r12,%rsi,8), %xmm11
mulpd %xmm8, %xmm11
addpd %xmm9, %xmm11
addpd %xmm10, %xmm11
addpd (%r8,%rsi,8), %xmm11
movapd %xmm11, (%r8,%rsi,8)
addq $0x2, %rsi
jmp 0x220cb
movq -0x70(%rsp), %rsi
movapd -0x8(%rsi,%r15,8), %xmm12
movq -0x78(%rsp), %rsi
movapd -0x10(%rsi,%r15,8), %xmm10
movq -0x80(%rsp), %rsi
movapd -0x18(%rsi,%r15,8), %xmm9
movq -0x48(%rsp), %r11
movq -0x28(%rsp), %rsi
leaq 0x18(%r11,%rsi), %rsi
movq %rsi, -0x70(%rsp)
leaq 0x8(%r11,%rbx), %rsi
movq %rsi, -0x78(%rsp)
movq -0x30(%rsp), %rsi
leaq 0x10(%r11,%rsi), %rsi
movq %rsi, -0x80(%rsp)
movq -0x20(%rsp), %rsi
leaq (%r11,%rsi), %rbx
addq $0x10, %rbx
movq %r15, %r11
cmpq -0x8(%rsp), %r11
jge 0x22294
movq -0x70(%rsp), %rsi
movapd -0x10(%rsi,%r11,8), %xmm11
movapd -0x10(%rbx,%r11,8), %xmm14
mulpd %xmm5, %xmm14
addpd (%r8,%r11,8), %xmm14
movapd (%rbx,%r11,8), %xmm13
mulpd %xmm5, %xmm13
addpd 0x10(%r8,%r11,8), %xmm13
shufpd $0x1, %xmm11, %xmm12    # xmm12 = xmm12[1],xmm11[0]
mulpd %xmm6, %xmm12
addpd %xmm14, %xmm12
mulpd %xmm7, %xmm10
addpd %xmm12, %xmm10
movapd (%rsi,%r11,8), %xmm12
shufpd $0x1, %xmm12, %xmm11    # xmm11 = xmm11[1],xmm12[0]
mulpd %xmm8, %xmm9
addpd %xmm10, %xmm9
movq -0x80(%rsp), %rsi
movapd (%rsi,%r11,8), %xmm10
movapd -0x10(%rsi,%r11,8), %xmm14
mulpd %xmm7, %xmm14
movq -0x78(%rsp), %rsi
movapd -0x10(%rsi,%r11,8), %xmm15
mulpd %xmm8, %xmm15
movapd %xmm9, (%r8,%r11,8)
movapd (%rsi,%r11,8), %xmm9
mulpd %xmm6, %xmm11
addpd %xmm13, %xmm11
addpd %xmm14, %xmm11
addpd %xmm15, %xmm11
movapd %xmm11, 0x10(%r8,%r11,8)
addq $0x4, %r11
jmp 0x22184
movapd (%r13,%r11,8), %xmm9
mulpd %xmm5, %xmm9
movupd (%rbp,%r11,8), %xmm10
mulpd %xmm6, %xmm10
addpd %xmm9, %xmm10
movupd (%r14,%r11,8), %xmm9
mulpd %xmm7, %xmm9
movupd (%r12,%r11,8), %xmm11
mulpd %xmm8, %xmm11
addpd %xmm9, %xmm11
addpd %xmm10, %xmm11
addpd (%r8,%r11,8), %xmm11
movapd %xmm11, (%r8,%r11,8)
addq $0x2, %r11
cmpq %r10, %r11
jl 0x22247
movq %r10, %rsi
cmpq %rdi, %rsi
jge 0x222f9
movsd (%r13,%rsi,8), %xmm5
mulsd %xmm1, %xmm5
addsd (%r8,%rsi,8), %xmm5
movsd %xmm5, (%r8,%rsi,8)
movsd (%rbp,%rsi,8), %xmm6
mulsd %xmm2, %xmm6
addsd %xmm5, %xmm6
movsd %xmm6, (%r8,%rsi,8)
movsd (%r14,%rsi,8), %xmm5
mulsd %xmm3, %xmm5
addsd %xmm6, %xmm5
movsd %xmm5, (%r8,%rsi,8)
movsd (%r12,%rsi,8), %xmm6
mulsd %xmm4, %xmm6
addsd %xmm5, %xmm6
movsd %xmm6, (%r8,%rsi,8)
incq %rsi
jmp 0x2229c
movq -0x18(%rsp), %r13
addq $0x4, %r13
movq -0x10(%rsp), %rbp
addq $0x20, %rbp
addq $0x20, -0x50(%rsp)
addq $0x20, -0x58(%rsp)
addq $0x20, -0x60(%rsp)
movq -0x38(%rsp), %rsi
movq -0x40(%rsp), %r14
jmp 0x21f4c
movq %r15, %rsi
cmpq %r10, %rsi
jge 0x22299
movupd (%r13,%rsi,8), %xmm9
mulpd %xmm5, %xmm9
movupd (%rbp,%rsi,8), %xmm10
mulpd %xmm6, %xmm10
addpd %xmm9, %xmm10
movupd (%r14,%rsi,8), %xmm9
mulpd %xmm7, %xmm9
movupd (%r12,%rsi,8), %xmm11
mulpd %xmm8, %xmm11
addpd %xmm9, %xmm11
addpd %xmm10, %xmm11
addpd (%r8,%rsi,8), %xmm11
movapd %xmm11, (%r8,%rsi,8)
addq $0x2, %rsi
jmp 0x2232f
movq %r15, %rsi
cmpq %r10, %rsi
jge 0x22299
movapd (%r13,%rsi,8), %xmm9
mulpd %xmm5, %xmm9
movupd (%rbp,%rsi,8), %xmm10
mulpd %xmm6, %xmm10
addpd %xmm9, %xmm10
movapd (%r14,%rsi,8), %xmm9
mulpd %xmm7, %xmm9
movupd (%r12,%rsi,8), %xmm11
mulpd %xmm8, %xmm11
addpd %xmm9, %xmm11
addpd %xmm10, %xmm11
addpd (%r8,%rsi,8), %xmm11
movapd %xmm11, (%r8,%rsi,8)
addq $0x2, %rsi
jmp 0x2238a
movq %r14, %rbx
leaq (,%rdx,8), %r14
cmpq %rsi, %rdx
jge 0x224cc
movq (%rcx), %r12
movq 0x8(%rcx), %r11
imulq %rdx, %r11
movsd (%r12,%r11,8), %xmm1
mulsd %xmm0, %xmm1
movapd %xmm1, %xmm2
movq (%r9), %r12
movq 0x8(%r9), %r11
movq %r11, %r13
imulq %rdx, %r13
leaq (%r12,%r13,8), %r13
imulq %r14, %r11
addq %r12, %r11
xorl %r12d, %r12d
cmpq %r12, %rax
je 0x2244b
movsd (%r11,%r12,8), %xmm3
mulsd %xmm1, %xmm3
addsd (%r8,%r12,8), %xmm3
movsd %xmm3, (%r8,%r12,8)
incq %r12
jmp 0x2242b
unpcklpd %xmm2, %xmm2            # xmm2 = xmm2[0,0]
leal (%r13,%r15,8), %ebp
movq %r15, %r12
testb $0xf, %bpl
je 0x22498
cmpq %r10, %r12
jge 0x2249d
movupd (%r11,%r12,8), %xmm3
mulpd %xmm2, %xmm3
addpd (%r8,%r12,8), %xmm3
movapd %xmm3, (%r8,%r12,8)
addq $0x2, %r12
jmp 0x2245d
movapd (%r11,%r12,8), %xmm3
mulpd %xmm2, %xmm3
addpd (%r8,%r12,8), %xmm3
movapd %xmm3, (%r8,%r12,8)
addq $0x2, %r12
cmpq %r10, %r12
jl 0x2247e
movq %r10, %r12
cmpq %rdi, %r12
jge 0x224c0
movsd (%r11,%r12,8), %xmm2
mulsd %xmm1, %xmm2
addsd (%r8,%r12,8), %xmm2
movsd %xmm2, (%r8,%r12,8)
incq %r12
jmp 0x224a0
incq %rdx
addq $0x8, %r14
jmp 0x223ed
movl $0x0, %r14d
movq %rbx, %rsi
movl $0x0, %edx
testq %rbx, %rbx
jne 0x223e2
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq