slasrt2/dlasrt2 always use insertion sort

Post here if you want to report a bug to the LAPACK team

slasrt2/dlasrt2 always use insertion sort

Postby quenteridean » Tue May 03, 2016 2:48 pm

slasrt2/dlasrt2
line 110 as referenced in
http://www.netlib.org/scalapack/explore-html/d6/da6/dlasrt2_8f_source.html
http://www.netlib.org/scalapack/explore-html/d9/d96/slasrt2_8f_source.html
is
Code: Select all
IF( ENDD-START.GT.0 ) THEN

which is almost always true.

It should be
Code: Select all
IF( ENDD-START.LE.SELECT .AND. ENDD-START.GT.0 ) THEN

(same as line 105 of dlapst.f) so that insertion sort is only performed for small arrays.
quenteridean
 
Posts: 1
Joined: Tue May 03, 2016 2:37 pm

Return to Bug report

Who is online

Users browsing this forum: No registered users and 1 guest