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

Diff for /rpl/lapack/lapack/dlarrv.f between versions 1.8 and 1.9

version 1.8, 2010/12/21 13:53:32 version 1.9, 2011/07/22 07:38:07
Line 4 Line 4
      $                   IBLOCK, INDEXW, GERS, Z, LDZ, ISUPPZ,       $                   IBLOCK, INDEXW, GERS, Z, LDZ, ISUPPZ,
      $                   WORK, IWORK, INFO )       $                   WORK, IWORK, INFO )
 *  *
 *  -- LAPACK auxiliary routine (version 3.2.2) --  *  -- LAPACK auxiliary 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..--
 *     June 2010  *  -- April 2011                                                      --
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            DOL, DOU, INFO, LDZ, M, N        INTEGER            DOL, DOU, INFO, LDZ, M, N
Line 25 Line 25
 *  =======  *  =======
 *  *
 *  DLARRV computes the eigenvectors of the tridiagonal matrix  *  DLARRV computes the eigenvectors of the tridiagonal matrix
 *  T = L D L^T given L, D and APPROXIMATIONS to the eigenvalues of L D L^T.  *  T = L D L**T given L, D and APPROXIMATIONS to the eigenvalues of L D L**T.
 *  The input eigenvalues should have been computed by DLARRE.  *  The input eigenvalues should have been computed by DLARRE.
 *  *
 *  Arguments  *  Arguments
Line 332 Line 332
 *        high relative accuracy is required for the computation of the  *        high relative accuracy is required for the computation of the
 *        corresponding eigenvectors.  *        corresponding eigenvectors.
          CALL DCOPY( IM, W( WBEGIN ), 1,           CALL DCOPY( IM, W( WBEGIN ), 1,
      &                   WORK( WBEGIN ), 1 )       $                   WORK( WBEGIN ), 1 )
   
 *        We store in W the eigenvalue approximations w.r.t. the original  *        We store in W the eigenvalue approximations w.r.t. the original
 *        matrix T.  *        matrix T.
Line 431 Line 431
                   Q = INDEXW( WBEGIN-1+OLDLST )                    Q = INDEXW( WBEGIN-1+OLDLST )
 *                 Offset for the arrays WORK, WGAP and WERR, i.e., the P-OFFSET  *                 Offset for the arrays WORK, WGAP and WERR, i.e., the P-OFFSET
 *                 through the Q-OFFSET elements of these arrays are to be used.  *                 through the Q-OFFSET elements of these arrays are to be used.
 C                  OFFSET = P-OLDFST  *                  OFFSET = P-OLDFST
                   OFFSET = INDEXW( WBEGIN ) - 1                    OFFSET = INDEXW( WBEGIN ) - 1
 *                 perform limited bisection (if necessary) to get approximate  *                 perform limited bisection (if necessary) to get approximate
 *                 eigenvalues to the precision needed.  *                 eigenvalues to the precision needed.
Line 570  C                  OFFSET = P-OLDFST Line 570  C                  OFFSET = P-OLDFST
 *                    Compute RRR of child cluster.  *                    Compute RRR of child cluster.
 *                    Note that the new RRR is stored in Z  *                    Note that the new RRR is stored in Z
 *  *
 C                    DLARRF needs LWORK = 2*N  *                    DLARRF needs LWORK = 2*N
                      CALL DLARRF( IN, D( IBEGIN ), L( IBEGIN ),                       CALL DLARRF( IN, D( IBEGIN ), L( IBEGIN ),
      $                         WORK(INDLD+IBEGIN-1),       $                         WORK(INDLD+IBEGIN-1),
      $                         NEWFST, NEWLST, WORK(WBEGIN),       $                         NEWFST, NEWLST, WORK(WBEGIN),
Removed from v.1.8  
changed lines
  Added in v.1.9

CVSweb interface <joel.bertrand@systella.fr>