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

Diff for /rpl/lapack/lapack/dlaqr5.f between versions 1.6 and 1.7

version 1.6, 2010/08/13 21:03:50 version 1.7, 2010/12/21 13:48:05
Line 2 Line 2
      $                   SR, SI, H, LDH, ILOZ, IHIZ, Z, LDZ, V, LDV, U,       $                   SR, SI, H, LDH, ILOZ, IHIZ, Z, LDZ, V, LDV, U,
      $                   LDU, NV, WV, LDWV, NH, WH, LDWH )       $                   LDU, NV, WV, LDWV, NH, WH, LDWH )
 *  *
 *  -- LAPACK auxiliary routine (version 3.2) --  *  -- LAPACK auxiliary routine (version 3.3.0) --
 *  -- 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 2006  *     November 2010
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            IHIZ, ILOZ, KACC22, KBOT, KTOP, LDH, LDU, LDV,        INTEGER            IHIZ, ILOZ, KACC22, KBOT, KTOP, LDH, LDU, LDV,
Line 453 Line 453
 *           ==== Special case: 2-by-2 reflection (if needed) ====  *           ==== Special case: 2-by-2 reflection (if needed) ====
 *  *
             K = KRCOL + 3*( M22-1 )              K = KRCOL + 3*( M22-1 )
             IF( BMP22 .AND. ( V( 1, M22 ).NE.ZERO ) ) THEN              IF( BMP22 ) THEN
                DO 100 J = JTOP, MIN( KBOT, K+3 )                 IF ( V( 1, M22 ).NE.ZERO ) THEN
                   REFSUM = V( 1, M22 )*( H( J, K+1 )+V( 2, M22 )*                    DO 100 J = JTOP, MIN( KBOT, K+3 )
      $                     H( J, K+2 ) )                       REFSUM = V( 1, M22 )*( H( J, K+1 )+V( 2, M22 )*
                   H( J, K+1 ) = H( J, K+1 ) - REFSUM       $                        H( J, K+2 ) )
                   H( J, K+2 ) = H( J, K+2 ) - REFSUM*V( 2, M22 )                       H( J, K+1 ) = H( J, K+1 ) - REFSUM
   100          CONTINUE                       H( J, K+2 ) = H( J, K+2 ) - REFSUM*V( 2, M22 )
 *    100             CONTINUE
                IF( ACCUM ) THEN  *
                   KMS = K - INCOL                    IF( ACCUM ) THEN
                   DO 110 J = MAX( 1, KTOP-INCOL ), KDU                       KMS = K - INCOL
                      REFSUM = V( 1, M22 )*( U( J, KMS+1 )+V( 2, M22 )*                       DO 110 J = MAX( 1, KTOP-INCOL ), KDU
      $                        U( J, KMS+2 ) )                          REFSUM = V( 1, M22 )*( U( J, KMS+1 )+
                      U( J, KMS+1 ) = U( J, KMS+1 ) - REFSUM       $                           V( 2, M22 )*U( J, KMS+2 ) )
                      U( J, KMS+2 ) = U( J, KMS+2 ) - REFSUM*V( 2, M22 )                          U( J, KMS+1 ) = U( J, KMS+1 ) - REFSUM
                           U( J, KMS+2 ) = U( J, KMS+2 ) -
        $                                  REFSUM*V( 2, M22 )
   110             CONTINUE    110             CONTINUE
                ELSE IF( WANTZ ) THEN                    ELSE IF( WANTZ ) THEN
                   DO 120 J = ILOZ, IHIZ                       DO 120 J = ILOZ, IHIZ
                      REFSUM = V( 1, M22 )*( Z( J, K+1 )+V( 2, M22 )*                          REFSUM = V( 1, M22 )*( Z( J, K+1 )+V( 2, M22 )*
      $                        Z( J, K+2 ) )       $                           Z( J, K+2 ) )
                      Z( J, K+1 ) = Z( J, K+1 ) - REFSUM                          Z( J, K+1 ) = Z( J, K+1 ) - REFSUM
                      Z( J, K+2 ) = Z( J, K+2 ) - REFSUM*V( 2, M22 )                          Z( J, K+2 ) = Z( J, K+2 ) - REFSUM*V( 2, M22 )
   120             CONTINUE    120                CONTINUE
                     END IF
                END IF                 END IF
             END IF              END IF
 *  *
Removed from v.1.6  
changed lines
  Added in v.1.7

CVSweb interface <joel.bertrand@systella.fr>