[BACK]Return to dhgeqz.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/dhgeqz.f between versions 1.20 and 1.21

version 1.20, 2018/05/29 07:17:54 version 1.21, 2023/08/07 08:38:52
Line 282 Line 282
 *> \author Univ. of Colorado Denver  *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.  *> \author NAG Ltd.
 *  *
 *> \date June 2016  
 *  
 *> \ingroup doubleGEcomputational  *> \ingroup doubleGEcomputational
 *  *
 *> \par Further Details:  *> \par Further Details:
Line 304 Line 302
      $                   ALPHAR, ALPHAI, BETA, Q, LDQ, Z, LDZ, WORK,       $                   ALPHAR, ALPHAI, BETA, Q, LDQ, Z, LDZ, WORK,
      $                   LWORK, INFO )       $                   LWORK, INFO )
 *  *
 *  -- LAPACK computational routine (version 3.7.0) --  *  -- LAPACK computational 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 2016  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          COMPQ, COMPZ, JOB        CHARACTER          COMPQ, COMPZ, JOB
Line 531 Line 528
 *  *
             GO TO 80              GO TO 80
          ELSE           ELSE
             IF( ABS( H( ILAST, ILAST-1 ) ).LE.ATOL ) THEN              IF( ABS( H( ILAST, ILAST-1 ) ).LE.MAX( SAFMIN, ULP*( 
        $         ABS( H( ILAST, ILAST ) ) + ABS( H( ILAST-1, ILAST-1 ) ) 
        $         ) ) ) THEN
                H( ILAST, ILAST-1 ) = ZERO                 H( ILAST, ILAST-1 ) = ZERO
                GO TO 80                 GO TO 80
             END IF              END IF
Line 551 Line 550
             IF( J.EQ.ILO ) THEN              IF( J.EQ.ILO ) THEN
                ILAZRO = .TRUE.                 ILAZRO = .TRUE.
             ELSE              ELSE
                IF( ABS( H( J, J-1 ) ).LE.ATOL ) THEN                 IF( ABS( H( J, J-1 ) ).LE.MAX( SAFMIN, ULP*( 
        $         ABS( H( J, J ) ) + ABS( H( J-1, J-1 ) ) 
        $         ) ) ) THEN
                   H( J, J-1 ) = ZERO                    H( J, J-1 ) = ZERO
                   ILAZRO = .TRUE.                    ILAZRO = .TRUE.
                ELSE                 ELSE
Removed from v.1.20  
changed lines
  Added in v.1.21

CVSweb interface <joel.bertrand@systella.fr>