version 1.12, 2017/06/17 10:53:52
|
version 1.17, 2023/08/07 08:38:52
|
Line 54
|
Line 54
|
*> 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 65
|
Line 65
|
*> \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] UPLO |
*> \param[in] UPLO |
Line 153
|
Line 153
|
*> |
*> |
*> \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 DPOTRS. |
*> On entry, the solution matrix X, as computed by DPOTRS. |
*> On exit, the improved solution matrix Y. |
*> On exit, the improved solution matrix Y. |
*> \endverbatim |
*> \endverbatim |
Line 186
|
Line 185
|
*> |
*> |
*> \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 233
|
Line 231
|
*> |
*> |
*> \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 249
|
Line 246
|
*> 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 348
|
Line 345
|
*> 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 377
|
Line 374
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date December 2016 |
|
* |
|
*> \ingroup doublePOcomputational |
*> \ingroup doublePOcomputational |
* |
* |
* ===================================================================== |
* ===================================================================== |
Line 390
|
Line 385
|
$ RTHRESH, DZ_UB, IGNORE_CWISE, |
$ RTHRESH, DZ_UB, IGNORE_CWISE, |
$ INFO ) |
$ 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 679
|
Line 673
|
END DO |
END DO |
* |
* |
RETURN |
RETURN |
|
* |
|
* End of DLA_PORFSX_EXTENDED |
|
* |
END |
END |