Method MultTransAdd

parallel/parallel_matrices.cpp:580–597

  }

  void ParallelMatrix :: MultTransAdd (double s, const BaseVector & x, BaseVector & y) const
  {
    // const auto & xpar = dynamic_cast_ParallelBaseVector(x);
    // auto & ypar = dynamic_cast_ParallelBaseVector(y);
    // if (op & char(1))
    if (ColType(op)==CUMULATED)
      x.Distribute();
    else
      x.Cumulate();
    
    //if (op & char(2))
    if (RowType(op)==CUMULATED)
      y.Distribute();
    else
      y.Cumulate();
    // mat->MultTransAdd (s, *xpar.GetLocalVector(), *ypar.GetLocalVector());
    mat->MultTransAdd (s, *x.GetLocalVector(), *y.GetLocalVector());
  }

  shared_ptr<BaseMatrix> ParallelMatrix :: CreateMatrix () const
  {