diff options
| author | julie <julielangou@users.noreply.github.com> | 2011-04-02 11:08:56 +0000 |
|---|---|---|
| committer | julie <julielangou@users.noreply.github.com> | 2011-04-02 11:08:56 +0000 |
| commit | f2953573ede24d7f8c01fdb18de48f65f00a9943 (patch) | |
| tree | 53172aa9083b9aa1abe2d6c130f7c173d8d8725b /SRC/strevc.f | |
| parent | 53b71f5605f83d116ab6bcf477bfb6d2ca757de1 (diff) | |
| download | lapack-f2953573ede24d7f8c01fdb18de48f65f00a9943.tar.gz lapack-f2953573ede24d7f8c01fdb18de48f65f00a9943.tar.bz2 lapack-f2953573ede24d7f8c01fdb18de48f65f00a9943.zip | |
First pass to homgenize notation for transpose (**T) and conjugate transpose (**H)
Corresponds to bug0024
Diffstat (limited to 'SRC/strevc.f')
| -rw-r--r-- | SRC/strevc.f | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/SRC/strevc.f b/SRC/strevc.f index 1ffb629d..82fe2fb4 100644 --- a/SRC/strevc.f +++ b/SRC/strevc.f @@ -27,9 +27,9 @@ * The right eigenvector x and the left eigenvector y of T corresponding * to an eigenvalue w are defined by: * -* T*x = w*x, (y**H)*T = w*(y**H) +* T*x = w*x, (y**T)*T = w*(y**T) * -* where y**H denotes the conjugate transpose of y. +* where y**T denotes the transpose of y. * The eigenvalues are not input to this routine, but are read directly * from the diagonal blocks of T. * @@ -652,7 +652,7 @@ 160 CONTINUE * * Solve the quasi-triangular system: -* (T(KI+1:N,KI+1:N) - WR)'*X = SCALE*WORK +* (T(KI+1:N,KI+1:N) - WR)**T*X = SCALE*WORK * VMAX = ONE VCRIT = BIGNUM @@ -689,7 +689,7 @@ $ SDOT( J-KI-1, T( KI+1, J ), 1, $ WORK( KI+1+N ), 1 ) * -* Solve (T(J,J)-WR)'*X = WORK +* Solve (T(J,J)-WR)**T*X = WORK * CALL SLALN2( .FALSE., 1, 1, SMIN, ONE, T( J, J ), $ LDT, ONE, ONE, WORK( J+N ), N, WR, @@ -779,7 +779,7 @@ * Complex left eigenvector. * * Initial solve: -* ((T(KI,KI) T(KI,KI+1) )' - (WR - I* WI))*X = 0. +* ((T(KI,KI) T(KI,KI+1) )**T - (WR - I* WI))*X = 0. * ((T(KI+1,KI) T(KI+1,KI+1)) ) * IF( ABS( T( KI, KI+1 ) ).GE.ABS( T( KI+1, KI ) ) ) THEN |