SUBROUTINE DERRAB( NUNIT ) IMPLICIT NONE * * -- LAPACK test routine (version 3.1.1) -- * Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. * January 2007 * * .. Scalar Arguments .. INTEGER NUNIT * .. * * Purpose * ======= * * DERRAB tests the error exits for DSGESV. * * Arguments * ========= * * NUNIT (input) INTEGER * The unit number for output. * * ===================================================================== * * .. Parameters .. INTEGER NMAX PARAMETER ( NMAX = 4 ) * .. * .. Local Scalars .. INTEGER I, INFO, ITER, J * .. * .. Local Arrays .. INTEGER IP( NMAX ) DOUBLE PRECISION A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ), $ C( NMAX ), R( NMAX ), R1( NMAX ), R2( NMAX ), $ W( 2*NMAX ), X( NMAX ) DOUBLE PRECISION WORK(1) REAL SWORK(1) * .. * .. External Subroutines .. EXTERNAL CHKXER, DSGESV * .. * .. 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 IP( J ) = J 20 CONTINUE OK = .TRUE. * SRNAMT = 'DSGESV' INFOT = 1 CALL DSGESV(-1,0,A,1,IP,B,1,X,1,WORK,SWORK,ITER,INFO) CALL CHKXER( 'DSGESV', INFOT, NOUT, LERR, OK ) INFOT = 2 CALL DSGESV(0,-1,A,1,IP,B,1,X,1,WORK,SWORK,ITER,INFO) CALL CHKXER( 'DSGESV', INFOT, NOUT, LERR, OK ) INFOT = 4 CALL DSGESV(2,1,A,1,IP,B,2,X,2,WORK,SWORK,ITER,INFO) CALL CHKXER( 'DSGESV', INFOT, NOUT, LERR, OK ) INFOT = 7 CALL DSGESV(2,1,A,2,IP,B,1,X,2,WORK,SWORK,ITER,INFO) CALL CHKXER( 'DSGESV', INFOT, NOUT, LERR, OK ) INFOT = 9 CALL DSGESV(2,1,A,2,IP,B,2,X,1,WORK,SWORK,ITER,INFO) CALL CHKXER( 'DSGESV', INFOT, NOUT, LERR, OK ) * * Print a summary line. * IF( OK ) THEN WRITE( NOUT, FMT = 9999 )'DSGESV' ELSE WRITE( NOUT, FMT = 9998 )'DSGESV' 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 DERRAB * END