--- rpl/lapack/lapack/dlarrr.f 2011/11/21 22:19:34 1.9 +++ rpl/lapack/lapack/dlarrr.f 2018/05/29 06:55:19 1.17 @@ -1,33 +1,33 @@ -*> \brief \b DLARRR +*> \brief \b DLARRR performs tests to decide whether the symmetric tridiagonal matrix T warrants expensive computations which guarantee high relative accuracy in the eigenvalues. * * =========== 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 DLARRR + dependencies -*> -*> [TGZ] -*> -*> [ZIP] -*> +*> Download DLARRR + dependencies +*> +*> [TGZ] +*> +*> [ZIP] +*> *> [TXT] -*> \endhtmlonly +*> \endhtmlonly * * Definition: * =========== * * SUBROUTINE DLARRR( N, D, E, INFO ) -* +* * .. Scalar Arguments .. * INTEGER N, INFO * .. * .. Array Arguments .. * DOUBLE PRECISION D( * ), E( * ) * .. -* -* +* +* * *> \par Purpose: * ============= @@ -73,14 +73,14 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * -*> \date November 2011 +*> \date June 2017 * -*> \ingroup auxOTHERauxiliary +*> \ingroup OTHERauxiliary * *> \par Contributors: * ================== @@ -94,10 +94,10 @@ * ===================================================================== SUBROUTINE DLARRR( N, D, E, 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, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* November 2011 +* June 2017 * * .. Scalar Arguments .. INTEGER N, INFO @@ -130,6 +130,13 @@ * .. * .. Executable Statements .. * +* Quick return if possible +* + IF( N.LE.0 ) THEN + INFO = 0 + RETURN + END IF +* * As a default, do NOT go for relative-accuracy preserving computations. INFO = 1