--- rpl/lapack/lapack/dlaed6.f 2012/08/22 09:48:16 1.12 +++ rpl/lapack/lapack/dlaed6.f 2015/11/26 11:44:17 1.16 @@ -1,4 +1,4 @@ -*> \brief \b DLAED6 +*> \brief \b DLAED6 used by sstedc. Computes one Newton step in solution of the secular equation. * * =========== DOCUMENTATION =========== * @@ -115,7 +115,7 @@ *> \author Univ. of Colorado Denver *> \author NAG Ltd. * -*> \date April 2012 +*> \date November 2015 * *> \ingroup auxOTHERcomputational * @@ -140,10 +140,10 @@ * ===================================================================== SUBROUTINE DLAED6( KNITER, ORGATI, RHO, D, Z, FINIT, TAU, INFO ) * -* -- LAPACK computational routine (version 3.4.1) -- +* -- LAPACK computational routine (version 3.6.0) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* April 2012 +* November 2015 * * .. Scalar Arguments .. LOGICAL ORGATI @@ -381,13 +381,14 @@ DF = DF + TEMP2 DDF = DDF + TEMP3 ELSE - GO TO 60 + GO TO 60 END IF 40 CONTINUE F = FINIT + TAU*FC ERRETM = EIGHT*( ABS( FINIT )+ABS( TAU )*ERRETM ) + $ ABS( TAU )*DF - IF( ABS( F ).LE.EPS*ERRETM ) + IF( ( ABS( F ).LE.FOUR*EPS*ERRETM ) .OR. + $ ( (UBD-LBD).LE.FOUR*EPS*ABS(TAU) ) ) $ GO TO 60 IF( F .LE. ZERO )THEN LBD = TAU