Files
lapack/TESTING/LIN/zerrac.f
Julie e18d437924 Making LAPACK's code eternal... no more version and date in source files.
GitHub is now enabling us to track accurately version and date.
No need for this anymore.
2021-03-25 10:16:58 -07:00

147 lines
3.8 KiB
FortranFixed

*> \brief \b ZERRAC
*
* =========== DOCUMENTATION ===========
*
* Online html documentation available at
* http://www.netlib.org/lapack/explore-html/
*
* Definition:
* ===========
*
* SUBROUTINE ZERRAC( NUNIT )
*
* .. Scalar Arguments ..
* INTEGER NUNIT
* ..
*
*
*> \par Purpose:
* =============
*>
*> \verbatim
*>
*> ZERRPX tests the error exits for ZCPOSV.
*> \endverbatim
*
* Arguments:
* ==========
*
*> \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 complex16_lin
*
* =====================================================================
SUBROUTINE ZERRAC( 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 ..
INTEGER NUNIT
* ..
*
* =====================================================================
*
* .. Parameters ..
INTEGER NMAX
PARAMETER ( NMAX = 4 )
* ..
* .. Local Scalars ..
INTEGER I, INFO, ITER, J
* ..
* .. Local Arrays ..
COMPLEX*16 A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
$ C( NMAX ), R( NMAX ), R1( NMAX ), R2( NMAX ),
$ W( 2*NMAX ), X( NMAX )
DOUBLE PRECISION RWORK( NMAX )
COMPLEX*16 WORK(NMAX*NMAX)
COMPLEX SWORK(NMAX*NMAX)
* ..
* .. External Subroutines ..
EXTERNAL CHKXER, ZCPOSV
* ..
* .. Scalars in Common ..
LOGICAL LERR, OK
CHARACTER*32 SRNAMT
INTEGER INFOT, NOUT
* ..
* .. Common blocks ..
COMMON / INFOC / INFOT, NOUT, OK, LERR
COMMON / SRNAMC / SRNAMT
* ..
* .. Intrinsic Functions ..
INTRINSIC DBLE
* ..
* .. Executable Statements ..
*
NOUT = NUNIT
WRITE( NOUT, FMT = * )
*
* Set the variables to innocuous values.
*
DO 20 J = 1, NMAX
DO 10 I = 1, NMAX
A( I, J ) = 1.D0 / DBLE( I+J )
AF( I, J ) = 1.D0 / DBLE( I+J )
10 CONTINUE
B( J ) = 0.D0
R1( J ) = 0.D0
R2( J ) = 0.D0
W( J ) = 0.D0
X( J ) = 0.D0
C( J ) = 0.D0
R( J ) = 0.D0
20 CONTINUE
OK = .TRUE.
*
SRNAMT = 'ZCPOSV'
INFOT = 1
CALL ZCPOSV('/',0,0,A,1,B,1,X,1,WORK,SWORK,RWORK,ITER,INFO)
CALL CHKXER( 'ZCPOSV', INFOT, NOUT, LERR, OK )
INFOT = 2
CALL ZCPOSV('U',-1,0,A,1,B,1,X,1,WORK,SWORK,RWORK,ITER,INFO)
CALL CHKXER( 'ZCPOSV', INFOT, NOUT, LERR, OK )
INFOT = 3
CALL ZCPOSV('U',0,-1,A,1,B,1,X,1,WORK,SWORK,RWORK,ITER,INFO)
CALL CHKXER( 'ZCPOSV', INFOT, NOUT, LERR, OK )
INFOT = 5
CALL ZCPOSV('U',2,1,A,1,B,2,X,2,WORK,SWORK,RWORK,ITER,INFO)
CALL CHKXER( 'ZCPOSV', INFOT, NOUT, LERR, OK )
INFOT = 7
CALL ZCPOSV('U',2,1,A,2,B,1,X,2,WORK,SWORK,RWORK,ITER,INFO)
CALL CHKXER( 'ZCPOSV', INFOT, NOUT, LERR, OK )
INFOT = 9
CALL ZCPOSV('U',2,1,A,2,B,2,X,1,WORK,SWORK,RWORK,ITER,INFO)
CALL CHKXER( 'ZCPOSV', INFOT, NOUT, LERR, OK )
*
* Print a summary line.
*
IF( OK ) THEN
WRITE( NOUT, FMT = 9999 )'ZCPOSV'
ELSE
WRITE( NOUT, FMT = 9998 )'ZCPOSV'
END IF
*
9999 FORMAT( 1X, A6, ' drivers passed the tests of the error exits' )
9998 FORMAT( ' *** ', A6, ' drivers failed the tests of the error ',
$ 'exits ***' )
*
RETURN
*
* End of ZERRAC
*
END