$ ITEMP, ITEMP2, IUPLO, IZERO, J, K, KL, KU, LDA,
$ LWORK, MODE, N, NB, NERRS, NFAIL, NIMAT, NRHS,
$ NRUN, NT
- REAL ALPHA, ANORM, CNDNUM, CONST, LAM_MAX, LAM_MIN,
- $ RCOND, RCONDC, STEMP
+ REAL ALPHA, ANORM, CNDNUM, CONST, SING_MAX,
+ $ SING_MIN, RCOND, RCONDC, STEMP
* ..
* .. Local Arrays ..
CHARACTER UPLOS( 2 )
$ CHETRI_ROOK, CHETRS_ROOK, XLAENV
* ..
* .. Intrinsic Functions ..
- INTRINSIC MAX, MIN, SQRT
+ INTRINSIC CONJG, MAX, MIN, SQRT
* ..
* .. Scalars in Common ..
LOGICAL LERR, OK
*
BLOCK( 1, 1 ) = AFAC( ( K-2 )*LDA+K-1 )
BLOCK( 1, 2 ) = AFAC( (K-1)*LDA+K-1 )
- BLOCK( 2, 1 ) = BLOCK( 1, 2 )
+ BLOCK( 2, 1 ) = CONJG( BLOCK( 1, 2 ) )
BLOCK( 2, 2 ) = AFAC( (K-1)*LDA+K )
*
CALL CGESVD( 'N', 'N', 2, 2, BLOCK, 2, RWORK,
$ WORK, 6, RWORK( 3 ), INFO )
*
*
- LAM_MAX = RWORK( 1 )
- LAM_MIN = RWORK( 2 )
+ SING_MAX = RWORK( 1 )
+ SING_MIN = RWORK( 2 )
*
- STEMP = LAM_MAX / LAM_MIN
+ STEMP = SING_MAX / SING_MIN
*
* STEMP should be bounded by CONST
*
*
BLOCK( 1, 1 ) = AFAC( ( K-1 )*LDA+K )
BLOCK( 2, 1 ) = AFAC( ( K-1 )*LDA+K+1 )
- BLOCK( 1, 2 ) = BLOCK( 2, 1 )
+ BLOCK( 1, 2 ) = CONJG( BLOCK( 2, 1 ) )
BLOCK( 2, 2 ) = AFAC( K*LDA+K+1 )
*
CALL CGESVD( 'N', 'N', 2, 2, BLOCK, 2, RWORK,
$ CDUMMY, 1, CDUMMY, 1,
$ WORK, 6, RWORK(3), INFO )
*
- LAM_MAX = RWORK( 1 )
- LAM_MIN = RWORK( 2 )
+ SING_MAX = RWORK( 1 )
+ SING_MIN = RWORK( 2 )
*
- STEMP = LAM_MAX / LAM_MIN
+ STEMP = SING_MAX / SING_MIN
*
* STEMP should be bounded by CONST
*
$ ITEMP, ITEMP2, IUPLO, IZERO, J, K, KL, KU, LDA,
$ LWORK, MODE, N, NB, NERRS, NFAIL, NIMAT, NRHS,
$ NRUN, NT
- REAL ALPHA, ANORM, CNDNUM, CONST, LAM_MAX, LAM_MIN,
- $ RCOND, RCONDC, STEMP
+ REAL ALPHA, ANORM, CNDNUM, CONST, SING_MAX,
+ $ SING_MIN, RCOND, RCONDC, STEMP
* ..
* .. Local Arrays ..
CHARACTER UPLOS( 2 )
$ WORK, 6, RWORK( 3 ), INFO )
*
*
- LAM_MAX = RWORK( 1 )
- LAM_MIN = RWORK( 2 )
+ SING_MAX = RWORK( 1 )
+ SING_MIN = RWORK( 2 )
*
- STEMP = LAM_MAX / LAM_MIN
+ STEMP = SING_MAX / SING_MIN
*
* STEMP should be bounded by CONST
*
$ CDUMMY, 1, CDUMMY, 1,
$ WORK, 6, RWORK(3), INFO )
*
- LAM_MAX = RWORK( 1 )
- LAM_MIN = RWORK( 2 )
+ SING_MAX = RWORK( 1 )
+ SING_MIN = RWORK( 2 )
*
- STEMP = LAM_MAX / LAM_MIN
+ STEMP = SING_MAX / SING_MIN
*
* STEMP should be bounded by CONST
*
$ ITEMP, IUPLO, IZERO, J, K, KL, KU, LDA, LWORK,
$ MODE, N, NB, NERRS, NFAIL, NIMAT, NRHS, NRUN,
$ NT
- DOUBLE PRECISION ALPHA, ANORM, CNDNUM, CONST, DTEMP, LAM_MAX,
- $ LAM_MIN, RCOND, RCONDC
+ DOUBLE PRECISION ALPHA, ANORM, CNDNUM, CONST, DTEMP, SING_MAX,
+ $ SING_MIN, RCOND, RCONDC
* ..
* .. Local Arrays ..
CHARACTER UPLOS( 2 )
$ DDUMMY, 1, DDUMMY, 1,
$ WORK, 10, INFO )
*
- LAM_MAX = RWORK( 1 )
- LAM_MIN = RWORK( 2 )
+ SING_MAX = RWORK( 1 )
+ SING_MIN = RWORK( 2 )
*
- DTEMP = LAM_MAX / LAM_MIN
+ DTEMP = SING_MAX / SING_MIN
*
* DTEMP should be bounded by CONST
*
$ WORK, 10, INFO )
*
*
- LAM_MAX = RWORK( 1 )
- LAM_MIN = RWORK( 2 )
+ SING_MAX = RWORK( 1 )
+ SING_MIN = RWORK( 2 )
*
- DTEMP = LAM_MAX / LAM_MIN
+ DTEMP = SING_MAX / SING_MIN
*
* DTEMP should be bounded by CONST
*
$ ITEMP, IUPLO, IZERO, J, K, KL, KU, LDA, LWORK,
$ MODE, N, NB, NERRS, NFAIL, NIMAT, NRHS, NRUN,
$ NT
- REAL ALPHA, ANORM, CNDNUM, CONST, LAM_MAX, LAM_MIN,
- $ RCOND, RCONDC, STEMP
+ REAL ALPHA, ANORM, CNDNUM, CONST, SING_MAX,
+ $ SING_MIN, RCOND, RCONDC, STEMP
* ..
* .. Local Arrays ..
CHARACTER UPLOS( 2 )
$ WORK, 10, INFO )
*
*
- LAM_MAX = RWORK( 1 )
- LAM_MIN = RWORK( 2 )
+ SING_MAX = RWORK( 1 )
+ SING_MIN = RWORK( 2 )
*
- STEMP = LAM_MAX / LAM_MIN
+ STEMP = SING_MAX / SING_MIN
*
* STEMP should be bounded by CONST
*
$ WORK, 10, INFO )
*
*
- LAM_MAX = RWORK( 1 )
- LAM_MIN = RWORK( 2 )
+ SING_MAX = RWORK( 1 )
+ SING_MIN = RWORK( 2 )
*
- STEMP = LAM_MAX / LAM_MIN
+ STEMP = SING_MAX / SING_MIN
*
* STEMP should be bounded by CONST
*
$ ITEMP, ITEMP2, IUPLO, IZERO, J, K, KL, KU, LDA,
$ LWORK, MODE, N, NB, NERRS, NFAIL, NIMAT, NRHS,
$ NRUN, NT
- DOUBLE PRECISION ALPHA, ANORM, CNDNUM, CONST, LAM_MAX, LAM_MIN,
- $ RCOND, RCONDC, DTEMP
+ DOUBLE PRECISION ALPHA, ANORM, CNDNUM, CONST, SING_MAX,
+ $ SING_MIN, RCOND, RCONDC, DTEMP
* ..
* .. Local Arrays ..
CHARACTER UPLOS( 2 )
$ ZHETRI_ROOK, ZHETRS_ROOK, XLAENV
* ..
* .. Intrinsic Functions ..
- INTRINSIC MAX, MIN, SQRT
+ INTRINSIC CONJG, MAX, MIN, SQRT
* ..
* .. Scalars in Common ..
LOGICAL LERR, OK
*
BLOCK( 1, 1 ) = AFAC( ( K-2 )*LDA+K-1 )
BLOCK( 1, 2 ) = AFAC( (K-1)*LDA+K-1 )
- BLOCK( 2, 1 ) = BLOCK( 1, 2 )
+ BLOCK( 2, 1 ) = CONJG( BLOCK( 1, 2 ) )
BLOCK( 2, 2 ) = AFAC( (K-1)*LDA+K )
*
CALL ZGESVD( 'N', 'N', 2, 2, BLOCK, 2, RWORK,
$ WORK, 6, RWORK( 3 ), INFO )
*
*
- LAM_MAX = RWORK( 1 )
- LAM_MIN = RWORK( 2 )
+ SING_MAX = RWORK( 1 )
+ SING_MIN = RWORK( 2 )
*
- DTEMP = LAM_MAX / LAM_MIN
+ DTEMP = SING_MAX / SING_MIN
*
* DTEMP should be bounded by CONST
*
*
BLOCK( 1, 1 ) = AFAC( ( K-1 )*LDA+K )
BLOCK( 2, 1 ) = AFAC( ( K-1 )*LDA+K+1 )
- BLOCK( 1, 2 ) = BLOCK( 2, 1 )
+ BLOCK( 1, 2 ) = CONJG( BLOCK( 2, 1 ) )
BLOCK( 2, 2 ) = AFAC( K*LDA+K+1 )
*
CALL ZGESVD( 'N', 'N', 2, 2, BLOCK, 2, RWORK,
$ ZDUMMY, 1, ZDUMMY, 1,
$ WORK, 6, RWORK(3), INFO )
*
- LAM_MAX = RWORK( 1 )
- LAM_MIN = RWORK( 2 )
+ SING_MAX = RWORK( 1 )
+ SING_MIN = RWORK( 2 )
*
- DTEMP = LAM_MAX / LAM_MIN
+ DTEMP = SING_MAX / SING_MIN
*
* DTEMP should be bounded by CONST
*
$ ITEMP, ITEMP2, IUPLO, IZERO, J, K, KL, KU, LDA,
$ LWORK, MODE, N, NB, NERRS, NFAIL, NIMAT, NRHS,
$ NRUN, NT
- DOUBLE PRECISION ALPHA, ANORM, CNDNUM, CONST, DTEMP, LAM_MAX,
- $ LAM_MIN, RCOND, RCONDC
+ DOUBLE PRECISION ALPHA, ANORM, CNDNUM, CONST, DTEMP, SING_MAX,
+ $ SING_MIN, RCOND, RCONDC
* ..
* .. Local Arrays ..
CHARACTER UPLOS( 2 )
$ WORK, 6, RWORK( 3 ), INFO )
*
*
- LAM_MAX = RWORK( 1 )
- LAM_MIN = RWORK( 2 )
+ SING_MAX = RWORK( 1 )
+ SING_MIN = RWORK( 2 )
*
- DTEMP = LAM_MAX / LAM_MIN
+ DTEMP = SING_MAX / SING_MIN
*
* DTEMP should be bounded by CONST
*
$ ZDUMMY, 1, ZDUMMY, 1,
$ WORK, 6, RWORK(3), INFO )
*
- LAM_MAX = RWORK( 1 )
- LAM_MIN = RWORK( 2 )
+ SING_MAX = RWORK( 1 )
+ SING_MIN = RWORK( 2 )
*
- DTEMP = LAM_MAX / LAM_MIN
+ DTEMP = SING_MAX / SING_MIN
*
* DTEMP should be bounded by CONST
*