--- rpl/lapack/lapack/dgesvxx.f 2011/11/21 20:42:52 1.5
+++ rpl/lapack/lapack/dgesvxx.f 2023/08/07 08:38:50 1.16
@@ -2,18 +2,18 @@
*
* =========== DOCUMENTATION ===========
*
-* Online html documentation available at
-* http://www.netlib.org/lapack/explore-html/
+* Online html documentation available at
+* http://www.netlib.org/lapack/explore-html/
*
*> \htmlonly
-*> Download DGESVXX + dependencies
-*>
-*> [TGZ]
-*>
-*> [ZIP]
-*>
+*> Download DGESVXX + dependencies
+*>
+*> [TGZ]
+*>
+*> [ZIP]
+*>
*> [TXT]
-*> \endhtmlonly
+*> \endhtmlonly
*
* Definition:
* ===========
@@ -23,7 +23,7 @@
* BERR, N_ERR_BNDS, ERR_BNDS_NORM,
* ERR_BNDS_COMP, NPARAMS, PARAMS, WORK, IWORK,
* INFO )
-*
+*
* .. Scalar Arguments ..
* CHARACTER EQUED, FACT, TRANS
* INTEGER INFO, LDA, LDAF, LDB, LDX, N, NRHS, NPARAMS,
@@ -38,7 +38,7 @@
* $ ERR_BNDS_NORM( NRHS, * ),
* $ ERR_BNDS_COMP( NRHS, * )
* ..
-*
+*
*
*> \par Purpose:
* =============
@@ -183,7 +183,7 @@
*>
*> \param[in,out] AF
*> \verbatim
-*> AF is or output) DOUBLE PRECISION array, dimension (LDAF,N)
+*> AF is DOUBLE PRECISION array, dimension (LDAF,N)
*> If FACT = 'F', then AF is an input argument and on entry
*> contains the factors L and U from the factorization
*> A = P*L*U as computed by DGETRF. If EQUED .ne. 'N', then
@@ -207,7 +207,7 @@
*>
*> \param[in,out] IPIV
*> \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
*> contains the pivot indices from the factorization A = P*L*U
*> as computed by DGETRF; row i of the matrix was interchanged
@@ -224,7 +224,7 @@
*>
*> \param[in,out] EQUED
*> \verbatim
-*> EQUED is or output) CHARACTER*1
+*> EQUED is CHARACTER*1
*> Specifies the form of equilibration that was done.
*> = 'N': No equilibration (always true if FACT = 'N').
*> = 'R': Row equilibration, i.e., A has been premultiplied by
@@ -239,7 +239,7 @@
*>
*> \param[in,out] R
*> \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
*> multiplied on the left by diag(R); if EQUED = 'N' or 'C', R
*> is not accessed. R is an input argument if FACT = 'F';
@@ -257,7 +257,7 @@
*>
*> \param[in,out] C
*> \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
*> multiplied on the right by diag(C); if EQUED = 'N' or 'R', C
*> is not accessed. C is an input argument if FACT = 'F';
@@ -411,7 +411,7 @@
*> information as described below. There currently are up to three
*> pieces of information returned for each right-hand side. If
*> 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 index in ERR_BNDS_COMP(i,:) corresponds to the ith
@@ -447,14 +447,14 @@
*> \param[in] NPARAMS
*> \verbatim
*> 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.
*> \endverbatim
*>
*> \param[in,out] PARAMS
*> \verbatim
-*> PARAMS is / output) DOUBLE PRECISION array, dimension (NPARAMS)
-*> Specifies algorithm parameters. If an entry is .LT. 0.0, then
+*> PARAMS is DOUBLE PRECISION array, dimension (NPARAMS)
+*> Specifies algorithm parameters. If an entry is < 0.0, then
*> that entry will be filled with default value used for that
*> parameter. Only positions up to NPARAMS are accessed; defaults
*> are used for higher-numbered parameters.
@@ -462,9 +462,9 @@
*> PARAMS(LA_LINRX_ITREF_I = 1) : Whether to perform iterative
*> refinement or not.
*> 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.
-*> = 1.0 : Use the extra-precise refinement algorithm.
+*> = 1.0: Use the extra-precise refinement algorithm.
*> (other values are reserved for future use)
*>
*> PARAMS(LA_LINRX_ITHRESH_I = 2) : Maximum number of residual
@@ -524,12 +524,10 @@
* Authors:
* ========
*
-*> \author Univ. of Tennessee
-*> \author Univ. of California Berkeley
-*> \author Univ. of Colorado Denver
-*> \author NAG Ltd.
-*
-*> \date November 2011
+*> \author Univ. of Tennessee
+*> \author Univ. of California Berkeley
+*> \author Univ. of Colorado Denver
+*> \author NAG Ltd.
*
*> \ingroup doubleGEsolve
*
@@ -540,10 +538,9 @@
$ ERR_BNDS_COMP, NPARAMS, PARAMS, WORK, IWORK,
$ INFO )
*
-* -- LAPACK driver routine (version 3.4.0) --
+* -- LAPACK driver routine --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
*
* .. Scalar Arguments ..
CHARACTER EQUED, FACT, TRANS
@@ -581,9 +578,9 @@
$ SMLNUM
* ..
* .. External Functions ..
- EXTERNAL LSAME, DLAMCH, DLA_RPVGRW
+ EXTERNAL LSAME, DLAMCH, DLA_GERPVGRW
LOGICAL LSAME
- DOUBLE PRECISION DLAMCH, DLA_RPVGRW
+ DOUBLE PRECISION DLAMCH, DLA_GERPVGRW
* ..
* .. External Subroutines ..
EXTERNAL DGEEQUB, DGETRF, DGETRS, DLACPY, DLAQGE,
@@ -732,14 +729,14 @@
* Compute the reciprocal pivot growth factor of the
* leading rank-deficient INFO columns of A.
*
- RPVGRW = DLA_RPVGRW( N, INFO, A, LDA, AF, LDAF )
+ RPVGRW = DLA_GERPVGRW( N, INFO, A, LDA, AF, LDAF )
RETURN
END IF
END IF
*
* Compute the reciprocal pivot growth factor RPVGRW.
*
- RPVGRW = DLA_RPVGRW( N, N, A, LDA, AF, LDAF )
+ RPVGRW = DLA_GERPVGRW( N, N, A, LDA, AF, LDAF )
*
* Compute the solution matrix X.
*