Skip to content

Commit a37f229

Browse files
authored
Replace the partially defined constant INDIBL with 1
1 parent 7ac31fa commit a37f229

File tree

10 files changed

+70
-80
lines changed

10 files changed

+70
-80
lines changed

SRC/chbgvx.f

Lines changed: 7 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -327,7 +327,7 @@ SUBROUTINE CHBGVX( JOBZ, RANGE, UPLO, N, KA, KB, AB, LDAB, BB,
327327
* .. Local Scalars ..
328328
LOGICAL ALLEIG, INDEIG, TEST, UPPER, VALEIG, WANTZ
329329
CHARACTER ORDER, VECT
330-
INTEGER I, IINFO, INDD, INDE, INDEE, INDIBL, INDISP,
330+
INTEGER I, IINFO, INDD, INDE, INDEE, INDISP,
331331
$ INDIWK, INDRWK, INDWRK, ITMP1, J, JJ, NSPLIT
332332
REAL TMP1
333333
* ..
@@ -470,17 +470,16 @@ SUBROUTINE CHBGVX( JOBZ, RANGE, UPLO, N, KA, KB, AB, LDAB, BB,
470470
ELSE
471471
ORDER = 'E'
472472
END IF
473-
INDIBL = 1
474-
INDISP = INDIBL + N
473+
INDISP = 1 + N
475474
INDIWK = INDISP + N
476475
CALL SSTEBZ( RANGE, ORDER, N, VL, VU, IL, IU, ABSTOL,
477476
$ RWORK( INDD ), RWORK( INDE ), M, NSPLIT, W,
478-
$ IWORK( INDIBL ), IWORK( INDISP ), RWORK( INDRWK ),
477+
$ IWORK( 1 ), IWORK( INDISP ), RWORK( INDRWK ),
479478
$ IWORK( INDIWK ), INFO )
480479
*
481480
IF( WANTZ ) THEN
482481
CALL CSTEIN( N, RWORK( INDD ), RWORK( INDE ), M, W,
483-
$ IWORK( INDIBL ), IWORK( INDISP ), Z, LDZ,
482+
$ IWORK( 1 ), IWORK( INDISP ), Z, LDZ,
484483
$ RWORK( INDRWK ), IWORK( INDIWK ), IFAIL, INFO )
485484
*
486485
* Apply unitary matrix used in reduction to tridiagonal
@@ -510,11 +509,11 @@ SUBROUTINE CHBGVX( JOBZ, RANGE, UPLO, N, KA, KB, AB, LDAB, BB,
510509
40 CONTINUE
511510
*
512511
IF( I.NE.0 ) THEN
513-
ITMP1 = IWORK( INDIBL+I-1 )
512+
ITMP1 = IWORK( 1 + I-1 )
514513
W( I ) = W( J )
515-
IWORK( INDIBL+I-1 ) = IWORK( INDIBL+J-1 )
514+
IWORK( 1 + I-1 ) = IWORK( 1 + J-1 )
516515
W( J ) = TMP1
517-
IWORK( INDIBL+J-1 ) = ITMP1
516+
IWORK( 1 + J-1 ) = ITMP1
518517
CALL CSWAP( N, Z( 1, I ), 1, Z( 1, J ), 1 )
519518
IF( INFO.NE.0 ) THEN
520519
ITMP1 = IFAIL( I )

SRC/chpevx.f

Lines changed: 7 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -264,7 +264,7 @@ SUBROUTINE CHPEVX( JOBZ, RANGE, UPLO, N, AP, VL, VU, IL, IU,
264264
* .. Local Scalars ..
265265
LOGICAL ALLEIG, INDEIG, TEST, VALEIG, WANTZ
266266
CHARACTER ORDER
267-
INTEGER I, IINFO, IMAX, INDD, INDE, INDEE, INDIBL,
267+
INTEGER I, IINFO, IMAX, INDD, INDE, INDEE,
268268
$ INDISP, INDIWK, INDRWK, INDTAU, INDWRK, ISCALE,
269269
$ ITMP1, J, JJ, NSPLIT
270270
REAL ABSTLL, ANRM, BIGNUM, EPS, RMAX, RMIN, SAFMIN,
@@ -434,17 +434,16 @@ SUBROUTINE CHPEVX( JOBZ, RANGE, UPLO, N, AP, VL, VU, IL, IU,
434434
ELSE
435435
ORDER = 'E'
436436
END IF
437-
INDIBL = 1
438-
INDISP = INDIBL + N
437+
INDISP = 1 + N
439438
INDIWK = INDISP + N
440439
CALL SSTEBZ( RANGE, ORDER, N, VLL, VUU, IL, IU, ABSTLL,
441440
$ RWORK( INDD ), RWORK( INDE ), M, NSPLIT, W,
442-
$ IWORK( INDIBL ), IWORK( INDISP ), RWORK( INDRWK ),
441+
$ IWORK( 1 ), IWORK( INDISP ), RWORK( INDRWK ),
443442
$ IWORK( INDIWK ), INFO )
444443
*
445444
IF( WANTZ ) THEN
446445
CALL CSTEIN( N, RWORK( INDD ), RWORK( INDE ), M, W,
447-
$ IWORK( INDIBL ), IWORK( INDISP ), Z, LDZ,
446+
$ IWORK( 1 ), IWORK( INDISP ), Z, LDZ,
448447
$ RWORK( INDRWK ), IWORK( INDIWK ), IFAIL, INFO )
449448
*
450449
* Apply unitary matrix used in reduction to tridiagonal
@@ -482,11 +481,11 @@ SUBROUTINE CHPEVX( JOBZ, RANGE, UPLO, N, AP, VL, VU, IL, IU,
482481
30 CONTINUE
483482
*
484483
IF( I.NE.0 ) THEN
485-
ITMP1 = IWORK( INDIBL+I-1 )
484+
ITMP1 = IWORK( 1 + I-1 )
486485
W( I ) = W( J )
487-
IWORK( INDIBL+I-1 ) = IWORK( INDIBL+J-1 )
486+
IWORK( 1 + I-1 ) = IWORK( 1 + J-1 )
488487
W( J ) = TMP1
489-
IWORK( INDIBL+J-1 ) = ITMP1
488+
IWORK( 1 + J-1 ) = ITMP1
490489
CALL CSWAP( N, Z( 1, I ), 1, Z( 1, J ), 1 )
491490
IF( INFO.NE.0 ) THEN
492491
ITMP1 = IFAIL( I )

SRC/dsbgvx.f

Lines changed: 7 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -317,7 +317,7 @@ SUBROUTINE DSBGVX( JOBZ, RANGE, UPLO, N, KA, KB, AB, LDAB, BB,
317317
* .. Local Scalars ..
318318
LOGICAL ALLEIG, INDEIG, TEST, UPPER, VALEIG, WANTZ
319319
CHARACTER ORDER, VECT
320-
INTEGER I, IINFO, INDD, INDE, INDEE, INDIBL, INDISP,
320+
INTEGER I, IINFO, INDD, INDE, INDEE, INDISP,
321321
$ INDIWO, INDWRK, ITMP1, J, JJ, NSPLIT
322322
DOUBLE PRECISION TMP1
323323
* ..
@@ -457,17 +457,16 @@ SUBROUTINE DSBGVX( JOBZ, RANGE, UPLO, N, KA, KB, AB, LDAB, BB,
457457
ELSE
458458
ORDER = 'E'
459459
END IF
460-
INDIBL = 1
461-
INDISP = INDIBL + N
460+
INDISP = 1 + N
462461
INDIWO = INDISP + N
463462
CALL DSTEBZ( RANGE, ORDER, N, VL, VU, IL, IU, ABSTOL,
464463
$ WORK( INDD ), WORK( INDE ), M, NSPLIT, W,
465-
$ IWORK( INDIBL ), IWORK( INDISP ), WORK( INDWRK ),
464+
$ IWORK( 1 ), IWORK( INDISP ), WORK( INDWRK ),
466465
$ IWORK( INDIWO ), INFO )
467466
*
468467
IF( WANTZ ) THEN
469468
CALL DSTEIN( N, WORK( INDD ), WORK( INDE ), M, W,
470-
$ IWORK( INDIBL ), IWORK( INDISP ), Z, LDZ,
469+
$ IWORK( 1 ), IWORK( INDISP ), Z, LDZ,
471470
$ WORK( INDWRK ), IWORK( INDIWO ), IFAIL, INFO )
472471
*
473472
* Apply transformation matrix used in reduction to tridiagonal
@@ -497,11 +496,11 @@ SUBROUTINE DSBGVX( JOBZ, RANGE, UPLO, N, KA, KB, AB, LDAB, BB,
497496
40 CONTINUE
498497
*
499498
IF( I.NE.0 ) THEN
500-
ITMP1 = IWORK( INDIBL+I-1 )
499+
ITMP1 = IWORK( 1 + I-1 )
501500
W( I ) = W( J )
502-
IWORK( INDIBL+I-1 ) = IWORK( INDIBL+J-1 )
501+
IWORK( 1 + I-1 ) = IWORK( 1 + J-1 )
503502
W( J ) = TMP1
504-
IWORK( INDIBL+J-1 ) = ITMP1
503+
IWORK( 1 + J-1 ) = ITMP1
505504
CALL DSWAP( N, Z( 1, I ), 1, Z( 1, J ), 1 )
506505
IF( INFO.NE.0 ) THEN
507506
ITMP1 = IFAIL( I )

SRC/dspevx.f

Lines changed: 7 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -255,7 +255,7 @@ SUBROUTINE DSPEVX( JOBZ, RANGE, UPLO, N, AP, VL, VU, IL, IU,
255255
* .. Local Scalars ..
256256
LOGICAL ALLEIG, INDEIG, TEST, VALEIG, WANTZ
257257
CHARACTER ORDER
258-
INTEGER I, IINFO, IMAX, INDD, INDE, INDEE, INDIBL,
258+
INTEGER I, IINFO, IMAX, INDD, INDE, INDEE,
259259
$ INDISP, INDIWO, INDTAU, INDWRK, ISCALE, ITMP1,
260260
$ J, JJ, NSPLIT
261261
DOUBLE PRECISION ABSTLL, ANRM, BIGNUM, EPS, RMAX, RMIN, SAFMIN,
@@ -424,17 +424,16 @@ SUBROUTINE DSPEVX( JOBZ, RANGE, UPLO, N, AP, VL, VU, IL, IU,
424424
ELSE
425425
ORDER = 'E'
426426
END IF
427-
INDIBL = 1
428-
INDISP = INDIBL + N
427+
INDISP = 1 + N
429428
INDIWO = INDISP + N
430429
CALL DSTEBZ( RANGE, ORDER, N, VLL, VUU, IL, IU, ABSTLL,
431430
$ WORK( INDD ), WORK( INDE ), M, NSPLIT, W,
432-
$ IWORK( INDIBL ), IWORK( INDISP ), WORK( INDWRK ),
431+
$ IWORK( 1 ), IWORK( INDISP ), WORK( INDWRK ),
433432
$ IWORK( INDIWO ), INFO )
434433
*
435434
IF( WANTZ ) THEN
436435
CALL DSTEIN( N, WORK( INDD ), WORK( INDE ), M, W,
437-
$ IWORK( INDIBL ), IWORK( INDISP ), Z, LDZ,
436+
$ IWORK( 1 ), IWORK( INDISP ), Z, LDZ,
438437
$ WORK( INDWRK ), IWORK( INDIWO ), IFAIL, INFO )
439438
*
440439
* Apply orthogonal matrix used in reduction to tridiagonal
@@ -471,11 +470,11 @@ SUBROUTINE DSPEVX( JOBZ, RANGE, UPLO, N, AP, VL, VU, IL, IU,
471470
30 CONTINUE
472471
*
473472
IF( I.NE.0 ) THEN
474-
ITMP1 = IWORK( INDIBL+I-1 )
473+
ITMP1 = IWORK( 1 + I-1 )
475474
W( I ) = W( J )
476-
IWORK( INDIBL+I-1 ) = IWORK( INDIBL+J-1 )
475+
IWORK( 1 + I-1 ) = IWORK( 1 + J-1 )
477476
W( J ) = TMP1
478-
IWORK( INDIBL+J-1 ) = ITMP1
477+
IWORK( 1 + J-1 ) = ITMP1
479478
CALL DSWAP( N, Z( 1, I ), 1, Z( 1, J ), 1 )
480479
IF( INFO.NE.0 ) THEN
481480
ITMP1 = IFAIL( I )

SRC/dstevx.f

Lines changed: 7 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -248,7 +248,7 @@ SUBROUTINE DSTEVX( JOBZ, RANGE, N, D, E, VL, VU, IL, IU, ABSTOL,
248248
* .. Local Scalars ..
249249
LOGICAL ALLEIG, INDEIG, TEST, VALEIG, WANTZ
250250
CHARACTER ORDER
251-
INTEGER I, IMAX, INDIBL, INDISP, INDIWO, INDWRK,
251+
INTEGER I, IMAX, INDISP, INDIWO, INDWRK,
252252
$ ISCALE, ITMP1, J, JJ, NSPLIT
253253
DOUBLE PRECISION BIGNUM, EPS, RMAX, RMIN, SAFMIN, SIGMA, SMLNUM,
254254
$ TMP1, TNRM, VLL, VUU
@@ -399,15 +399,14 @@ SUBROUTINE DSTEVX( JOBZ, RANGE, N, D, E, VL, VU, IL, IU, ABSTOL,
399399
ORDER = 'E'
400400
END IF
401401
INDWRK = 1
402-
INDIBL = 1
403-
INDISP = INDIBL + N
402+
INDISP = 1 + N
404403
INDIWO = INDISP + N
405404
CALL DSTEBZ( RANGE, ORDER, N, VLL, VUU, IL, IU, ABSTOL, D, E, M,
406-
$ NSPLIT, W, IWORK( INDIBL ), IWORK( INDISP ),
405+
$ NSPLIT, W, IWORK( 1 ), IWORK( INDISP ),
407406
$ WORK( INDWRK ), IWORK( INDIWO ), INFO )
408407
*
409408
IF( WANTZ ) THEN
410-
CALL DSTEIN( N, D, E, M, W, IWORK( INDIBL ), IWORK( INDISP ),
409+
CALL DSTEIN( N, D, E, M, W, IWORK( 1 ), IWORK( INDISP ),
411410
$ Z, LDZ, WORK( INDWRK ), IWORK( INDIWO ), IFAIL,
412411
$ INFO )
413412
END IF
@@ -439,11 +438,11 @@ SUBROUTINE DSTEVX( JOBZ, RANGE, N, D, E, VL, VU, IL, IU, ABSTOL,
439438
30 CONTINUE
440439
*
441440
IF( I.NE.0 ) THEN
442-
ITMP1 = IWORK( INDIBL+I-1 )
441+
ITMP1 = IWORK( 1 + I-1 )
443442
W( I ) = W( J )
444-
IWORK( INDIBL+I-1 ) = IWORK( INDIBL+J-1 )
443+
IWORK( 1 + I-1 ) = IWORK( 1 + J-1 )
445444
W( J ) = TMP1
446-
IWORK( INDIBL+J-1 ) = ITMP1
445+
IWORK( 1 + J-1 ) = ITMP1
447446
CALL DSWAP( N, Z( 1, I ), 1, Z( 1, J ), 1 )
448447
IF( INFO.NE.0 ) THEN
449448
ITMP1 = IFAIL( I )

SRC/ssbgvx.f

Lines changed: 7 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -317,7 +317,7 @@ SUBROUTINE SSBGVX( JOBZ, RANGE, UPLO, N, KA, KB, AB, LDAB, BB,
317317
* .. Local Scalars ..
318318
LOGICAL ALLEIG, INDEIG, TEST, UPPER, VALEIG, WANTZ
319319
CHARACTER ORDER, VECT
320-
INTEGER I, IINFO, INDD, INDE, INDEE, INDIBL, INDISP,
320+
INTEGER I, IINFO, INDD, INDE, INDEE, INDISP,
321321
$ INDIWO, INDWRK, ITMP1, J, JJ, NSPLIT
322322
REAL TMP1
323323
* ..
@@ -457,17 +457,16 @@ SUBROUTINE SSBGVX( JOBZ, RANGE, UPLO, N, KA, KB, AB, LDAB, BB,
457457
ELSE
458458
ORDER = 'E'
459459
END IF
460-
INDIBL = 1
461-
INDISP = INDIBL + N
460+
INDISP = 1 + N
462461
INDIWO = INDISP + N
463462
CALL SSTEBZ( RANGE, ORDER, N, VL, VU, IL, IU, ABSTOL,
464463
$ WORK( INDD ), WORK( INDE ), M, NSPLIT, W,
465-
$ IWORK( INDIBL ), IWORK( INDISP ), WORK( INDWRK ),
464+
$ IWORK( 1 ), IWORK( INDISP ), WORK( INDWRK ),
466465
$ IWORK( INDIWO ), INFO )
467466
*
468467
IF( WANTZ ) THEN
469468
CALL SSTEIN( N, WORK( INDD ), WORK( INDE ), M, W,
470-
$ IWORK( INDIBL ), IWORK( INDISP ), Z, LDZ,
469+
$ IWORK( 1 ), IWORK( INDISP ), Z, LDZ,
471470
$ WORK( INDWRK ), IWORK( INDIWO ), IFAIL, INFO )
472471
*
473472
* Apply transformation matrix used in reduction to tridiagonal
@@ -497,11 +496,11 @@ SUBROUTINE SSBGVX( JOBZ, RANGE, UPLO, N, KA, KB, AB, LDAB, BB,
497496
40 CONTINUE
498497
*
499498
IF( I.NE.0 ) THEN
500-
ITMP1 = IWORK( INDIBL+I-1 )
499+
ITMP1 = IWORK( 1 + I-1 )
501500
W( I ) = W( J )
502-
IWORK( INDIBL+I-1 ) = IWORK( INDIBL+J-1 )
501+
IWORK( 1 + I-1 ) = IWORK( 1 + J-1 )
503502
W( J ) = TMP1
504-
IWORK( INDIBL+J-1 ) = ITMP1
503+
IWORK( 1 + J-1 ) = ITMP1
505504
CALL SSWAP( N, Z( 1, I ), 1, Z( 1, J ), 1 )
506505
IF( INFO.NE.0 ) THEN
507506
ITMP1 = IFAIL( I )

SRC/sspevx.f

Lines changed: 7 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -255,7 +255,7 @@ SUBROUTINE SSPEVX( JOBZ, RANGE, UPLO, N, AP, VL, VU, IL, IU,
255255
* .. Local Scalars ..
256256
LOGICAL ALLEIG, INDEIG, TEST, VALEIG, WANTZ
257257
CHARACTER ORDER
258-
INTEGER I, IINFO, IMAX, INDD, INDE, INDEE, INDIBL,
258+
INTEGER I, IINFO, IMAX, INDD, INDE, INDEE,
259259
$ INDISP, INDIWO, INDTAU, INDWRK, ISCALE, ITMP1,
260260
$ J, JJ, NSPLIT
261261
REAL ABSTLL, ANRM, BIGNUM, EPS, RMAX, RMIN, SAFMIN,
@@ -424,17 +424,16 @@ SUBROUTINE SSPEVX( JOBZ, RANGE, UPLO, N, AP, VL, VU, IL, IU,
424424
ELSE
425425
ORDER = 'E'
426426
END IF
427-
INDIBL = 1
428-
INDISP = INDIBL + N
427+
INDISP = 1 + N
429428
INDIWO = INDISP + N
430429
CALL SSTEBZ( RANGE, ORDER, N, VLL, VUU, IL, IU, ABSTLL,
431430
$ WORK( INDD ), WORK( INDE ), M, NSPLIT, W,
432-
$ IWORK( INDIBL ), IWORK( INDISP ), WORK( INDWRK ),
431+
$ IWORK( 1 ), IWORK( INDISP ), WORK( INDWRK ),
433432
$ IWORK( INDIWO ), INFO )
434433
*
435434
IF( WANTZ ) THEN
436435
CALL SSTEIN( N, WORK( INDD ), WORK( INDE ), M, W,
437-
$ IWORK( INDIBL ), IWORK( INDISP ), Z, LDZ,
436+
$ IWORK( 1 ), IWORK( INDISP ), Z, LDZ,
438437
$ WORK( INDWRK ), IWORK( INDIWO ), IFAIL, INFO )
439438
*
440439
* Apply orthogonal matrix used in reduction to tridiagonal
@@ -471,11 +470,11 @@ SUBROUTINE SSPEVX( JOBZ, RANGE, UPLO, N, AP, VL, VU, IL, IU,
471470
30 CONTINUE
472471
*
473472
IF( I.NE.0 ) THEN
474-
ITMP1 = IWORK( INDIBL+I-1 )
473+
ITMP1 = IWORK( 1 + I-1 )
475474
W( I ) = W( J )
476-
IWORK( INDIBL+I-1 ) = IWORK( INDIBL+J-1 )
475+
IWORK( 1 + I-1 ) = IWORK( 1 + J-1 )
477476
W( J ) = TMP1
478-
IWORK( INDIBL+J-1 ) = ITMP1
477+
IWORK( 1 + J-1 ) = ITMP1
479478
CALL SSWAP( N, Z( 1, I ), 1, Z( 1, J ), 1 )
480479
IF( INFO.NE.0 ) THEN
481480
ITMP1 = IFAIL( I )

SRC/sstevx.f

Lines changed: 7 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -248,7 +248,7 @@ SUBROUTINE SSTEVX( JOBZ, RANGE, N, D, E, VL, VU, IL, IU, ABSTOL,
248248
* .. Local Scalars ..
249249
LOGICAL ALLEIG, INDEIG, TEST, VALEIG, WANTZ
250250
CHARACTER ORDER
251-
INTEGER I, IMAX, INDIBL, INDISP, INDIWO, INDWRK,
251+
INTEGER I, IMAX, INDISP, INDIWO, INDWRK,
252252
$ ISCALE, ITMP1, J, JJ, NSPLIT
253253
REAL BIGNUM, EPS, RMAX, RMIN, SAFMIN, SIGMA, SMLNUM,
254254
$ TMP1, TNRM, VLL, VUU
@@ -399,15 +399,14 @@ SUBROUTINE SSTEVX( JOBZ, RANGE, N, D, E, VL, VU, IL, IU, ABSTOL,
399399
ORDER = 'E'
400400
END IF
401401
INDWRK = 1
402-
INDIBL = 1
403-
INDISP = INDIBL + N
402+
INDISP = 1 + N
404403
INDIWO = INDISP + N
405404
CALL SSTEBZ( RANGE, ORDER, N, VLL, VUU, IL, IU, ABSTOL, D, E, M,
406-
$ NSPLIT, W, IWORK( INDIBL ), IWORK( INDISP ),
405+
$ NSPLIT, W, IWORK( 1 ), IWORK( INDISP ),
407406
$ WORK( INDWRK ), IWORK( INDIWO ), INFO )
408407
*
409408
IF( WANTZ ) THEN
410-
CALL SSTEIN( N, D, E, M, W, IWORK( INDIBL ), IWORK( INDISP ),
409+
CALL SSTEIN( N, D, E, M, W, IWORK( 1 ), IWORK( INDISP ),
411410
$ Z, LDZ, WORK( INDWRK ), IWORK( INDIWO ), IFAIL,
412411
$ INFO )
413412
END IF
@@ -439,11 +438,11 @@ SUBROUTINE SSTEVX( JOBZ, RANGE, N, D, E, VL, VU, IL, IU, ABSTOL,
439438
30 CONTINUE
440439
*
441440
IF( I.NE.0 ) THEN
442-
ITMP1 = IWORK( INDIBL+I-1 )
441+
ITMP1 = IWORK( 1 + I-1 )
443442
W( I ) = W( J )
444-
IWORK( INDIBL+I-1 ) = IWORK( INDIBL+J-1 )
443+
IWORK( 1 + I-1 ) = IWORK( 1 + J-1 )
445444
W( J ) = TMP1
446-
IWORK( INDIBL+J-1 ) = ITMP1
445+
IWORK( 1 + J-1 ) = ITMP1
447446
CALL SSWAP( N, Z( 1, I ), 1, Z( 1, J ), 1 )
448447
IF( INFO.NE.0 ) THEN
449448
ITMP1 = IFAIL( I )

0 commit comments

Comments
 (0)