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