Diff for /rpl/lapack/lapack/dlarrf.f between versions 1.10 and 1.21

version 1.10, 2011/11/21 22:19:34 version 1.21, 2018/05/29 06:55:19
Line 1 Line 1
 *> \brief \b DLARRF  *> \brief \b DLARRF finds a new relatively robust representation such that at least one of the eigenvalues is relatively isolated.
 *  *
 *  =========== 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 DLARRF + dependencies   *> Download DLARRF + dependencies
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dlarrf.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dlarrf.f">
 *> [TGZ]</a>   *> [TGZ]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dlarrf.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dlarrf.f">
 *> [ZIP]</a>   *> [ZIP]</a>
 *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dlarrf.f">   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dlarrf.f">
 *> [TXT]</a>  *> [TXT]</a>
 *> \endhtmlonly   *> \endhtmlonly
 *  *
 *  Definition:  *  Definition:
 *  ===========  *  ===========
Line 22 Line 22
 *                          W, WGAP, WERR,  *                          W, WGAP, WERR,
 *                          SPDIAM, CLGAPL, CLGAPR, PIVMIN, SIGMA,  *                          SPDIAM, CLGAPL, CLGAPR, PIVMIN, SIGMA,
 *                          DPLUS, LPLUS, WORK, INFO )  *                          DPLUS, LPLUS, WORK, INFO )
 *   *
 *       .. Scalar Arguments ..  *       .. Scalar Arguments ..
 *       INTEGER            CLSTRT, CLEND, INFO, N  *       INTEGER            CLSTRT, CLEND, INFO, N
 *       DOUBLE PRECISION   CLGAPL, CLGAPR, PIVMIN, SIGMA, SPDIAM  *       DOUBLE PRECISION   CLGAPL, CLGAPR, PIVMIN, SIGMA, SPDIAM
Line 31 Line 31
 *       DOUBLE PRECISION   D( * ), DPLUS( * ), L( * ), LD( * ),  *       DOUBLE PRECISION   D( * ), DPLUS( * ), L( * ), LD( * ),
 *      $          LPLUS( * ), W( * ), WGAP( * ), WERR( * ), WORK( * )  *      $          LPLUS( * ), W( * ), WGAP( * ), WERR( * ), WORK( * )
 *       ..  *       ..
 *    *
 *  *
 *> \par Purpose:  *> \par Purpose:
 *  =============  *  =============
Line 51 Line 51
 *> \param[in] N  *> \param[in] N
 *> \verbatim  *> \verbatim
 *>          N is INTEGER  *>          N is INTEGER
 *>          The order of the matrix (subblock, if the matrix splitted).  *>          The order of the matrix (subblock, if the matrix split).
 *> \endverbatim  *> \endverbatim
 *>  *>
 *> \param[in] D  *> \param[in] D
Line 169 Line 169
 *  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  *> \date June 2016
 *  *
 *> \ingroup auxOTHERauxiliary  *> \ingroup OTHERauxiliary
 *  *
 *> \par Contributors:  *> \par Contributors:
 *  ==================  *  ==================
Line 193 Line 193
      $                   SPDIAM, CLGAPL, CLGAPR, PIVMIN, SIGMA,       $                   SPDIAM, CLGAPL, CLGAPR, PIVMIN, SIGMA,
      $                   DPLUS, LPLUS, WORK, INFO )       $                   DPLUS, LPLUS, WORK, INFO )
 *  *
 *  -- LAPACK auxiliary routine (version 3.4.0) --  *  -- LAPACK auxiliary routine (version 3.7.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  *     June 2016
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            CLSTRT, CLEND, INFO, N        INTEGER            CLSTRT, CLEND, INFO, N
Line 210 Line 210
 *  =====================================================================  *  =====================================================================
 *  *
 *     .. Parameters ..  *     .. Parameters ..
       DOUBLE PRECISION   FOUR, MAXGROWTH1, MAXGROWTH2, ONE, QUART, TWO,        DOUBLE PRECISION   FOUR, MAXGROWTH1, MAXGROWTH2, ONE, QUART, TWO
      $                   ZERO        PARAMETER          ( ONE = 1.0D0, TWO = 2.0D0, FOUR = 4.0D0,
       PARAMETER          ( ZERO = 0.0D0, ONE = 1.0D0, TWO = 2.0D0,       $                     QUART = 0.25D0,
      $                     FOUR = 4.0D0, QUART = 0.25D0,  
      $                     MAXGROWTH1 = 8.D0,       $                     MAXGROWTH1 = 8.D0,
      $                     MAXGROWTH2 = 8.D0 )       $                     MAXGROWTH2 = 8.D0 )
 *     ..  *     ..
Line 240 Line 239
 *     .. Executable Statements ..  *     .. Executable Statements ..
 *  *
       INFO = 0        INFO = 0
   *
   *     Quick return if possible
   *
         IF( N.LE.0 ) THEN
            RETURN
         END IF
   *
       FACT = DBLE(2**KTRYMAX)        FACT = DBLE(2**KTRYMAX)
       EPS = DLAMCH( 'Precision' )        EPS = DLAMCH( 'Precision' )
       SHIFT = 0        SHIFT = 0
Line 258 Line 264
   
 *     Decide whether the code should accept the best among all  *     Decide whether the code should accept the best among all
 *     representations despite large element growth or signal INFO=1  *     representations despite large element growth or signal INFO=1
       NOFAIL = .TRUE.  *     Setting NOFAIL to .FALSE. for quick fix for bug 113
         NOFAIL = .FALSE.
 *  *
   
 *     Compute the average gap length of the cluster  *     Compute the average gap length of the cluster

Removed from v.1.10  
changed lines
  Added in v.1.21


CVSweb interface <joel.bertrand@systella.fr>