version 1.13, 2017/06/17 11:06:21
|
version 1.17, 2023/08/07 08:38:52
|
Line 55
|
Line 55
|
*> 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 ERR_BNDS_NORM |
*> componentwise error bound if possible. See comments for ERR_BNDS_NORM |
*> and ERR_BNDS_COMP for details of the error bounds. Note that this |
*> and ERR_BNDS_COMP 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 |
*> ERR_BNDS_NORM and ERR_BNDS_COMP. |
*> ERR_BNDS_NORM and ERR_BNDS_COMP. |
*> \endverbatim |
*> \endverbatim |
* |
* |
Line 66
|
Line 66
|
*> \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 177
|
Line 177
|
*> |
*> |
*> \param[in,out] Y |
*> \param[in,out] Y |
*> \verbatim |
*> \verbatim |
*> Y is DOUBLE PRECISION array, dimension |
*> Y is DOUBLE PRECISION array, dimension (LDY,NRHS) |
*> (LDY,NRHS) |
|
*> On entry, the solution matrix X, as computed by DGBTRS. |
*> On entry, the solution matrix X, as computed by DGBTRS. |
*> On exit, the improved solution matrix Y. |
*> On exit, the improved solution matrix Y. |
*> \endverbatim |
*> \endverbatim |
Line 210
|
Line 209
|
*> |
*> |
*> \param[in,out] ERR_BNDS_NORM |
*> \param[in,out] ERR_BNDS_NORM |
*> \verbatim |
*> \verbatim |
*> ERR_BNDS_NORM is DOUBLE PRECISION array, dimension |
*> ERR_BNDS_NORM 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 257
|
Line 255
|
*> |
*> |
*> \param[in,out] ERR_BNDS_COMP |
*> \param[in,out] ERR_BNDS_COMP |
*> \verbatim |
*> \verbatim |
*> ERR_BNDS_COMP is DOUBLE PRECISION array, dimension |
*> ERR_BNDS_COMP 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 273
|
Line 270
|
*> 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 |
*> ERR_BNDS_COMP is not accessed. If N_ERR_BNDS .LT. 3, then at most |
*> ERR_BNDS_COMP 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 ERR_BNDS_COMP(i,:) corresponds to the ith |
*> The first index in ERR_BNDS_COMP(i,:) corresponds to the ith |
Line 372
|
Line 369
|
*> 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 401
|
Line 398
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date December 2016 |
|
* |
|
*> \ingroup doubleGBcomputational |
*> \ingroup doubleGBcomputational |
* |
* |
* ===================================================================== |
* ===================================================================== |
Line 414
|
Line 409
|
$ 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, LDAB, LDAFB, LDB, LDY, N, KL, KU, NRHS, |
INTEGER INFO, LDAB, LDAFB, LDB, LDY, N, KL, KU, NRHS, |
Line 707
|
Line 701
|
END DO |
END DO |
* |
* |
RETURN |
RETURN |
|
* |
|
* End of DLA_GBRFSX_EXTENDED |
|
* |
END |
END |