derrvx.f
来自「famous linear algebra library (LAPACK) p」· F 代码 · 共 653 行 · 第 1/2 页
F
653 行
SUBROUTINE DERRVX( 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
* =======
*
* DERRVX tests the error exits for the DOUBLE PRECISION 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
DOUBLE PRECISION RCOND
* ..
* .. Local Arrays ..
INTEGER IP( NMAX ), IW( NMAX )
DOUBLE PRECISION 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, DGBSV, DGBSVX, DGESV, DGESVX, DGTSV,
$ DGTSVX, DPBSV, DPBSVX, DPOSV, DPOSVX, DPPSV,
$ DPPSVX, DPTSV, DPTSVX, DSPSV, DSPSVX, DSYSV,
$ DSYSVX
* ..
* .. 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 DBLE
* ..
* .. 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.D0 / DBLE( I+J )
AF( I, J ) = 1.D0 / DBLE( I+J )
10 CONTINUE
B( J ) = 0.D0
R1( J ) = 0.D0
R2( J ) = 0.D0
W( J ) = 0.D0
X( J ) = 0.D0
C( J ) = 0.D0
R( J ) = 0.D0
IP( J ) = J
20 CONTINUE
EQ = ' '
OK = .TRUE.
*
IF( LSAMEN( 2, C2, 'GE' ) ) THEN
*
* DGESV
*
SRNAMT = 'DGESV '
INFOT = 1
CALL DGESV( -1, 0, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'DGESV ', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DGESV( 0, -1, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'DGESV ', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL DGESV( 2, 1, A, 1, IP, B, 2, INFO )
CALL CHKXER( 'DGESV ', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL DGESV( 2, 1, A, 2, IP, B, 1, INFO )
CALL CHKXER( 'DGESV ', INFOT, NOUT, LERR, OK )
*
* DGESVX
*
SRNAMT = 'DGESVX'
INFOT = 1
CALL DGESVX( '/', 'N', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'DGESVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DGESVX( 'N', '/', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'DGESVX', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DGESVX( 'N', 'N', -1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'DGESVX', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL DGESVX( 'N', 'N', 0, -1, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'DGESVX', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL DGESVX( 'N', 'N', 2, 1, A, 1, AF, 2, IP, EQ, R, C, B, 2,
$ X, 2, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'DGESVX', INFOT, NOUT, LERR, OK )
INFOT = 8
CALL DGESVX( 'N', 'N', 2, 1, A, 2, AF, 1, IP, EQ, R, C, B, 2,
$ X, 2, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'DGESVX', INFOT, NOUT, LERR, OK )
INFOT = 10
EQ = '/'
CALL DGESVX( 'F', 'N', 0, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'DGESVX', INFOT, NOUT, LERR, OK )
INFOT = 11
EQ = 'R'
CALL DGESVX( 'F', 'N', 1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'DGESVX', INFOT, NOUT, LERR, OK )
INFOT = 12
EQ = 'C'
CALL DGESVX( 'F', 'N', 1, 0, A, 1, AF, 1, IP, EQ, R, C, B, 1,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'DGESVX', INFOT, NOUT, LERR, OK )
INFOT = 14
CALL DGESVX( 'N', 'N', 2, 1, A, 2, AF, 2, IP, EQ, R, C, B, 1,
$ X, 2, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'DGESVX', INFOT, NOUT, LERR, OK )
INFOT = 16
CALL DGESVX( 'N', 'N', 2, 1, A, 2, AF, 2, IP, EQ, R, C, B, 2,
$ X, 1, RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'DGESVX', INFOT, NOUT, LERR, OK )
*
ELSE IF( LSAMEN( 2, C2, 'GB' ) ) THEN
*
* DGBSV
*
SRNAMT = 'DGBSV '
INFOT = 1
CALL DGBSV( -1, 0, 0, 0, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'DGBSV ', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DGBSV( 1, -1, 0, 0, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'DGBSV ', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DGBSV( 1, 0, -1, 0, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'DGBSV ', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL DGBSV( 0, 0, 0, -1, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'DGBSV ', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL DGBSV( 1, 1, 1, 0, A, 3, IP, B, 1, INFO )
CALL CHKXER( 'DGBSV ', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DGBSV( 2, 0, 0, 0, A, 1, IP, B, 1, INFO )
CALL CHKXER( 'DGBSV ', INFOT, NOUT, LERR, OK )
*
* DGBSVX
*
SRNAMT = 'DGBSVX'
INFOT = 1
CALL DGBSVX( '/', '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( 'DGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DGBSVX( '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( 'DGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DGBSVX( '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( 'DGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL DGBSVX( '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( 'DGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL DGBSVX( '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( 'DGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL DGBSVX( '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( 'DGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 8
CALL DGBSVX( '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( 'DGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 10
CALL DGBSVX( '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( 'DGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 12
EQ = '/'
CALL DGBSVX( '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( 'DGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 13
EQ = 'R'
CALL DGBSVX( '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( 'DGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 14
EQ = 'C'
CALL DGBSVX( '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( 'DGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 16
CALL DGBSVX( '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( 'DGBSVX', INFOT, NOUT, LERR, OK )
INFOT = 18
CALL DGBSVX( '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( 'DGBSVX', INFOT, NOUT, LERR, OK )
*
ELSE IF( LSAMEN( 2, C2, 'GT' ) ) THEN
*
* DGTSV
*
SRNAMT = 'DGTSV '
INFOT = 1
CALL DGTSV( -1, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1,
$ INFO )
CALL CHKXER( 'DGTSV ', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DGTSV( 0, -1, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1,
$ INFO )
CALL CHKXER( 'DGTSV ', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL DGTSV( 2, 0, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B, 1, INFO )
CALL CHKXER( 'DGTSV ', INFOT, NOUT, LERR, OK )
*
* DGTSVX
*
SRNAMT = 'DGTSVX'
INFOT = 1
CALL DGTSVX( '/', '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( 'DGTSVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DGTSVX( '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( 'DGTSVX', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DGTSVX( '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( 'DGTSVX', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL DGTSVX( '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( 'DGTSVX', INFOT, NOUT, LERR, OK )
INFOT = 14
CALL DGTSVX( '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( 'DGTSVX', INFOT, NOUT, LERR, OK )
INFOT = 16
CALL DGTSVX( '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( 'DGTSVX', INFOT, NOUT, LERR, OK )
*
ELSE IF( LSAMEN( 2, C2, 'PO' ) ) THEN
*
* DPOSV
*
SRNAMT = 'DPOSV '
INFOT = 1
CALL DPOSV( '/', 0, 0, A, 1, B, 1, INFO )
CALL CHKXER( 'DPOSV ', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DPOSV( 'U', -1, 0, A, 1, B, 1, INFO )
CALL CHKXER( 'DPOSV ', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DPOSV( 'U', 0, -1, A, 1, B, 1, INFO )
CALL CHKXER( 'DPOSV ', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL DPOSV( 'U', 2, 0, A, 1, B, 2, INFO )
CALL CHKXER( 'DPOSV ', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL DPOSV( 'U', 2, 0, A, 2, B, 1, INFO )
CALL CHKXER( 'DPOSV ', INFOT, NOUT, LERR, OK )
*
* DPOSVX
*
SRNAMT = 'DPOSVX'
INFOT = 1
CALL DPOSVX( '/', 'U', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
$ RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'DPOSVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DPOSVX( 'N', '/', 0, 0, A, 1, AF, 1, EQ, C, B, 1, X, 1,
$ RCOND, R1, R2, W, IW, INFO )
CALL CHKXER( 'DPOSVX', INFOT, NOUT, LERR, OK )
INFOT = 3
⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?