Re: [ublas] lu factorize cannot be compiled on sparse complex matrix?

From: Gunter Winkler (guwi17_at_[hidden])
Date: 2007-09-10 03:07:09

pinacle wrote:
> Hi Gunter
>
> In the reply to Karl, I proposed the solution as changing the LINE 147 of
> lu.cpp from
>   project (mci, range (i + 1, size1)) *= value_type (1) / m (i, i); 
>
> to
>   project (mci, range (i + 1, size1)) *= value_type (1);
>   project (mci, range (i + 1, size1)) /= m (i, i);
>
> At least this works for me now.
>   
Ok. In this case you can simply drop the multiplication by 1. Does my
proposol work too? There are still some processors available where
division is slower than multiplication.

mfg
Gunter