Skip to content

Commit ff5dc3e

Browse files
martin-frbgmartin-frbg
authored
Change loop ordering to improve performance (Reference-LAPACK PR 1023)
1 parent d26960a commit ff5dc3e

4 files changed

Lines changed: 20 additions & 28 deletions

File tree

‎lapack-netlib/SRC/cgetc2.f‎

Lines changed: 5 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -5,15 +5,13 @@
55
* Online html documentation available at
66
* http://www.netlib.org/lapack/explore-html/
77
*
8-
*> \htmlonly
98
*> Download CGETC2 + dependencies
109
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/cgetc2.f">
1110
*> [TGZ]</a>
1211
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/cgetc2.f">
1312
*> [ZIP]</a>
1413
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/cgetc2.f">
1514
*> [TXT]</a>
16-
*> \endhtmlonly
1715
*
1816
* Definition:
1917
* ===========
@@ -98,7 +96,7 @@
9896
*> \author Univ. of Colorado Denver
9997
*> \author NAG Ltd.
10098
*
101-
*> \ingroup complexGEauxiliary
99+
*> \ingroup getc2
102100
*
103101
*> \par Contributors:
104102
* ==================
@@ -108,6 +106,7 @@
108106
*
109107
* =====================================================================
110108
SUBROUTINE CGETC2( N, A, LDA, IPIV, JPIV, INFO )
109+
IMPLICIT NONE
111110
*
112111
* -- LAPACK auxiliary routine --
113112
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -132,7 +131,7 @@ SUBROUTINE CGETC2( N, A, LDA, IPIV, JPIV, INFO )
132131
REAL BIGNUM, EPS, SMIN, SMLNUM, XMAX
133132
* ..
134133
* .. External Subroutines ..
135-
EXTERNAL CGERU, CSWAP, SLABAD
134+
EXTERNAL CGERU, CSWAP
136135
* ..
137136
* .. External Functions ..
138137
REAL SLAMCH
@@ -155,7 +154,6 @@ SUBROUTINE CGETC2( N, A, LDA, IPIV, JPIV, INFO )
155154
EPS = SLAMCH( 'P' )
156155
SMLNUM = SLAMCH( 'S' ) / EPS
157156
BIGNUM = ONE / SMLNUM
158-
CALL SLABAD( SMLNUM, BIGNUM )
159157
*
160158
* Handle the case N=1 by itself
161159
*
@@ -177,8 +175,8 @@ SUBROUTINE CGETC2( N, A, LDA, IPIV, JPIV, INFO )
177175
* Find max element in matrix A
178176
*
179177
XMAX = ZERO
180-
DO 20 IP = I, N
181-
DO 10 JP = I, N
178+
DO 20 JP = I, N
179+
DO 10 IP = I, N
182180
IF( ABS( A( IP, JP ) ).GE.XMAX ) THEN
183181
XMAX = ABS( A( IP, JP ) )
184182
IPV = IP

‎lapack-netlib/SRC/dgetc2.f‎

Lines changed: 5 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -5,15 +5,13 @@
55
* Online html documentation available at
66
* http://www.netlib.org/lapack/explore-html/
77
*
8-
*> \htmlonly
98
*> Download DGETC2 + dependencies
109
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dgetc2.f">
1110
*> [TGZ]</a>
1211
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dgetc2.f">
1312
*> [ZIP]</a>
1413
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dgetc2.f">
1514
*> [TXT]</a>
16-
*> \endhtmlonly
1715
*
1816
* Definition:
1917
* ===========
@@ -98,7 +96,7 @@
9896
*> \author Univ. of Colorado Denver
9997
*> \author NAG Ltd.
10098
*
101-
*> \ingroup doubleGEauxiliary
99+
*> \ingroup getc2
102100
*
103101
*> \par Contributors:
104102
* ==================
@@ -108,6 +106,7 @@
108106
*
109107
* =====================================================================
110108
SUBROUTINE DGETC2( N, A, LDA, IPIV, JPIV, INFO )
109+
IMPLICIT NONE
111110
*
112111
* -- LAPACK auxiliary routine --
113112
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -132,7 +131,7 @@ SUBROUTINE DGETC2( N, A, LDA, IPIV, JPIV, INFO )
132131
DOUBLE PRECISION BIGNUM, EPS, SMIN, SMLNUM, XMAX
133132
* ..
134133
* .. External Subroutines ..
135-
EXTERNAL DGER, DSWAP, DLABAD
134+
EXTERNAL DGER, DSWAP
136135
* ..
137136
* .. External Functions ..
138137
DOUBLE PRECISION DLAMCH
@@ -155,7 +154,6 @@ SUBROUTINE DGETC2( N, A, LDA, IPIV, JPIV, INFO )
155154
EPS = DLAMCH( 'P' )
156155
SMLNUM = DLAMCH( 'S' ) / EPS
157156
BIGNUM = ONE / SMLNUM
158-
CALL DLABAD( SMLNUM, BIGNUM )
159157
*
160158
* Handle the case N=1 by itself
161159
*
@@ -177,8 +175,8 @@ SUBROUTINE DGETC2( N, A, LDA, IPIV, JPIV, INFO )
177175
* Find max element in matrix A
178176
*
179177
XMAX = ZERO
180-
DO 20 IP = I, N
181-
DO 10 JP = I, N
178+
DO 20 JP = I, N
179+
DO 10 IP = I, N
182180
IF( ABS( A( IP, JP ) ).GE.XMAX ) THEN
183181
XMAX = ABS( A( IP, JP ) )
184182
IPV = IP

‎lapack-netlib/SRC/sgetc2.f‎

Lines changed: 5 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -5,15 +5,13 @@
55
* Online html documentation available at
66
* http://www.netlib.org/lapack/explore-html/
77
*
8-
*> \htmlonly
98
*> Download SGETC2 + dependencies
109
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/sgetc2.f">
1110
*> [TGZ]</a>
1211
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/sgetc2.f">
1312
*> [ZIP]</a>
1413
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/sgetc2.f">
1514
*> [TXT]</a>
16-
*> \endhtmlonly
1715
*
1816
* Definition:
1917
* ===========
@@ -98,7 +96,7 @@
9896
*> \author Univ. of Colorado Denver
9997
*> \author NAG Ltd.
10098
*
101-
*> \ingroup realGEauxiliary
99+
*> \ingroup getc2
102100
*
103101
*> \par Contributors:
104102
* ==================
@@ -108,6 +106,7 @@
108106
*
109107
* =====================================================================
110108
SUBROUTINE SGETC2( N, A, LDA, IPIV, JPIV, INFO )
109+
IMPLICIT NONE
111110
*
112111
* -- LAPACK auxiliary routine --
113112
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -132,7 +131,7 @@ SUBROUTINE SGETC2( N, A, LDA, IPIV, JPIV, INFO )
132131
REAL BIGNUM, EPS, SMIN, SMLNUM, XMAX
133132
* ..
134133
* .. External Subroutines ..
135-
EXTERNAL SGER, SLABAD, SSWAP
134+
EXTERNAL SGER, SSWAP
136135
* ..
137136
* .. External Functions ..
138137
REAL SLAMCH
@@ -155,7 +154,6 @@ SUBROUTINE SGETC2( N, A, LDA, IPIV, JPIV, INFO )
155154
EPS = SLAMCH( 'P' )
156155
SMLNUM = SLAMCH( 'S' ) / EPS
157156
BIGNUM = ONE / SMLNUM
158-
CALL SLABAD( SMLNUM, BIGNUM )
159157
*
160158
* Handle the case N=1 by itself
161159
*
@@ -177,8 +175,8 @@ SUBROUTINE SGETC2( N, A, LDA, IPIV, JPIV, INFO )
177175
* Find max element in matrix A
178176
*
179177
XMAX = ZERO
180-
DO 20 IP = I, N
181-
DO 10 JP = I, N
178+
DO 20 JP = I, N
179+
DO 10 IP = I, N
182180
IF( ABS( A( IP, JP ) ).GE.XMAX ) THEN
183181
XMAX = ABS( A( IP, JP ) )
184182
IPV = IP

‎lapack-netlib/SRC/zgetc2.f‎

Lines changed: 5 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -5,15 +5,13 @@
55
* Online html documentation available at
66
* http://www.netlib.org/lapack/explore-html/
77
*
8-
*> \htmlonly
98
*> Download ZGETC2 + dependencies
109
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zgetc2.f">
1110
*> [TGZ]</a>
1211
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zgetc2.f">
1312
*> [ZIP]</a>
1413
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zgetc2.f">
1514
*> [TXT]</a>
16-
*> \endhtmlonly
1715
*
1816
* Definition:
1917
* ===========
@@ -98,7 +96,7 @@
9896
*> \author Univ. of Colorado Denver
9997
*> \author NAG Ltd.
10098
*
101-
*> \ingroup complex16GEauxiliary
99+
*> \ingroup getc2
102100
*
103101
*> \par Contributors:
104102
* ==================
@@ -108,6 +106,7 @@
108106
*
109107
* =====================================================================
110108
SUBROUTINE ZGETC2( N, A, LDA, IPIV, JPIV, INFO )
109+
IMPLICIT NONE
111110
*
112111
* -- LAPACK auxiliary routine --
113112
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -132,7 +131,7 @@ SUBROUTINE ZGETC2( N, A, LDA, IPIV, JPIV, INFO )
132131
DOUBLE PRECISION BIGNUM, EPS, SMIN, SMLNUM, XMAX
133132
* ..
134133
* .. External Subroutines ..
135-
EXTERNAL ZGERU, ZSWAP, DLABAD
134+
EXTERNAL ZGERU, ZSWAP
136135
* ..
137136
* .. External Functions ..
138137
DOUBLE PRECISION DLAMCH
@@ -155,7 +154,6 @@ SUBROUTINE ZGETC2( N, A, LDA, IPIV, JPIV, INFO )
155154
EPS = DLAMCH( 'P' )
156155
SMLNUM = DLAMCH( 'S' ) / EPS
157156
BIGNUM = ONE / SMLNUM
158-
CALL DLABAD( SMLNUM, BIGNUM )
159157
*
160158
* Handle the case N=1 by itself
161159
*
@@ -177,8 +175,8 @@ SUBROUTINE ZGETC2( N, A, LDA, IPIV, JPIV, INFO )
177175
* Find max element in matrix A
178176
*
179177
XMAX = ZERO
180-
DO 20 IP = I, N
181-
DO 10 JP = I, N
178+
DO 20 JP = I, N
179+
DO 10 IP = I, N
182180
IF( ABS( A( IP, JP ) ).GE.XMAX ) THEN
183181
XMAX = ABS( A( IP, JP ) )
184182
IPV = IP

0 commit comments

Comments
 (0)