Method MultTransAdd

linalg/special_matrix.cpp:76–83

  }
  
  void PermutationMatrix :: MultTransAdd (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(ind[i]) += s * fvx(i);
  }