version 1.14, 2018/05/29 07:18:32
|
version 1.15, 2020/05/21 21:46:10
|
Line 45
|
Line 45
|
*> |
*> |
*> ZPOSVXX uses the Cholesky factorization A = U**T*U or A = L*L**T |
*> ZPOSVXX uses the Cholesky factorization A = U**T*U or A = L*L**T |
*> to compute the solution to a complex*16 system of linear equations |
*> to compute the solution to a complex*16 system of linear equations |
*> A * X = B, where A is an N-by-N symmetric positive definite matrix |
*> A * X = B, where A is an N-by-N Hermitian positive definite matrix |
*> and X and B are N-by-NRHS matrices. |
*> and X and B are N-by-NRHS matrices. |
*> |
*> |
*> If requested, both normwise and maximum componentwise error bounds |
*> If requested, both normwise and maximum componentwise error bounds |
Line 157
|
Line 157
|
*> \param[in,out] A |
*> \param[in,out] A |
*> \verbatim |
*> \verbatim |
*> A is COMPLEX*16 array, dimension (LDA,N) |
*> A is COMPLEX*16 array, dimension (LDA,N) |
*> On entry, the symmetric matrix A, except if FACT = 'F' and EQUED = |
*> On entry, the Hermitian matrix A, except if FACT = 'F' and EQUED = |
*> 'Y', then A must contain the equilibrated matrix |
*> 'Y', then A must contain the equilibrated matrix |
*> diag(S)*A*diag(S). If UPLO = 'U', the leading N-by-N upper |
*> diag(S)*A*diag(S). If UPLO = 'U', the leading N-by-N upper |
*> triangular part of A contains the upper triangular part of the |
*> triangular part of A contains the upper triangular part of the |
Line 365
|
Line 365
|
*> 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 401
|
Line 401
|
*> \param[in] NPARAMS |
*> \param[in] NPARAMS |
*> \verbatim |
*> \verbatim |
*> NPARAMS is INTEGER |
*> NPARAMS is INTEGER |
*> Specifies the number of parameters set in PARAMS. If .LE. 0, the |
*> Specifies the number of parameters set in PARAMS. If <= 0, the |
*> PARAMS array is never referenced and default values are used. |
*> PARAMS array is never referenced and default values are used. |
*> \endverbatim |
*> \endverbatim |
*> |
*> |
*> \param[in,out] PARAMS |
*> \param[in,out] PARAMS |
*> \verbatim |
*> \verbatim |
*> PARAMS is DOUBLE PRECISION array, dimension NPARAMS |
*> PARAMS is DOUBLE PRECISION array, dimension NPARAMS |
*> Specifies algorithm parameters. If an entry is .LT. 0.0, then |
*> Specifies algorithm parameters. If an entry is < 0.0, then |
*> that entry will be filled with default value used for that |
*> that entry will be filled with default value used for that |
*> parameter. Only positions up to NPARAMS are accessed; defaults |
*> parameter. Only positions up to NPARAMS are accessed; defaults |
*> are used for higher-numbered parameters. |
*> are used for higher-numbered parameters. |
Line 416
|
Line 416
|
*> PARAMS(LA_LINRX_ITREF_I = 1) : Whether to perform iterative |
*> PARAMS(LA_LINRX_ITREF_I = 1) : Whether to perform iterative |
*> refinement or not. |
*> refinement or not. |
*> Default: 1.0D+0 |
*> Default: 1.0D+0 |
*> = 0.0 : No refinement is performed, and no error bounds are |
*> = 0.0: No refinement is performed, and no error bounds are |
*> computed. |
*> computed. |
*> = 1.0 : Use the extra-precise refinement algorithm. |
*> = 1.0: Use the extra-precise refinement algorithm. |
*> (other values are reserved for future use) |
*> (other values are reserved for future use) |
*> |
*> |
*> PARAMS(LA_LINRX_ITHRESH_I = 2) : Maximum number of residual |
*> PARAMS(LA_LINRX_ITHRESH_I = 2) : Maximum number of residual |