Method MultTransAdd

linalg/sparsematrix_impl.hpp:345–375

  template <class TM, class TV_ROW, class TV_COL>
  void SparseMatrix<TM,TV_ROW,TV_COL> ::
  MultTransAdd (double s, const BaseVector & x, BaseVector & y) const
  {
    static Timer timer ("SparseMatrix::MultTransAdd");
    RegionTimer reg (timer);

    FlatVector<TVY> fx = x.FV<TVY>();
    FlatVector<TVX> fy = y.FV<TVX>();
    
    for (int i = 0; i < this->Height(); i++)
      AddRowTransToVector (i, s*fx(i), fy);

    timer.AddFlops (this->NZE());


    /*
    static Timer t("SparseMatrix::MultTransAdd"); RegionTimer reg(t);
    t.AddFlops (this->NZE()*sizeof(TV_ROW)*sizeof(TV_COL)/sqr(sizeof(double)));

    ParallelForRange
      (balance, [&] (IntRange myrange)
       {
         FlatVector<TVY> fx = x.FV<TVY>(); 
         FlatVector<TVX> fy = y.FV<TVX>(); 

         for (auto i : myrange)
           AddRowTransToVectorAtomic (i, s*fx(i), fy);
       });
    */
    
  }


  template <class TM, class TV_ROW, class TV_COL>