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

Diff for /rpl/lapack/lapack/ztrsen.f between versions 1.7 and 1.8

version 1.7, 2010/12/21 13:53:57 version 1.8, 2011/07/22 07:38:21
Line 1 Line 1
       SUBROUTINE ZTRSEN( JOB, COMPQ, SELECT, N, T, LDT, Q, LDQ, W, M, S,        SUBROUTINE ZTRSEN( JOB, COMPQ, SELECT, N, T, LDT, Q, LDQ, W, M, S,
      $                   SEP, WORK, LWORK, INFO )       $                   SEP, WORK, LWORK, INFO )
 *  *
 *  -- LAPACK routine (version 3.2) --  *  -- LAPACK routine (version 3.3.1) --
 *  -- 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  *  -- April 2011                                                      --
 *  *
 *     Modified to call ZLACN2 in place of ZLACON, 10 Feb 03, SJH.  *     Modified to call ZLACN2 in place of ZLACON, 10 Feb 03, SJH.
 *  *
Line 119 Line 119
 *  transformation Z to move them to the top left corner of T. In other  *  transformation Z to move them to the top left corner of T. In other
 *  words, the selected eigenvalues are the eigenvalues of T11 in:  *  words, the selected eigenvalues are the eigenvalues of T11 in:
 *  *
 *                Z'*T*Z = ( T11 T12 ) n1  *          Z**H * T * Z = ( T11 T12 ) n1
 *                         (  0  T22 ) n2  *                         (  0  T22 ) n2
 *                            n1  n2  *                            n1  n2
 *  *
 *  where N = n1+n2 and Z' means the conjugate transpose of Z. The first  *  where N = n1+n2. The first
 *  n1 columns of Z span the specified invariant subspace of T.  *  n1 columns of Z span the specified invariant subspace of T.
 *  *
 *  If T has been obtained from the Schur factorization of a matrix  *  If T has been obtained from the Schur factorization of a matrix
 *  A = Q*T*Q', then the reordered Schur factorization of A is given by  *  A = Q*T*Q**H, then the reordered Schur factorization of A is given by
 *  A = (Q*Z)*(Z'*T*Z)*(Q*Z)', and the first n1 columns of Q*Z span the  *  A = (Q*Z)*(Z**H*T*Z)*(Q*Z)**H, and the first n1 columns of Q*Z span the
 *  corresponding invariant subspace of A.  *  corresponding invariant subspace of A.
 *  *
 *  The reciprocal condition number of the average of the eigenvalues of  *  The reciprocal condition number of the average of the eigenvalues of
Line 331 Line 331
      $                      IERR )       $                      IERR )
             ELSE              ELSE
 *  *
 *              Solve T11'*R - R*T22' = scale*X.  *              Solve T11**H*R - R*T22**H = scale*X.
 *  *
                CALL ZTRSYL( 'C', 'C', -1, N1, N2, T, LDT,                 CALL ZTRSYL( 'C', 'C', -1, N1, N2, T, LDT,
      $                      T( N1+1, N1+1 ), LDT, WORK, N1, SCALE,       $                      T( N1+1, N1+1 ), LDT, WORK, N1, SCALE,
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>