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

Diff for /rpl/lapack/lapack/zstedc.f between versions 1.12 and 1.13

version 1.12, 2014/01/27 09:28:42 version 1.13, 2015/11/26 11:44:26
Line 199 Line 199
 *> \author Univ. of Colorado Denver   *> \author Univ. of Colorado Denver 
 *> \author NAG Ltd.   *> \author NAG Ltd. 
 *  *
 *> \date November 2011  *> \date November 2015
 *  *
 *> \ingroup complex16OTHERcomputational  *> \ingroup complex16OTHERcomputational
 *  *
Line 213 Line 213
       SUBROUTINE ZSTEDC( COMPZ, N, D, E, Z, LDZ, WORK, LWORK, RWORK,        SUBROUTINE ZSTEDC( COMPZ, N, D, E, Z, LDZ, WORK, LWORK, RWORK,
      $                   LRWORK, IWORK, LIWORK, INFO )       $                   LRWORK, IWORK, LIWORK, INFO )
 *  *
 *  -- LAPACK computational routine (version 3.4.0) --  *  -- LAPACK computational routine (version 3.6.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 2011  *     November 2015
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          COMPZ        CHARACTER          COMPZ
Line 454 Line 454
 *  *
 *        endwhile  *        endwhile
 *  *
 *        If the problem split any number of times, then the eigenvalues  *
 *        will not be properly ordered.  Here we permute the eigenvalues  *        Use Selection Sort to minimize swaps of eigenvectors
 *        (and the associated eigenvectors) into ascending order.  *
 *           DO 60 II = 2, N
          IF( M.NE.N ) THEN             I = II - 1
 *             K = I
 *           Use Selection Sort to minimize swaps of eigenvectors             P = D( I )
 *             DO 50 J = II, N
             DO 60 II = 2, N                IF( D( J ).LT.P ) THEN
                I = II - 1                   K = J
                K = I                   P = D( J )
                P = D( I )                END IF
                DO 50 J = II, N     50      CONTINUE
                   IF( D( J ).LT.P ) THEN             IF( K.NE.I ) THEN
                      K = J                D( K ) = D( I )
                      P = D( J )                D( I ) = P
                   END IF                CALL ZSWAP( N, Z( 1, I ), 1, Z( 1, K ), 1 )
    50          CONTINUE             END IF
                IF( K.NE.I ) THEN     60    CONTINUE
                   D( K ) = D( I )  
                   D( I ) = P  
                   CALL ZSWAP( N, Z( 1, I ), 1, Z( 1, K ), 1 )  
                END IF  
    60       CONTINUE  
          END IF  
       END IF        END IF
 *  *
    70 CONTINUE     70 CONTINUE
Removed from v.1.12  
changed lines
  Added in v.1.13

CVSweb interface <joel.bertrand@systella.fr>