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

From: pinacle (pinacle2000_at_[hidden])
Date: 2007-09-10 18:30:25

Your way to change "value_type (1) / m (i, i)" to "value_type (1) /
value_type ( m (i, i) )" also works.

Best,
Feng

-----Original Message-----
From: ublas-bounces_at_[hidden] [mailto:ublas-bounces_at_[hidden]]
On Behalf Of Gunter Winkler
Sent: Monday, September 10, 2007 12:07 AM
To: ublas mailing list
Subject: Re: [ublas] lu factorize cannot be compiled on sparse complex
matrix?

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

_______________________________________________
ublas mailing list
ublas_at_[hidden]
http://lists.boost.org/mailman/listinfo.cgi/ublas