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

version 1.8, 2010/12/21 13:53:24 version 1.9, 2011/07/22 07:38:04
Line 1 Line 1
       SUBROUTINE DBDSDC( UPLO, COMPQ, N, D, E, U, LDU, VT, LDVT, Q, IQ,        SUBROUTINE DBDSDC( UPLO, COMPQ, N, D, E, U, LDU, VT, LDVT, Q, IQ,
      $                   WORK, IWORK, INFO )       $                   WORK, IWORK, INFO )
 *  *
 *  -- LAPACK routine (version 3.2.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..--
 *     June 2010  *  -- April 2011                                                      --
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          COMPQ, UPLO        CHARACTER          COMPQ, UPLO
Line 77 Line 77
 *  *
 *  VT      (output) DOUBLE PRECISION array, dimension (LDVT,N)  *  VT      (output) DOUBLE PRECISION array, dimension (LDVT,N)
 *          If  COMPQ = 'I', then:  *          If  COMPQ = 'I', then:
 *             On exit, if INFO = 0, VT' contains the right singular  *             On exit, if INFO = 0, VT**T contains the right singular
 *             vectors of the bidiagonal matrix.  *             vectors of the bidiagonal matrix.
 *          For other values of COMPQ, VT is not referenced.  *          For other values of COMPQ, VT is not referenced.
 *  *
Line 287 Line 287
       CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, N, 1, D, N, IERR )        CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, N, 1, D, N, IERR )
       CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, NM1, 1, E, NM1, IERR )        CALL DLASCL( 'G', 0, 0, ORGNRM, ONE, NM1, 1, E, NM1, IERR )
 *  *
       EPS = DLAMCH( 'Epsilon' )        EPS = (0.9D+0)*DLAMCH( 'Epsilon' )
 *  *
       MLVL = INT( LOG( DBLE( N ) / DBLE( SMLSIZ+1 ) ) / LOG( TWO ) ) + 1        MLVL = INT( LOG( DBLE( N ) / DBLE( SMLSIZ+1 ) ) / LOG( TWO ) ) + 1
       SMLSZP = SMLSIZ + 1        SMLSZP = SMLSIZ + 1

Removed from v.1.8  
changed lines
  Added in v.1.9


CVSweb interface <joel.bertrand@systella.fr>