derrgg.f
来自「famous linear algebra library (LAPACK) p」· F 代码 · 共 954 行 · 第 1/3 页
F
954 行
CALL DGGEV( 'N', '/', 1, A, 1, B, 1, R1, R2, R3, Q, 1, U, 1, W,
$ 1, INFO )
CALL CHKXER( 'DGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DGGEV( 'V', 'V', -1, A, 1, B, 1, R1, R2, R3, Q, 1, U, 1,
$ W, 1, INFO )
CALL CHKXER( 'DGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL DGGEV( 'V', 'V', 1, A, 0, B, 1, R1, R2, R3, Q, 1, U, 1, W,
$ 1, INFO )
CALL CHKXER( 'DGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL DGGEV( 'V', 'V', 1, A, 1, B, 0, R1, R2, R3, Q, 1, U, 1, W,
$ 1, INFO )
CALL CHKXER( 'DGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 12
CALL DGGEV( 'N', 'V', 1, A, 1, B, 1, R1, R2, R3, Q, 0, U, 1, W,
$ 1, INFO )
CALL CHKXER( 'DGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 12
CALL DGGEV( 'V', 'V', 2, A, 2, B, 2, R1, R2, R3, Q, 1, U, 2, W,
$ 1, INFO )
CALL CHKXER( 'DGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 14
CALL DGGEV( 'V', 'N', 2, A, 2, B, 2, R1, R2, R3, Q, 2, U, 0, W,
$ 1, INFO )
CALL CHKXER( 'DGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 14
CALL DGGEV( 'V', 'V', 2, A, 2, B, 2, R1, R2, R3, Q, 2, U, 1, W,
$ 1, INFO )
CALL CHKXER( 'DGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 16
CALL DGGEV( 'V', 'V', 1, A, 1, B, 1, R1, R2, R3, Q, 1, U, 1, W,
$ 1, INFO )
CALL CHKXER( 'DGGEV ', INFOT, NOUT, LERR, OK )
NT = NT + 10
*
* DGGEVX
*
SRNAMT = 'DGGEVX'
INFOT = 1
CALL DGGEVX( '/', 'N', 'N', 'N', 1, A, 1, B, 1, R1, R2, R3, Q,
$ 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ IW, BW, INFO )
CALL CHKXER( 'DGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DGGEVX( 'N', '/', 'N', 'N', 1, A, 1, B, 1, R1, R2, R3, Q,
$ 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ IW, BW, INFO )
CALL CHKXER( 'DGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DGGEVX( 'N', 'N', '/', 'N', 1, A, 1, B, 1, R1, R2, R3, Q,
$ 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ IW, BW, INFO )
CALL CHKXER( 'DGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL DGGEVX( 'N', 'N', 'N', '/', 1, A, 1, B, 1, R1, R2, R3, Q,
$ 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ IW, BW, INFO )
CALL CHKXER( 'DGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL DGGEVX( 'N', 'N', 'N', 'N', -1, A, 1, B, 1, R1, R2, R3, Q,
$ 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ IW, BW, INFO )
CALL CHKXER( 'DGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL DGGEVX( 'N', 'N', 'N', 'N', 1, A, 0, B, 1, R1, R2, R3, Q,
$ 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ IW, BW, INFO )
CALL CHKXER( 'DGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DGGEVX( 'N', 'N', 'N', 'N', 1, A, 1, B, 0, R1, R2, R3, Q,
$ 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ IW, BW, INFO )
CALL CHKXER( 'DGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 14
CALL DGGEVX( 'N', 'N', 'N', 'N', 1, A, 1, B, 1, R1, R2, R3, Q,
$ 0, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ IW, BW, INFO )
CALL CHKXER( 'DGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 14
CALL DGGEVX( 'N', 'V', 'N', 'N', 2, A, 2, B, 2, R1, R2, R3, Q,
$ 1, U, 2, 1, 2, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ IW, BW, INFO )
CALL CHKXER( 'DGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 16
CALL DGGEVX( 'N', 'N', 'N', 'N', 1, A, 1, B, 1, R1, R2, R3, Q,
$ 1, U, 0, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ IW, BW, INFO )
CALL CHKXER( 'DGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 16
CALL DGGEVX( 'N', 'N', 'V', 'N', 2, A, 2, B, 2, R1, R2, R3, Q,
$ 2, U, 1, 1, 2, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ IW, BW, INFO )
CALL CHKXER( 'DGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 26
CALL DGGEVX( 'N', 'N', 'V', 'N', 2, A, 2, B, 2, R1, R2, R3, Q,
$ 2, U, 2, 1, 2, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ IW, BW, INFO )
CALL CHKXER( 'DGGEVX', INFOT, NOUT, LERR, OK )
NT = NT + 12
*
* DTGEXC
*
SRNAMT = 'DTGEXC'
INFOT = 3
CALL DTGEXC( .TRUE., .TRUE., -1, A, 1, B, 1, Q, 1, Z, 1, IFST,
$ ILST, W, 1, INFO )
CALL CHKXER( 'DTGEXC', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL DTGEXC( .TRUE., .TRUE., 1, A, 0, B, 1, Q, 1, Z, 1, IFST,
$ ILST, W, 1, INFO )
CALL CHKXER( 'DTGEXC', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL DTGEXC( .TRUE., .TRUE., 1, A, 1, B, 0, Q, 1, Z, 1, IFST,
$ ILST, W, 1, INFO )
CALL CHKXER( 'DTGEXC', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DTGEXC( .FALSE., .TRUE., 1, A, 1, B, 1, Q, 0, Z, 1, IFST,
$ ILST, W, 1, INFO )
CALL CHKXER( 'DTGEXC', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DTGEXC( .TRUE., .TRUE., 1, A, 1, B, 1, Q, 0, Z, 1, IFST,
$ ILST, W, 1, INFO )
CALL CHKXER( 'DTGEXC', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL DTGEXC( .TRUE., .FALSE., 1, A, 1, B, 1, Q, 1, Z, 0, IFST,
$ ILST, W, 1, INFO )
CALL CHKXER( 'DTGEXC', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL DTGEXC( .TRUE., .TRUE., 1, A, 1, B, 1, Q, 1, Z, 0, IFST,
$ ILST, W, 1, INFO )
CALL CHKXER( 'DTGEXC', INFOT, NOUT, LERR, OK )
INFOT = 15
CALL DTGEXC( .TRUE., .TRUE., 1, A, 1, B, 1, Q, 1, Z, 1, IFST,
$ ILST, W, 0, INFO )
CALL CHKXER( 'DTGEXC', INFOT, NOUT, LERR, OK )
NT = NT + 8
*
* DTGSEN
*
SRNAMT = 'DTGSEN'
INFOT = 1
CALL DTGSEN( -1, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, R1, R2,
$ R3, Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'DTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL DTGSEN( 1, .TRUE., .TRUE., SEL, -1, A, 1, B, 1, R1, R2,
$ R3, Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'DTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL DTGSEN( 1, .TRUE., .TRUE., SEL, 1, A, 0, B, 1, R1, R2, R3,
$ Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'DTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL DTGSEN( 1, .TRUE., .TRUE., SEL, 1, A, 1, B, 0, R1, R2, R3,
$ Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'DTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 14
CALL DTGSEN( 1, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, R1, R2, R3,
$ Q, 0, Z, 1, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'DTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 16
CALL DTGSEN( 1, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, R1, R2, R3,
$ Q, 1, Z, 0, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'DTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 22
CALL DTGSEN( 0, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, R1, R2, R3,
$ Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'DTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 22
CALL DTGSEN( 1, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, R1, R2, R3,
$ Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'DTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 22
CALL DTGSEN( 2, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, R1, R2, R3,
$ Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'DTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 24
CALL DTGSEN( 0, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, R1, R2, R3,
$ Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 20, IW, 0,
$ INFO )
CALL CHKXER( 'DTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 24
CALL DTGSEN( 1, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, R1, R2, R3,
$ Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 20, IW, 0,
$ INFO )
CALL CHKXER( 'DTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 24
CALL DTGSEN( 2, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, R1, R2, R3,
$ Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 20, IW, 1,
$ INFO )
CALL CHKXER( 'DTGSEN', INFOT, NOUT, LERR, OK )
NT = NT + 12
*
* DTGSNA
*
SRNAMT = 'DTGSNA'
INFOT = 1
CALL DTGSNA( '/', 'A', SEL, 1, A, 1, B, 1, Q, 1, U, 1, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'DTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DTGSNA( 'B', '/', SEL, 1, A, 1, B, 1, Q, 1, U, 1, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'DTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL DTGSNA( 'B', 'A', SEL, -1, A, 1, B, 1, Q, 1, U, 1, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'DTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL DTGSNA( 'B', 'A', SEL, 1, A, 0, B, 1, Q, 1, U, 1, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'DTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 8
CALL DTGSNA( 'B', 'A', SEL, 1, A, 1, B, 0, Q, 1, U, 1, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'DTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 10
CALL DTGSNA( 'E', 'A', SEL, 1, A, 1, B, 1, Q, 0, U, 1, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'DTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 12
CALL DTGSNA( 'E', 'A', SEL, 1, A, 1, B, 1, Q, 1, U, 0, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'DTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 15
CALL DTGSNA( 'E', 'A', SEL, 1, A, 1, B, 1, Q, 1, U, 1, R1, R2,
$ 0, M, W, 1, IW, INFO )
CALL CHKXER( 'DTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 18
CALL DTGSNA( 'E', 'A', SEL, 1, A, 1, B, 1, Q, 1, U, 1, R1, R2,
$ 1, M, W, 0, IW, INFO )
CALL CHKXER( 'DTGSNA', INFOT, NOUT, LERR, OK )
NT = NT + 9
*
* DTGSYL
*
SRNAMT = 'DTGSYL'
INFOT = 1
CALL DTGSYL( '/', 0, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'DTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL DTGSYL( 'N', -1, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'DTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL DTGSYL( 'N', 0, 0, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'DTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL DTGSYL( 'N', 0, 1, 0, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'DTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL DTGSYL( 'N', 0, 1, 1, A, 0, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'DTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 8
CALL DTGSYL( 'N', 0, 1, 1, A, 1, B, 0, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'DTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 10
CALL DTGSYL( 'N', 0, 1, 1, A, 1, B, 1, Q, 0, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'DTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 12
CALL DTGSYL( 'N', 0, 1, 1, A, 1, B, 1, Q, 1, U, 0, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'DTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 14
CALL DTGSYL( 'N', 0, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 0, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'DTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 16
CALL DTGSYL( 'N', 0, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 0,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'DTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 20
CALL DTGSYL( 'N', 1, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'DTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 20
CALL DTGSYL( 'N', 2, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'DTGSYL', INFOT, NOUT, LERR, OK )
NT = NT + 12
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 DERRGG
*
END
⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?