sgesv in 1.1 is slow...

[dan]

[jpeinado]

[cjest]

Here comes the result out of my machine, working on complex precision

>> A = rand(2048,2048,'single')+1i*rand(2048,2048,'single');
>> B = rand(2048,64,'single')+1i*rand(2048,64,'single');
>> tic; A\B; toc
Elapsed time is 1.594184 seconds.
>> tic; x = culasv(A,B ); toc;
Elapsed time is 0.338316 seconds.
>> B = rand(2048,2048,'single')+1i*rand(2048,2048,'single');
>> tic; A\B; toc;
Elapsed time is 2.029847 seconds.
>> tic; x = culasv(A,B ); toc;
Elapsed time is 0.770844 seconds.
>> B = rand(2048,5000,'single')+1i*rand(2048,5000,'single');
>> tic; A\B; toc;
Elapsed time is 3.651626 seconds.
>> tic; x = culasv(A,B ); toc;
Elapsed time is 1.451895 seconds.

Note: culasv is "Cula Solver" using "Cgesv".
Speedup is obtained even where size of B >= size of A.


BR/CJ

[jpeinado]

[Boxed Cylon]

[jpeinado]

[john]

[Boxed Cylon]

[jpeinado]

[cjest]

[jpeinado]

[jpeinado]

[john]

Should be very soon. We have expanded the scope of work to have some far-reaching speedups across many of the CULA routines including gels, getrf, posv, and more.

[jpeinado]

[john]

1.3 is released. I'm anxious to see the new results from this thread!