Re: [ublas] Solving linear equations in a loop (Optimization issues)

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

Sourabh schrieb:
> On Fri, 2 Mar 2007, Gunter Winkler wrote:
>   
> The program (see below), after some modifications is still taking around 
> 8 minutes for the arguments 1 200000 100. 
>
> Even if I remove solve, and addition methods, the statement
>  lembda1.assign (ident-lembda) is causing the program to take 1 minute.
>   
1) You can avoid the assign if you fill lembda with the negative values, 
such that you have to solve with (I+lemda).
2) you should try the current CVS head of boost.
3) use inplace_solve_lower<true>(A, x, row_major_tag()); from my examples.

mfg
Gunter