zchkrfp.f

Go to the documentation of this file.

*> \brief \b ZCHKRFP
*
*  =========== DOCUMENTATION ===========
*
* Online html documentation available at 
*            http://www.netlib.org/lapack/explore-html/ 
*
*  Definition:
*  ===========
*
*       PROGRAM ZCHKRFP
* 
*
*> \par Purpose:
*  =============
*>
*> \verbatim
*>
*> ZCHKRFP is the main test program for the COMPLEX*16 linear equation
*> routines with RFP storage format
*>
*> \endverbatim
*
*  Arguments:
*  ==========
*
*> \verbatim
*>  MAXIN   INTEGER
*>          The number of different values that can be used for each of
*>          M, N, or NB
*>
*>  MAXRHS  INTEGER
*>          The maximum number of right hand sides
*>
*>  NTYPES  INTEGER
*>
*>  NMAX    INTEGER
*>          The maximum allowable value for N.
*>
*>  NIN     INTEGER
*>          The unit number for input
*>
*>  NOUT    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. 
*
*> \date April 2012
*
*> \ingroup complex16_lin
*
*  =====================================================================
      PROGRAM zchkrfp
*
*  -- LAPACK test routine (version 3.4.1) --
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*     April 2012
*
*  =====================================================================
*
*     .. Parameters ..
      INTEGER            MAXIN
      parameter                ( maxin = 12 )
      INTEGER            NMAX
      parameter                ( nmax =  50 )
      INTEGER            MAXRHS
      parameter                ( maxrhs = 16 )
      INTEGER            NTYPES
      parameter                ( ntypes = 9 )
      INTEGER            NIN, NOUT
      parameter                ( nin = 5, nout = 6 )
*     ..
*     .. Local Scalars ..
      LOGICAL            FATAL, TSTERR
      INTEGER            VERS_MAJOR, VERS_MINOR, VERS_PATCH
      INTEGER            I, NN, NNS, NNT
      DOUBLE PRECISION   EPS, S1, S2, THRESH

*     ..
*     .. Local Arrays ..
      INTEGER            NVAL( maxin ), NSVAL( maxin ), NTVAL( ntypes )
      COMPLEX*16         WORKA( nmax, nmax )
      COMPLEX*16         WORKASAV( nmax, nmax )
      COMPLEX*16         WORKB( nmax, maxrhs )
      COMPLEX*16         WORKXACT( nmax, maxrhs )
      COMPLEX*16         WORKBSAV( nmax, maxrhs )
      COMPLEX*16         WORKX( nmax, maxrhs )
      COMPLEX*16         WORKAFAC( nmax, nmax )
      COMPLEX*16         WORKAINV( nmax, nmax )
      COMPLEX*16         WORKARF( (nmax*(nmax+1))/2 )
      COMPLEX*16         WORKAP( (nmax*(nmax+1))/2 )
      COMPLEX*16         WORKARFINV( (nmax*(nmax+1))/2 )
      COMPLEX*16         Z_WORK_ZLATMS( 3 * nmax )
      COMPLEX*16         Z_WORK_ZPOT02( nmax, maxrhs )
      COMPLEX*16         Z_WORK_ZPOT03( nmax, nmax )
      DOUBLE PRECISION   D_WORK_ZLATMS( nmax )
      DOUBLE PRECISION   D_WORK_ZLANHE( nmax )
      DOUBLE PRECISION   D_WORK_ZPOT01( nmax )
      DOUBLE PRECISION   D_WORK_ZPOT02( nmax )
      DOUBLE PRECISION   D_WORK_ZPOT03( nmax )
*     ..
*     .. External Functions ..
      DOUBLE PRECISION   DLAMCH, DSECND
      EXTERNAL           dlamch, dsecnd
*     ..
*     .. External Subroutines ..
      EXTERNAL           ilaver, zdrvrfp, zdrvrf1, zdrvrf2, zdrvrf3,
     +                   zdrvrf4
*     ..
*     .. Executable Statements ..
*
      s1 = dsecnd( )
      fatal = .false.
*
*     Read a dummy line.
*
      READ( nin, fmt = * )
*
*     Report LAPACK version tag (e.g. LAPACK-3.2.0)
*
      CALL ilaver( vers_major, vers_minor, vers_patch )
      WRITE( nout, fmt = 9994 ) vers_major, vers_minor, vers_patch
*
*     Read the values of N
*
      READ( nin, fmt = * )nn
      IF( nn.LT.1 ) THEN
         WRITE( nout, fmt = 9996 )' NN ', nn, 1
         nn = 0
         fatal = .true.
      ELSE IF( nn.GT.maxin ) THEN
         WRITE( nout, fmt = 9995 )' NN ', nn, maxin
         nn = 0
         fatal = .true.
      END IF
      READ( nin, fmt = * )( nval( i ), i = 1, nn )
      DO 10 i = 1, nn
         IF( nval( i ).LT.0 ) THEN
            WRITE( nout, fmt = 9996 )' M  ', nval( i ), 0
            fatal = .true.
         ELSE IF( nval( i ).GT.nmax ) THEN
            WRITE( nout, fmt = 9995 )' M  ', nval( i ), nmax
            fatal = .true.
         END IF
   10 CONTINUE
      IF( nn.GT.0 )
     $   WRITE( nout, fmt = 9993 )'N   ', ( nval( i ), i = 1, nn )
*
*     Read the values of NRHS
*
      READ( nin, fmt = * )nns
      IF( nns.LT.1 ) THEN
         WRITE( nout, fmt = 9996 )' NNS', nns, 1
         nns = 0
         fatal = .true.
      ELSE IF( nns.GT.maxin ) THEN
         WRITE( nout, fmt = 9995 )' NNS', nns, maxin
         nns = 0
         fatal = .true.
      END IF
      READ( nin, fmt = * )( nsval( i ), i = 1, nns )
      DO 30 i = 1, nns
         IF( nsval( i ).LT.0 ) THEN
            WRITE( nout, fmt = 9996 )'NRHS', nsval( i ), 0
            fatal = .true.
         ELSE IF( nsval( i ).GT.maxrhs ) THEN
            WRITE( nout, fmt = 9995 )'NRHS', nsval( i ), maxrhs
            fatal = .true.
         END IF
   30 CONTINUE
      IF( nns.GT.0 )
     $   WRITE( nout, fmt = 9993 )'NRHS', ( nsval( i ), i = 1, nns )
*
*     Read the matrix types
*
      READ( nin, fmt = * )nnt
      IF( nnt.LT.1 ) THEN
         WRITE( nout, fmt = 9996 )' NMA', nnt, 1
         nnt = 0
         fatal = .true.
      ELSE IF( nnt.GT.ntypes ) THEN
         WRITE( nout, fmt = 9995 )' NMA', nnt, ntypes
         nnt = 0
         fatal = .true.
      END IF
      READ( nin, fmt = * )( ntval( i ), i = 1, nnt )
      DO 320 i = 1, nnt
         IF( ntval( i ).LT.0 ) THEN
            WRITE( nout, fmt = 9996 )'TYPE', ntval( i ), 0
            fatal = .true.
         ELSE IF( ntval( i ).GT.ntypes ) THEN
            WRITE( nout, fmt = 9995 )'TYPE', ntval( i ), ntypes
            fatal = .true.
         END IF
  320 CONTINUE
      IF( nnt.GT.0 )
     $   WRITE( nout, fmt = 9993 )'TYPE', ( ntval( i ), i = 1, nnt )
*
*     Read the threshold value for the test ratios.
*
      READ( nin, fmt = * )thresh
      WRITE( nout, fmt = 9992 )thresh
*
*     Read the flag that indicates whether to test the error exits.
*
      READ( nin, fmt = * )tsterr
*
      IF( fatal ) THEN
         WRITE( nout, fmt = 9999 )
         stop
      END IF
*
      IF( fatal ) THEN
         WRITE( nout, fmt = 9999 )
         stop
      END IF
*
*     Calculate and print the machine dependent constants.
*
      eps = dlamch( 'Underflow threshold' )
      WRITE( nout, fmt = 9991 )'underflow', eps
      eps = dlamch( 'Overflow threshold' )
      WRITE( nout, fmt = 9991 )'overflow ', eps
      eps = dlamch( 'Epsilon' )
      WRITE( nout, fmt = 9991 )'precision', eps
      WRITE( nout, fmt = * )
*
*     Test the error exit of:
*
      IF( tsterr )
     $   CALL zerrrfp( nout )
*
*    Test the routines: zpftrf, zpftri, zpftrs (as in ZDRVPO).
*    This also tests the routines: ztfsm, ztftri, ztfttr, ztrttf.
*
      CALL zdrvrfp( nout, nn, nval, nns, nsval, nnt, ntval, thresh,
     $              worka, workasav, workafac, workainv, workb,
     $              workbsav, workxact, workx, workarf, workarfinv,
     $              z_work_zlatms, z_work_zpot02,
     $              z_work_zpot03, d_work_zlatms, d_work_zlanhe,
     $              d_work_zpot01, d_work_zpot02, d_work_zpot03 )
*
*    Test the routine: zlanhf
*
      CALL zdrvrf1( nout, nn, nval, thresh, worka, nmax, workarf,
     +              d_work_zlanhe )
*
*    Test the convertion routines:
*       zhfttp, ztpthf, ztfttr, ztrttf, ztrttp and ztpttr.
*
      CALL zdrvrf2( nout, nn, nval, worka, nmax, workarf,
     +              workap, workasav )
*
*    Test the routine: ztfsm
*
      CALL zdrvrf3( nout, nn, nval, thresh, worka, nmax, workarf,
     +              workainv, workafac, d_work_zlanhe,
     +              z_work_zpot03, z_work_zpot02 )
 
*
*    Test the routine: zhfrk
*
      CALL zdrvrf4( nout, nn, nval, thresh, worka, workafac, nmax,
     +              workarf, workainv, nmax,d_work_zlanhe)
*
      CLOSE ( nin )
      s2 = dsecnd( )
      WRITE( nout, fmt = 9998 )
      WRITE( nout, fmt = 9997 )s2 - s1
*
 9999 FORMAT( / ' Execution not attempted due to input errors' )
 9998 FORMAT( / ' End of tests' )
 9997 FORMAT( ' Total time used = ', f12.2, ' seconds', / )
 9996 FORMAT( ' !! Invalid input value: ', a4, '=', i6, '; must be >=',
     $      i6 )
 9995 FORMAT( ' !! Invalid input value: ', a4, '=', i6, '; must be <=',
     $      i6 )
 9994 FORMAT( /  ' Tests of the COMPLEX*16 LAPACK RFP routines ',
     $      / ' LAPACK VERSION ', i1, '.', i1, '.', i1,
     $      / / ' The following parameter values will be used:' )
 9993 FORMAT( 4x, a4, ':  ', 10i6, / 11x, 10i6 )
 9992 FORMAT( / ' Routines pass computational tests if test ratio is ',
     $      'less than', f8.2, / )
 9991 FORMAT( ' Relative machine ', a, ' is taken to be', d16.6 )
*
*     End of ZCHKRFP
*
      END