version 1.12, 2014/01/27 09:28:42
|
version 1.14, 2016/08/27 15:35:06
|
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 |