Diff for /rpl/lapack/lapack/dlaqr1.f between versions 1.17 and 1.20

version 1.17, 2018/05/29 06:55:17 version 1.20, 2023/08/07 08:38:55
Line 69 Line 69
 *> \verbatim  *> \verbatim
 *>          LDH is INTEGER  *>          LDH is INTEGER
 *>              The leading dimension of H as declared in  *>              The leading dimension of H as declared in
 *>              the calling procedure.  LDH.GE.N  *>              the calling procedure.  LDH >= N
 *> \endverbatim  *> \endverbatim
 *>  *>
 *> \param[in] SR1  *> \param[in] SR1
Line 108 Line 108
 *> \author Univ. of Colorado Denver  *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.  *> \author NAG Ltd.
 *  *
 *> \date June 2017  
 *  
 *> \ingroup doubleOTHERauxiliary  *> \ingroup doubleOTHERauxiliary
 *  *
 *> \par Contributors:  *> \par Contributors:
Line 121 Line 119
 *  =====================================================================  *  =====================================================================
       SUBROUTINE DLAQR1( N, H, LDH, SR1, SI1, SR2, SI2, V )        SUBROUTINE DLAQR1( N, H, LDH, SR1, SI1, SR2, SI2, V )
 *  *
 *  -- LAPACK auxiliary routine (version 3.7.1) --  *  -- 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..--
 *     June 2017  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       DOUBLE PRECISION   SI1, SI2, SR1, SR2        DOUBLE PRECISION   SI1, SI2, SR1, SR2
Line 147 Line 144
       INTRINSIC          ABS        INTRINSIC          ABS
 *     ..  *     ..
 *     .. Executable Statements ..  *     .. Executable Statements ..
   *
   *     Quick return if possible
   *
         IF( N.NE.2 .AND. N.NE.3 ) THEN
            RETURN
         END IF
   *
       IF( N.EQ.2 ) THEN        IF( N.EQ.2 ) THEN
          S = ABS( H( 1, 1 )-SR2 ) + ABS( SI2 ) + ABS( H( 2, 1 ) )           S = ABS( H( 1, 1 )-SR2 ) + ABS( SI2 ) + ABS( H( 2, 1 ) )
          IF( S.EQ.ZERO ) THEN           IF( S.EQ.ZERO ) THEN

Removed from v.1.17  
changed lines
  Added in v.1.20


CVSweb interface <joel.bertrand@systella.fr>