cerrgg.f
来自「famous linear algebra library (LAPACK) p」· F 代码 · 共 957 行 · 第 1/3 页
F
957 行
INFOT = 24
CALL CGGESX( 'V', 'V', 'S', CLCTSX, 'V', 1, A, 1, B, 1, SDIM,
$ ALPHA, BETA, Q, 1, U, 1, RCE, RCV, W, 32, RW, IW,
$ 0, BW, INFO )
CALL CHKXER( 'CGGESX', INFOT, NOUT, LERR, OK )
NT = NT + 13
*
* CGGEV
*
SRNAMT = 'CGGEV '
INFOT = 1
CALL CGGEV( '/', 'N', 1, A, 1, B, 1, ALPHA, BETA, Q, 1, U, 1,
$ W, 1, RW, INFO )
CALL CHKXER( 'CGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL CGGEV( 'N', '/', 1, A, 1, B, 1, ALPHA, BETA, Q, 1, U, 1,
$ W, 1, RW, INFO )
CALL CHKXER( 'CGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL CGGEV( 'V', 'V', -1, A, 1, B, 1, ALPHA, BETA, Q, 1, U, 1,
$ W, 1, RW, INFO )
CALL CHKXER( 'CGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL CGGEV( 'V', 'V', 1, A, 0, B, 1, ALPHA, BETA, Q, 1, U, 1,
$ W, 1, RW, INFO )
CALL CHKXER( 'CGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL CGGEV( 'V', 'V', 1, A, 1, B, 0, ALPHA, BETA, Q, 1, U, 1,
$ W, 1, RW, INFO )
CALL CHKXER( 'CGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL CGGEV( 'N', 'V', 1, A, 1, B, 1, ALPHA, BETA, Q, 0, U, 1,
$ W, 1, RW, INFO )
CALL CHKXER( 'CGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL CGGEV( 'V', 'V', 2, A, 2, B, 2, ALPHA, BETA, Q, 1, U, 2,
$ W, 1, RW, INFO )
CALL CHKXER( 'CGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 13
CALL CGGEV( 'V', 'N', 2, A, 2, B, 2, ALPHA, BETA, Q, 2, U, 0,
$ W, 1, RW, INFO )
CALL CHKXER( 'CGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 13
CALL CGGEV( 'V', 'V', 2, A, 2, B, 2, ALPHA, BETA, Q, 2, U, 1,
$ W, 1, RW, INFO )
CALL CHKXER( 'CGGEV ', INFOT, NOUT, LERR, OK )
INFOT = 15
CALL CGGEV( 'V', 'V', 1, A, 1, B, 1, ALPHA, BETA, Q, 1, U, 1,
$ W, 1, RW, INFO )
CALL CHKXER( 'CGGEV ', INFOT, NOUT, LERR, OK )
NT = NT + 10
*
* CGGEVX
*
SRNAMT = 'CGGEVX'
INFOT = 1
CALL CGGEVX( '/', 'N', 'N', 'N', 1, A, 1, B, 1, ALPHA, BETA, Q,
$ 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ RW, IW, BW, INFO )
CALL CHKXER( 'CGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL CGGEVX( 'N', '/', 'N', 'N', 1, A, 1, B, 1, ALPHA, BETA, Q,
$ 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ RW, IW, BW, INFO )
CALL CHKXER( 'CGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL CGGEVX( 'N', 'N', '/', 'N', 1, A, 1, B, 1, ALPHA, BETA, Q,
$ 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ RW, IW, BW, INFO )
CALL CHKXER( 'CGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL CGGEVX( 'N', 'N', 'N', '/', 1, A, 1, B, 1, ALPHA, BETA, Q,
$ 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ RW, IW, BW, INFO )
CALL CHKXER( 'CGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL CGGEVX( 'N', 'N', 'N', 'N', -1, A, 1, B, 1, ALPHA, BETA,
$ Q, 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W,
$ 1, RW, IW, BW, INFO )
CALL CHKXER( 'CGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL CGGEVX( 'N', 'N', 'N', 'N', 1, A, 0, B, 1, ALPHA, BETA, Q,
$ 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ RW, IW, BW, INFO )
CALL CHKXER( 'CGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL CGGEVX( 'N', 'N', 'N', 'N', 1, A, 1, B, 0, ALPHA, BETA, Q,
$ 1, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ RW, IW, BW, INFO )
CALL CHKXER( 'CGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 13
CALL CGGEVX( 'N', 'N', 'N', 'N', 1, A, 1, B, 1, ALPHA, BETA, Q,
$ 0, U, 1, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ RW, IW, BW, INFO )
CALL CHKXER( 'CGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 13
CALL CGGEVX( 'N', 'V', 'N', 'N', 2, A, 2, B, 2, ALPHA, BETA, Q,
$ 1, U, 2, 1, 2, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ RW, IW, BW, INFO )
CALL CHKXER( 'CGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 15
CALL CGGEVX( 'N', 'N', 'N', 'N', 1, A, 1, B, 1, ALPHA, BETA, Q,
$ 1, U, 0, 1, 1, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ RW, IW, BW, INFO )
CALL CHKXER( 'CGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 15
CALL CGGEVX( 'N', 'N', 'V', 'N', 2, A, 2, B, 2, ALPHA, BETA, Q,
$ 2, U, 1, 1, 2, LS, RS, ANRM, BNRM, RCE, RCV, W, 1,
$ RW, IW, BW, INFO )
CALL CHKXER( 'CGGEVX', INFOT, NOUT, LERR, OK )
INFOT = 25
CALL CGGEVX( 'N', 'N', 'V', 'N', 2, A, 2, B, 2, ALPHA, BETA, Q,
$ 2, U, 2, 1, 2, LS, RS, ANRM, BNRM, RCE, RCV, W, 0,
$ RW, IW, BW, INFO )
CALL CHKXER( 'CGGEVX', INFOT, NOUT, LERR, OK )
NT = NT + 12
*
* CTGEXC
*
SRNAMT = 'CTGEXC'
INFOT = 3
CALL CTGEXC( .TRUE., .TRUE., -1, A, 1, B, 1, Q, 1, Z, 1, IFST,
$ ILST, INFO )
CALL CHKXER( 'CTGEXC', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL CTGEXC( .TRUE., .TRUE., 1, A, 0, B, 1, Q, 1, Z, 1, IFST,
$ ILST, INFO )
CALL CHKXER( 'CTGEXC', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL CTGEXC( .TRUE., .TRUE., 1, A, 1, B, 0, Q, 1, Z, 1, IFST,
$ ILST, INFO )
CALL CHKXER( 'CTGEXC', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL CTGEXC( .FALSE., .TRUE., 1, A, 1, B, 1, Q, 0, Z, 1, IFST,
$ ILST, INFO )
CALL CHKXER( 'CTGEXC', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL CTGEXC( .TRUE., .TRUE., 1, A, 1, B, 1, Q, 0, Z, 1, IFST,
$ ILST, INFO )
CALL CHKXER( 'CTGEXC', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL CTGEXC( .TRUE., .FALSE., 1, A, 1, B, 1, Q, 1, Z, 0, IFST,
$ ILST, INFO )
CALL CHKXER( 'CTGEXC', INFOT, NOUT, LERR, OK )
INFOT = 11
CALL CTGEXC( .TRUE., .TRUE., 1, A, 1, B, 1, Q, 1, Z, 0, IFST,
$ ILST, INFO )
CALL CHKXER( 'CTGEXC', INFOT, NOUT, LERR, OK )
NT = NT + 7
*
* CTGSEN
*
SRNAMT = 'CTGSEN'
INFOT = 1
CALL CTGSEN( -1, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, ALPHA,
$ BETA, Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'CTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL CTGSEN( 1, .TRUE., .TRUE., SEL, -1, A, 1, B, 1, ALPHA,
$ BETA, Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'CTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL CTGSEN( 1, .TRUE., .TRUE., SEL, 1, A, 0, B, 1, ALPHA,
$ BETA, Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'CTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL CTGSEN( 1, .TRUE., .TRUE., SEL, 1, A, 1, B, 0, ALPHA,
$ BETA, Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'CTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 13
CALL CTGSEN( 1, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, ALPHA,
$ BETA, Q, 0, Z, 1, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'CTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 15
CALL CTGSEN( 1, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, ALPHA,
$ BETA, Q, 1, Z, 0, M, TOLA, TOLB, RCV, W, 1, IW, 1,
$ INFO )
CALL CHKXER( 'CTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 21
CALL CTGSEN( 3, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, ALPHA,
$ BETA, Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, -5, IW,
$ 1, INFO )
CALL CHKXER( 'CTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 23
CALL CTGSEN( 0, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, ALPHA,
$ BETA, Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 20, IW,
$ 0, INFO )
CALL CHKXER( 'CTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 23
CALL CTGSEN( 1, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, ALPHA,
$ BETA, Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 20, IW,
$ 0, INFO )
CALL CHKXER( 'CTGSEN', INFOT, NOUT, LERR, OK )
INFOT = 23
CALL CTGSEN( 5, .TRUE., .TRUE., SEL, 1, A, 1, B, 1, ALPHA,
$ BETA, Q, 1, Z, 1, M, TOLA, TOLB, RCV, W, 20, IW,
$ 1, INFO )
CALL CHKXER( 'CTGSEN', INFOT, NOUT, LERR, OK )
NT = NT + 11
*
* CTGSNA
*
SRNAMT = 'CTGSNA'
INFOT = 1
CALL CTGSNA( '/', 'A', SEL, 1, A, 1, B, 1, Q, 1, U, 1, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'CTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL CTGSNA( 'B', '/', SEL, 1, A, 1, B, 1, Q, 1, U, 1, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'CTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL CTGSNA( 'B', 'A', SEL, -1, A, 1, B, 1, Q, 1, U, 1, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'CTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL CTGSNA( 'B', 'A', SEL, 1, A, 0, B, 1, Q, 1, U, 1, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'CTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 8
CALL CTGSNA( 'B', 'A', SEL, 1, A, 1, B, 0, Q, 1, U, 1, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'CTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 10
CALL CTGSNA( 'E', 'A', SEL, 1, A, 1, B, 1, Q, 0, U, 1, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'CTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 12
CALL CTGSNA( 'E', 'A', SEL, 1, A, 1, B, 1, Q, 1, U, 0, R1, R2,
$ 1, M, W, 1, IW, INFO )
CALL CHKXER( 'CTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 15
CALL CTGSNA( 'E', 'A', SEL, 1, A, 1, B, 1, Q, 1, U, 1, R1, R2,
$ 0, M, W, 1, IW, INFO )
CALL CHKXER( 'CTGSNA', INFOT, NOUT, LERR, OK )
INFOT = 18
CALL CTGSNA( 'E', 'A', SEL, 1, A, 1, B, 1, Q, 1, U, 1, R1, R2,
$ 1, M, W, 0, IW, INFO )
CALL CHKXER( 'CTGSNA', INFOT, NOUT, LERR, OK )
NT = NT + 9
*
* CTGSYL
*
SRNAMT = 'CTGSYL'
INFOT = 1
CALL CTGSYL( '/', 0, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'CTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL CTGSYL( 'N', -1, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'CTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL CTGSYL( 'N', 0, 0, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'CTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 4
CALL CTGSYL( 'N', 0, 1, 0, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'CTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 6
CALL CTGSYL( 'N', 0, 1, 1, A, 0, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'CTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 8
CALL CTGSYL( 'N', 0, 1, 1, A, 1, B, 0, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'CTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 10
CALL CTGSYL( 'N', 0, 1, 1, A, 1, B, 1, Q, 0, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'CTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 12
CALL CTGSYL( 'N', 0, 1, 1, A, 1, B, 1, Q, 1, U, 0, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'CTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 14
CALL CTGSYL( 'N', 0, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 0, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'CTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 16
CALL CTGSYL( 'N', 0, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 0,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'CTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 20
CALL CTGSYL( 'N', 1, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'CTGSYL', INFOT, NOUT, LERR, OK )
INFOT = 20
CALL CTGSYL( 'N', 2, 1, 1, A, 1, B, 1, Q, 1, U, 1, V, 1, Z, 1,
$ SCALE, DIF, W, 1, IW, INFO )
CALL CHKXER( 'CTGSYL', 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 CERRGG
*
END
⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?