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

Diff for /rpl/lapack/lapack/dsyrfsx.f between versions 1.6 and 1.16

version 1.6, 2011/11/21 22:19:41 version 1.16, 2023/08/07 08:39:09
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 DSYRFSX + dependencies   *> Download DSYRFSX + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dsyrfsx.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dsyrfsx.f">
 *> [TGZ]</a>   *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dsyrfsx.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dsyrfsx.f">
 *> [ZIP]</a>   *> [ZIP]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dsyrfsx.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dsyrfsx.f">
 *> [TXT]</a>  *> [TXT]</a>
 *> \endhtmlonly   *> \endhtmlonly
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
Line 22 Line 22
 *                           S, B, LDB, X, LDX, RCOND, BERR, N_ERR_BNDS,  *                           S, B, LDB, X, LDX, RCOND, BERR, N_ERR_BNDS,
 *                           ERR_BNDS_NORM, ERR_BNDS_COMP, NPARAMS, PARAMS,  *                           ERR_BNDS_NORM, ERR_BNDS_COMP, NPARAMS, PARAMS,
 *                           WORK, IWORK, INFO )  *                           WORK, IWORK, INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       CHARACTER          UPLO, EQUED  *       CHARACTER          UPLO, EQUED
 *       INTEGER            INFO, LDA, LDAF, LDB, LDX, N, NRHS, NPARAMS,  *       INTEGER            INFO, LDA, LDAF, LDB, LDX, N, NRHS, NPARAMS,
Line 37 Line 37
 *      $                   ERR_BNDS_NORM( NRHS, * ),  *      $                   ERR_BNDS_NORM( NRHS, * ),
 *      $                   ERR_BNDS_COMP( NRHS, * )  *      $                   ERR_BNDS_COMP( NRHS, * )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 142 Line 142
 *>  *>
 *> \param[in,out] S  *> \param[in,out] S
 *> \verbatim  *> \verbatim
 *>          S is or output) DOUBLE PRECISION array, dimension (N)  *>          S is DOUBLE PRECISION array, dimension (N)
 *>     The scale factors for A.  If EQUED = 'Y', A is multiplied on  *>     The scale factors for A.  If EQUED = 'Y', A is multiplied on
 *>     the left and right by diag(S).  S is an input argument if FACT =  *>     the left and right by diag(S).  S is an input argument if FACT =
 *>     'F'; otherwise, S is an output argument.  If FACT = 'F' and EQUED  *>     'F'; otherwise, S is an output argument.  If FACT = 'F' and EQUED
Line 271 Line 271
 *>     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 307 Line 307
 *> \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 / 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 < 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 322 Line 322
 *>       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 double-precision refinement algorithm,  *>            = 1.0:  Use the double-precision refinement algorithm,
 *>                    possibly with doubled-single computations if the  *>                    possibly with doubled-single computations if the
 *>                    compilation environment does not support DOUBLE  *>                    compilation environment does not support DOUBLE
 *>                    PRECISION.  *>                    PRECISION.
Line 387 Line 387
 *  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 doubleSYcomputational  *> \ingroup doubleSYcomputational
 *  *
Line 402 Line 400
      $                    ERR_BNDS_NORM, ERR_BNDS_COMP, NPARAMS, PARAMS,       $                    ERR_BNDS_NORM, ERR_BNDS_COMP, NPARAMS, PARAMS,
      $                    WORK, IWORK, INFO )       $                    WORK, IWORK, INFO )
 *  *
 *  -- LAPACK computational routine (version 3.4.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..--
 *     November 2011  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          UPLO, EQUED        CHARACTER          UPLO, EQUED
Line 462 Line 459
       INTRINSIC          MAX, SQRT        INTRINSIC          MAX, SQRT
 *     ..  *     ..
 *     .. External Functions ..  *     .. External Functions ..
       EXTERNAL           LSAME, BLAS_FPINFO_X, ILATRANS, ILAPREC        EXTERNAL           LSAME, ILAPREC
       EXTERNAL           DLAMCH, DLANSY, DLA_SYRCOND        EXTERNAL           DLAMCH, DLANSY, DLA_SYRCOND
       DOUBLE PRECISION   DLAMCH, DLANSY, DLA_SYRCOND        DOUBLE PRECISION   DLAMCH, DLANSY, DLA_SYRCOND
       LOGICAL            LSAME        LOGICAL            LSAME
       INTEGER            BLAS_FPINFO_X        INTEGER            ILAPREC
       INTEGER            ILATRANS, ILAPREC  
 *     ..  *     ..
 *     .. Executable Statements ..  *     .. Executable Statements ..
 *  *
Line 534 Line 530
       ELSE IF( LDAF.LT.MAX( 1, N ) ) THEN        ELSE IF( LDAF.LT.MAX( 1, N ) ) THEN
         INFO = -8          INFO = -8
       ELSE IF( LDB.LT.MAX( 1, N ) ) THEN        ELSE IF( LDB.LT.MAX( 1, N ) ) THEN
         INFO = -11          INFO = -12
       ELSE IF( LDX.LT.MAX( 1, N ) ) THEN        ELSE IF( LDX.LT.MAX( 1, N ) ) THEN
         INFO = -13          INFO = -14
       END IF        END IF
       IF( INFO.NE.0 ) THEN        IF( INFO.NE.0 ) THEN
         CALL XERBLA( 'DSYRFSX', -INFO )          CALL XERBLA( 'DSYRFSX', -INFO )
Line 677 Line 673
             IF ( RCOND_TMP .LT. ILLRCOND_THRESH ) THEN              IF ( RCOND_TMP .LT. ILLRCOND_THRESH ) THEN
                ERR_BNDS_COMP( J, LA_LINRX_ERR_I ) = 1.0D+0                 ERR_BNDS_COMP( J, LA_LINRX_ERR_I ) = 1.0D+0
                ERR_BNDS_COMP( J, LA_LINRX_TRUST_I ) = 0.0D+0                 ERR_BNDS_COMP( J, LA_LINRX_TRUST_I ) = 0.0D+0
                IF ( PARAMS( LA_LINRX_CWISE_I ) .EQ. 1.0D+0                 IF ( .NOT. IGNORE_CWISE
      $              .AND. INFO.LT.N + J ) INFO = N + J       $              .AND. INFO.LT.N + J ) INFO = N + J
             ELSE IF ( ERR_BNDS_COMP( J, LA_LINRX_ERR_I )              ELSE IF ( ERR_BNDS_COMP( J, LA_LINRX_ERR_I )
      $              .LT. ERR_LBND ) THEN       $              .LT. ERR_LBND ) THEN
Removed from v.1.6  
changed lines
  Added in v.1.16

CVSweb interface <joel.bertrand@systella.fr>