serrst.f
来自「famous linear algebra library (LAPACK) p」· F 代码 · 共 887 行 · 第 1/3 页
F
887 行
INFOT = 8
CALL SSYEVX( 'N', 'V', 'U', 1, A, 1, 0.0, 0.0, 0, 0, 0.0, M, X,
$ Z, 1, W, 8, IW, I3, INFO )
CALL CHKXER( 'SSYEVX', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL SSYEVX( 'N', 'I', 'U', 1, A, 1, 0.0, 0.0, 0, 0, 0.0, M, X,
$ Z, 1, W, 8, IW, I3, INFO )
CALL CHKXER( 'SSYEVX', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL SSYEVX( 'N', 'I', 'U', 1, A, 1, 0.0, 0.0, 2, 1, 0.0, M, X,
$ Z, 1, W, 8, IW, I3, INFO )
CALL CHKXER( 'SSYEVX', INFOT, NOUT, LERR, OK )
INFOT = 10
CALL SSYEVX( 'N', 'I', 'U', 2, A, 2, 0.0, 0.0, 2, 1, 0.0, M, X,
$ Z, 1, W, 16, IW, I3, INFO )
CALL CHKXER( 'SSYEVX', INFOT, NOUT, LERR, OK )
INFOT = 10
CALL SSYEVX( 'N', 'I', 'U', 1, A, 1, 0.0, 0.0, 1, 2, 0.0, M, X,
$ Z, 1, W, 8, IW, I3, INFO )
CALL CHKXER( 'SSYEVX', INFOT, NOUT, LERR, OK )
INFOT = 15
CALL SSYEVX( 'V', 'A', 'U', 2, A, 2, 0.0, 0.0, 0, 0, 0.0, M, X,
$ Z, 1, W, 16, IW, I3, INFO )
CALL CHKXER( 'SSYEVX', INFOT, NOUT, LERR, OK )
INFOT = 17
CALL SSYEVX( 'V', 'A', 'U', 1, A, 1, 0.0, 0.0, 0, 0, 0.0, M, X,
$ Z, 1, W, 0, IW, I3, INFO )
CALL CHKXER( 'SSYEVX', INFOT, NOUT, LERR, OK )
NT = NT + 12
*
* SSPEVD
*
SRNAMT = 'SSPEVD'
INFOT = 1
CALL SSPEVD( '/', 'U', 0, A, X, Z, 1, W, 1, IW, 1, INFO )
CALL CHKXER( 'SSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SSPEVD( 'N', '/', 0, A, X, Z, 1, W, 1, IW, 1, INFO )
CALL CHKXER( 'SSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL SSPEVD( 'N', 'U', -1, A, X, Z, 1, W, 1, IW, 1, INFO )
CALL CHKXER( 'SSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL SSPEVD( 'V', 'U', 2, A, X, Z, 1, W, 23, IW, 12, INFO )
CALL CHKXER( 'SSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL SSPEVD( 'N', 'U', 1, A, X, Z, 1, W, 0, IW, 1, INFO )
CALL CHKXER( 'SSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL SSPEVD( 'N', 'U', 2, A, X, Z, 1, W, 3, IW, 1, INFO )
CALL CHKXER( 'SSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL SSPEVD( 'V', 'U', 2, A, X, Z, 2, W, 16, IW, 12, INFO )
CALL CHKXER( 'SSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL SSPEVD( 'N', 'U', 1, A, X, Z, 1, W, 1, IW, 0, INFO )
CALL CHKXER( 'SSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL SSPEVD( 'N', 'U', 2, A, X, Z, 1, W, 4, IW, 0, INFO )
CALL CHKXER( 'SSPEVD', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL SSPEVD( 'V', 'U', 2, A, X, Z, 2, W, 23, IW, 11, INFO )
CALL CHKXER( 'SSPEVD', INFOT, NOUT, LERR, OK )
NT = NT + 10
*
* SSPEV
*
SRNAMT = 'SSPEV '
INFOT = 1
CALL SSPEV( '/', 'U', 0, A, W, Z, 1, X, INFO )
CALL CHKXER( 'SSPEV ', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SSPEV( 'N', '/', 0, A, W, Z, 1, X, INFO )
CALL CHKXER( 'SSPEV ', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL SSPEV( 'N', 'U', -1, A, W, Z, 1, X, INFO )
CALL CHKXER( 'SSPEV ', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL SSPEV( 'V', 'U', 2, A, W, Z, 1, X, INFO )
CALL CHKXER( 'SSPEV ', INFOT, NOUT, LERR, OK )
NT = NT + 4
*
* SSPEVX
*
SRNAMT = 'SSPEVX'
INFOT = 1
CALL SSPEVX( '/', 'A', 'U', 0, A, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
$ 1, W, IW, I3, INFO )
CALL CHKXER( 'SSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SSPEVX( 'N', '/', 'U', 0, A, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
$ 1, W, IW, I3, INFO )
CALL CHKXER( 'SSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL SSPEVX( 'N', 'A', '/', 0, A, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
$ 1, W, IW, I3, INFO )
INFOT = 4
CALL SSPEVX( 'N', 'A', 'U', -1, A, 0.0, 0.0, 0, 0, 0.0, M, X,
$ Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'SSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL SSPEVX( 'N', 'V', 'U', 1, A, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
$ 1, W, IW, I3, INFO )
CALL CHKXER( 'SSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 8
CALL SSPEVX( 'N', 'I', 'U', 1, A, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
$ 1, W, IW, I3, INFO )
CALL CHKXER( 'SSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 8
CALL SSPEVX( 'N', 'I', 'U', 1, A, 0.0, 0.0, 2, 1, 0.0, M, X, Z,
$ 1, W, IW, I3, INFO )
CALL CHKXER( 'SSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL SSPEVX( 'N', 'I', 'U', 2, A, 0.0, 0.0, 2, 1, 0.0, M, X, Z,
$ 1, W, IW, I3, INFO )
CALL CHKXER( 'SSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL SSPEVX( 'N', 'I', 'U', 1, A, 0.0, 0.0, 1, 2, 0.0, M, X, Z,
$ 1, W, IW, I3, INFO )
CALL CHKXER( 'SSPEVX', INFOT, NOUT, LERR, OK )
INFOT = 14
CALL SSPEVX( 'V', 'A', 'U', 2, A, 0.0, 0.0, 0, 0, 0.0, M, X, Z,
$ 1, W, IW, I3, INFO )
CALL CHKXER( 'SSPEVX', INFOT, NOUT, LERR, OK )
NT = NT + 10
*
* Test error exits for the SB path.
*
ELSE IF( LSAMEN( 2, C2, 'SB' ) ) THEN
*
* SSBTRD
*
SRNAMT = 'SSBTRD'
INFOT = 1
CALL SSBTRD( '/', 'U', 0, 0, A, 1, D, E, Z, 1, W, INFO )
CALL CHKXER( 'SSBTRD', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SSBTRD( 'N', '/', 0, 0, A, 1, D, E, Z, 1, W, INFO )
CALL CHKXER( 'SSBTRD', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL SSBTRD( 'N', 'U', -1, 0, A, 1, D, E, Z, 1, W, INFO )
CALL CHKXER( 'SSBTRD', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL SSBTRD( 'N', 'U', 0, -1, A, 1, D, E, Z, 1, W, INFO )
CALL CHKXER( 'SSBTRD', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL SSBTRD( 'N', 'U', 1, 1, A, 1, D, E, Z, 1, W, INFO )
CALL CHKXER( 'SSBTRD', INFOT, NOUT, LERR, OK )
INFOT = 10
CALL SSBTRD( 'V', 'U', 2, 0, A, 1, D, E, Z, 1, W, INFO )
CALL CHKXER( 'SSBTRD', INFOT, NOUT, LERR, OK )
NT = NT + 6
*
* SSBEVD
*
SRNAMT = 'SSBEVD'
INFOT = 1
CALL SSBEVD( '/', 'U', 0, 0, A, 1, X, Z, 1, W, 1, IW, 1, INFO )
CALL CHKXER( 'SSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SSBEVD( 'N', '/', 0, 0, A, 1, X, Z, 1, W, 1, IW, 1, INFO )
CALL CHKXER( 'SSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL SSBEVD( 'N', 'U', -1, 0, A, 1, X, Z, 1, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'SSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL SSBEVD( 'N', 'U', 0, -1, A, 1, X, Z, 1, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'SSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL SSBEVD( 'N', 'U', 2, 1, A, 1, X, Z, 1, W, 4, IW, 1, INFO )
CALL CHKXER( 'SSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL SSBEVD( 'V', 'U', 2, 1, A, 2, X, Z, 1, W, 25, IW, 12,
$ INFO )
CALL CHKXER( 'SSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL SSBEVD( 'N', 'U', 1, 0, A, 1, X, Z, 1, W, 0, IW, 1, INFO )
CALL CHKXER( 'SSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL SSBEVD( 'N', 'U', 2, 0, A, 1, X, Z, 1, W, 3, IW, 1, INFO )
CALL CHKXER( 'SSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL SSBEVD( 'V', 'U', 2, 0, A, 1, X, Z, 2, W, 18, IW, 12,
$ INFO )
CALL CHKXER( 'SSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 13
CALL SSBEVD( 'N', 'U', 1, 0, A, 1, X, Z, 1, W, 1, IW, 0, INFO )
CALL CHKXER( 'SSBEVD', INFOT, NOUT, LERR, OK )
INFOT = 13
CALL SSBEVD( 'V', 'U', 2, 0, A, 1, X, Z, 2, W, 25, IW, 11,
$ INFO )
CALL CHKXER( 'SSBEVD', INFOT, NOUT, LERR, OK )
NT = NT + 11
*
* SSBEV
*
SRNAMT = 'SSBEV '
INFOT = 1
CALL SSBEV( '/', 'U', 0, 0, A, 1, X, Z, 1, W, INFO )
CALL CHKXER( 'SSBEV ', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SSBEV( 'N', '/', 0, 0, A, 1, X, Z, 1, W, INFO )
CALL CHKXER( 'SSBEV ', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL SSBEV( 'N', 'U', -1, 0, A, 1, X, Z, 1, W, INFO )
CALL CHKXER( 'SSBEV ', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL SSBEV( 'N', 'U', 0, -1, A, 1, X, Z, 1, W, INFO )
CALL CHKXER( 'SSBEV ', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL SSBEV( 'N', 'U', 2, 1, A, 1, X, Z, 1, W, INFO )
CALL CHKXER( 'SSBEV ', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL SSBEV( 'V', 'U', 2, 0, A, 1, X, Z, 1, W, INFO )
CALL CHKXER( 'SSBEV ', INFOT, NOUT, LERR, OK )
NT = NT + 6
*
* SSBEVX
*
SRNAMT = 'SSBEVX'
INFOT = 1
CALL SSBEVX( '/', 'A', 'U', 0, 0, A, 1, Q, 1, 0.0, 0.0, 0, 0,
$ 0.0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'SSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SSBEVX( 'N', '/', 'U', 0, 0, A, 1, Q, 1, 0.0, 0.0, 0, 0,
$ 0.0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'SSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL SSBEVX( 'N', 'A', '/', 0, 0, A, 1, Q, 1, 0.0, 0.0, 0, 0,
$ 0.0, M, X, Z, 1, W, IW, I3, INFO )
INFOT = 4
CALL SSBEVX( 'N', 'A', 'U', -1, 0, A, 1, Q, 1, 0.0, 0.0, 0, 0,
$ 0.0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'SSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL SSBEVX( 'N', 'A', 'U', 0, -1, A, 1, Q, 1, 0.0, 0.0, 0, 0,
$ 0.0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'SSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL SSBEVX( 'N', 'A', 'U', 2, 1, A, 1, Q, 1, 0.0, 0.0, 0, 0,
$ 0.0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'SSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL SSBEVX( 'V', 'A', 'U', 2, 0, A, 1, Q, 1, 0.0, 0.0, 0, 0,
$ 0.0, M, X, Z, 2, W, IW, I3, INFO )
CALL CHKXER( 'SSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL SSBEVX( 'N', 'V', 'U', 1, 0, A, 1, Q, 1, 0.0, 0.0, 0, 0,
$ 0.0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'SSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 12
CALL SSBEVX( 'N', 'I', 'U', 1, 0, A, 1, Q, 1, 0.0, 0.0, 0, 0,
$ 0.0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'SSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 12
CALL SSBEVX( 'N', 'I', 'U', 1, 0, A, 1, Q, 1, 0.0, 0.0, 2, 1,
$ 0.0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'SSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 13
CALL SSBEVX( 'N', 'I', 'U', 2, 0, A, 1, Q, 1, 0.0, 0.0, 2, 1,
$ 0.0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'SSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 13
CALL SSBEVX( 'N', 'I', 'U', 1, 0, A, 1, Q, 1, 0.0, 0.0, 1, 2,
$ 0.0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'SSBEVX', INFOT, NOUT, LERR, OK )
INFOT = 18
CALL SSBEVX( 'V', 'A', 'U', 2, 0, A, 1, Q, 2, 0.0, 0.0, 0, 0,
$ 0.0, M, X, Z, 1, W, IW, I3, INFO )
CALL CHKXER( 'SSBEVX', 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 SERRST
*
END
⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?