DGETC2 is slow

Post here if you have a question about LAPACK or ScaLAPACK algorithm or data format

DGETC2 is slow

Postby yshl0000 » Sat Jan 26, 2013 1:05 am

Hello, all.

I tested DGETC2 subroutine that is LU factorization with complete pivoting of general matrix,
and it is slower than my subroutine.
In dgetc2.f, the subroutine find maximum element by row major order.

Code: Select all
         DO 20 IP = I, N
            DO 10 JP = I, N
               IF( ABS( A( IP, JP ) ).GE.XMAX ) THEN
                  XMAX = ABS( A( IP, JP ) )
                  IPV = IP
                  JPV = JP
               END IF
   10       CONTINUE
   20    CONTINUE

I think DGETC2 will be faster if it find by column order and use IDAMAX in BLAS.
Are there some reasons for not using BLAS subroutine ?

Thanks in advance.
Posts: 1
Joined: Fri Jan 25, 2013 11:55 pm

Return to Algorithm / Data

Who is online

Users browsing this forum: No registered users and 1 guest