Browse Source

Handle corner cases of LWORK (Reference-LAPACK PR 942)

tags/v0.3.26
Martin Kroeker GitHub 2 years ago
parent
commit
8ce44c18a0
No known key found for this signature in database GPG Key ID: 4AEE18F83AFDEB23
5 changed files with 17 additions and 17 deletions
  1. +6
    -6
      lapack-netlib/TESTING/EIG/cerrst.f
  2. +1
    -1
      lapack-netlib/TESTING/EIG/chkxer.f
  3. +2
    -2
      lapack-netlib/TESTING/EIG/derrst.f
  4. +2
    -2
      lapack-netlib/TESTING/EIG/serrst.f
  5. +6
    -6
      lapack-netlib/TESTING/EIG/zerrst.f

+ 6
- 6
lapack-netlib/TESTING/EIG/cerrst.f View File

@@ -748,17 +748,17 @@
CALL CHKXER( 'CHEEVR', INFOT, NOUT, LERR, OK )
INFOT = 18
CALL CHEEVR( 'V', 'I', 'U', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
$ M, R, Z, 1, IW, Q, 2*N-1, RW, 24*N, IW( 2*N+1 ),
$ M, R, Z, 1, IW, Q, 0, RW, 24*N, IW( 2*N+1 ),
$ 10*N, INFO )
CALL CHKXER( 'CHEEVR', INFOT, NOUT, LERR, OK )
INFOT = 20
CALL CHEEVR( 'V', 'I', 'U', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
$ M, R, Z, 1, IW, Q, 2*N, RW, 24*N-1, IW( 2*N-1 ),
$ M, R, Z, 1, IW, Q, 2*N, RW, 0, IW( 2*N-1 ),
$ 10*N, INFO )
CALL CHKXER( 'CHEEVR', INFOT, NOUT, LERR, OK )
INFOT = 22
CALL CHEEVR( 'V', 'I', 'U', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
$ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW, 10*N-1,
$ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW, 0,
$ INFO )
CALL CHKXER( 'CHEEVR', INFOT, NOUT, LERR, OK )
NT = NT + 12
@@ -830,19 +830,19 @@
INFOT = 18
CALL CHEEVR_2STAGE( 'N', 'I', 'U', 1, A, 1,
$ 0.0, 0.0, 1, 1, 0.0,
$ M, R, Z, 1, IW, Q, 2*N-1, RW, 24*N, IW( 2*N+1 ),
$ M, R, Z, 1, IW, Q, 0, RW, 24*N, IW( 2*N+1 ),
$ 10*N, INFO )
CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK )
INFOT = 20
CALL CHEEVR_2STAGE( 'N', 'I', 'U', 1, A, 1,
$ 0.0, 0.0, 1, 1, 0.0,
$ M, R, Z, 1, IW, Q, 26*N, RW, 24*N-1, IW( 2*N-1 ),
$ M, R, Z, 1, IW, Q, 26*N, RW, 0, IW( 2*N-1 ),
$ 10*N, INFO )
CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK )
INFOT = 22
CALL CHEEVR_2STAGE( 'N', 'I', 'U', 1, A, 1,
$ 0.0, 0.0, 1, 1, 0.0,
$ M, R, Z, 1, IW, Q, 26*N, RW, 24*N, IW, 10*N-1,
$ M, R, Z, 1, IW, Q, 26*N, RW, 24*N, IW, 0,
$ INFO )
CALL CHKXER( 'CHEEVR_2STAGE', INFOT, NOUT, LERR, OK )
NT = NT + 13


+ 1
- 1
lapack-netlib/TESTING/EIG/chkxer.f View File

@@ -61,7 +61,7 @@
RETURN
*
9999 FORMAT( ' *** Illegal value of parameter number ', I2,
$ ' not detected by ', A6, ' ***' )
$ ' not detected by ', A, ' ***' )
*
* End of CHKXER
*


+ 2
- 2
lapack-netlib/TESTING/EIG/derrst.f View File

@@ -735,12 +735,12 @@
CALL CHKXER( 'DSYEVR', INFOT, NOUT, LERR, OK )
INFOT = 18
CALL DSYEVR( 'V', 'I', 'U', 1, A, 1, 0.0D0, 0.0D0, 1, 1, 0.0D0,
$ M, R, Z, 1, IW, Q, 26*N-1, IW( 2*N+1 ), 10*N,
$ M, R, Z, 1, IW, Q, 0, IW( 2*N+1 ), 10*N,
$ INFO )
CALL CHKXER( 'DSYEVR', INFOT, NOUT, LERR, OK )
INFOT = 20
CALL DSYEVR( 'V', 'I', 'U', 1, A, 1, 0.0D0, 0.0D0, 1, 1, 0.0D0,
$ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N-1,
$ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 0,
$ INFO )
CALL CHKXER( 'DSYEVR', INFOT, NOUT, LERR, OK )
NT = NT + 11


+ 2
- 2
lapack-netlib/TESTING/EIG/serrst.f View File

@@ -733,12 +733,12 @@
CALL CHKXER( 'SSYEVR', INFOT, NOUT, LERR, OK )
INFOT = 18
CALL SSYEVR( 'V', 'I', 'U', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
$ M, R, Z, 1, IW, Q, 26*N-1, IW( 2*N+1 ), 10*N,
$ M, R, Z, 1, IW, Q, 0, IW( 2*N+1 ), 10*N,
$ INFO )
CALL CHKXER( 'SSYEVR', INFOT, NOUT, LERR, OK )
INFOT = 20
CALL SSYEVR( 'V', 'I', 'U', 1, A, 1, 0.0E0, 0.0E0, 1, 1, 0.0,
$ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 10*N-1,
$ M, R, Z, 1, IW, Q, 26*N, IW( 2*N+1 ), 0,
$ INFO )
CALL CHKXER( 'SSYEVR', INFOT, NOUT, LERR, OK )
NT = NT + 11


+ 6
- 6
lapack-netlib/TESTING/EIG/zerrst.f View File

@@ -748,17 +748,17 @@
CALL CHKXER( 'ZHEEVR', INFOT, NOUT, LERR, OK )
INFOT = 18
CALL ZHEEVR( 'V', 'I', 'U', 1, A, 1, 0.0D0, 0.0D0, 1, 1, 0.0D0,
$ M, R, Z, 1, IW, Q, 2*N-1, RW, 24*N, IW( 2*N+1 ),
$ M, R, Z, 1, IW, Q, 0, RW, 24*N, IW( 2*N+1 ),
$ 10*N, INFO )
CALL CHKXER( 'ZHEEVR', INFOT, NOUT, LERR, OK )
INFOT = 20
CALL ZHEEVR( 'V', 'I', 'U', 1, A, 1, 0.0D0, 0.0D0, 1, 1, 0.0D0,
$ M, R, Z, 1, IW, Q, 2*N, RW, 24*N-1, IW( 2*N-1 ),
$ M, R, Z, 1, IW, Q, 2*N, RW, 0, IW( 2*N-1 ),
$ 10*N, INFO )
CALL CHKXER( 'ZHEEVR', INFOT, NOUT, LERR, OK )
INFOT = 22
CALL ZHEEVR( 'V', 'I', 'U', 1, A, 1, 0.0D0, 0.0D0, 1, 1, 0.0D0,
$ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW, 10*N-1,
$ M, R, Z, 1, IW, Q, 2*N, RW, 24*N, IW, 0,
$ INFO )
CALL CHKXER( 'ZHEEVR', INFOT, NOUT, LERR, OK )
NT = NT + 12
@@ -830,19 +830,19 @@
INFOT = 18
CALL ZHEEVR_2STAGE( 'N', 'I', 'U', 1, A, 1,
$ 0.0D0, 0.0D0, 1, 1, 0.0D0,
$ M, R, Z, 1, IW, Q, 2*N-1, RW, 24*N, IW( 2*N+1 ),
$ M, R, Z, 1, IW, Q, 0, RW, 24*N, IW( 2*N+1 ),
$ 10*N, INFO )
CALL CHKXER( 'ZHEEVR_2STAGE', INFOT, NOUT, LERR, OK )
INFOT = 20
CALL ZHEEVR_2STAGE( 'N', 'I', 'U', 1, A, 1,
$ 0.0D0, 0.0D0, 1, 1, 0.0D0,
$ M, R, Z, 1, IW, Q, 26*N, RW, 24*N-1, IW( 2*N-1 ),
$ M, R, Z, 1, IW, Q, 26*N, RW, 0, IW( 2*N-1 ),
$ 10*N, INFO )
CALL CHKXER( 'ZHEEVR_2STAGE', INFOT, NOUT, LERR, OK )
INFOT = 22
CALL ZHEEVR_2STAGE( 'N', 'I', 'U', 1, A, 1,
$ 0.0D0, 0.0D0, 1, 1, 0.0D0,
$ M, R, Z, 1, IW, Q, 26*N, RW, 24*N, IW, 10*N-1,
$ M, R, Z, 1, IW, Q, 26*N, RW, 24*N, IW, 0,
$ INFO )
CALL CHKXER( 'ZHEEVR_2STAGE', INFOT, NOUT, LERR, OK )
NT = NT + 13


Loading…
Cancel
Save