serrvx.f
来自「famous linear algebra library (LAPACK) p」· F 代码 · 共 653 行 · 第 1/2 页
F
653 行
SUBROUTINE SERRVX( PATH, NUNIT )
*
* -- LAPACK test routine (version 3.1.1) --
* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
* January 2007
*
* .. Scalar Arguments ..
CHARACTER*3 PATH
INTEGER NUNIT
* ..
*
* Purpose
* =======
*
* SERRVX tests the error exits for the REAL driver routines
* for solving linear systems of equations.
*
* Arguments
* =========
*
* PATH (input) CHARACTER*3
* The LAPACK path name for the routines to be tested.
*
* NUNIT (input) INTEGER
* The unit number for output.
*
* =====================================================================
*
* .. Parameters ..
INTEGER NMAX
PARAMETER ( NMAX = 4 )
* ..
* .. Local Scalars ..
CHARACTER EQ
CHARACTER*2 C2
INTEGER I, INFO, J
REAL RCOND
* ..
* .. Local Arrays ..
INTEGER IP( NMAX ), IW( NMAX )
REAL A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
$ C( NMAX ), R( NMAX ), R1( NMAX ), R2( NMAX ),
$ W( 2*NMAX ), X( NMAX )
* ..
* .. External Functions ..
LOGICAL LSAMEN
EXTERNAL LSAMEN
* ..
* .. External Subroutines ..
EXTERNAL CHKXER, SGBSV, SGBSVX, SGESV, SGESVX, SGTSV,
$ SGTSVX, SPBSV, SPBSVX, SPOSV, SPOSVX, SPPSV,
$ SPPSVX, SPTSV, SPTSVX, SSPSV, SSPSVX, SSYSV,
$ SSYSVX
* ..
* .. Scalars in Common ..
LOGICAL LERR, OK
CHARACTER*6 SRNAMT
INTEGER INFOT, NOUT
* ..
* .. Common blocks ..
COMMON / INFOC / INFOT, NOUT, OK, LERR
COMMON / SRNAMC / SRNAMT
* ..
* .. Intrinsic Functions ..
INTRINSIC REAL
* ..
* .. Executable Statements ..
*
NOUT = NUNIT
WRITE( NOUT, FMT = * )
C2 = PATH( 2: 3 )
*
* Set the variables to innocuous values.
*
DO 20 J = 1, NMAX
DO 10 I = 1, NMAX
A( I, J ) = 1. / REAL( I+J )
AF( I, J ) = 1. / REAL( I+J )
10 CONTINUE
B( J ) = 0.
R1( J ) = 0.
R2( J ) = 0.
W( J ) = 0.
X( J ) = 0.
C( J ) = 0.
R( J ) = 0.
IP( J ) = J
20 CONTINUE
EQ = ' '
OK = .TRUE.
*
IF( LSAMEN( 2, C2, 'GE' ) ) THEN
*
* SGESV
*
SRNAMT = 'SGESV '
INFOT = 1
CALL SGESV( -1, 0, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'SGESV ', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SGESV( 0, -1, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'SGESV ', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL SGESV( 2, 1, A, 1, IP, B, 2, INFO )
CALL CHKXER( 'SGESV ', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL SGESV( 2, 1, A, 2, IP, B, 1, INFO )
CALL CHKXER( 'SGESV ', INFOT, NOUT, LERR, OK )
*
* SGESVX
*
SRNAMT = 'SGESVX'
INFOT = 1
CALL SGESVX( '/', 'N', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGESVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SGESVX( 'N', '/', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGESVX', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL SGESVX( 'N', 'N', -1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGESVX', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL SGESVX( 'N', 'N', 0, -1, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGESVX', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL SGESVX( 'N', 'N', 2, 1, A, 1, AF, 2, IP, EQ, R, C, B, 2,
$ X, 2, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGESVX', INFOT, NOUT, LERR, OK )
INFOT = 8
CALL SGESVX( 'N', 'N', 2, 1, A, 2, AF, 1, IP, EQ, R, C, B, 2,
$ X, 2, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGESVX', INFOT, NOUT, LERR, OK )
INFOT = 10
EQ = '/'
CALL SGESVX( 'F', 'N', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGESVX', INFOT, NOUT, LERR, OK )
INFOT = 11
EQ = 'R'
CALL SGESVX( 'F', 'N', 1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGESVX', INFOT, NOUT, LERR, OK )
INFOT = 12
EQ = 'C'
CALL SGESVX( 'F', 'N', 1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGESVX', INFOT, NOUT, LERR, OK )
INFOT = 14
CALL SGESVX( 'N', 'N', 2, 1, A, 2, AF, 2, IP, EQ, R, C, B, 1,
$ X, 2, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGESVX', INFOT, NOUT, LERR, OK )
INFOT = 16
CALL SGESVX( 'N', 'N', 2, 1, A, 2, AF, 2, IP, EQ, R, C, B, 2,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGESVX', INFOT, NOUT, LERR, OK )
*
ELSE IF( LSAMEN( 2, C2, 'GB' ) ) THEN
*
* SGBSV
*
SRNAMT = 'SGBSV '
INFOT = 1
CALL SGBSV( -1, 0, 0, 0, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'SGBSV ', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SGBSV( 1, -1, 0, 0, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'SGBSV ', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL SGBSV( 1, 0, -1, 0, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'SGBSV ', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL SGBSV( 0, 0, 0, -1, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'SGBSV ', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL SGBSV( 1, 1, 1, 0, A, 3, IP, B, 1, INFO )
CALL CHKXER( 'SGBSV ', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL SGBSV( 2, 0, 0, 0, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'SGBSV ', INFOT, NOUT, LERR, OK )
*
* SGBSVX
*
SRNAMT = 'SGBSVX'
INFOT = 1
CALL SGBSVX( '/', 'N', 0, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
$ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SGBSVX( 'N', '/', 0, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
$ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL SGBSVX( 'N', 'N', -1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
$ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL SGBSVX( 'N', 'N', 1, -1, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
$ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL SGBSVX( 'N', 'N', 1, 0, -1, 0, A, 1, AF, 1, IP, EQ, R, C,
$ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL SGBSVX( 'N', 'N', 0, 0, 0, -1, A, 1, AF, 1, IP, EQ, R, C,
$ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 8
CALL SGBSVX( 'N', 'N', 1, 1, 1, 0, A, 2, AF, 4, IP, EQ, R, C,
$ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 10
CALL SGBSVX( 'N', 'N', 1, 1, 1, 0, A, 3, AF, 3, IP, EQ, R, C,
$ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 12
EQ = '/'
CALL SGBSVX( 'F', 'N', 0, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
$ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 13
EQ = 'R'
CALL SGBSVX( 'F', 'N', 1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
$ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 14
EQ = 'C'
CALL SGBSVX( 'F', 'N', 1, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
$ B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 16
CALL SGBSVX( 'N', 'N', 2, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
$ B, 1, X, 2, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 18
CALL SGBSVX( 'N', 'N', 2, 0, 0, 0, A, 1, AF, 1, IP, EQ, R, C,
$ B, 2, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGBSVX', INFOT, NOUT, LERR, OK )
*
ELSE IF( LSAMEN( 2, C2, 'GT' ) ) THEN
*
* SGTSV
*
SRNAMT = 'SGTSV '
INFOT = 1
CALL SGTSV( -1, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1,
$ INFO )
CALL CHKXER( 'SGTSV ', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SGTSV( 0, -1, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1,
$ INFO )
CALL CHKXER( 'SGTSV ', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL SGTSV( 2, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1, INFO )
CALL CHKXER( 'SGTSV ', INFOT, NOUT, LERR, OK )
*
* SGTSVX
*
SRNAMT = 'SGTSVX'
INFOT = 1
CALL SGTSVX( '/', 'N', 0, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
$ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
$ IP, B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGTSVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SGTSVX( 'N', '/', 0, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
$ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
$ IP, B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGTSVX', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL SGTSVX( 'N', 'N', -1, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
$ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
$ IP, B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGTSVX', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL SGTSVX( 'N', 'N', 0, -1, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
$ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
$ IP, B, 1, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGTSVX', INFOT, NOUT, LERR, OK )
INFOT = 14
CALL SGTSVX( 'N', 'N', 2, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
$ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
$ IP, B, 1, X, 2, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGTSVX', INFOT, NOUT, LERR, OK )
INFOT = 16
CALL SGTSVX( 'N', 'N', 2, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
$ AF( 1, 1 ), AF( 1, 2 ), AF( 1, 3 ), AF( 1, 4 ),
$ IP, B, 2, X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SGTSVX', INFOT, NOUT, LERR, OK )
*
ELSE IF( LSAMEN( 2, C2, 'PO' ) ) THEN
*
* SPOSV
*
SRNAMT = 'SPOSV '
INFOT = 1
CALL SPOSV( '/', 0, 0, A, 1, B, 1, INFO )
CALL CHKXER( 'SPOSV ', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SPOSV( 'U', -1, 0, A, 1, B, 1, INFO )
CALL CHKXER( 'SPOSV ', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL SPOSV( 'U', 0, -1, A, 1, B, 1, INFO )
CALL CHKXER( 'SPOSV ', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL SPOSV( 'U', 2, 0, A, 1, B, 2, INFO )
CALL CHKXER( 'SPOSV ', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL SPOSV( 'U', 2, 0, A, 2, B, 1, INFO )
CALL CHKXER( 'SPOSV ', INFOT, NOUT, LERR, OK )
*
* SPOSVX
*
SRNAMT = 'SPOSVX'
INFOT = 1
CALL SPOSVX( '/', 'U', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
$ RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SPOSVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL SPOSVX( 'N', '/', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
$ RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'SPOSVX', INFOT, NOUT, LERR, OK )
INFOT = 3
⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?