derrst.f
来自「famous linear algebra library (LAPACK) p」· F 代码 · 共 889 行 · 第 1/3 页
F
889 行
INFOT = 8
CALL DSYEVX( 'N', 'V', 'U', 1, A, 1, 0.0D0, 0.0D0, 0, 0, 0.0D0,
$ M, X, Z, 1, W, 8, IW, I3, INFO )
CALL CHKXER( 'DSYEVX', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DSYEVX( 'N', 'I', 'U', 1, A, 1, 0.0D0, 0.0D0, 0, 0, 0.0D0,
$ M, X, Z, 1, W, 8, IW, I3, INFO )
CALL CHKXER( 'DSYEVX', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DSYEVX( 'N', 'I', 'U', 1, A, 1, 0.0D0, 0.0D0, 2, 1, 0.0D0,
$ M, X, Z, 1, W, 8, IW, I3, INFO )
CALL CHKXER( 'DSYEVX', INFOT, NOUT, LERR, OK )
INFOT = 10
CALL DSYEVX( 'N', 'I', 'U', 2, A, 2, 0.0D0, 0.0D0, 2, 1, 0.0D0,
$ M, X, Z, 1, W, 16, IW, I3, INFO )
CALL CHKXER( 'DSYEVX', INFOT, NOUT, LERR, OK )
INFOT = 10
CALL DSYEVX( 'N', 'I', 'U', 1, A, 1, 0.0D0, 0.0D0, 1, 2, 0.0D0,
$ M, X, Z, 1, W, 8, IW, I3, INFO )
CALL CHKXER( 'DSYEVX', INFOT, NOUT, LERR, OK )
INFOT = 15
CALL DSYEVX( 'V', 'A', 'U', 2, A, 2, 0.0D0, 0.0D0, 0, 0, 0.0D0,
$ M, X, Z, 1, W, 16, IW, I3, INFO )
CALL CHKXER( 'DSYEVX', INFOT, NOUT, LERR, OK )
INFOT = 17
CALL DSYEVX( 'V', 'A', 'U', 1, A, 1, 0.0D0, 0.0D0, 0, 0, 0.0D0,
$ M, X, Z, 1, W, 0, IW, I3, INFO )
CALL CHKXER( 'DSYEVX', INFOT, NOUT, LERR, OK )
NT = NT + 12
*
* DSPEVD
*
SRNAMT = 'DSPEVD'
INFOT = 1
CALL DSPEVD( '/', 'U', 0, A, X, Z, 1, W, 1, IW, 1, INFO )
CALL CHKXER( 'DSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DSPEVD( 'N', '/', 0, A, X, Z, 1, W, 1, IW, 1, INFO )
CALL CHKXER( 'DSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DSPEVD( 'N', 'U', -1, A, X, Z, 1, W, 1, IW, 1, INFO )
CALL CHKXER( 'DSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL DSPEVD( 'V', 'U', 2, A, X, Z, 1, W, 23, IW, 12, INFO )
CALL CHKXER( 'DSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DSPEVD( 'N', 'U', 1, A, X, Z, 1, W, 0, IW, 1, INFO )
CALL CHKXER( 'DSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DSPEVD( 'N', 'U', 2, A, X, Z, 1, W, 3, IW, 1, INFO )
CALL CHKXER( 'DSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DSPEVD( 'V', 'U', 2, A, X, Z, 2, W, 16, IW, 12, INFO )
CALL CHKXER( 'DSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL DSPEVD( 'N', 'U', 1, A, X, Z, 1, W, 1, IW, 0, INFO )
CALL CHKXER( 'DSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL DSPEVD( 'N', 'U', 2, A, X, Z, 1, W, 4, IW, 0, INFO )
CALL CHKXER( 'DSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL DSPEVD( 'V', 'U', 2, A, X, Z, 2, W, 23, IW, 11, INFO )
CALL CHKXER( 'DSPEVD', INFOT, NOUT, LERR, OK )
NT = NT + 10
*
* DSPEV
*
SRNAMT = 'DSPEV '
INFOT = 1
CALL DSPEV( '/', 'U', 0, A, W, Z, 1, X, INFO )
CALL CHKXER( 'DSPEV ', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DSPEV( 'N', '/', 0, A, W, Z, 1, X, INFO )
CALL CHKXER( 'DSPEV ', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DSPEV( 'N', 'U', -1, A, W, Z, 1, X, INFO )
CALL CHKXER( 'DSPEV ', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL DSPEV( 'V', 'U', 2, A, W, Z, 1, X, INFO )
CALL CHKXER( 'DSPEV ', INFOT, NOUT, LERR, OK )
NT = NT + 4
*
* DSPEVX
*
SRNAMT = 'DSPEVX'
INFOT = 1
CALL DSPEVX( '/', 'A', 'U', 0, A, 0.0D0, 0.0D0, 0, 0, 0.0D0, M,
$ X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DSPEVX( 'N', '/', 'U', 0, A, 0.0D0, 0.0D0, 0, 0, 0.0D0, M,
$ X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DSPEVX( 'N', 'A', '/', 0, A, 0.0D0, 0.0D0, 0, 0, 0.0D0, M,
$ X, Z, 1, W, IW, I3, INFO )
INFOT = 4
CALL DSPEVX( 'N', 'A', 'U', -1, A, 0.0D0, 0.0D0, 0, 0, 0.0D0,
$ M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL DSPEVX( 'N', 'V', 'U', 1, A, 0.0D0, 0.0D0, 0, 0, 0.0D0, M,
$ X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 8
CALL DSPEVX( 'N', 'I', 'U', 1, A, 0.0D0, 0.0D0, 0, 0, 0.0D0, M,
$ X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 8
CALL DSPEVX( 'N', 'I', 'U', 1, A, 0.0D0, 0.0D0, 2, 1, 0.0D0, M,
$ X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DSPEVX( 'N', 'I', 'U', 2, A, 0.0D0, 0.0D0, 2, 1, 0.0D0, M,
$ X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DSPEVX( 'N', 'I', 'U', 1, A, 0.0D0, 0.0D0, 1, 2, 0.0D0, M,
$ X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 14
CALL DSPEVX( 'V', 'A', 'U', 2, A, 0.0D0, 0.0D0, 0, 0, 0.0D0, M,
$ X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSPEVX', INFOT, NOUT, LERR, OK )
NT = NT + 10
*
* Test error exits for the SB path.
*
ELSE IF( LSAMEN( 2, C2, 'SB' ) ) THEN
*
* DSBTRD
*
SRNAMT = 'DSBTRD'
INFOT = 1
CALL DSBTRD( '/', 'U', 0, 0, A, 1, D, E, Z, 1, W, INFO )
CALL CHKXER( 'DSBTRD', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DSBTRD( 'N', '/', 0, 0, A, 1, D, E, Z, 1, W, INFO )
CALL CHKXER( 'DSBTRD', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DSBTRD( 'N', 'U', -1, 0, A, 1, D, E, Z, 1, W, INFO )
CALL CHKXER( 'DSBTRD', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL DSBTRD( 'N', 'U', 0, -1, A, 1, D, E, Z, 1, W, INFO )
CALL CHKXER( 'DSBTRD', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL DSBTRD( 'N', 'U', 1, 1, A, 1, D, E, Z, 1, W, INFO )
CALL CHKXER( 'DSBTRD', INFOT, NOUT, LERR, OK )
INFOT = 10
CALL DSBTRD( 'V', 'U', 2, 0, A, 1, D, E, Z, 1, W, INFO )
CALL CHKXER( 'DSBTRD', INFOT, NOUT, LERR, OK )
NT = NT + 6
*
* DSBEVD
*
SRNAMT = 'DSBEVD'
INFOT = 1
CALL DSBEVD( '/', 'U', 0, 0, A, 1, X, Z, 1, W, 1, IW, 1, INFO )
CALL CHKXER( 'DSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DSBEVD( 'N', '/', 0, 0, A, 1, X, Z, 1, W, 1, IW, 1, INFO )
CALL CHKXER( 'DSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DSBEVD( 'N', 'U', -1, 0, A, 1, X, Z, 1, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'DSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL DSBEVD( 'N', 'U', 0, -1, A, 1, X, Z, 1, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'DSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL DSBEVD( 'N', 'U', 2, 1, A, 1, X, Z, 1, W, 4, IW, 1, INFO )
CALL CHKXER( 'DSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DSBEVD( 'V', 'U', 2, 1, A, 2, X, Z, 1, W, 25, IW, 12,
$ INFO )
CALL CHKXER( 'DSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL DSBEVD( 'N', 'U', 1, 0, A, 1, X, Z, 1, W, 0, IW, 1, INFO )
CALL CHKXER( 'DSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL DSBEVD( 'N', 'U', 2, 0, A, 1, X, Z, 1, W, 3, IW, 1, INFO )
CALL CHKXER( 'DSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL DSBEVD( 'V', 'U', 2, 0, A, 1, X, Z, 2, W, 18, IW, 12,
$ INFO )
CALL CHKXER( 'DSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 13
CALL DSBEVD( 'N', 'U', 1, 0, A, 1, X, Z, 1, W, 1, IW, 0, INFO )
CALL CHKXER( 'DSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 13
CALL DSBEVD( 'V', 'U', 2, 0, A, 1, X, Z, 2, W, 25, IW, 11,
$ INFO )
CALL CHKXER( 'DSBEVD', INFOT, NOUT, LERR, OK )
NT = NT + 11
*
* DSBEV
*
SRNAMT = 'DSBEV '
INFOT = 1
CALL DSBEV( '/', 'U', 0, 0, A, 1, X, Z, 1, W, INFO )
CALL CHKXER( 'DSBEV ', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DSBEV( 'N', '/', 0, 0, A, 1, X, Z, 1, W, INFO )
CALL CHKXER( 'DSBEV ', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DSBEV( 'N', 'U', -1, 0, A, 1, X, Z, 1, W, INFO )
CALL CHKXER( 'DSBEV ', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL DSBEV( 'N', 'U', 0, -1, A, 1, X, Z, 1, W, INFO )
CALL CHKXER( 'DSBEV ', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL DSBEV( 'N', 'U', 2, 1, A, 1, X, Z, 1, W, INFO )
CALL CHKXER( 'DSBEV ', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DSBEV( 'V', 'U', 2, 0, A, 1, X, Z, 1, W, INFO )
CALL CHKXER( 'DSBEV ', INFOT, NOUT, LERR, OK )
NT = NT + 6
*
* DSBEVX
*
SRNAMT = 'DSBEVX'
INFOT = 1
CALL DSBEVX( '/', 'A', 'U', 0, 0, A, 1, Q, 1, 0.0D0, 0.0D0, 0,
$ 0, 0.0D0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DSBEVX( 'N', '/', 'U', 0, 0, A, 1, Q, 1, 0.0D0, 0.0D0, 0,
$ 0, 0.0D0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DSBEVX( 'N', 'A', '/', 0, 0, A, 1, Q, 1, 0.0D0, 0.0D0, 0,
$ 0, 0.0D0, M, X, Z, 1, W, IW, I3, INFO )
INFOT = 4
CALL DSBEVX( 'N', 'A', 'U', -1, 0, A, 1, Q, 1, 0.0D0, 0.0D0, 0,
$ 0, 0.0D0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL DSBEVX( 'N', 'A', 'U', 0, -1, A, 1, Q, 1, 0.0D0, 0.0D0, 0,
$ 0, 0.0D0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL DSBEVX( 'N', 'A', 'U', 2, 1, A, 1, Q, 1, 0.0D0, 0.0D0, 0,
$ 0, 0.0D0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DSBEVX( 'V', 'A', 'U', 2, 0, A, 1, Q, 1, 0.0D0, 0.0D0, 0,
$ 0, 0.0D0, M, X, Z, 2, W, IW, I3, INFO )
CALL CHKXER( 'DSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL DSBEVX( 'N', 'V', 'U', 1, 0, A, 1, Q, 1, 0.0D0, 0.0D0, 0,
$ 0, 0.0D0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 12
CALL DSBEVX( 'N', 'I', 'U', 1, 0, A, 1, Q, 1, 0.0D0, 0.0D0, 0,
$ 0, 0.0D0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 12
CALL DSBEVX( 'N', 'I', 'U', 1, 0, A, 1, Q, 1, 0.0D0, 0.0D0, 2,
$ 1, 0.0D0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 13
CALL DSBEVX( 'N', 'I', 'U', 2, 0, A, 1, Q, 1, 0.0D0, 0.0D0, 2,
$ 1, 0.0D0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 13
CALL DSBEVX( 'N', 'I', 'U', 1, 0, A, 1, Q, 1, 0.0D0, 0.0D0, 1,
$ 2, 0.0D0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 18
CALL DSBEVX( 'V', 'A', 'U', 2, 0, A, 1, Q, 2, 0.0D0, 0.0D0, 0,
$ 0, 0.0D0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'DSBEVX', INFOT, NOUT, LERR, OK )
NT = NT + 13
END IF
*
* Print a summary line.
*
IF( OK ) THEN
WRITE( NOUT, FMT = 9999 )PATH, NT
ELSE
WRITE( NOUT, FMT = 9998 )PATH
END IF
*
9999 FORMAT( 1X, A3, ' routines passed the tests of the error exits',
$ ' (', I3, ' tests done)' )
9998 FORMAT( ' *** ', A3, ' routines failed the tests of the error ',
$ 'exits ***' )
*
RETURN
*
* End of DERRST
*
END
⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?