Re: [ublas] Iteration over non zero elements

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 
at
http://www.cise.ufl.edu/research/sparse/SIAM06/

mfg
Gunter