From: Gunter Winkler (guwi17_at_[hidden])
Date: 2006-11-09 03:38:51
On Wednesday 08 November 2006 21:16, Fred wrote:
> I'm building one simple version of Gauss Elimination, but I would like to
> iterate only over non zero elements.
This happens by design. If you want to iterate over non_zeros you have to use
const_iterator, the mutable iterator runs over all elements because it has to
give access to all assignable elements.
Did you have a look at my cholesky and incomplete cholesky examples? Maybe
btw. If you are interested in sparse LU decomposition you should have a look