version 1.5, 2011/11/21 20:43:14
|
version 1.9, 2012/12/14 12:30:30
|
Line 1
|
Line 1
|
*> \brief \b ZLA_SYRFSX_EXTENDED |
*> \brief \b ZLA_SYRFSX_EXTENDED improves the computed solution to a system of linear equations for symmetric indefinite matrices by performing extra-precise iterative refinement and provides error bounds and backward error estimates for the solution. |
* |
* |
* =========== DOCUMENTATION =========== |
* =========== DOCUMENTATION =========== |
* |
* |
Line 373
|
Line 373
|
*> \verbatim |
*> \verbatim |
*> INFO is INTEGER |
*> INFO is INTEGER |
*> = 0: Successful exit. |
*> = 0: Successful exit. |
*> < 0: if INFO = -i, the ith argument to ZSYTRS had an illegal |
*> < 0: if INFO = -i, the ith argument to ZLA_HERFSX_EXTENDED had an illegal |
*> value |
*> value |
*> \endverbatim |
*> \endverbatim |
* |
* |
Line 385
|
Line 385
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date November 2011 |
*> \date September 2012 |
* |
* |
*> \ingroup complex16SYcomputational |
*> \ingroup complex16SYcomputational |
* |
* |
Line 398
|
Line 398
|
$ RTHRESH, DZ_UB, IGNORE_CWISE, |
$ RTHRESH, DZ_UB, IGNORE_CWISE, |
$ INFO ) |
$ INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.4.0) -- |
* -- LAPACK computational routine (version 3.4.2) -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* November 2011 |
* September 2012 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
INTEGER INFO, LDA, LDAF, LDB, LDY, N, NRHS, PREC_TYPE, |
INTEGER INFO, LDA, LDAF, LDB, LDY, N, NRHS, PREC_TYPE, |
Line 428
|
Line 428
|
$ DZRAT, PREVNORMDX, PREV_DZ_Z, DXRATMAX, |
$ DZRAT, PREVNORMDX, PREV_DZ_Z, DXRATMAX, |
$ DZRATMAX, DX_X, DZ_Z, FINAL_DX_X, FINAL_DZ_Z, |
$ DZRATMAX, DX_X, DZ_Z, FINAL_DX_X, FINAL_DZ_Z, |
$ EPS, HUGEVAL, INCR_THRESH |
$ EPS, HUGEVAL, INCR_THRESH |
LOGICAL INCR_PREC |
LOGICAL INCR_PREC, UPPER |
COMPLEX*16 ZDUM |
COMPLEX*16 ZDUM |
* .. |
* .. |
* .. Parameters .. |
* .. Parameters .. |
Line 479
|
Line 479
|
* .. |
* .. |
* .. Executable Statements .. |
* .. Executable Statements .. |
* |
* |
IF ( INFO.NE.0 ) RETURN |
INFO = 0 |
|
UPPER = LSAME( UPLO, 'U' ) |
|
IF( .NOT.UPPER .AND. .NOT.LSAME( UPLO, 'L' ) ) THEN |
|
INFO = -2 |
|
ELSE IF( N.LT.0 ) THEN |
|
INFO = -3 |
|
ELSE IF( NRHS.LT.0 ) THEN |
|
INFO = -4 |
|
ELSE IF( LDA.LT.MAX( 1, N ) ) THEN |
|
INFO = -6 |
|
ELSE IF( LDAF.LT.MAX( 1, N ) ) THEN |
|
INFO = -8 |
|
ELSE IF( LDB.LT.MAX( 1, N ) ) THEN |
|
INFO = -13 |
|
ELSE IF( LDY.LT.MAX( 1, N ) ) THEN |
|
INFO = -15 |
|
END IF |
|
IF( INFO.NE.0 ) THEN |
|
CALL XERBLA( 'ZLA_HERFSX_EXTENDED', -INFO ) |
|
RETURN |
|
END IF |
EPS = DLAMCH( 'Epsilon' ) |
EPS = DLAMCH( 'Epsilon' ) |
HUGEVAL = DLAMCH( 'Overflow' ) |
HUGEVAL = DLAMCH( 'Overflow' ) |
* Force HUGEVAL to Inf |
* Force HUGEVAL to Inf |