version 1.12, 2017/06/17 11:06:51
|
version 1.16, 2023/08/07 08:39:27
|
Line 53
|
Line 53
|
*> In addition to normwise error bound, the code provides maximum |
*> In addition to normwise error bound, the code provides maximum |
*> componentwise error bound if possible. See comments for ERRS_N |
*> componentwise error bound if possible. See comments for ERRS_N |
*> and ERRS_C for details of the error bounds. Note that this |
*> and ERRS_C for details of the error bounds. Note that this |
*> subroutine is only resonsible for setting the second fields of |
*> subroutine is only responsible for setting the second fields of |
*> ERRS_N and ERRS_C. |
*> ERRS_N and ERRS_C. |
*> \endverbatim |
*> \endverbatim |
* |
* |
Line 64
|
Line 64
|
*> \verbatim |
*> \verbatim |
*> PREC_TYPE is INTEGER |
*> PREC_TYPE is INTEGER |
*> Specifies the intermediate precision to be used in refinement. |
*> Specifies the intermediate precision to be used in refinement. |
*> The value is defined by ILAPREC(P) where P is a CHARACTER and |
*> The value is defined by ILAPREC(P) where P is a CHARACTER and P |
*> P = 'S': Single |
*> = 'S': Single |
*> = 'D': Double |
*> = 'D': Double |
*> = 'I': Indigenous |
*> = 'I': Indigenous |
*> = 'X', 'E': Extra |
*> = 'X' or 'E': Extra |
*> \endverbatim |
*> \endverbatim |
*> |
*> |
*> \param[in] TRANS_TYPE |
*> \param[in] TRANS_TYPE |
*> \verbatim |
*> \verbatim |
*> TRANS_TYPE is INTEGER |
*> TRANS_TYPE is INTEGER |
*> Specifies the transposition operation on A. |
*> Specifies the transposition operation on A. |
*> The value is defined by ILATRANS(T) where T is a CHARACTER and |
*> The value is defined by ILATRANS(T) where T is a CHARACTER and T |
*> T = 'N': No transpose |
*> = 'N': No transpose |
*> = 'T': Transpose |
*> = 'T': Transpose |
*> = 'C': Conjugate transpose |
*> = 'C': Conjugate transpose |
*> \endverbatim |
*> \endverbatim |
Line 195
|
Line 195
|
*> |
*> |
*> \param[in,out] ERRS_N |
*> \param[in,out] ERRS_N |
*> \verbatim |
*> \verbatim |
*> ERRS_N is DOUBLE PRECISION array, dimension |
*> ERRS_N is DOUBLE PRECISION array, dimension (NRHS, N_ERR_BNDS) |
*> (NRHS, N_ERR_BNDS) |
|
*> For each right-hand side, this array contains information about |
*> For each right-hand side, this array contains information about |
*> various error bounds and condition numbers corresponding to the |
*> various error bounds and condition numbers corresponding to the |
*> normwise relative error, which is defined as follows: |
*> normwise relative error, which is defined as follows: |
Line 242
|
Line 241
|
*> |
*> |
*> \param[in,out] ERRS_C |
*> \param[in,out] ERRS_C |
*> \verbatim |
*> \verbatim |
*> ERRS_C is DOUBLE PRECISION array, dimension |
*> ERRS_C is DOUBLE PRECISION array, dimension (NRHS, N_ERR_BNDS) |
*> (NRHS, N_ERR_BNDS) |
|
*> For each right-hand side, this array contains information about |
*> For each right-hand side, this array contains information about |
*> various error bounds and condition numbers corresponding to the |
*> various error bounds and condition numbers corresponding to the |
*> componentwise relative error, which is defined as follows: |
*> componentwise relative error, which is defined as follows: |
Line 258
|
Line 256
|
*> information as described below. There currently are up to three |
*> information as described below. There currently are up to three |
*> pieces of information returned for each right-hand side. If |
*> pieces of information returned for each right-hand side. If |
*> componentwise accuracy is not requested (PARAMS(3) = 0.0), then |
*> componentwise accuracy is not requested (PARAMS(3) = 0.0), then |
*> ERRS_C is not accessed. If N_ERR_BNDS .LT. 3, then at most |
*> ERRS_C is not accessed. If N_ERR_BNDS < 3, then at most |
*> the first (:,N_ERR_BNDS) entries are returned. |
*> the first (:,N_ERR_BNDS) entries are returned. |
*> |
*> |
*> The first index in ERRS_C(i,:) corresponds to the ith |
*> The first index in ERRS_C(i,:) corresponds to the ith |
Line 357
|
Line 355
|
*> DZ_UB is DOUBLE PRECISION |
*> DZ_UB is DOUBLE PRECISION |
*> Determines when to start considering componentwise convergence. |
*> Determines when to start considering componentwise convergence. |
*> Componentwise convergence is only considered after each component |
*> Componentwise convergence is only considered after each component |
*> of the solution Y is stable, which we definte as the relative |
*> of the solution Y is stable, which we define as the relative |
*> change in each component being less than DZ_UB. The default value |
*> change in each component being less than DZ_UB. The default value |
*> is 0.25, requiring the first bit to be stable. See LAWN 165 for |
*> is 0.25, requiring the first bit to be stable. See LAWN 165 for |
*> more details. |
*> more details. |
Line 386
|
Line 384
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date December 2016 |
|
* |
|
*> \ingroup complex16GEcomputational |
*> \ingroup complex16GEcomputational |
* |
* |
* ===================================================================== |
* ===================================================================== |
Line 398
|
Line 394
|
$ Y_TAIL, RCOND, ITHRESH, RTHRESH, |
$ Y_TAIL, RCOND, ITHRESH, RTHRESH, |
$ DZ_UB, IGNORE_CWISE, INFO ) |
$ DZ_UB, IGNORE_CWISE, INFO ) |
* |
* |
* -- LAPACK computational routine (version 3.7.0) -- |
* -- LAPACK computational routine -- |
* -- 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..-- |
* December 2016 |
|
* |
* |
* .. 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 696
|
Line 691
|
END DO |
END DO |
* |
* |
RETURN |
RETURN |
|
* |
|
* End of ZLA_GERFSX_EXTENDED |
|
* |
END |
END |