[BACK]Return to dsgesv.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/dsgesv.f between versions 1.11 and 1.22

version 1.11, 2012/07/31 11:06:37 version 1.22, 2023/08/07 08:39:06
Line 2 Line 2
 *  *
 *  =========== DOCUMENTATION ===========  *  =========== DOCUMENTATION ===========
 *  *
 * Online html documentation available at   * Online html documentation available at
 *            http://www.netlib.org/lapack/explore-html/   *            http://www.netlib.org/lapack/explore-html/
 *  *
 *> \htmlonly  *> \htmlonly
 *> Download DSGESV + dependencies   *> Download DSGESV + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dsgesv.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dsgesv.f">
 *> [TGZ]</a>   *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dsgesv.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dsgesv.f">
 *> [ZIP]</a>   *> [ZIP]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dsgesv.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dsgesv.f">
 *> [TXT]</a>  *> [TXT]</a>
 *> \endhtmlonly   *> \endhtmlonly
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
 *  *
 *       SUBROUTINE DSGESV( N, NRHS, A, LDA, IPIV, B, LDB, X, LDX, WORK,  *       SUBROUTINE DSGESV( N, NRHS, A, LDA, IPIV, B, LDB, X, LDX, WORK,
 *                          SWORK, ITER, INFO )  *                          SWORK, ITER, INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       INTEGER            INFO, ITER, LDA, LDB, LDX, N, NRHS  *       INTEGER            INFO, ITER, LDA, LDB, LDX, N, NRHS
 *       ..  *       ..
Line 30 Line 30
 *       DOUBLE PRECISION   A( LDA, * ), B( LDB, * ), WORK( N, * ),  *       DOUBLE PRECISION   A( LDA, * ), B( LDB, * ), WORK( N, * ),
 *      $                   X( LDX, * )  *      $                   X( LDX, * )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 92 Line 92
 *>          dimension (LDA,N)  *>          dimension (LDA,N)
 *>          On entry, the N-by-N coefficient matrix A.  *>          On entry, the N-by-N coefficient matrix A.
 *>          On exit, if iterative refinement has been successfully used  *>          On exit, if iterative refinement has been successfully used
 *>          (INFO.EQ.0 and ITER.GE.0, see description below), then A is  *>          (INFO = 0 and ITER >= 0, see description below), then A is
 *>          unchanged, if double precision factorization has been used  *>          unchanged, if double precision factorization has been used
 *>          (INFO.EQ.0 and ITER.LT.0, see description below), then the  *>          (INFO = 0 and ITER < 0, see description below), then the
 *>          array A contains the factors L and U from the factorization  *>          array A contains the factors L and U from the factorization
 *>          A = P*L*U; the unit diagonal elements of L are not stored.  *>          A = P*L*U; the unit diagonal elements of L are not stored.
 *> \endverbatim  *> \endverbatim
Line 111 Line 111
 *>          The pivot indices that define the permutation matrix P;  *>          The pivot indices that define the permutation matrix P;
 *>          row i of the matrix was interchanged with row IPIV(i).  *>          row i of the matrix was interchanged with row IPIV(i).
 *>          Corresponds either to the single precision factorization  *>          Corresponds either to the single precision factorization
 *>          (if INFO.EQ.0 and ITER.GE.0) or the double precision  *>          (if INFO = 0 and ITER >= 0) or the double precision
 *>          factorization (if INFO.EQ.0 and ITER.LT.0).  *>          factorization (if INFO = 0 and ITER < 0).
 *> \endverbatim  *> \endverbatim
 *>  *>
 *> \param[in] B  *> \param[in] B
Line 164 Line 164
 *>               -3 : failure of SGETRF  *>               -3 : failure of SGETRF
 *>               -31: stop the iterative refinement after the 30th  *>               -31: stop the iterative refinement after the 30th
 *>                    iterations  *>                    iterations
 *>          > 0: iterative refinement has been sucessfully used.  *>          > 0: iterative refinement has been successfully used.
 *>               Returns the number of iterations  *>               Returns the number of iterations
 *> \endverbatim  *> \endverbatim
 *>  *>
Line 182 Line 182
 *  Authors:  *  Authors:
 *  ========  *  ========
 *  *
 *> \author Univ. of Tennessee   *> \author Univ. of Tennessee
 *> \author Univ. of California Berkeley   *> \author Univ. of California Berkeley
 *> \author Univ. of Colorado Denver   *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.   *> \author NAG Ltd.
 *  
 *> \date November 2011  
 *  *
 *> \ingroup doubleGEsolve  *> \ingroup doubleGEsolve
 *  *
Line 195 Line 193
       SUBROUTINE DSGESV( N, NRHS, A, LDA, IPIV, B, LDB, X, LDX, WORK,        SUBROUTINE DSGESV( N, NRHS, A, LDA, IPIV, B, LDB, X, LDX, WORK,
      $                   SWORK, ITER, INFO )       $                   SWORK, ITER, INFO )
 *  *
 *  -- LAPACK driver routine (version 3.4.0) --  *  -- LAPACK driver 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..--
 *     November 2011  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            INFO, ITER, LDA, LDB, LDX, N, NRHS        INTEGER            INFO, ITER, LDA, LDB, LDX, N, NRHS
Line 230 Line 227
       DOUBLE PRECISION   ANRM, CTE, EPS, RNRM, XNRM        DOUBLE PRECISION   ANRM, CTE, EPS, RNRM, XNRM
 *  *
 *     .. External Subroutines ..  *     .. External Subroutines ..
       EXTERNAL           DAXPY, DGEMM, DLACPY, DLAG2S, SLAG2D, SGETRF,        EXTERNAL           DAXPY, DGEMM, DLACPY, DLAG2S, DGETRF, DGETRS,
      $                   SGETRS, XERBLA       $                   SGETRF, SGETRS, SLAG2D, XERBLA
 *     ..  *     ..
 *     .. External Functions ..  *     .. External Functions ..
       INTEGER            IDAMAX        INTEGER            IDAMAX
Line 406 Line 403
    30 CONTINUE     30 CONTINUE
 *  *
 *     If we are at this place of the code, this is because we have  *     If we are at this place of the code, this is because we have
 *     performed ITER=ITERMAX iterations and never satisified the  *     performed ITER=ITERMAX iterations and never satisfied the
 *     stopping criterion, set up the ITER flag accordingly and follow up  *     stopping criterion, set up the ITER flag accordingly and follow up
 *     on double precision routine.  *     on double precision routine.
 *  *
Line 428 Line 425
 *  *
       RETURN        RETURN
 *  *
 *     End of DSGESV.  *     End of DSGESV
 *  *
       END        END
Removed from v.1.11  
changed lines
  Added in v.1.22

CVSweb interface <joel.bertrand@systella.fr>