Diff for /rpl/lapack/lapack/dtrsen.f between versions 1.1 and 1.8

version 1.1, 2010/01/26 15:22:45 version 1.8, 2011/07/22 07:38:13
Line 1 Line 1
       SUBROUTINE DTRSEN( JOB, COMPQ, SELECT, N, T, LDT, Q, LDQ, WR, WI,        SUBROUTINE DTRSEN( JOB, COMPQ, SELECT, N, T, LDT, Q, LDQ, WR, WI,
      $                   M, S, SEP, WORK, LWORK, IWORK, LIWORK, INFO )       $                   M, S, SEP, WORK, LWORK, IWORK, LIWORK, 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                                                      --
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          COMPQ, JOB        CHARACTER          COMPQ, JOB
Line 156 Line 156
 *  In other words, the selected eigenvalues are the eigenvalues of T11  *  In other words, the selected eigenvalues are the eigenvalues of T11
 *  in:  *  in:
 *  *
 *                Z'*T*Z = ( T11 T12 ) n1  *          Z**T * T * Z = ( T11 T12 ) n1
 *                         (  0  T22 ) n2  *                         (  0  T22 ) n2
 *                            n1  n2  *                            n1  n2
 *  *
 *  where N = n1+n2 and Z' means the transpose of Z. The first n1 columns  *  where N = n1+n2 and Z**T means the transpose of Z. The first n1 columns
 *  of Z span the specified invariant subspace of T.  *  of Z span the specified invariant subspace of T.
 *  *
 *  If T has been obtained from the real Schur factorization of a matrix  *  If T has been obtained from the real Schur factorization of a matrix
 *  A = Q*T*Q', then the reordered real Schur factorization of A is given  *  A = Q*T*Q**T, then the reordered real Schur factorization of A is given
 *  by A = (Q*Z)*(Z'*T*Z)*(Q*Z)', and the first n1 columns of Q*Z span  *  by A = (Q*Z)*(Z**T*T*Z)*(Q*Z)**T, and the first n1 columns of Q*Z span
 *  the corresponding invariant subspace of A.  *  the 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 422 Line 422
      $                      IERR )       $                      IERR )
             ELSE              ELSE
 *  *
 *              Solve  T11'*R - R*T22' = scale*X.  *              Solve T11**T*R - R*T22**T = scale*X.
 *  *
                CALL DTRSYL( 'T', 'T', -1, N1, N2, T, LDT,                 CALL DTRSYL( 'T', 'T', -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.1  
changed lines
  Added in v.1.8


CVSweb interface <joel.bertrand@systella.fr>