Method MultAdd

linalg/special_matrix.cpp:67–74

  }

  void PermutationMatrix :: MultAdd (double s, const BaseVector & x, BaseVector & y) const
  {
    auto fvx = x.FV<double>();
    auto fvy = y.FV<double>();
    for (size_t i = 0; i < ind.Size(); i++)
      if(ind[i] != size_t(-1))
        fvy(i) += s * fvx(ind[i]);
  }
  
  void PermutationMatrix :: MultTransAdd (double s, const BaseVector & x, BaseVector & y) const
  {