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

From: Gunter Winkler (guwi17_at_[hidden])
Date: 2007-03-03 10:20:22

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 ?