Exemplo n.º 1
0
void SparseMatrixTemplate_RM<T>::eraseZeros(T zeroTol)
{
  for(size_t i=0;i<rows.size();i++) {
    RowT temp;
    temp.n = n;
    typename RowT::const_iterator j;
    for(j=rows[i].begin();j!=rows[i].end();j++)
      if(Abs(j->second) > Abs(zeroTol)) temp.insert(j->first,j->second);
    rows[i].entries.clear();
    rows[i] = temp;
  }
}