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

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

version 1.6, 2010/08/13 21:04:10 version 1.7, 2010/12/21 13:48:06
Line 2 Line 2
      $                   H, LDH, ILOZ, IHIZ, Z, LDZ, V, LDV, U, LDU, NV,       $                   H, LDH, ILOZ, IHIZ, Z, LDZ, V, LDV, U, LDU, NV,
      $                   WV, LDWV, NH, WH, LDWH )       $                   WV, LDWV, NH, WH, LDWH )
 *  *
 *  -- LAPACK auxiliary routine (version 3.2) --  *  -- LAPACK auxiliary routine (version 3.3.0) --
 *     Univ. of Tennessee, Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..  *     Univ. of Tennessee, 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 441 Line 441
 *           ==== 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 90 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 90 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 ) -                       H( J, K+1 ) = H( J, K+1 ) - REFSUM
      $                          REFSUM*DCONJG( V( 2, M22 ) )                       H( J, K+2 ) = H( J, K+2 ) -
    90          CONTINUE  
 *  
                IF( ACCUM ) THEN  
                   KMS = K - INCOL  
                   DO 100 J = MAX( 1, KTOP-INCOL ), KDU  
                      REFSUM = V( 1, M22 )*( U( J, KMS+1 )+V( 2, M22 )*  
      $                        U( J, KMS+2 ) )  
                      U( J, KMS+1 ) = U( J, KMS+1 ) - REFSUM  
                      U( J, KMS+2 ) = U( J, KMS+2 ) -  
      $                               REFSUM*DCONJG( V( 2, M22 ) )  
   100             CONTINUE  
                ELSE IF( WANTZ ) THEN  
                   DO 110 J = ILOZ, IHIZ  
                      REFSUM = V( 1, M22 )*( Z( J, K+1 )+V( 2, M22 )*  
      $                        Z( J, K+2 ) )  
                      Z( J, K+1 ) = Z( J, K+1 ) - REFSUM  
                      Z( J, K+2 ) = Z( J, K+2 ) -  
      $                             REFSUM*DCONJG( V( 2, M22 ) )       $                             REFSUM*DCONJG( V( 2, M22 ) )
   110             CONTINUE     90             CONTINUE
   *
                     IF( ACCUM ) THEN
                        KMS = K - INCOL
                        DO 100 J = MAX( 1, KTOP-INCOL ), KDU
                           REFSUM = V( 1, M22 )*( U( J, KMS+1 )+
        $                           V( 2, M22 )*U( J, KMS+2 ) )
                           U( J, KMS+1 ) = U( J, KMS+1 ) - REFSUM
                           U( J, KMS+2 ) = U( J, KMS+2 ) -
        $                                  REFSUM*DCONJG( V( 2, M22 ) )
     100                CONTINUE
                     ELSE IF( WANTZ ) THEN
                        DO 110 J = ILOZ, IHIZ
                           REFSUM = V( 1, M22 )*( Z( J, K+1 )+V( 2, M22 )*
        $                           Z( J, K+2 ) )
                           Z( J, K+1 ) = Z( J, K+1 ) - REFSUM
                           Z( J, K+2 ) = Z( J, K+2 ) -
        $                                REFSUM*DCONJG( V( 2, M22 ) )
     110                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>