// File: /workspace/github_cpp_makefile_star2_O3/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 $0x108, %rsp            # imm = 0x108
movsd %xmm1, 0x48(%rsp)
movsd %xmm0, 0x40(%rsp)
movq %r9, 0x78(%rsp)
movq %r8, 0x70(%rsp)
movq %rcx, 0x30(%rsp)
movq %rdx, %r12
movq %rsi, %r13
movq %rdi, %rbp
testb $0xf, (%rdx)
je 0x10c29
movl $0x3, %eax
movq %rax, 0x60(%rsp)
movq $0x0, 0x8(%rsp)
jmp 0x10c65
xorl %eax, %eax
movq 0x30(%rsp), %rcx
testb $0xf, (%rcx)
setne %al
movq %rax, %rcx
negq %rcx
testq %r13, %r13
sete %dl
cmpq %rbp, %rcx
sete %cl
orb %dl, %al
orb %cl, %al
movzbl %al, %eax
leaq (%rax,%rax,2), %rcx
movq %rcx, 0x60(%rsp)
xorl %ecx, %ecx
testb %al, %al
cmoveq %r13, %rcx
movq %rcx, 0x8(%rsp)
leaq 0x3(%rbp), %rax
testq %rbp, %rbp
cmovnsq %rbp, %rax
movq %rax, 0xb8(%rsp)
andq $-0x4, %rax
movq %rax, 0x38(%rsp)
cmpq $0x4, %rbp
jl 0x11490
movq 0x8(%rsp), %rax
movq %rax, %rdx
movq 0x60(%rsp), %rcx
subq %rcx, %rdx
movq %rdx, 0xc8(%rsp)
shll $0x4, %ecx
movq %rcx, 0xc0(%rsp)
movq %r13, %rcx
subq %rax, %rcx
movq %rcx, 0xd8(%rsp)
shlq $0x4, %rax
movq %rax, 0xd0(%rsp)
movl $0x30, %ebx
movl $0x20, %r15d
movl $0x10, %eax
movq %rax, 0x28(%rsp)
movq $0x0, 0x20(%rsp)
movapd 0xc781(%rip), %xmm13    # 0x1d470
movapd 0xc758(%rip), %xmm14    # 0x1d450
movq $0x0, 0x68(%rsp)
movq %r12, 0xf0(%rsp)
movq %r13, 0xe8(%rsp)
movq %rbp, 0xe0(%rsp)
xorpd %xmm0, %xmm0
movapd %xmm0, 0x50(%rsp)
movapd %xmm0, 0xa0(%rsp)
movapd %xmm0, 0x90(%rsp)
movapd %xmm0, 0x80(%rsp)
movq (%r12), %rcx
movq 0x8(%r12), %rax
xorpd %xmm1, %xmm1
xorpd %xmm2, %xmm2
xorpd %xmm3, %xmm3
cmpq $0x0, 0x8(%rsp)
jle 0x111b0
movq 0x30(%rsp), %rdx
movq (%rdx), %rdx
movq 0x60(%rsp), %rsi
leaq 0xc75d(%rip), %rdi      # 0x1d4d0
movslq (%rdi,%rsi,4), %rsi
addq %rdi, %rsi
jmpq *%rsi
movq %rax, %rsi
imulq %rbx, %rsi
addq %rcx, %rsi
movq %rax, %rdi
imulq %r15, %rdi
addq %rcx, %rdi
movq %rax, %r8
imulq 0x28(%rsp), %r8
addq %rcx, %r8
movq %rax, %r9
imulq 0x20(%rsp), %r9
addq %rcx, %r9
movq 0xc0(%rsp), %r10
movq 0xc8(%rsp), %r11
movapd (%rdx,%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]
movapd (%r8,%r10), %xmm7
movapd %xmm7, %xmm8
unpcklpd %xmm7, %xmm8            # xmm8 = xmm8[0],xmm7[0]
mulpd %xmm6, %xmm8
movapd (%rdi,%r10), %xmm9
movapd %xmm9, %xmm10
unpcklpd %xmm9, %xmm10           # xmm10 = xmm10[0],xmm9[0]
mulpd %xmm6, %xmm10
movapd (%rsi,%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, %xmm0
xorpd %xmm13, %xmm8
unpckhpd %xmm7, %xmm7            # xmm7 = xmm7[1,1]
mulpd %xmm6, %xmm7
addpd %xmm8, %xmm7
addpd %xmm7, %xmm1
xorpd %xmm13, %xmm10
unpckhpd %xmm9, %xmm9            # xmm9 = xmm9[1,1]
mulpd %xmm6, %xmm9
addpd %xmm10, %xmm9
addpd %xmm9, %xmm2
xorpd %xmm13, %xmm12
unpckhpd %xmm11, %xmm11          # xmm11 = xmm11[1,1]
mulpd %xmm6, %xmm11
addpd %xmm12, %xmm11
addpd %xmm11, %xmm3
addq $0x10, %r10
decq %r11
jne 0x10db8
jmp 0x1117b
movq %rax, %rsi
imulq 0x20(%rsp), %rsi
addq %rcx, %rsi
movq %rax, %rdi
imulq %rbx, %rdi
addq %rcx, %rdi
movq %rax, %r8
imulq %r15, %r8
addq %rcx, %r8
movq %rax, %r9
imulq 0x28(%rsp), %r9
addq %rcx, %r9
xorl %r10d, %r10d
movq 0x8(%rsp), %r11
movapd (%rdx,%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]
movupd (%r9,%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 (%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, %xmm0
xorpd %xmm13, %xmm8
unpckhpd %xmm7, %xmm7            # xmm7 = xmm7[1,1]
mulpd %xmm6, %xmm7
addpd %xmm8, %xmm7
addpd %xmm7, %xmm1
xorpd %xmm13, %xmm10
unpckhpd %xmm9, %xmm9            # xmm9 = xmm9[1,1]
mulpd %xmm6, %xmm9
addpd %xmm10, %xmm9
addpd %xmm9, %xmm2
xorpd %xmm13, %xmm12
unpckhpd %xmm11, %xmm11          # xmm11 = xmm11[1,1]
mulpd %xmm6, %xmm11
addpd %xmm12, %xmm11
addpd %xmm11, %xmm3
addq $0x10, %r10
decq %r11
jne 0x10eb7
jmp 0x1117b
movq %rax, %rsi
imulq %rbx, %rsi
addq %rcx, %rsi
movq %rax, %rdi
imulq %r15, %rdi
addq %rcx, %rdi
movq %rax, %r8
imulq 0x28(%rsp), %r8
addq %rcx, %r8
movq %rax, %r9
imulq 0x20(%rsp), %r9
addq %rcx, %r9
xorl %r10d, %r10d
movq 0x8(%rsp), %r11
movapd (%rdx,%r10), %xmm6
movupd (%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 (%r8,%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 (%rsi,%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, %xmm0
xorpd %xmm13, %xmm8
unpckhpd %xmm7, %xmm7            # xmm7 = xmm7[1,1]
mulpd %xmm6, %xmm7
addpd %xmm8, %xmm7
addpd %xmm7, %xmm1
xorpd %xmm13, %xmm10
unpckhpd %xmm9, %xmm9            # xmm9 = xmm9[1,1]
mulpd %xmm6, %xmm9
addpd %xmm10, %xmm9
addpd %xmm9, %xmm2
xorpd %xmm13, %xmm12
unpckhpd %xmm11, %xmm11          # xmm11 = xmm11[1,1]
mulpd %xmm6, %xmm11
addpd %xmm12, %xmm11
addpd %xmm11, %xmm3
addq $0x10, %r10
decq %r11
jne 0x10fb6
jmp 0x1117b
movq %rax, %rsi
imulq %r15, %rsi
addq %rcx, %rsi
movq %rax, %rdi
imulq 0x20(%rsp), %rdi
addq %rcx, %rdi
movq %rax, %r8
imulq %rbx, %r8
addq %rcx, %r8
movq %rax, %r9
imulq 0x28(%rsp), %r9
addq %rcx, %r9
xorl %r10d, %r10d
movq 0x8(%rsp), %r11
movapd (%rdx,%r10), %xmm6
movapd (%rdi,%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 (%r9,%r10), %xmm7
movapd %xmm7, %xmm8
unpcklpd %xmm7, %xmm8            # xmm8 = xmm8[0],xmm7[0]
mulpd %xmm6, %xmm8
movapd (%rsi,%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, %xmm0
xorpd %xmm13, %xmm8
unpckhpd %xmm7, %xmm7            # xmm7 = xmm7[1,1]
mulpd %xmm6, %xmm7
addpd %xmm8, %xmm7
addpd %xmm7, %xmm1
xorpd %xmm13, %xmm10
unpckhpd %xmm9, %xmm9            # xmm9 = xmm9[1,1]
mulpd %xmm6, %xmm9
addpd %xmm10, %xmm9
addpd %xmm9, %xmm2
xorpd %xmm13, %xmm12
unpckhpd %xmm11, %xmm11          # xmm11 = xmm11[1,1]
mulpd %xmm6, %xmm11
addpd %xmm12, %xmm11
addpd %xmm11, %xmm3
addq $0x10, %r10
decq %r11
jne 0x110b5
xorpd %xmm4, %xmm4
addpd %xmm4, %xmm0
movapd %xmm0, 0x50(%rsp)
addpd %xmm4, %xmm1
movapd %xmm1, 0xa0(%rsp)
addpd %xmm4, %xmm2
movapd %xmm2, 0x90(%rsp)
addpd %xmm4, %xmm3
movapd %xmm3, 0x80(%rsp)
cmpq 0x8(%rsp), %r13
jle 0x112e2
movq 0x30(%rsp), %rdx
movq (%rdx), %rdx
movq %rax, %rsi
imulq %rbx, %rsi
addq %rcx, %rsi
movq %rax, %rdi
imulq %r15, %rdi
addq %rcx, %rdi
movq %rax, %r8
imulq 0x28(%rsp), %r8
addq %rcx, %r8
imulq 0x20(%rsp), %rax
addq %rcx, %rax
movq 0xd0(%rsp), %rcx
movq 0xd8(%rsp), %r9
movupd (%rax,%rcx), %xmm6
movapd %xmm6, %xmm7
xorpd %xmm14, %xmm7
movsd (%rdx,%rcx), %xmm4
movsd 0x8(%rdx,%rcx), %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, %xmm0
movapd %xmm0, 0x50(%rsp)
movupd (%r8,%rcx), %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, 0xa0(%rsp)
movupd (%rdi,%rcx), %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 (%rsi,%rcx), %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, %xmm3
movapd %xmm3, 0x80(%rsp)
addq $0x10, %rcx
decq %r9
jne 0x111fc
movq %r15, 0xf8(%rsp)
movq 0x68(%rsp), %r14
movq %r14, %r15
orq $0x1, %r15
movq %rbx, 0x100(%rsp)
movq %r14, %rbx
orq $0x2, %rbx
movq %r14, %r12
orq $0x3, %r12
movsd 0x40(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x48(%rsp), %xmm0
movsd %xmm0, 0x18(%rsp)
leaq 0x10(%rsp), %rdi
leaq 0x50(%rsp), %rsi
callq 0xff34
movq %r14, %rax
movq 0x78(%rsp), %rbp
imulq %rbp, %rax
shlq $0x4, %rax
movq 0x70(%rsp), %r13
movupd (%r13,%rax), %xmm0
addpd 0x10(%rsp), %xmm0
movupd %xmm0, (%r13,%rax)
movsd 0x40(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x48(%rsp), %xmm0
movsd %xmm0, 0x18(%rsp)
leaq 0x10(%rsp), %rdi
leaq 0xa0(%rsp), %rsi
callq 0xff34
imulq %rbp, %r15
shlq $0x4, %r15
movupd (%r13,%r15), %xmm0
addpd 0x10(%rsp), %xmm0
movupd %xmm0, (%r13,%r15)
movq 0xf8(%rsp), %r15
movsd 0x40(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x48(%rsp), %xmm0
movsd %xmm0, 0x18(%rsp)
leaq 0x10(%rsp), %rdi
leaq 0x90(%rsp), %rsi
callq 0xff34
imulq %rbp, %rbx
shlq $0x4, %rbx
movupd (%r13,%rbx), %xmm0
addpd 0x10(%rsp), %xmm0
movupd %xmm0, (%r13,%rbx)
movq 0x100(%rsp), %rbx
movsd 0x40(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x48(%rsp), %xmm0
movsd %xmm0, 0x18(%rsp)
leaq 0x10(%rsp), %rdi
leaq 0x80(%rsp), %rsi
callq 0xff34
movapd 0xc025(%rip), %xmm14    # 0x1d450
imulq %rbp, %r12
shlq $0x4, %r12
movupd (%r13,%r12), %xmm0
addpd 0x10(%rsp), %xmm0
movupd %xmm0, (%r13,%r12)
addq $0x4, %r14
addq $0x40, %rbx
addq $0x40, %r15
addq $0x40, 0x28(%rsp)
addq $0x40, 0x20(%rsp)
movq %r14, 0x68(%rsp)
cmpq 0x38(%rsp), %r14
movq 0xf0(%rsp), %r12
movq 0xe8(%rsp), %r13
movq 0xe0(%rsp), %rbp
movapd 0xbfe6(%rip), %xmm13    # 0x1d470
jl 0x10d19
cmpq %rbp, 0x38(%rsp)
jge 0x1169a
movq 0xb8(%rsp), %rbx
sarq $0x2, %rbx
shlq $0x6, %rbx
movq %r13, %r14
movq 0x8(%rsp), %r15
subq %r15, %r14
shlq $0x4, %r15
movapd 0xbfae(%rip), %xmm5     # 0x1d470
movapd 0xbf86(%rip), %xmm6     # 0x1d450
xorpd %xmm0, %xmm0
movapd %xmm0, 0x50(%rsp)
movq (%r12), %rax
movq 0x8(%r12), %rcx
cmpq $0x0, 0x8(%rsp)
jle 0x11558
movq %rcx, %rsi
imulq 0x38(%rsp), %rsi
shlq $0x4, %rsi
addq %rax, %rsi
movapd 0x50(%rsp), %xmm0
movq 0x30(%rsp), %rdx
movq (%rdx), %rdx
testb $0xf, %sil
je 0x1155e
movq %rcx, %rsi
imulq %rbx, %rsi
addq %rax, %rsi
xorl %edi, %edi
movq 0x8(%rsp), %r8
movapd %xmm0, %xmm1
movapd (%rdx,%rdi), %xmm2
movupd (%rsi,%rdi), %xmm3
movapd %xmm3, %xmm4
unpcklpd %xmm3, %xmm4            # xmm4 = xmm4[0],xmm3[0]
mulpd %xmm2, %xmm4
xorpd %xmm5, %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
decq %r8
jne 0x1151e
jmp 0x115ab
xorpd %xmm2, %xmm2
jmp 0x115bd
movq %rcx, %rsi
imulq %rbx, %rsi
addq %rax, %rsi
xorl %edi, %edi
movq 0x8(%rsp), %r8
movapd %xmm0, %xmm1
movapd (%rdx,%rdi), %xmm2
movapd (%rsi,%rdi), %xmm3
movapd %xmm3, %xmm4
unpcklpd %xmm3, %xmm4            # xmm4 = xmm4[0],xmm3[0]
mulpd %xmm2, %xmm4
xorpd %xmm5, %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
decq %r8
jne 0x11573
movq %xmm0, %xmm2            # xmm2 = xmm0[0],zero
addpd %xmm1, %xmm2
movapd %xmm2, 0x50(%rsp)
unpckhpd %xmm2, %xmm2            # xmm2 = xmm2[1,1]
cmpq 0x8(%rsp), %r13
jle 0x11627
movapd 0x50(%rsp), %xmm0
unpcklpd %xmm2, %xmm0            # xmm0 = xmm0[0],xmm2[0]
movq 0x30(%rsp), %rdx
movq (%rdx), %rdx
imulq %rbx, %rcx
addq %rcx, %rax
movq %r15, %rcx
movq %r14, %rsi
movupd (%rax,%rcx), %xmm1
movsd (%rdx,%rcx), %xmm2
movapd %xmm1, %xmm3
unpcklpd %xmm2, %xmm2            # xmm2 = xmm2[0,0]
mulpd %xmm1, %xmm2
xorpd %xmm6, %xmm1
movsd 0x8(%rdx,%rcx), %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, 0x50(%rsp)
addq $0x10, %rcx
decq %rsi
jne 0x115e3
movsd 0x40(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x48(%rsp), %xmm0
movsd %xmm0, 0x18(%rsp)
leaq 0x10(%rsp), %rdi
leaq 0x50(%rsp), %rsi
callq 0xff34
movapd 0xbdfa(%rip), %xmm6     # 0x1d450
movq 0x38(%rsp), %rdx
movq %rdx, %rax
imulq 0x78(%rsp), %rax
shlq $0x4, %rax
movq 0x70(%rsp), %rcx
movupd (%rcx,%rax), %xmm0
addpd 0x10(%rsp), %xmm0
movupd %xmm0, (%rcx,%rax)
incq %rdx
addq $0x10, %rbx
movq %rdx, 0x38(%rsp)
cmpq %rbp, %rdx
movapd 0xbddc(%rip), %xmm5     # 0x1d470
jne 0x114ca
addq $0x108, %rsp            # imm = 0x108
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq