serrec.f

Go to the documentation of this file.

*> \brief \b SERREC
*
*  =========== DOCUMENTATION ===========
*
* Online html documentation available at
*            http://www.netlib.org/lapack/explore-html/
*
*  Definition:
*  ===========
*
*       SUBROUTINE SERREC( PATH, NUNIT )
*
*       .. Scalar Arguments ..
*       CHARACTER*3        PATH
*       INTEGER            NUNIT
*       ..
*
*
*> \par Purpose:
*  =============
*>
*> \verbatim
*>
*> SERREC tests the error exits for the routines for eigen- condition
*> estimation for REAL matrices:
*>    STRSYL, STREXC, STRSNA and STRSEN.
*> \endverbatim
*
*  Arguments:
*  ==========
*
*> \param[in] PATH
*> \verbatim
*>          PATH is CHARACTER*3
*>          The LAPACK path name for the routines to be tested.
*> \endverbatim
*>
*> \param[in] NUNIT
*> \verbatim
*>          NUNIT is INTEGER
*>          The unit number for output.
*> \endverbatim
*
*  Authors:
*  ========
*
*> \author Univ. of Tennessee
*> \author Univ. of California Berkeley
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
*> \ingroup single_eig
*
*  =====================================================================
      SUBROUTINE serrec( PATH, NUNIT )
*
*  -- LAPACK test routine --
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*
*     .. Scalar Arguments ..
      CHARACTER*3        PATH
      INTEGER            NUNIT
*     ..
*
*  =====================================================================
*
*     .. Parameters ..
      INTEGER            NMAX
      REAL               ONE, ZERO
      parameter( nmax = 4, one = 1.0e0, zero = 0.0e0 )
*     ..
*     .. Local Scalars ..
      INTEGER            I, IFST, ILST, INFO, J, M, NT
      REAL               SCALE
*     ..
*     .. Local Arrays ..
      LOGICAL            SEL( NMAX )
      INTEGER            IWORK( NMAX )
      REAL               A( NMAX, NMAX ), B( NMAX, NMAX ),
     $                   C( NMAX, NMAX ), S( NMAX ), SEP( NMAX ),
     $                   WI( NMAX ), WORK( NMAX ), WR( NMAX )
*     ..
*     .. External Subroutines ..
      EXTERNAL           chkxer, strexc, strsen, strsna, strsyl
*     ..
*     .. Scalars in Common ..
      LOGICAL            LERR, OK
      CHARACTER*32       SRNAMT
      INTEGER            INFOT, NOUT
*     ..
*     .. Common blocks ..
      COMMON             / infoc / infot, nout, ok, lerr
      COMMON             / srnamc / srnamt
*     ..
*     .. Executable Statements ..
*
      nout = nunit
      ok = .true.
      nt = 0
*
*     Initialize A, B and SEL
*
      DO 20 j = 1, nmax
         DO 10 i = 1, nmax
            a( i, j ) = zero
            b( i, j ) = zero
   10    CONTINUE
   20 CONTINUE
      DO 30 i = 1, nmax
         a( i, i ) = one
         sel( i ) = .true.
   30 CONTINUE
*
*     Test STRSYL
*
      srnamt = 'STRSYL'
      infot = 1
      CALL strsyl( 'X', 'N', 1, 0, 0, a, 1, b, 1, c, 1, scale, info )
      CALL chkxer( 'STRSYL', infot, nout, lerr, ok )
      infot = 2
      CALL strsyl( 'N', 'X', 1, 0, 0, a, 1, b, 1, c, 1, scale, info )
      CALL chkxer( 'STRSYL', infot, nout, lerr, ok )
      infot = 3
      CALL strsyl( 'N', 'N', 0, 0, 0, a, 1, b, 1, c, 1, scale, info )
      CALL chkxer( 'STRSYL', infot, nout, lerr, ok )
      infot = 4
      CALL strsyl( 'N', 'N', 1, -1, 0, a, 1, b, 1, c, 1, scale, info )
      CALL chkxer( 'STRSYL', infot, nout, lerr, ok )
      infot = 5
      CALL strsyl( 'n', 'n', 1, 0, -1, A, 1, B, 1, C, 1, SCALE, INFO )
      CALL CHKXER( 'strsyl', INFOT, NOUT, LERR, OK )
      INFOT = 7
      CALL STRSYL( 'n', 'n', 1, 2, 0, A, 1, B, 1, C, 2, SCALE, INFO )
      CALL CHKXER( 'strsyl', INFOT, NOUT, LERR, OK )
      INFOT = 9
      CALL STRSYL( 'n', 'n', 1, 0, 2, A, 1, B, 1, C, 1, SCALE, INFO )
      CALL CHKXER( 'strsyl', INFOT, NOUT, LERR, OK )
      INFOT = 11
      CALL STRSYL( 'n', 'n', 1, 2, 0, A, 2, B, 1, C, 1, SCALE, INFO )
      CALL CHKXER( 'strsyl', infot, nout, lerr, ok )
      nt = nt + 8
*
*     Test STREXC
*
      srnamt = 'STREXC'
      ifst = 1
      ilst = 1
      infot = 1
      CALL strexc( 'X', 1, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'STREXC', infot, nout, lerr, ok )
      infot = 2
      CALL strexc( 'N', -1, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'STREXC', infot, nout, lerr, ok )
      infot = 4
      ilst = 2
      CALL strexc( 'N', 2, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'STREXC', infot, nout, lerr, ok )
      infot = 6
      CALL strexc( 'V', 2, a, 2, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'STREXC', infot, nout, lerr, ok )
      infot = 7
      ifst = 0
      ilst = 1
      CALL strexc( 'V', 1, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'STREXC', infot, nout, lerr, ok )
      infot = 7
      ifst = 2
      CALL strexc( 'V', 1, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'STREXC', infot, nout, lerr, ok )
      infot = 8
      ifst = 1
      ilst = 0
      CALL strexc( 'V', 1, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'STREXC', infot, nout, lerr, ok )
      infot = 8
      ilst = 2
      CALL strexc( 'V', 1, a, 1, b, 1, ifst, ilst, work, info )
      CALL chkxer( 'STREXC', infot, nout, lerr, ok )
      nt = nt + 8
*
*     Test STRSNA
*
      srnamt = 'STRSNA'
      infot = 1
      CALL strsna( 'X', 'A', sel, 0, a, 1, b, 1, c, 1, s, sep, 1, m,
     $             work, 1, iwork, info )
      CALL chkxer( 'STRSNA', infot, nout, lerr, ok )
      infot = 2
      CALL strsna( 'b', 'x', SEL, 0, A, 1, B, 1, C, 1, S, SEP, 1, M,
     $             WORK, 1, IWORK, INFO )
      CALL CHKXER( 'strsna', INFOT, NOUT, LERR, OK )
      INFOT = 4
      CALL STRSNA( 'b', 'a', SEL, -1, A, 1, B, 1, C, 1, S, SEP, 1, M,
     $             WORK, 1, IWORK, INFO )
      CALL CHKXER( 'strsna', INFOT, NOUT, LERR, OK )
      INFOT = 6
      CALL STRSNA( 'v', 'a', SEL, 2, A, 1, B, 1, C, 1, S, SEP, 2, M,
     $             WORK, 2, IWORK, INFO )
      CALL CHKXER( 'strsna', INFOT, NOUT, LERR, OK )
      INFOT = 8
      CALL STRSNA( 'b', 'a', SEL, 2, A, 2, B, 1, C, 2, S, SEP, 2, M,
     $             WORK, 2, IWORK, INFO )
      CALL CHKXER( 'strsna', INFOT, NOUT, LERR, OK )
      INFOT = 10
      CALL STRSNA( 'b', 'a', SEL, 2, A, 2, B, 2, C, 1, S, SEP, 2, M,
     $             WORK, 2, IWORK, INFO )
      CALL CHKXER( 'strsna', INFOT, NOUT, LERR, OK )
      INFOT = 13
      CALL STRSNA( 'b', 'a', SEL, 1, A, 1, B, 1, C, 1, S, SEP, 0, M,
     $             WORK, 1, IWORK, INFO )
      CALL CHKXER( 'strsna', INFOT, NOUT, LERR, OK )
      INFOT = 13
      CALL STRSNA( 'b', 's', sel, 2, a, 2, b, 2, c, 2, s, sep, 1, m,
     $             work, 2, iwork, info )
      CALL chkxer( 'STRSNA', infot, nout, lerr, ok )
      infot = 16
      CALL strsna( 'B', 'A', sel, 2, a, 2, b, 2, c, 2, s, sep, 2, m,
     $             work, 1, iwork, info )
      CALL chkxer( 'STRSNA', infot, nout, lerr, ok )
      nt = nt + 9
*
*     Test STRSEN
*
      sel( 1 ) = .false.
      srnamt = 'STRSEN'
      infot = 1
      CALL strsen( 'X', 'N', sel, 0, a, 1, b, 1, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 1, iwork, 1, info )
      CALL chkxer( 'STRSEN', infot, nout, lerr, ok )
      infot = 2
      CALL strsen( 'N', 'X', sel, 0, a, 1, b, 1, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 1, iwork, 1, info )
      CALL chkxer( 'STRSEN', infot, nout, lerr, ok )
      infot = 4
      CALL strsen( 'N', 'N', sel, -1, a, 1, b, 1, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 1, iwork, 1, info )
      CALL chkxer( 'STRSEN', infot, nout, lerr, ok )
      infot = 6
      CALL strsen( 'N', 'N', sel, 2, a, 1, b, 1, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 2, iwork, 1, info )
      CALL chkxer( 'STRSEN', infot, nout, lerr, ok )
      infot = 8
      CALL strsen( 'N', 'V', sel, 2, a, 2, b, 1, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 1, iwork, 1, info )
      CALL chkxer( 'STRSEN', infot, nout, lerr, ok )
      infot = 15
      CALL strsen( 'N', 'V', sel, 2, a, 2, b, 2, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 0, iwork, 1, info )
      CALL chkxer( 'STRSEN', infot, nout, lerr, ok )
      infot = 15
      CALL strsen( 'E', 'V', sel, 3, a, 3, b, 3, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 1, iwork, 1, info )
      CALL chkxer( 'STRSEN', infot, nout, lerr, ok )
      infot = 15
      CALL strsen( 'V', 'V', sel, 3, a, 3, b, 3, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 3, iwork, 2, info )
      CALL chkxer( 'STRSEN', infot, nout, lerr, ok )
      infot = 17
      CALL strsen( 'E', 'V', sel, 2, a, 2, b, 2, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 1, iwork, 0, info )
      CALL chkxer( 'STRSEN', infot, nout, lerr, ok )
      infot = 17
      CALL strsen( 'V', 'V', sel, 3, a, 3, b, 3, wr, wi, m, s( 1 ),
     $             sep( 1 ), work, 4, iwork, 1, info )
      CALL chkxer( 'STRSEN', infot, nout, lerr, ok )
      nt = nt + 10
*
*     Print a summary line.
*
      IF( ok ) THEN
         WRITE( nout, fmt = 9999 )path, nt
      ELSE
         WRITE( nout, fmt = 9998 )path
      END IF
*
      RETURN
 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 ex',
     $      'its ***' )
*
*     End of SERREC
*
      END