Boost logo

Ublas :

From: Gunter Winkler (guwi17_at_[hidden])
Date: 2006-11-09 03:38:51

On Wednesday 08 November 2006 21:16, Fred wrote:
> Hi,
> 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
this helps.

btw. If you are interested in sparse LU decomposition you should have a look