zchkunhr_col()

subroutine zchkunhr_col	(	double precision	THRESH,
		logical	TSTERR,
		integer	NM,
		integer, dimension( * )	MVAL,
		integer	NN,
		integer, dimension( * )	NVAL,
		integer	NNB,
		integer, dimension( * )	NBVAL,
		integer	NOUT
	)

ZCHKUNHR_COL

Purpose:

 ZCHKUNHR_COL tests ZUNHR_COL using ZLATSQR and ZGEMQRT. Therefore, ZLATSQR
 (used in ZGEQR) and ZGEMQRT (used in ZGEMQR) have to be tested
 before this test.

Parameters

[in]	THRESH	THRESH is DOUBLE PRECISION The threshold value for the test ratios. A result is included in the output file if RESULT >= THRESH. To have every test ratio printed, use THRESH = 0.
[in]	TSTERR	TSTERR is LOGICAL Flag that indicates whether error exits are to be tested.
[in]	NM	NM is INTEGER The number of values of M contained in the vector MVAL.
[in]	MVAL	MVAL is INTEGER array, dimension (NM) The values of the matrix row dimension M.
[in]	NN	NN is INTEGER The number of values of N contained in the vector NVAL.
[in]	NVAL	NVAL is INTEGER array, dimension (NN) The values of the matrix column dimension N.
[in]	NNB	NNB is INTEGER The number of values of NB contained in the vector NBVAL.
[in]	NBVAL	NBVAL is INTEGER array, dimension (NBVAL) The values of the blocksize NB.
[in]	NOUT	NOUT is INTEGER The unit number for output.

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.

Date

November 2019

Definition at line 107 of file zchkunhr_col.f.

      IMPLICIT NONE
*
*  -- LAPACK test routine (version 3.7.0) --
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*     December 2016
*
*     .. Scalar Arguments ..
      LOGICAL            TSTERR
      INTEGER            NM, NN, NNB, NOUT
      DOUBLE PRECISION   THRESH
*     ..
*     .. Array Arguments ..
      INTEGER            MVAL( * ), NBVAL( * ), NVAL( * )
*     ..
*
*  =====================================================================
*
*     .. Parameters ..
      INTEGER            NTESTS
      parameter( ntests = 6 )
*     ..
*     .. Local Scalars ..
      CHARACTER(LEN=3)   PATH
      INTEGER            I, IMB1, INB1, INB2, J, T, M, N, MB1, NB1,
     $                   NB2, NFAIL, NERRS, NRUN
*
*     .. Local Arrays ..
      DOUBLE PRECISION   RESULT( NTESTS )
*     ..
*     .. External Subroutines ..
      EXTERNAL           alahd, alasum, zerrunhr_col, zunhr_col01
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC  max, min
*     ..
*     .. Scalars in Common ..
      LOGICAL            LERR, OK
      CHARACTER(LEN=32)  SRNAMT
      INTEGER            INFOT, NUNIT
*     ..
*     .. Common blocks ..
      COMMON             / infoc / infot, nunit, ok, lerr
      COMMON             / srnamc / srnamt
*     ..
*     .. Executable Statements ..
*
*     Initialize constants
*
      path( 1: 1 ) = 'Z'
      path( 2: 3 ) = 'HH'
      nrun = 0
      nfail = 0
      nerrs = 0
*
*     Test the error exits
*
      IF( tsterr ) CALL zerrunhr_col( path, nout )
      infot = 0
*
*     Do for each value of M in MVAL.
*
      DO i = 1, nm
         m = mval( i )
*
*        Do for each value of N in NVAL.
*
         DO j = 1, nn
            n = nval( j )
*
*           Only for M >= N
*
            IF ( min( m, n ).GT.0 .AND. m.GE.n ) THEN
*
*              Do for each possible value of MB1
*
               DO imb1 = 1, nnb
                  mb1 = nbval( imb1 )
*
*                 Only for MB1 > N
*
                  IF ( mb1.GT.n ) THEN
*
*                    Do for each possible value of NB1
*
                     DO inb1 = 1, nnb
                        nb1 = nbval( inb1 )
*
*                       Do for each possible value of NB2
*
                        DO inb2 = 1, nnb
                           nb2 = nbval( inb2 )
*
                           IF( nb1.GT.0 .AND. nb2.GT.0 ) THEN
*
*                             Test ZUNHR_COL
*
                              CALL zunhr_col01( m, n, mb1, nb1, nb2,
     $                                          result )
*
*                             Print information about the tests that did
*                             not pass the threshold.
*
                              DO t = 1, ntests
                                 IF( result( t ).GE.thresh ) THEN
                                    IF( nfail.EQ.0 .AND. nerrs.EQ.0 )
     $                              CALL alahd( nout, path )
                                    WRITE( nout, fmt = 9999 ) m, n, mb1,
     $                                     nb1, nb2, t, result( t )
                                    nfail = nfail + 1
                                 END IF
                              END DO
                              nrun = nrun + ntests
                           END IF
                        END DO
                     END DO
                  END IF
                END DO
            END IF
         END DO
      END DO
*
*     Print a summary of the results.
*
      CALL alasum( path, nout, nfail, nrun, nerrs )
*
 9999 FORMAT( 'M=', i5, ', N=', i5, ', MB1=', i5,
     $        ', NB1=', i5, ', NB2=', i5,' test(', i2, ')=', g12.5 )
      RETURN
*
*     End of ZCHKUNHR_COL
*

Here is the call graph for this function:

Here is the caller graph for this function: