Re: [ublas] ublas double precision issues

Subject: Re: [ublas] ublas double precision issues
From: Umut Tabak (u.tabak_at_[hidden])
Date: 2011-06-13 11:47:16

On 06/13/2011 02:11 PM, Umut Tabak wrote:
> Dear all,
>
> I have implemented a lanczos solver in ublas, however there is a 
> subtle question from the point of view of accuracy related to floating 
> point operations.
>
As a follow up to my question, the algorithm does some 
orthogonalizations, using the inner_prod function, this is where the 
accuracy is being lost. Could you comment on this also?

-- 
If I have a thousand ideas and only one turns out to be good,
I am satisfied.
Alfred Nobel