version 1.10, 2011/11/21 22:19:34
|
version 1.21, 2023/08/07 08:38:57
|
Line 1
|
Line 1
|
*> \brief \b DLARRJ |
*> \brief \b DLARRJ performs refinement of the initial estimates of the eigenvalues of the matrix T. |
* |
* |
* =========== 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 DLARRJ + dependencies |
*> Download DLARRJ + dependencies |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dlarrj.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/dlarrj.f"> |
*> [TGZ]</a> |
*> [TGZ]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dlarrj.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/dlarrj.f"> |
*> [ZIP]</a> |
*> [ZIP]</a> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dlarrj.f"> |
*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/dlarrj.f"> |
*> [TXT]</a> |
*> [TXT]</a> |
*> \endhtmlonly |
*> \endhtmlonly |
* |
* |
* Definition: |
* Definition: |
* =========== |
* =========== |
Line 21
|
Line 21
|
* SUBROUTINE DLARRJ( N, D, E2, IFIRST, ILAST, |
* SUBROUTINE DLARRJ( N, D, E2, IFIRST, ILAST, |
* RTOL, OFFSET, W, WERR, WORK, IWORK, |
* RTOL, OFFSET, W, WERR, WORK, IWORK, |
* PIVMIN, SPDIAM, INFO ) |
* PIVMIN, SPDIAM, INFO ) |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
* INTEGER IFIRST, ILAST, INFO, N, OFFSET |
* INTEGER IFIRST, ILAST, INFO, N, OFFSET |
* DOUBLE PRECISION PIVMIN, RTOL, SPDIAM |
* DOUBLE PRECISION PIVMIN, RTOL, SPDIAM |
Line 31
|
Line 31
|
* DOUBLE PRECISION D( * ), E2( * ), W( * ), |
* DOUBLE PRECISION D( * ), E2( * ), W( * ), |
* $ WERR( * ), WORK( * ) |
* $ WERR( * ), WORK( * ) |
* .. |
* .. |
* |
* |
* |
* |
*> \par Purpose: |
*> \par Purpose: |
* ============= |
* ============= |
Line 85
|
Line 85
|
*> RTOL is DOUBLE PRECISION |
*> RTOL is DOUBLE PRECISION |
*> Tolerance for the convergence of the bisection intervals. |
*> Tolerance for the convergence of the bisection intervals. |
*> An interval [LEFT,RIGHT] has converged if |
*> An interval [LEFT,RIGHT] has converged if |
*> RIGHT-LEFT.LT.RTOL*MAX(|LEFT|,|RIGHT|). |
*> RIGHT-LEFT < RTOL*MAX(|LEFT|,|RIGHT|). |
*> \endverbatim |
*> \endverbatim |
*> |
*> |
*> \param[in] OFFSET |
*> \param[in] OFFSET |
Line 145
|
Line 145
|
* 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 auxOTHERauxiliary |
*> \ingroup OTHERauxiliary |
* |
* |
*> \par Contributors: |
*> \par Contributors: |
* ================== |
* ================== |
Line 168
|
Line 166
|
$ RTOL, OFFSET, W, WERR, WORK, IWORK, |
$ RTOL, OFFSET, W, WERR, WORK, IWORK, |
$ PIVMIN, SPDIAM, INFO ) |
$ PIVMIN, SPDIAM, INFO ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.4.0) -- |
* -- LAPACK auxiliary 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 IFIRST, ILAST, INFO, N, OFFSET |
INTEGER IFIRST, ILAST, INFO, N, OFFSET |
Line 204
|
Line 201
|
* |
* |
INFO = 0 |
INFO = 0 |
* |
* |
|
* Quick return if possible |
|
* |
|
IF( N.LE.0 ) THEN |
|
RETURN |
|
END IF |
|
* |
MAXITR = INT( ( LOG( SPDIAM+PIVMIN )-LOG( PIVMIN ) ) / |
MAXITR = INT( ( LOG( SPDIAM+PIVMIN )-LOG( PIVMIN ) ) / |
$ LOG( TWO ) ) + 2 |
$ LOG( TWO ) ) + 2 |
* |
* |