Boost logo

Ublas :

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