// File: /workspace/github_cpp_makefile_star2_O2/Pyboon[P]sage/./Eigen/src/Core/products/GeneralMatrixVector.h
// Function: Eigen::internal::general_matrix_vector_product<long, std::complex<double>, Eigen::internal::const_blas_data_mapper<std::complex<double>, long, 1>, 1, false, std::complex<double>, Eigen::internal::const_blas_data_mapper<std::complex<double>, long, 0>, true, 0>::run(long, long, Eigen::internal::const_blas_data_mapper<std::complex<double>, long, 1> const&, Eigen::internal::const_blas_data_mapper<std::complex<double>, long, 0> const&, std::complex<double>*, long, std::complex<double>)
EIGEN_DONT_INLINE void general_matrix_vector_product<Index,LhsScalar,LhsMapper,RowMajor,ConjugateLhs,RhsScalar,RhsMapper,ConjugateRhs,Version>::run(
  Index rows, Index cols,
  const LhsMapper& lhs,
  const RhsMapper& rhs,
  ResScalar* res, Index resIncr,
  ResScalar alpha)
{
  eigen_internal_assert(rhs.stride()==1);

  #ifdef _EIGEN_ACCUMULATE_PACKETS
  #error _EIGEN_ACCUMULATE_PACKETS has already been defined
  #endif

  #define _EIGEN_ACCUMULATE_PACKETS(Alignment0,Alignment13,Alignment2) {\
    RhsPacket b = rhs.getVectorMapper(j, 0).template load<RhsPacket, Aligned>(0);  \
    ptmp0 = pcj.pmadd(lhs0.template load<LhsPacket, Alignment0>(j), b, ptmp0); \
    ptmp1 = pcj.pmadd(lhs1.template load<LhsPacket, Alignment13>(j), b, ptmp1); \
    ptmp2 = pcj.pmadd(lhs2.template load<LhsPacket, Alignment2>(j), b, ptmp2); \
    ptmp3 = pcj.pmadd(lhs3.template load<LhsPacket, Alignment13>(j), b, ptmp3); }

  conj_helper<LhsScalar,RhsScalar,ConjugateLhs,ConjugateRhs> cj;
  conj_helper<LhsPacket,RhsPacket,ConjugateLhs,ConjugateRhs> pcj;

  typedef typename LhsMapper::VectorMapper LhsScalars;

  enum { AllAligned=0, EvenAligned=1, FirstAligned=2, NoneAligned=3 };
  const Index rowsAtOnce = 4;
  const Index peels = 2;
  const Index RhsPacketAlignedMask = RhsPacketSize-1;
  const Index LhsPacketAlignedMask = LhsPacketSize-1;
  const Index depth = cols;
  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
  // if that's not the case then vectorization is discarded, see below.
  Index alignedStart = rhs.firstAligned(depth);
  Index alignedSize = RhsPacketSize>1 ? alignedStart + ((depth-alignedStart) & ~RhsPacketAlignedMask) : 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(depth);
  const Index rhsAlignmentOffset = rhs.firstAligned(rows);

  // find how many rows do we have to skip to be aligned with rhs (if possible)
  Index skipRows = 0;
  // if the data cannot be aligned (TODO add some compile time tests when possible, e.g. for floats)
  if( (sizeof(LhsScalar)!=sizeof(RhsScalar)) ||
      (lhsAlignmentOffset < 0) || (lhsAlignmentOffset == depth) ||
      (rhsAlignmentOffset < 0) || (rhsAlignmentOffset == rows) )
  {
    alignedSize = 0;
    alignedStart = 0;
    alignmentPattern = NoneAligned;
  }
  else if(LhsPacketSize > 4)
  {
    // TODO: extend the code to support aligned loads whenever possible when LhsPacketSize > 4.
    alignmentPattern = NoneAligned;
  }
  else if (LhsPacketSize>1)
  {
  //    eigen_internal_assert(size_t(firstLhs+lhsAlignmentOffset)%sizeof(LhsPacket)==0  || depth<LhsPacketSize);

    while (skipRows<LhsPacketSize &&
           alignedStart != ((lhsAlignmentOffset + alignmentStep*skipRows)%LhsPacketSize))
      ++skipRows;
    if (skipRows==LhsPacketSize)
    {
      // nothing can be aligned, no need to skip any column
      alignmentPattern = NoneAligned;
      skipRows = 0;
    }
    else
    {
      skipRows = (std::min)(skipRows,Index(rows));
      // note that the skiped columns are processed later.
    }
    /*    eigen_internal_assert(  alignmentPattern==NoneAligned
                      || LhsPacketSize==1
                      || (skipRows + rowsAtOnce >= rows)
                      || LhsPacketSize > depth
                      || (size_t(firstLhs+alignedStart+lhsStride*skipRows)%sizeof(LhsPacket))==0);*/
  }
  else if(Vectorizable)
  {
    alignedStart = 0;
    alignedSize = depth;
    alignmentPattern = AllAligned;
  }

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

  Index rowBound = ((rows-skipRows)/rowsAtOnce)*rowsAtOnce + skipRows;
  for (Index i=skipRows; i<rowBound; i+=rowsAtOnce)
  {
    // FIXME: what is the purpose of this EIGEN_ALIGN_DEFAULT ??
    EIGEN_ALIGN_MAX ResScalar tmp0 = ResScalar(0);
    ResScalar tmp1 = ResScalar(0), tmp2 = ResScalar(0), tmp3 = ResScalar(0);

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

    if (Vectorizable)
    {
      /* explicit vectorization */
      ResPacket ptmp0 = pset1<ResPacket>(ResScalar(0)), ptmp1 = pset1<ResPacket>(ResScalar(0)),
                ptmp2 = pset1<ResPacket>(ResScalar(0)), ptmp3 = pset1<ResPacket>(ResScalar(0));

      // process initial unaligned coeffs
      // FIXME this loop get vectorized by the compiler !
      for (Index j=0; j<alignedStart; ++j)
      {
        RhsScalar b = rhs(j, 0);
        tmp0 += cj.pmul(lhs0(j),b); tmp1 += cj.pmul(lhs1(j),b);
        tmp2 += cj.pmul(lhs2(j),b); tmp3 += cj.pmul(lhs3(j),b);
      }

      if (alignedSize>alignedStart)
      {
        switch(alignmentPattern)
        {
          case AllAligned:
            for (Index j = alignedStart; j<alignedSize; j+=RhsPacketSize)
              _EIGEN_ACCUMULATE_PACKETS(Aligned,Aligned,Aligned);
            break;
          case EvenAligned:
            for (Index j = alignedStart; j<alignedSize; j+=RhsPacketSize)
              _EIGEN_ACCUMULATE_PACKETS(Aligned,Unaligned,Aligned);
            break;
          case FirstAligned:
          {
            Index j = alignedStart;
            if (peels>1)
            {
              /* Here we proccess 4 rows with with two peeled iterations to hide
               * the overhead of unaligned loads. Moreover unaligned loads are handled
               * using special shift/move operations between the two aligned packets
               * overlaping the desired unaligned packet. This is *much* more efficient
               * than basic unaligned loads.
               */
              LhsPacket A01, A02, A03, A11, A12, A13;
              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*RhsPacketSize)
              {
                RhsPacket b = rhs.getVectorMapper(j, 0).template load<RhsPacket, Aligned>(0);
                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);

                ptmp0 = pcj.pmadd(lhs0.template load<LhsPacket, Aligned>(j), b, ptmp0);
                ptmp1 = pcj.pmadd(A01, b, ptmp1);
                A01 = lhs1.template load<LhsPacket, Aligned>(j-1+2*LhsPacketSize);  palign<1>(A11,A01);
                ptmp2 = pcj.pmadd(A02, b, ptmp2);
                A02 = lhs2.template load<LhsPacket, Aligned>(j-2+2*LhsPacketSize);  palign<2>(A12,A02);
                ptmp3 = pcj.pmadd(A03, b, ptmp3);
                A03 = lhs3.template load<LhsPacket, Aligned>(j-3+2*LhsPacketSize);  palign<3>(A13,A03);

                b = rhs.getVectorMapper(j+RhsPacketSize, 0).template load<RhsPacket, Aligned>(0);
                ptmp0 = pcj.pmadd(lhs0.template load<LhsPacket, Aligned>(j+LhsPacketSize), b, ptmp0);
                ptmp1 = pcj.pmadd(A11, b, ptmp1);
                ptmp2 = pcj.pmadd(A12, b, ptmp2);
                ptmp3 = pcj.pmadd(A13, b, ptmp3);
              }
            }
            for (; j<alignedSize; j+=RhsPacketSize)
              _EIGEN_ACCUMULATE_PACKETS(Aligned,Unaligned,Unaligned);
            break;
          }
          default:
            for (Index j = alignedStart; j<alignedSize; j+=RhsPacketSize)
              _EIGEN_ACCUMULATE_PACKETS(Unaligned,Unaligned,Unaligned);
            break;
        }
        tmp0 += predux(ptmp0);
        tmp1 += predux(ptmp1);
        tmp2 += predux(ptmp2);
        tmp3 += predux(ptmp3);
      }
    } // end explicit vectorization

    // process remaining coeffs (or all if no explicit vectorization)
    // FIXME this loop get vectorized by the compiler !
    for (Index j=alignedSize; j<depth; ++j)
    {
      RhsScalar b = rhs(j, 0);
      tmp0 += cj.pmul(lhs0(j),b); tmp1 += cj.pmul(lhs1(j),b);
      tmp2 += cj.pmul(lhs2(j),b); tmp3 += cj.pmul(lhs3(j),b);
    }
    res[i*resIncr]            += alpha*tmp0;
    res[(i+offset1)*resIncr]  += alpha*tmp1;
    res[(i+2)*resIncr]        += alpha*tmp2;
    res[(i+offset3)*resIncr]  += alpha*tmp3;
  }

  // process remaining first and last rows (at most columnsAtOnce-1)
  Index end = rows;
  Index start = rowBound;
  do
  {
    for (Index i=start; i<end; ++i)
    {
      EIGEN_ALIGN_MAX ResScalar tmp0 = ResScalar(0);
      ResPacket ptmp0 = pset1<ResPacket>(tmp0);
      const LhsScalars lhs0 = lhs.getVectorMapper(i, 0);
      // process first unaligned result's coeffs
      // FIXME this loop get vectorized by the compiler !
      for (Index j=0; j<alignedStart; ++j)
        tmp0 += cj.pmul(lhs0(j), rhs(j, 0));

      if (alignedSize>alignedStart)
      {
        // process aligned rhs coeffs
        if (lhs0.template aligned<LhsPacket>(alignedStart))
          for (Index j = alignedStart;j<alignedSize;j+=RhsPacketSize)
            ptmp0 = pcj.pmadd(lhs0.template load<LhsPacket, Aligned>(j), rhs.getVectorMapper(j, 0).template load<RhsPacket, Aligned>(0), ptmp0);
        else
          for (Index j = alignedStart;j<alignedSize;j+=RhsPacketSize)
            ptmp0 = pcj.pmadd(lhs0.template load<LhsPacket, Unaligned>(j), rhs.getVectorMapper(j, 0).template load<RhsPacket, Aligned>(0), ptmp0);
        tmp0 += predux(ptmp0);
      }

      // process remaining scalars
      // FIXME this loop get vectorized by the compiler !
      for (Index j=alignedSize; j<depth; ++j)
        tmp0 += cj.pmul(lhs0(j), rhs(j, 0));
      res[i*resIncr] += alpha*tmp0;
    }
    if (skipRows)
    {
      start = 0;
      end = skipRows;
      skipRows = 0;
    }
    else
      break;
  } while(Vectorizable);

  #undef _EIGEN_ACCUMULATE_PACKETS
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0xe8, %rsp
movq %r9, 0x58(%rsp)
movq %r8, 0x50(%rsp)
movq %rcx, %rbx
movq %rdx, %r15
movq %rsi, %r13
movq %rdi, 0x38(%rsp)
movsd %xmm0, 0x40(%rsp)
movsd %xmm1, 0x48(%rsp)
testb $0xf, (%rdx)
je 0xdd40
pushq $0x3
popq %rsi
movq $0x0, 0x8(%rsp)
jmp 0xdd74
xorl %eax, %eax
testb $0xf, (%rbx)
setne %al
movq %rax, %rcx
negq %rcx
testq %r13, %r13
sete %dl
cmpq 0x38(%rsp), %rcx
sete %cl
orb %dl, %al
orb %cl, %al
movzbl %al, %eax
leaq (%rax,%rax,2), %rsi
xorl %ecx, %ecx
testb %al, %al
cmoveq %r13, %rcx
movq %rcx, 0x8(%rsp)
pushq $0x4
popq %rcx
movq 0x38(%rsp), %rax
cqto
idivq %rcx
movq %rax, 0x18(%rsp)
leaq (,%rax,4), %r9
movq 0x8(%rsp), %rcx
movq %rcx, %rax
subq %rsi, %rax
movq %rax, 0xa8(%rsp)
movq %rsi, 0xb0(%rsp)
shll $0x4, %esi
movq %rsi, 0xa0(%rsp)
movq %rcx, %rax
shlq $0x4, %rax
movq %rax, 0xb8(%rsp)
pushq $0x10
popq %rbp
pushq $0x20
popq %r10
pushq $0x30
popq %r11
movq $0x0, 0x10(%rsp)
movapd 0xf670(%rip), %xmm13    # 0x1d450
movapd 0xf647(%rip), %xmm14    # 0x1d430
xorl %r12d, %r12d
movq %rbx, 0xd0(%rsp)
movq %r15, 0xc8(%rsp)
movq %r9, 0xc0(%rsp)
cmpq %r9, %r12
jge 0xe546
xorpd %xmm0, %xmm0
movapd %xmm0, 0x20(%rsp)
movapd %xmm0, 0x90(%rsp)
movapd %xmm0, 0x80(%rsp)
movapd %xmm0, 0x70(%rsp)
movq (%r15), %rdx
movq 0x8(%r15), %rax
movq (%rbx), %rcx
cmpq $0x0, 0x8(%rsp)
movq %r10, 0x68(%rsp)
movq %r11, 0x60(%rsp)
jle 0xdf79
movq 0xb0(%rsp), %rsi
leaq 0xf622(%rip), %rdi      # 0x1d480
movslq (%rdi,%rsi,4), %rsi
addq %rdi, %rsi
jmpq *%rsi
movq %rax, %rsi
imulq 0x10(%rsp), %rsi
addq %rdx, %rsi
movq %rax, %rdi
imulq %rbp, %rdi
addq %rdx, %rdi
movq %rax, %r8
imulq %r10, %r8
addq %rdx, %r8
movq %rax, %r9
imulq %r11, %r9
addq %rdx, %r9
movq 0xa0(%rsp), %r10
movq 0xa8(%rsp), %r11
xorpd %xmm1, %xmm1
xorpd %xmm2, %xmm2
xorpd %xmm3, %xmm3
subq $0x1, %r11
jb 0xe2a8
movapd (%rcx,%r10), %xmm6
movapd (%rsi,%r10), %xmm4
movapd %xmm4, %xmm5
unpcklpd %xmm4, %xmm5            # xmm5 = xmm5[0],xmm4[0]
mulpd %xmm6, %xmm5
xorpd %xmm13, %xmm5
unpckhpd %xmm4, %xmm4            # xmm4 = xmm4[1,1]
movapd (%rdi,%r10), %xmm7
movapd %xmm7, %xmm8
unpcklpd %xmm7, %xmm8            # xmm8 = xmm8[0],xmm7[0]
mulpd %xmm6, %xmm8
movapd (%r8,%r10), %xmm9
movapd %xmm9, %xmm10
unpcklpd %xmm9, %xmm10           # xmm10 = xmm10[0],xmm9[0]
mulpd %xmm6, %xmm10
movapd (%r9,%r10), %xmm11
movapd %xmm11, %xmm12
unpcklpd %xmm11, %xmm12          # xmm12 = xmm12[0],xmm11[0]
mulpd %xmm6, %xmm12
shufpd $0x1, %xmm6, %xmm6      # xmm6 = xmm6[1,0]
mulpd %xmm6, %xmm4
addpd %xmm5, %xmm4
addpd %xmm4, %xmm3
xorpd %xmm13, %xmm8
unpckhpd %xmm7, %xmm7            # xmm7 = xmm7[1,1]
mulpd %xmm6, %xmm7
addpd %xmm8, %xmm7
addpd %xmm7, %xmm2
xorpd %xmm13, %xmm10
unpckhpd %xmm9, %xmm9            # xmm9 = xmm9[1,1]
mulpd %xmm6, %xmm9
addpd %xmm10, %xmm9
addpd %xmm9, %xmm1
xorpd %xmm13, %xmm12
unpckhpd %xmm11, %xmm11          # xmm11 = xmm11[1,1]
mulpd %xmm6, %xmm11
addpd %xmm12, %xmm11
addpd %xmm11, %xmm0
addq $0x10, %r10
jmp 0xdead
xorpd %xmm3, %xmm3
xorpd %xmm2, %xmm2
xorpd %xmm1, %xmm1
jmp 0xe2e4
movq %rax, %rsi
imulq %rbp, %rsi
addq %rdx, %rsi
movq %rax, %rdi
imulq %r10, %rdi
addq %rdx, %rdi
movq %rax, %r8
imulq %r11, %r8
addq %rdx, %r8
movq %rax, %r9
imulq 0x10(%rsp), %r9
addq %rdx, %r9
xorl %r10d, %r10d
movq 0x8(%rsp), %r11
xorpd %xmm1, %xmm1
xorpd %xmm2, %xmm2
xorpd %xmm3, %xmm3
subq $0x1, %r11
jb 0xe2a8
movapd (%rcx,%r10), %xmm6
movapd (%r9,%r10), %xmm4
movapd %xmm4, %xmm5
unpcklpd %xmm4, %xmm5            # xmm5 = xmm5[0],xmm4[0]
mulpd %xmm6, %xmm5
xorpd %xmm13, %xmm5
unpckhpd %xmm4, %xmm4            # xmm4 = xmm4[1,1]
movupd (%rsi,%r10), %xmm7
movapd %xmm7, %xmm8
unpcklpd %xmm7, %xmm8            # xmm8 = xmm8[0],xmm7[0]
mulpd %xmm6, %xmm8
movupd (%rdi,%r10), %xmm9
movapd %xmm9, %xmm10
unpcklpd %xmm9, %xmm10           # xmm10 = xmm10[0],xmm9[0]
mulpd %xmm6, %xmm10
movupd (%r8,%r10), %xmm11
movapd %xmm11, %xmm12
unpcklpd %xmm11, %xmm12          # xmm12 = xmm12[0],xmm11[0]
mulpd %xmm6, %xmm12
shufpd $0x1, %xmm6, %xmm6      # xmm6 = xmm6[1,0]
mulpd %xmm6, %xmm4
addpd %xmm5, %xmm4
addpd %xmm4, %xmm3
xorpd %xmm13, %xmm8
unpckhpd %xmm7, %xmm7            # xmm7 = xmm7[1,1]
mulpd %xmm6, %xmm7
addpd %xmm8, %xmm7
addpd %xmm7, %xmm2
xorpd %xmm13, %xmm10
unpckhpd %xmm9, %xmm9            # xmm9 = xmm9[1,1]
mulpd %xmm6, %xmm9
addpd %xmm10, %xmm9
addpd %xmm9, %xmm1
xorpd %xmm13, %xmm12
unpckhpd %xmm11, %xmm11          # xmm11 = xmm11[1,1]
mulpd %xmm6, %xmm11
addpd %xmm12, %xmm11
addpd %xmm11, %xmm0
addq $0x10, %r10
jmp 0xdfc8
movq %rax, %rsi
imulq 0x10(%rsp), %rsi
addq %rdx, %rsi
movq %rax, %rdi
imulq %rbp, %rdi
addq %rdx, %rdi
movq %rax, %r8
imulq %r10, %r8
addq %rdx, %r8
movq %rax, %r9
imulq %r11, %r9
addq %rdx, %r9
xorl %r10d, %r10d
movq 0x8(%rsp), %r11
xorpd %xmm1, %xmm1
xorpd %xmm2, %xmm2
xorpd %xmm3, %xmm3
subq $0x1, %r11
jb 0xe2a8
movapd (%rcx,%r10), %xmm6
movupd (%rsi,%r10), %xmm4
movapd %xmm4, %xmm5
unpcklpd %xmm4, %xmm5            # xmm5 = xmm5[0],xmm4[0]
mulpd %xmm6, %xmm5
xorpd %xmm13, %xmm5
unpckhpd %xmm4, %xmm4            # xmm4 = xmm4[1,1]
movupd (%rdi,%r10), %xmm7
movapd %xmm7, %xmm8
unpcklpd %xmm7, %xmm8            # xmm8 = xmm8[0],xmm7[0]
mulpd %xmm6, %xmm8
movupd (%r8,%r10), %xmm9
movapd %xmm9, %xmm10
unpcklpd %xmm9, %xmm10           # xmm10 = xmm10[0],xmm9[0]
mulpd %xmm6, %xmm10
movupd (%r9,%r10), %xmm11
movapd %xmm11, %xmm12
unpcklpd %xmm11, %xmm12          # xmm12 = xmm12[0],xmm11[0]
mulpd %xmm6, %xmm12
shufpd $0x1, %xmm6, %xmm6      # xmm6 = xmm6[1,0]
mulpd %xmm6, %xmm4
addpd %xmm5, %xmm4
addpd %xmm4, %xmm3
xorpd %xmm13, %xmm8
unpckhpd %xmm7, %xmm7            # xmm7 = xmm7[1,1]
mulpd %xmm6, %xmm7
addpd %xmm8, %xmm7
addpd %xmm7, %xmm2
xorpd %xmm13, %xmm10
unpckhpd %xmm9, %xmm9            # xmm9 = xmm9[1,1]
mulpd %xmm6, %xmm9
addpd %xmm10, %xmm9
addpd %xmm9, %xmm1
xorpd %xmm13, %xmm12
unpckhpd %xmm11, %xmm11          # xmm11 = xmm11[1,1]
mulpd %xmm6, %xmm11
addpd %xmm12, %xmm11
addpd %xmm11, %xmm0
addq $0x10, %r10
jmp 0xe0d2
movq %rax, %rsi
imulq %rbp, %rsi
addq %rdx, %rsi
movq %rax, %rdi
imulq %r11, %rdi
addq %rdx, %rdi
movq %rax, %r8
imulq 0x10(%rsp), %r8
addq %rdx, %r8
movq %rax, %r9
imulq %r10, %r9
addq %rdx, %r9
xorl %r10d, %r10d
movq 0x8(%rsp), %r11
xorpd %xmm1, %xmm1
xorpd %xmm2, %xmm2
xorpd %xmm3, %xmm3
subq $0x1, %r11
jb 0xe2a8
movapd (%rcx,%r10), %xmm6
movapd (%r8,%r10), %xmm4
movapd %xmm4, %xmm5
unpcklpd %xmm4, %xmm5            # xmm5 = xmm5[0],xmm4[0]
mulpd %xmm6, %xmm5
xorpd %xmm13, %xmm5
unpckhpd %xmm4, %xmm4            # xmm4 = xmm4[1,1]
movupd (%rsi,%r10), %xmm7
movapd %xmm7, %xmm8
unpcklpd %xmm7, %xmm8            # xmm8 = xmm8[0],xmm7[0]
mulpd %xmm6, %xmm8
movapd (%r9,%r10), %xmm9
movapd %xmm9, %xmm10
unpcklpd %xmm9, %xmm10           # xmm10 = xmm10[0],xmm9[0]
mulpd %xmm6, %xmm10
movupd (%rdi,%r10), %xmm11
movapd %xmm11, %xmm12
unpcklpd %xmm11, %xmm12          # xmm12 = xmm12[0],xmm11[0]
mulpd %xmm6, %xmm12
shufpd $0x1, %xmm6, %xmm6      # xmm6 = xmm6[1,0]
mulpd %xmm6, %xmm4
addpd %xmm5, %xmm4
addpd %xmm4, %xmm3
xorpd %xmm13, %xmm8
unpckhpd %xmm7, %xmm7            # xmm7 = xmm7[1,1]
mulpd %xmm6, %xmm7
addpd %xmm8, %xmm7
addpd %xmm7, %xmm2
xorpd %xmm13, %xmm10
unpckhpd %xmm9, %xmm9            # xmm9 = xmm9[1,1]
mulpd %xmm6, %xmm9
addpd %xmm10, %xmm9
addpd %xmm9, %xmm1
xorpd %xmm13, %xmm12
unpckhpd %xmm11, %xmm11          # xmm11 = xmm11[1,1]
mulpd %xmm6, %xmm11
addpd %xmm12, %xmm11
addpd %xmm11, %xmm0
addq $0x10, %r10
jmp 0xe1dc
xorpd %xmm4, %xmm4
addpd %xmm4, %xmm3
movapd %xmm3, 0x20(%rsp)
addpd %xmm4, %xmm2
movapd %xmm2, 0x90(%rsp)
addpd %xmm4, %xmm1
movapd %xmm1, 0x80(%rsp)
addpd %xmm4, %xmm0
movapd %xmm0, 0x70(%rsp)
movq 0x68(%rsp), %r10
movq 0x60(%rsp), %r11
movq %r12, %r15
orq $0x1, %r15
movq %r12, %rbx
orq $0x2, %rbx
movq %r12, %r14
orq $0x3, %r12
movq %rax, %rsi
imulq 0x10(%rsp), %rsi
addq %rdx, %rsi
movq %rax, %rdi
imulq %rbp, %rdi
addq %rdx, %rdi
movq %rax, %r8
imulq %r10, %r8
addq %rdx, %r8
imulq %r11, %rax
addq %rdx, %rax
movq 0xb8(%rsp), %rdx
movq 0x8(%rsp), %r9
cmpq %r13, %r9
jge 0xe418
movupd (%rsi,%rdx), %xmm6
movapd %xmm6, %xmm7
xorpd %xmm14, %xmm7
movsd (%rcx,%rdx), %xmm4
movsd 0x8(%rcx,%rdx), %xmm5
unpcklpd %xmm5, %xmm5            # xmm5 = xmm5[0,0]
movapd %xmm6, %xmm8
shufpd $0x1, %xmm7, %xmm8      # xmm8 = xmm8[1],xmm7[0]
mulpd %xmm5, %xmm8
unpcklpd %xmm4, %xmm4            # xmm4 = xmm4[0,0]
mulpd %xmm4, %xmm6
addpd %xmm8, %xmm6
addpd %xmm6, %xmm3
movapd %xmm3, 0x20(%rsp)
movupd (%rdi,%rdx), %xmm6
movapd %xmm6, %xmm7
xorpd %xmm14, %xmm7
movapd %xmm6, %xmm8
shufpd $0x1, %xmm7, %xmm8      # xmm8 = xmm8[1],xmm7[0]
mulpd %xmm5, %xmm8
mulpd %xmm4, %xmm6
addpd %xmm8, %xmm6
addpd %xmm6, %xmm2
movapd %xmm2, 0x90(%rsp)
movupd (%r8,%rdx), %xmm6
movapd %xmm6, %xmm7
xorpd %xmm14, %xmm7
movapd %xmm6, %xmm8
shufpd $0x1, %xmm7, %xmm8      # xmm8 = xmm8[1],xmm7[0]
mulpd %xmm5, %xmm8
mulpd %xmm4, %xmm6
addpd %xmm8, %xmm6
addpd %xmm6, %xmm1
movapd %xmm1, 0x80(%rsp)
movupd (%rax,%rdx), %xmm6
movapd %xmm6, %xmm7
mulpd %xmm6, %xmm4
xorpd %xmm14, %xmm6
shufpd $0x1, %xmm6, %xmm7      # xmm7 = xmm7[1],xmm6[0]
mulpd %xmm5, %xmm7
addpd %xmm7, %xmm4
addpd %xmm4, %xmm0
movapd %xmm0, 0x70(%rsp)
incq %r9
addq $0x10, %rdx
jmp 0xe32d
leaq 0x40(%rsp), %rdi
leaq 0x20(%rsp), %rsi
callq 0xab66
movq %r14, 0xd8(%rsp)
movq %r14, %rax
movq %rbp, 0xe0(%rsp)
movq 0x58(%rsp), %rbp
imulq %rbp, %rax
shlq $0x4, %rax
movq 0x50(%rsp), %r14
movupd (%r14,%rax), %xmm2
unpcklpd %xmm1, %xmm0            # xmm0 = xmm0[0],xmm1[0]
addpd %xmm2, %xmm0
movupd %xmm0, (%r14,%rax)
leaq 0x40(%rsp), %rdi
leaq 0x90(%rsp), %rsi
callq 0xab66
imulq %rbp, %r15
shlq $0x4, %r15
movupd (%r14,%r15), %xmm2
unpcklpd %xmm1, %xmm0            # xmm0 = xmm0[0],xmm1[0]
addpd %xmm2, %xmm0
movupd %xmm0, (%r14,%r15)
leaq 0x40(%rsp), %rdi
leaq 0x80(%rsp), %rsi
callq 0xab66
imulq %rbp, %rbx
shlq $0x4, %rbx
movupd (%r14,%rbx), %xmm2
unpcklpd %xmm1, %xmm0            # xmm0 = xmm0[0],xmm1[0]
addpd %xmm2, %xmm0
movupd %xmm0, (%r14,%rbx)
leaq 0x40(%rsp), %rdi
leaq 0x70(%rsp), %rsi
callq 0xab66
movapd 0xef5c(%rip), %xmm14    # 0x1d430
movq 0x60(%rsp), %r11
movq 0x68(%rsp), %r10
movapd 0xef69(%rip), %xmm13    # 0x1d450
imulq %rbp, %r12
movq 0xe0(%rsp), %rbp
shlq $0x4, %r12
movupd (%r14,%r12), %xmm2
unpcklpd %xmm1, %xmm0            # xmm0 = xmm0[0],xmm1[0]
addpd %xmm2, %xmm0
movupd %xmm0, (%r14,%r12)
movq 0xd8(%rsp), %r12
addq $0x4, %r12
addq $0x40, 0x10(%rsp)
addq $0x40, %rbp
addq $0x40, %r10
addq $0x40, %r11
movq 0xd0(%rsp), %rbx
movq 0xc8(%rsp), %r15
movq 0xc0(%rsp), %r9
jmp 0xde04
shlq $0x6, 0x18(%rsp)
movq 0x8(%rsp), %r12
shlq $0x4, %r12
movapd 0xeed3(%rip), %xmm5     # 0x1d430
leaq 0x40(%rsp), %r14
cmpq 0x38(%rsp), %r9
jge 0xe722
xorpd %xmm0, %xmm0
movapd %xmm0, 0x20(%rsp)
movq (%r15), %rax
movq 0x8(%r15), %rdx
cmpq $0x0, 0x8(%rsp)
jle 0xe5fa
movq %rdx, %rsi
imulq %r9, %rsi
shlq $0x4, %rsi
addq %rax, %rsi
movapd 0x20(%rsp), %xmm0
movq (%rbx), %rcx
testb $0xf, %sil
je 0xe603
movq %rdx, %rsi
imulq 0x18(%rsp), %rsi
addq %rax, %rsi
xorl %edi, %edi
movq 0x8(%rsp), %r8
movapd %xmm0, %xmm1
subq $0x1, %r8
jb 0xe656
movapd (%rcx,%rdi), %xmm2
movupd (%rsi,%rdi), %xmm3
movapd %xmm3, %xmm4
unpcklpd %xmm3, %xmm4            # xmm4 = xmm4[0],xmm3[0]
mulpd %xmm2, %xmm4
xorpd %xmm13, %xmm4
unpckhpd %xmm3, %xmm3            # xmm3 = xmm3[1,1]
shufpd $0x1, %xmm2, %xmm2      # xmm2 = xmm2[1,0]
mulpd %xmm3, %xmm2
addpd %xmm4, %xmm2
addpd %xmm2, %xmm1
addq $0x10, %rdi
jmp 0xe5ba
movq (%rbx), %rcx
xorpd %xmm2, %xmm2
jmp 0xe668
movq %rdx, %rsi
imulq 0x18(%rsp), %rsi
addq %rax, %rsi
xorl %edi, %edi
movq 0x8(%rsp), %r8
movapd %xmm0, %xmm1
subq $0x1, %r8
jb 0xe656
movapd (%rcx,%rdi), %xmm2
movapd (%rsi,%rdi), %xmm3
movapd %xmm3, %xmm4
unpcklpd %xmm3, %xmm4            # xmm4 = xmm4[0],xmm3[0]
mulpd %xmm2, %xmm4
xorpd %xmm13, %xmm4
unpckhpd %xmm3, %xmm3            # xmm3 = xmm3[1,1]
shufpd $0x1, %xmm2, %xmm2      # xmm2 = xmm2[1,0]
mulpd %xmm3, %xmm2
addpd %xmm4, %xmm2
addpd %xmm2, %xmm1
addq $0x10, %rdi
jmp 0xe61a
movq %xmm0, %xmm2            # xmm2 = xmm0[0],zero
addpd %xmm1, %xmm2
movapd %xmm2, 0x20(%rsp)
unpckhpd %xmm2, %xmm2            # xmm2 = xmm2[1,1]
movapd 0x20(%rsp), %xmm0
unpcklpd %xmm2, %xmm0            # xmm0 = xmm0[0],xmm2[0]
imulq 0x18(%rsp), %rdx
addq %rdx, %rax
movq %r12, %rdx
movq 0x8(%rsp), %rsi
cmpq %r13, %rsi
jge 0xe6cc
movupd (%rax,%rdx), %xmm1
movsd (%rcx,%rdx), %xmm2
movapd %xmm1, %xmm3
unpcklpd %xmm2, %xmm2            # xmm2 = xmm2[0,0]
mulpd %xmm1, %xmm2
xorpd %xmm5, %xmm1
movsd 0x8(%rcx,%rdx), %xmm4
unpcklpd %xmm4, %xmm4            # xmm4 = xmm4[0,0]
shufpd $0x1, %xmm1, %xmm3      # xmm3 = xmm3[1],xmm1[0]
mulpd %xmm4, %xmm3
addpd %xmm3, %xmm2
addpd %xmm2, %xmm0
movapd %xmm0, 0x20(%rsp)
incq %rsi
addq $0x10, %rdx
jmp 0xe683
movq %r14, %rdi
leaq 0x20(%rsp), %rsi
movq %r9, %rbp
callq 0xab66
movapd 0xed4c(%rip), %xmm5     # 0x1d430
movapd 0xed63(%rip), %xmm13    # 0x1d450
movq %rbp, %r9
movq %rbp, %rax
imulq 0x58(%rsp), %rax
shlq $0x4, %rax
movq 0x50(%rsp), %rcx
movupd (%rcx,%rax), %xmm2
unpcklpd %xmm1, %xmm0            # xmm0 = xmm0[0],xmm1[0]
addpd %xmm2, %xmm0
movupd %xmm0, (%rcx,%rax)
incq %r9
addq $0x10, 0x18(%rsp)
jmp 0xe562
addq $0xe8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq