version 1.6, 2011/11/21 22:19:46
|
version 1.8, 2012/08/22 09:48:30
|
Line 183
|
Line 183
|
*> |
*> |
*> \param[in,out] AF |
*> \param[in,out] AF |
*> \verbatim |
*> \verbatim |
*> AF is or output) COMPLEX*16 array, dimension (LDAF,N) |
*> AF is COMPLEX*16 array, dimension (LDAF,N) |
*> If FACT = 'F', then AF is an input argument and on entry |
*> If FACT = 'F', then AF is an input argument and on entry |
*> contains the factors L and U from the factorization |
*> contains the factors L and U from the factorization |
*> A = P*L*U as computed by ZGETRF. If EQUED .ne. 'N', then |
*> A = P*L*U as computed by ZGETRF. If EQUED .ne. 'N', then |
Line 207
|
Line 207
|
*> |
*> |
*> \param[in,out] IPIV |
*> \param[in,out] IPIV |
*> \verbatim |
*> \verbatim |
*> IPIV is or output) INTEGER array, dimension (N) |
*> IPIV is INTEGER array, dimension (N) |
*> If FACT = 'F', then IPIV is an input argument and on entry |
*> If FACT = 'F', then IPIV is an input argument and on entry |
*> contains the pivot indices from the factorization A = P*L*U |
*> contains the pivot indices from the factorization A = P*L*U |
*> as computed by ZGETRF; row i of the matrix was interchanged |
*> as computed by ZGETRF; row i of the matrix was interchanged |
Line 224
|
Line 224
|
*> |
*> |
*> \param[in,out] EQUED |
*> \param[in,out] EQUED |
*> \verbatim |
*> \verbatim |
*> EQUED is or output) CHARACTER*1 |
*> EQUED is CHARACTER*1 |
*> Specifies the form of equilibration that was done. |
*> Specifies the form of equilibration that was done. |
*> = 'N': No equilibration (always true if FACT = 'N'). |
*> = 'N': No equilibration (always true if FACT = 'N'). |
*> = 'R': Row equilibration, i.e., A has been premultiplied by |
*> = 'R': Row equilibration, i.e., A has been premultiplied by |
Line 239
|
Line 239
|
*> |
*> |
*> \param[in,out] R |
*> \param[in,out] R |
*> \verbatim |
*> \verbatim |
*> R is or output) DOUBLE PRECISION array, dimension (N) |
*> R is DOUBLE PRECISION array, dimension (N) |
*> The row scale factors for A. If EQUED = 'R' or 'B', A is |
*> The row scale factors for A. If EQUED = 'R' or 'B', A is |
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R |
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R |
*> is not accessed. R is an input argument if FACT = 'F'; |
*> is not accessed. R is an input argument if FACT = 'F'; |
Line 257
|
Line 257
|
*> |
*> |
*> \param[in,out] C |
*> \param[in,out] C |
*> \verbatim |
*> \verbatim |
*> C is or output) DOUBLE PRECISION array, dimension (N) |
*> C is DOUBLE PRECISION array, dimension (N) |
*> The column scale factors for A. If EQUED = 'C' or 'B', A is |
*> The column scale factors for A. If EQUED = 'C' or 'B', A is |
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C |
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C |
*> is not accessed. C is an input argument if FACT = 'F'; |
*> is not accessed. C is an input argument if FACT = 'F'; |
Line 453
|
Line 453
|
*> |
*> |
*> \param[in,out] PARAMS |
*> \param[in,out] PARAMS |
*> \verbatim |
*> \verbatim |
*> PARAMS is / output) 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 .LT. 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 |
Line 529
|
Line 529
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date November 2011 |
*> \date April 2012 |
* |
* |
*> \ingroup complex16GEsolve |
*> \ingroup complex16GEsolve |
* |
* |
Line 540
|
Line 540
|
$ ERR_BNDS_COMP, NPARAMS, PARAMS, WORK, RWORK, |
$ ERR_BNDS_COMP, NPARAMS, PARAMS, WORK, RWORK, |
$ INFO ) |
$ INFO ) |
* |
* |
* -- LAPACK driver routine (version 3.4.0) -- |
* -- LAPACK driver routine (version 3.4.1) -- |
* -- 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 |
* April 2012 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER EQUED, FACT, TRANS |
CHARACTER EQUED, FACT, TRANS |
Line 581
|
Line 581
|
$ ROWCND, SMLNUM |
$ ROWCND, SMLNUM |
* .. |
* .. |
* .. External Functions .. |
* .. External Functions .. |
EXTERNAL LSAME, DLAMCH, ZLA_RPVGRW |
EXTERNAL LSAME, DLAMCH, ZLA_GERPVGRW |
LOGICAL LSAME |
LOGICAL LSAME |
DOUBLE PRECISION DLAMCH, ZLA_RPVGRW |
DOUBLE PRECISION DLAMCH, ZLA_GERPVGRW |
* .. |
* .. |
* .. External Subroutines .. |
* .. External Subroutines .. |
EXTERNAL ZGEEQUB, ZGETRF, ZGETRS, ZLACPY, ZLAQGE, |
EXTERNAL ZGEEQUB, ZGETRF, ZGETRS, ZLACPY, ZLAQGE, |
Line 732
|
Line 732
|
* Compute the reciprocal pivot growth factor of the |
* Compute the reciprocal pivot growth factor of the |
* leading rank-deficient INFO columns of A. |
* leading rank-deficient INFO columns of A. |
* |
* |
RPVGRW = ZLA_RPVGRW( N, INFO, A, LDA, AF, LDAF ) |
RPVGRW = ZLA_GERPVGRW( N, INFO, A, LDA, AF, LDAF ) |
RETURN |
RETURN |
END IF |
END IF |
END IF |
END IF |
* |
* |
* Compute the reciprocal pivot growth factor RPVGRW. |
* Compute the reciprocal pivot growth factor RPVGRW. |
* |
* |
RPVGRW = ZLA_RPVGRW( N, N, A, LDA, AF, LDAF ) |
RPVGRW = ZLA_GERPVGRW( N, N, A, LDA, AF, LDAF ) |
* |
* |
* Compute the solution matrix X. |
* Compute the solution matrix X. |
* |
* |