// File: /workspace/github_cpp_makefile_star2_O3/Luthaf[P]zeopp/./Eigen/src/Core/products/GeneralMatrixVector.h
// Function: Eigen::internal::general_matrix_vector_product<long, double, 0, false, double, false, 0>::run(long, long, double const*, long, double const*, long, double*, long, double)
EIGEN_DONT_INLINE void general_matrix_vector_product<Index,LhsScalar,ColMajor,ConjugateLhs,RhsScalar,ConjugateRhs,Version>::run(
  Index rows, Index cols,
  const LhsScalar* lhs, Index lhsStride,
  const RhsScalar* rhs, Index rhsIncr,
  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(A0,A13,A2) \
    pstore(&res[j], \
      padd(pload<ResPacket>(&res[j]), \
        padd( \
          padd(pcj.pmul(EIGEN_CAT(ploa , A0)<LhsPacket>(&lhs0[j]),    ptmp0), \
                  pcj.pmul(EIGEN_CAT(ploa , A13)<LhsPacket>(&lhs1[j]),   ptmp1)), \
          padd(pcj.pmul(EIGEN_CAT(ploa , A2)<LhsPacket>(&lhs2[j]),    ptmp2), \
                  pcj.pmul(EIGEN_CAT(ploa , A13)<LhsPacket>(&lhs3[j]),   ptmp3)) )))

  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;
  
  // 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_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 = internal::first_aligned(lhs,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( (size_t(lhs)%sizeof(LhsScalar)) || (size_t(res)%sizeof(ResScalar)) )
  {
    alignedSize = 0;
    alignedStart = 0;
  }
  else if (LhsPacketSize>1)
  {
    eigen_internal_assert(size_t(lhs+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(lhs+alignedStart+lhsStride*skipColumns)%sizeof(LhsPacket))==0);
  }
  else if(Vectorizable)
  {
    alignedStart = 0;
    alignedSize = size;
    alignmentPattern = AllAligned;
  }

  Index offset1 = (FirstAligned && alignmentStep==1?3:1);
  Index offset3 = (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*rhsIncr]),
              ptmp1 = pset1<RhsPacket>(alpha*rhs[(i+offset1)*rhsIncr]),
              ptmp2 = pset1<RhsPacket>(alpha*rhs[(i+2)*rhsIncr]),
              ptmp3 = pset1<RhsPacket>(alpha*rhs[(i+offset3)*rhsIncr]);

    // this helps a lot generating better binary code
    const LhsScalar *lhs0 = lhs + i*lhsStride,     *lhs1 = lhs + (i+offset1)*lhsStride,
                    *lhs2 = lhs + (i+2)*lhsStride, *lhs3 = lhs + (i+offset3)*lhsStride;

    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(d,d,d);
            break;
          case EvenAligned:
            for (Index j = alignedStart; j<alignedSize; j+=ResPacketSize)
              _EIGEN_ACCUMULATE_PACKETS(d,du,d);
            break;
          case FirstAligned:
          {
            Index j = alignedStart;
            if(peels>1)
            {
              LhsPacket A00, A01, A02, A03, A10, A11, A12, A13;
              ResPacket T0, T1;

              A01 = pload<LhsPacket>(&lhs1[alignedStart-1]);
              A02 = pload<LhsPacket>(&lhs2[alignedStart-2]);
              A03 = pload<LhsPacket>(&lhs3[alignedStart-3]);

              for (; j<peeledSize; j+=peels*ResPacketSize)
              {
                A11 = pload<LhsPacket>(&lhs1[j-1+LhsPacketSize]);  palign<1>(A01,A11);
                A12 = pload<LhsPacket>(&lhs2[j-2+LhsPacketSize]);  palign<2>(A02,A12);
                A13 = pload<LhsPacket>(&lhs3[j-3+LhsPacketSize]);  palign<3>(A03,A13);

                A00 = pload<LhsPacket>(&lhs0[j]);
                A10 = pload<LhsPacket>(&lhs0[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 = pload<LhsPacket>(&lhs1[j-1+2*LhsPacketSize]);  palign<1>(A11,A01);
                T0  = pcj.pmadd(A02, ptmp2, T0);
                A02 = pload<LhsPacket>(&lhs2[j-2+2*LhsPacketSize]);  palign<2>(A12,A02);
                T0  = pcj.pmadd(A03, ptmp3, T0);
                pstore(&res[j],T0);
                A03 = pload<LhsPacket>(&lhs3[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(d,du,du);
            break;
          }
          default:
            for (Index j = alignedStart; j<alignedSize; j+=ResPacketSize)
              _EIGEN_ACCUMULATE_PACKETS(du,du,du);
            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*rhsIncr]);
      const LhsScalar* lhs0 = lhs + k*lhsStride;

      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 ((size_t(lhs0+alignedStart)%sizeof(LhsPacket))==0)
          for (Index i = alignedStart;i<alignedSize;i+=ResPacketSize)
            pstore(&res[i], pcj.pmadd(pload<LhsPacket>(&lhs0[i]), ptmp0, pload<ResPacket>(&res[i])));
        else
          for (Index i = alignedStart;i<alignedSize;i+=ResPacketSize)
            pstore(&res[i], pcj.pmadd(ploadu<LhsPacket>(&lhs0[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:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x28, %rsp
movq %r9, -0x40(%rsp)
movq %r8, -0x48(%rsp)
movq %rcx, -0x60(%rsp)
movq %rdx, -0x68(%rsp)
movq 0x60(%rsp), %rax
movl %eax, %r13d
shrl $0x3, %r13d
andl $0x1, %r13d
cmpq %rdi, %r13
cmovgeq %rdi, %r13
testb $0x7, %al
cmovneq %rdi, %r13
movq %rdi, %r12
subq %r13, %r12
andq $-0x2, %r12
movl %ecx, %r8d
andl $0x1, %r8d
movl %eax, %r10d
orl %edx, %r10d
movl %edx, %r9d
shrl $0x3, %r9d
andl $0x1, %r9d
cmpq %rdi, %r9
cmovgeq %rdi, %r9
testb $0x7, %dl
cmovneq %rdi, %r9
testb $0x7, %r10b
je 0x10df26
xorl %ecx, %ecx
movq %r8, -0x58(%rsp)
xorl %r10d, %r10d
xorl %r11d, %r11d
jmp 0x10df82
leaq (%r12,%r13), %r10
movb $0x1, %bl
xorl %r11d, %r11d
movq %r11, %r14
imulq %r8, %r14
addq %r9, %r14
movq %r14, %r15
shrq $0x3f, %r15
addq %r14, %r15
andq $-0x2, %r15
subq %r15, %r14
cmpq %r14, %r13
je 0x10df70
movl $0x1, %r11d
testb $0x1, %bl
movl $0x0, %ebx
jne 0x10df2f
movl $0x3, %ecx
movq %rcx, -0x58(%rsp)
movq %r13, %r11
xorl %ecx, %ecx
jmp 0x10df82
cmpq %rsi, %r11
cmovgeq %rsi, %r11
movq %r11, %rcx
movq %r8, -0x58(%rsp)
movq %r13, %r11
xorl %r9d, %r9d
testq %r8, %r8
sete %bpl
movq %rsi, -0x28(%rsp)
movq %rsi, %r14
subq %rcx, %r14
leaq 0x3(%r14), %rbx
testq %r14, %r14
cmovnsq %r14, %rbx
andq $-0x4, %rbx
movq %rcx, -0x70(%rsp)
addq %rcx, %rbx
cmpq $0x4, %r14
jl 0x10e4d7
leaq (%r12,%r13), %rcx
addq $-0x5, %rcx
movq %rcx, -0x50(%rsp)
leaq 0x1(%r8,%r8), %r12
movb %bpl, %r9b
leaq 0x1(%r9,%r9), %rdx
movq -0x70(%rsp), %r9
movq %rdx, -0x10(%rsp)
leaq (%rdx,%r9), %r14
movq -0x60(%rsp), %r8
imulq %r8, %r14
movq -0x68(%rsp), %r15
leaq (%r15,%r14,8), %rsi
movq %r8, %rcx
shlq $0x5, %rcx
movq %rcx, -0x18(%rsp)
leaq 0x2(%r9), %rdx
imulq %r8, %rdx
leaq (%r15,%rdx,8), %rbp
movq %r12, -0x8(%rsp)
leaq (%r12,%r9), %rcx
imulq %r8, %rcx
leaq (%r15,%rcx,8), %r12
imulq %r8, %r9
leaq (%r15,%r9,8), %r13
movq %rdi, %r8
subq %r10, %r8
movq %r8, -0x20(%rsp)
leaq (,%r10,8), %r8
leaq (%r8,%r14,8), %r14
addq %r15, %r14
leaq (%r8,%rdx,8), %rdx
addq %r15, %rdx
movq %rdx, -0x38(%rsp)
leaq (%r8,%rcx,8), %rcx
addq %r15, %rcx
movq %rcx, -0x78(%rsp)
leaq (%rax,%r10,8), %rdx
leaq (%r8,%r9,8), %rcx
addq %r15, %rcx
movq -0x70(%rsp), %r8
movq %rsi, 0x10(%rsp)
movq %rcx, 0x20(%rsp)
movq %r14, -0x30(%rsp)
movq %r8, %rsi
movq -0x40(%rsp), %r15
imulq %r15, %rsi
movq -0x8(%rsp), %rcx
leaq (%r8,%rcx), %r14
movq %r14, -0x80(%rsp)
imulq %r15, %r14
leaq 0x2(%r8), %rcx
movq %rcx, 0x8(%rsp)
imulq %r15, %rcx
movq -0x10(%rsp), %r9
movq %r8, 0x18(%rsp)
addq %r8, %r9
movq %r9, (%rsp)
imulq %r15, %r9
movq -0x48(%rsp), %r15
movsd (%r15,%rsi,8), %xmm1
movq 0x10(%rsp), %rsi
mulsd %xmm0, %xmm1
movsd (%r15,%r14,8), %xmm2
mulsd %xmm0, %xmm2
movsd (%r15,%rcx,8), %xmm3
mulsd %xmm0, %xmm3
movsd (%r15,%r9,8), %xmm4
mulsd %xmm0, %xmm4
testq %r11, %r11
jle 0x10e14a
xorl %r8d, %r8d
movsd (%r13,%r8,8), %xmm5
mulsd %xmm1, %xmm5
addsd (%rax,%r8,8), %xmm5
movsd %xmm5, (%rax,%r8,8)
movsd (%r12,%r8,8), %xmm6
mulsd %xmm2, %xmm6
addsd %xmm5, %xmm6
movsd %xmm6, (%rax,%r8,8)
movsd (%rbp,%r8,8), %xmm5
mulsd %xmm3, %xmm5
addsd %xmm6, %xmm5
movsd %xmm5, (%rax,%r8,8)
movsd (%rsi,%r8,8), %xmm6
mulsd %xmm4, %xmm6
addsd %xmm5, %xmm6
movsd %xmm6, (%rax,%r8,8)
incq %r8
cmpq %r8, %r11
jne 0x10e0ee
cmpq %r11, %r10
movq -0x30(%rsp), %r14
movq -0x38(%rsp), %r15
movq -0x78(%rsp), %r9
jle 0x10e418
movapd %xmm1, %xmm5
unpcklpd %xmm1, %xmm5            # xmm5 = xmm5[0],xmm1[0]
movapd %xmm2, %xmm6
unpcklpd %xmm2, %xmm6            # xmm6 = xmm6[0],xmm2[0]
movapd %xmm3, %xmm7
unpcklpd %xmm3, %xmm7            # xmm7 = xmm7[0],xmm3[0]
movapd %xmm4, %xmm8
unpcklpd %xmm4, %xmm8            # xmm8 = xmm8[0],xmm4[0]
movq -0x58(%rsp), %rcx
leaq 0x49718(%rip), %r8      # 0x1578a8
movslq (%r8,%rcx,4), %rcx
addq %r8, %rcx
jmpq *%rcx
movq %r11, %r8
movapd (%r13,%r8,8), %xmm9
mulpd %xmm5, %xmm9
movapd (%r12,%r8,8), %xmm10
mulpd %xmm6, %xmm10
addpd %xmm9, %xmm10
movapd (%rbp,%r8,8), %xmm9
mulpd %xmm7, %xmm9
movapd (%rsi,%r8,8), %xmm11
mulpd %xmm8, %xmm11
addpd %xmm9, %xmm11
addpd %xmm10, %xmm11
addpd (%rax,%r8,8), %xmm11
movapd %xmm11, (%rax,%r8,8)
addq $0x2, %r8
cmpq %r10, %r8
jl 0x10e19c
jmp 0x10e418
movq %r11, %r8
cmpq -0x50(%rsp), %r11
jge 0x10e362
movq -0x60(%rsp), %r9
movq -0x80(%rsp), %r14
imulq %r9, %r14
movq 0x8(%rsp), %r8
imulq %r9, %r8
movq (%rsp), %rcx
imulq %r9, %rcx
movq -0x68(%rsp), %r9
leaq (%r9,%r14,8), %r14
movq %r14, -0x80(%rsp)
movq -0x30(%rsp), %r14
leaq (%r9,%r8,8), %r8
leaq (%r9,%rcx,8), %r9
movapd -0x18(%r9,%r11,8), %xmm9
movq -0x78(%rsp), %r9
movapd -0x10(%r8,%r11,8), %xmm11
movq -0x80(%rsp), %rcx
movapd -0x8(%rcx,%r11,8), %xmm12
movq -0x50(%rsp), %rcx
movq %r11, %r8
movapd 0x8(%r12,%r8,8), %xmm10
shufpd $0x1, %xmm10, %xmm12    # xmm12 = xmm12[1],xmm10[0]
movapd (%r13,%r8,8), %xmm14
mulpd %xmm5, %xmm14
addpd (%rax,%r8,8), %xmm14
movapd 0x10(%r13,%r8,8), %xmm13
mulpd %xmm5, %xmm13
addpd 0x10(%rax,%r8,8), %xmm13
mulpd %xmm6, %xmm12
addpd %xmm14, %xmm12
mulpd %xmm7, %xmm11
addpd %xmm12, %xmm11
movapd 0x18(%r12,%r8,8), %xmm12
shufpd $0x1, %xmm12, %xmm10    # xmm10 = xmm10[1],xmm12[0]
mulpd %xmm8, %xmm9
addpd %xmm11, %xmm9
movapd 0x10(%rbp,%r8,8), %xmm11
movapd (%rbp,%r8,8), %xmm14
mulpd %xmm7, %xmm14
movapd -0x8(%rsi,%r8,8), %xmm15
mulpd %xmm8, %xmm15
movapd %xmm9, (%rax,%r8,8)
movapd 0x8(%rsi,%r8,8), %xmm9
mulpd %xmm6, %xmm10
addpd %xmm13, %xmm10
addpd %xmm14, %xmm10
addpd %xmm15, %xmm10
movapd %xmm10, 0x10(%rax,%r8,8)
addq $0x4, %r8
cmpq %rcx, %r8
jl 0x10e262
jmp 0x10e362
movapd (%r13,%r8,8), %xmm9
mulpd %xmm5, %xmm9
movupd (%r12,%r8,8), %xmm10
mulpd %xmm6, %xmm10
addpd %xmm9, %xmm10
movupd (%rbp,%r8,8), %xmm9
mulpd %xmm7, %xmm9
movupd (%rsi,%r8,8), %xmm11
mulpd %xmm8, %xmm11
addpd %xmm9, %xmm11
addpd %xmm10, %xmm11
addpd (%rax,%r8,8), %xmm11
movapd %xmm11, (%rax,%r8,8)
addq $0x2, %r8
cmpq %r10, %r8
jl 0x10e315
jmp 0x10e418
movq %r11, %r8
movupd (%r13,%r8,8), %xmm9
mulpd %xmm5, %xmm9
movupd (%r12,%r8,8), %xmm10
mulpd %xmm6, %xmm10
addpd %xmm9, %xmm10
movupd (%rbp,%r8,8), %xmm9
mulpd %xmm7, %xmm9
movupd (%rsi,%r8,8), %xmm11
mulpd %xmm8, %xmm11
addpd %xmm9, %xmm11
addpd %xmm10, %xmm11
addpd (%rax,%r8,8), %xmm11
movapd %xmm11, (%rax,%r8,8)
addq $0x2, %r8
cmpq %r10, %r8
jl 0x10e36f
jmp 0x10e418
movq %r11, %r8
movapd (%r13,%r8,8), %xmm9
mulpd %xmm5, %xmm9
movupd (%r12,%r8,8), %xmm10
mulpd %xmm6, %xmm10
addpd %xmm9, %xmm10
movapd (%rbp,%r8,8), %xmm9
mulpd %xmm7, %xmm9
movupd (%rsi,%r8,8), %xmm11
mulpd %xmm8, %xmm11
addpd %xmm9, %xmm11
addpd %xmm10, %xmm11
addpd (%rax,%r8,8), %xmm11
movapd %xmm11, (%rax,%r8,8)
addq $0x2, %r8
cmpq %r10, %r8
jl 0x10e3c6
movq %r12, -0x78(%rsp)
movq %rbp, -0x80(%rsp)
cmpq %r10, %rdi
movq -0x20(%rsp), %rcx
movq 0x20(%rsp), %rbp
jle 0x10e48f
xorl %r8d, %r8d
movsd (%rbp,%r8,8), %xmm5
mulsd %xmm1, %xmm5
addsd (%rdx,%r8,8), %xmm5
movsd %xmm5, (%rdx,%r8,8)
movsd (%r9,%r8,8), %xmm6
mulsd %xmm2, %xmm6
addsd %xmm5, %xmm6
movsd %xmm6, (%rdx,%r8,8)
movsd (%r15,%r8,8), %xmm5
mulsd %xmm3, %xmm5
addsd %xmm6, %xmm5
movsd %xmm5, (%rdx,%r8,8)
movsd (%r14,%r8,8), %xmm6
mulsd %xmm4, %xmm6
addsd %xmm5, %xmm6
movsd %xmm6, (%rdx,%r8,8)
incq %r8
cmpq %r8, %rcx
jne 0x10e434
movq 0x18(%rsp), %r8
addq $0x4, %r8
movq -0x18(%rsp), %rcx
addq %rcx, %rsi
addq %rcx, -0x80(%rsp)
movq -0x78(%rsp), %r12
addq %rcx, %r12
addq %rcx, %r13
addq %rcx, %r14
addq %rcx, %r15
movq %r15, -0x38(%rsp)
addq %rcx, %r9
movq %r9, -0x78(%rsp)
addq %rcx, %rbp
cmpq %rbx, %r8
movq %rbp, %rcx
movq -0x80(%rsp), %rbp
jl 0x10e065
movq -0x60(%rsp), %r15
leaq (,%r15,8), %rcx
movq -0x40(%rsp), %r14
movq -0x48(%rsp), %rsi
movq -0x68(%rsp), %r12
movq -0x28(%rsp), %r13
movq -0x70(%rsp), %r8
movq %r8, %rdx
cmpq %r13, %rbx
jge 0x10e5ec
movq %rcx, %r8
imulq %rbx, %r8
addq %r12, %r8
movq %rbx, %r9
imulq %r14, %r9
movsd (%rsi,%r9,8), %xmm1
mulsd %xmm0, %xmm1
testq %r11, %r11
jle 0x10e54a
xorl %r9d, %r9d
movsd (%r8,%r9,8), %xmm2
mulsd %xmm1, %xmm2
addsd (%rax,%r9,8), %xmm2
movsd %xmm2, (%rax,%r9,8)
incq %r9
cmpq %r9, %r11
jne 0x10e52c
movq %rbx, %r9
imulq %r15, %r9
movapd %xmm1, %xmm2
unpcklpd %xmm1, %xmm2            # xmm2 = xmm2[0],xmm1[0]
leaq (%r12,%r9,8), %r9
leal (%r9,%r11,8), %r9d
testb $0xf, %r9b
je 0x10e590
movq %r11, %r9
cmpq %r10, %r11
jge 0x10e5b7
movupd (%r8,%r9,8), %xmm3
mulpd %xmm2, %xmm3
addpd (%rax,%r9,8), %xmm3
movapd %xmm3, (%rax,%r9,8)
addq $0x2, %r9
cmpq %r10, %r9
jl 0x10e56f
jmp 0x10e5b7
movq %r11, %r9
cmpq %r10, %r11
jge 0x10e5b7
movapd (%r8,%r9,8), %xmm3
mulpd %xmm2, %xmm3
addpd (%rax,%r9,8), %xmm3
movapd %xmm3, (%rax,%r9,8)
addq $0x2, %r9
cmpq %r10, %r9
jl 0x10e598
movq %r10, %r9
cmpq %r10, %rdi
jle 0x10e5dd
movsd (%r8,%r9,8), %xmm2
mulsd %xmm1, %xmm2
addsd (%rax,%r9,8), %xmm2
movsd %xmm2, (%rax,%r9,8)
incq %r9
cmpq %r9, %rdi
jne 0x10e5bf
incq %rbx
addq %rcx, %r8
cmpq %r13, %rbx
jne 0x10e513
movl $0x0, %r8d
movq %rdx, %r13
movl $0x0, %ebx
testq %rdx, %rdx
jne 0x10e4fd
addq $0x28, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq