|
|
Ublas : |
From: Paul C. Leopardi (leopardi_at_[hidden])
Date: 2005-07-10 18:49:06
Hi all, comments below. Best regards, Paul Leopardi
On Sun, 10 Jul 2005 05:38 pm, Gunter Winkler wrote:
> Zitat von Georg Baum <georg.baum_at_[hidden]>:
> > Hello,
> >
> > while playing a bit with the profiler I noticed that matrix_norm_inf
> > computes the norm by iterating over each element (i, j). The same is true
> > for matrix_norm_1 and matrix_norm_frobenius. This is of course
> > inefficient for sparse matrices.
> > I wonder why that is the case? I can't believe that such an inefficient
> > implementation is used without reason. Now I am wondering if I have to
> > expect more surprises like that? Or did I miss something?
>
> Maybe you are the first one who uses these norms for non-dense matrices.
GluCat uses these norms, notably in operator/= for matrix_multi<>. See
matrix_multi_imp.h