Merge pull request #1023 from Goddan-wq/lapack_issue1021

changing the order of loop to improve performance

Author: langou

SRC/cgetc2.f

@@ -176,8 +176,8 @@ SUBROUTINE CGETC2( N, A, LDA, IPIV, JPIV, INFO )
 *        Find max element in matrix A
 *
          XMAX = ZERO
-         DO 20 IP = I, N
-            DO 10 JP = I, N
+         DO 20 JP = I, N
+            DO 10 IP = I, N
                IF( ABS( A( IP, JP ) ).GE.XMAX ) THEN
                   XMAX = ABS( A( IP, JP ) )
                   IPV = IP

SRC/dgetc2.f

@@ -176,8 +176,8 @@ SUBROUTINE DGETC2( N, A, LDA, IPIV, JPIV, INFO )
 *        Find max element in matrix A
 *
          XMAX = ZERO
-         DO 20 IP = I, N
-            DO 10 JP = I, N
+         DO 20 JP = I, N
+            DO 10 IP = I, N
                IF( ABS( A( IP, JP ) ).GE.XMAX ) THEN
                   XMAX = ABS( A( IP, JP ) )
                   IPV = IP

SRC/sgetc2.f

@@ -176,8 +176,8 @@ SUBROUTINE SGETC2( N, A, LDA, IPIV, JPIV, INFO )
 *        Find max element in matrix A
 *
          XMAX = ZERO
-         DO 20 IP = I, N
-            DO 10 JP = I, N
+         DO 20 JP = I, N
+            DO 10 IP = I, N
                IF( ABS( A( IP, JP ) ).GE.XMAX ) THEN
                   XMAX = ABS( A( IP, JP ) )
                   IPV = IP

SRC/zgetc2.f

@@ -176,8 +176,8 @@ SUBROUTINE ZGETC2( N, A, LDA, IPIV, JPIV, INFO )
 *        Find max element in matrix A
 *
          XMAX = ZERO
-         DO 20 IP = I, N
-            DO 10 JP = I, N
+         DO 20 JP = I, N
+            DO 10 IP = I, N
                IF( ABS( A( IP, JP ) ).GE.XMAX ) THEN
                   XMAX = ABS( A( IP, JP ) )
                   IPV = IP