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

From: Sourabh (sourabh_at_[hidden])
Date: 2007-03-05 03:37:56

Thanks Gunter. 
Your inputs are proving to be of great help. 
Can I include your trisolve.hpp in my program ? 
Are there some licence issues ? Please let me know. 

On Sat, 3 Mar 2007, Gunter Winkler wrote:

> Am Samstag, 3. März 2007 14:39 schrieb Sourabh:
> > When I run this with arguments 200000 100 , the program takes 1.5
> > seconds on my machine. If I change solve to
> > inplace_solve_lower<true>(lembda1, omegas, row_major_tag());
> > it takes 35 seconds.
> 
> Did you define NDEBUG and use optimizations?
> 
> $ g++ -I $HOME/include -o sousolve sousolve.cpp -DNDEBUG -O2
> $ ./sousolve 200000 100
> 1.65 (solve)
> 0.99 (inplace solve)
> 
> mfg
> Gunter
> 
> PS: Did you read the documentation included in trisolve.hpp ?
> 

-- 
-- Sourabh