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