Add quick return for *LARR* routines

eugene.chereshnev · eugene.chereshnev · commit 16db9c24ea1d · 2017-03-31T13:06:11.000-07:00
diff --git a/SRC/cgejsv.f b/SRC/cgejsv.f
diff --git a/SRC/clarrv.f b/SRC/clarrv.f
@@ -348,6 +348,13 @@ SUBROUTINE CLARRV( N, VL, VU, D, L, PIVMIN,
 *     ..
 
       INFO = 0
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
+*
 *     The first N entries of WORK are reserved for the eigenvalues
       INDLD = N+1
       INDLLD= 2*N+1
diff --git a/SRC/dlarra.f b/SRC/dlarra.f
@@ -167,7 +167,13 @@ SUBROUTINE DLARRA( N, D, E, E2, SPLTOL, TNRM,
 *     .. Executable Statements ..
 *
       INFO = 0
-
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
+*
 *     Compute splitting points
       NSPLIT = 1
       IF(SPLTOL.LT.ZERO) THEN
diff --git a/SRC/dlarrb.f b/SRC/dlarrb.f
@@ -236,6 +236,12 @@ SUBROUTINE DLARRB( N, D, LLD, IFIRST, ILAST, RTOL1,
 *     .. Executable Statements ..
 *
       INFO = 0
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
 *
       MAXITR = INT( ( LOG( SPDIAM+PIVMIN )-LOG( PIVMIN ) ) /
      $           LOG( TWO ) ) + 2
diff --git a/SRC/dlarrc.f b/SRC/dlarrc.f
@@ -170,6 +170,13 @@ SUBROUTINE DLARRC( JOBT, N, VL, VU, D, E, PIVMIN,
 *     .. Executable Statements ..
 *
       INFO = 0
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
+*
       LCNT = 0
       RCNT = 0
       EIGCNT = 0
diff --git a/SRC/dlarrd.f b/SRC/dlarrd.f
@@ -385,6 +385,12 @@ SUBROUTINE DLARRD( RANGE, ORDER, N, VL, VU, IL, IU, GERS,
 *
       INFO = 0
 *
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
+*
 *     Decode RANGE
 *
       IF( LSAME( RANGE, 'A' ) ) THEN
diff --git a/SRC/dlarre.f b/SRC/dlarre.f
@@ -370,7 +370,12 @@ SUBROUTINE DLARRE( RANGE, N, VL, VU, IL, IU, D, E, E2,
 *
 
       INFO = 0
-
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
 *
 *     Decode RANGE
 *
diff --git a/SRC/dlarrf.f b/SRC/dlarrf.f
@@ -239,6 +239,13 @@ SUBROUTINE DLARRF( N, D, L, LD, CLSTRT, CLEND,
 *     .. Executable Statements ..
 *
       INFO = 0
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
+*
       FACT = DBLE(2**KTRYMAX)
       EPS = DLAMCH( 'Precision' )
       SHIFT = 0
diff --git a/SRC/dlarrj.f b/SRC/dlarrj.f
@@ -203,6 +203,12 @@ SUBROUTINE DLARRJ( N, D, E2, IFIRST, ILAST,
 *     .. Executable Statements ..
 *
       INFO = 0
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
 *
       MAXITR = INT( ( LOG( SPDIAM+PIVMIN )-LOG( PIVMIN ) ) /
      $           LOG( TWO ) ) + 2
diff --git a/SRC/dlarrk.f b/SRC/dlarrk.f
@@ -179,6 +179,13 @@ SUBROUTINE DLARRK( N, IW, GL, GU,
 *     ..
 *     .. Executable Statements ..
 *
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         INFO = 0
+         RETURN
+      END IF
+*
 *     Get machine constants
       EPS = DLAMCH( 'P' )
 
diff --git a/SRC/dlarrr.f b/SRC/dlarrr.f
@@ -130,6 +130,13 @@ SUBROUTINE DLARRR( N, D, E, INFO )
 *     ..
 *     .. Executable Statements ..
 *
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         INFO = 0
+         RETURN
+      END IF
+*
 *     As a default, do NOT go for relative-accuracy preserving computations.
       INFO = 1
 
diff --git a/SRC/dlarrv.f b/SRC/dlarrv.f
@@ -344,6 +344,13 @@ SUBROUTINE DLARRV( N, VL, VU, D, L, PIVMIN,
 *     ..
 
       INFO = 0
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
+*
 *     The first N entries of WORK are reserved for the eigenvalues
       INDLD = N+1
       INDLLD= 2*N+1
diff --git a/SRC/slarra.f b/SRC/slarra.f
@@ -167,7 +167,13 @@ SUBROUTINE SLARRA( N, D, E, E2, SPLTOL, TNRM,
 *     .. Executable Statements ..
 *
       INFO = 0
-
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
+*
 *     Compute splitting points
       NSPLIT = 1
       IF(SPLTOL.LT.ZERO) THEN
diff --git a/SRC/slarrb.f b/SRC/slarrb.f
@@ -236,6 +236,12 @@ SUBROUTINE SLARRB( N, D, LLD, IFIRST, ILAST, RTOL1,
 *     .. Executable Statements ..
 *
       INFO = 0
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
 *
       MAXITR = INT( ( LOG( SPDIAM+PIVMIN )-LOG( PIVMIN ) ) /
      $           LOG( TWO ) ) + 2
diff --git a/SRC/slarrc.f b/SRC/slarrc.f
@@ -170,6 +170,13 @@ SUBROUTINE SLARRC( JOBT, N, VL, VU, D, E, PIVMIN,
 *     .. Executable Statements ..
 *
       INFO = 0
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
+*
       LCNT = 0
       RCNT = 0
       EIGCNT = 0
diff --git a/SRC/slarrd.f b/SRC/slarrd.f
@@ -385,6 +385,12 @@ SUBROUTINE SLARRD( RANGE, ORDER, N, VL, VU, IL, IU, GERS,
 *
       INFO = 0
 *
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
+*
 *     Decode RANGE
 *
       IF( LSAME( RANGE, 'A' ) ) THEN
diff --git a/SRC/slarre.f b/SRC/slarre.f
@@ -370,7 +370,12 @@ SUBROUTINE SLARRE( RANGE, N, VL, VU, IL, IU, D, E, E2,
 *
 
       INFO = 0
-
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
 *
 *     Decode RANGE
 *
diff --git a/SRC/slarrf.f b/SRC/slarrf.f
@@ -239,6 +239,13 @@ SUBROUTINE SLARRF( N, D, L, LD, CLSTRT, CLEND,
 *     .. Executable Statements ..
 *
       INFO = 0
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
+*
       FACT = REAL(2**KTRYMAX)
       EPS = SLAMCH( 'Precision' )
       SHIFT = 0
diff --git a/SRC/slarrj.f b/SRC/slarrj.f
@@ -203,6 +203,12 @@ SUBROUTINE SLARRJ( N, D, E2, IFIRST, ILAST,
 *     .. Executable Statements ..
 *
       INFO = 0
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
 *
       MAXITR = INT( ( LOG( SPDIAM+PIVMIN )-LOG( PIVMIN ) ) /
      $           LOG( TWO ) ) + 2
diff --git a/SRC/slarrk.f b/SRC/slarrk.f
@@ -179,6 +179,13 @@ SUBROUTINE SLARRK( N, IW, GL, GU,
 *     ..
 *     .. Executable Statements ..
 *
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         INFO = 0
+         RETURN
+      END IF
+*
 *     Get machine constants
       EPS = SLAMCH( 'P' )
 
diff --git a/SRC/slarrr.f b/SRC/slarrr.f
@@ -130,6 +130,13 @@ SUBROUTINE SLARRR( N, D, E, INFO )
 *     ..
 *     .. Executable Statements ..
 *
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         INFO = 0
+         RETURN
+      END IF
+*
 *     As a default, do NOT go for relative-accuracy preserving computations.
       INFO = 1
 
diff --git a/SRC/slarrv.f b/SRC/slarrv.f
@@ -344,6 +344,13 @@ SUBROUTINE SLARRV( N, VL, VU, D, L, PIVMIN,
 *     ..
 
       INFO = 0
+*
+*     Quick return if possible
+*
+      IF( N.LE.0 ) THEN
+         RETURN
+      END IF
+*
 *     The first N entries of WORK are reserved for the eigenvalues
       INDLD = N+1
       INDLLD= 2*N+1
diff --git a/SRC/zgejsv.f b/SRC/zgejsv.f
diff --git a/SRC/zlarrv.f b/SRC/zlarrv.f

Original file line number	Diff line number	Diff line change
`@@ -236,6 +236,12 @@ SUBROUTINE DLARRB( N, D, LLD, IFIRST, ILAST, RTOL1,`
`236`	`236`	`* .. Executable Statements ..`
`237`	`237`	`*`
`238`	`238`	`INFO = 0`
	`239`	`+*`
	`240`	`+* Quick return if possible`
	`241`	`+*`
	`242`	`+ IF( N.LE.0 ) THEN`
	`243`	`+ RETURN`
	`244`	`+ END IF`
`239`	`245`	`*`
`240`	`246`	`MAXITR = INT( ( LOG( SPDIAM+PIVMIN )-LOG( PIVMIN ) ) /`
`241`	`247`	`$ LOG( TWO ) ) + 2`
Original file line number	Diff line number	Diff line change
`@@ -385,6 +385,12 @@ SUBROUTINE DLARRD( RANGE, ORDER, N, VL, VU, IL, IU, GERS,`
`385`	`385`	`*`
`386`	`386`	`INFO = 0`
`387`	`387`	`*`
	`388`	`+* Quick return if possible`
	`389`	`+*`
	`390`	`+ IF( N.LE.0 ) THEN`
	`391`	`+ RETURN`
	`392`	`+ END IF`
	`393`	`+*`
`388`	`394`	`* Decode RANGE`
`389`	`395`	`*`
`390`	`396`	`IF( LSAME( RANGE, 'A' ) ) THEN`
Original file line number	Diff line number	Diff line change
`@@ -370,7 +370,12 @@ SUBROUTINE DLARRE( RANGE, N, VL, VU, IL, IU, D, E, E2,`
`370`	`370`	`*`
`371`	`371`
`372`	`372`	`INFO = 0`
`373`		`-`
	`373`	`+*`
	`374`	`+* Quick return if possible`
	`375`	`+*`
	`376`	`+ IF( N.LE.0 ) THEN`
	`377`	`+ RETURN`
	`378`	`+ END IF`
`374`	`379`	`*`
`375`	`380`	`* Decode RANGE`
`376`	`381`	`*`
Original file line number	Diff line number	Diff line change
`@@ -203,6 +203,12 @@ SUBROUTINE DLARRJ( N, D, E2, IFIRST, ILAST,`
`203`	`203`	`* .. Executable Statements ..`
`204`	`204`	`*`
`205`	`205`	`INFO = 0`
	`206`	`+*`
	`207`	`+* Quick return if possible`
	`208`	`+*`
	`209`	`+ IF( N.LE.0 ) THEN`
	`210`	`+ RETURN`
	`211`	`+ END IF`
`206`	`212`	`*`
`207`	`213`	`MAXITR = INT( ( LOG( SPDIAM+PIVMIN )-LOG( PIVMIN ) ) /`
`208`	`214`	`$ LOG( TWO ) ) + 2`