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

Diff for /rpl/lapack/lapack/dlasd6.f between versions 1.1.1.1 and 1.10

version 1.1.1.1, 2010/01/26 15:22:45 version 1.10, 2011/07/22 07:38:08
Line 3 Line 3
      $                   LDGNUM, POLES, DIFL, DIFR, Z, K, C, S, WORK,       $                   LDGNUM, POLES, DIFL, DIFR, Z, K, C, S, WORK,
      $                   IWORK, INFO )       $                   IWORK, INFO )
 *  *
 *  -- LAPACK auxiliary routine (version 3.2) --  *  -- LAPACK auxiliary routine (version 3.3.0) --
 *  -- 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  *     November 2010
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            GIVPTR, ICOMPQ, INFO, K, LDGCOL, LDGNUM, NL,        INTEGER            GIVPTR, ICOMPQ, INFO, K, LDGCOL, LDGNUM, NL,
Line 34 Line 34
 *  *
 *  DLASD6 computes the SVD as follows:  *  DLASD6 computes the SVD as follows:
 *  *
 *                ( D1(in)  0    0     0 )  *                ( D1(in)    0    0       0 )
 *    B = U(in) * (   Z1'   a   Z2'    b ) * VT(in)  *    B = U(in) * (   Z1**T   a   Z2**T    b ) * VT(in)
 *                (   0     0   D2(in) 0 )  *                (   0       0   D2(in)   0 )
 *  *
 *      = U(out) * ( D(out) 0) * VT(out)  *      = U(out) * ( D(out) 0) * VT(out)
 *  *
 *  where Z' = (Z1' a Z2' b) = u' VT', and u is a vector of dimension M  *  where Z**T = (Z1**T a Z2**T b) = u**T VT**T, and u is a vector of dimension M
 *  with ALPHA and BETA in the NL+1 and NL+2 th entries and zeros  *  with ALPHA and BETA in the NL+1 and NL+2 th entries and zeros
 *  elsewhere; and the entry b is empty if SQRE = 0.  *  elsewhere; and the entry b is empty if SQRE = 0.
 *  *
Line 191 Line 191
 *  INFO   (output) INTEGER  *  INFO   (output) INTEGER
 *          = 0:  successful exit.  *          = 0:  successful exit.
 *          < 0:  if INFO = -i, the i-th argument had an illegal value.  *          < 0:  if INFO = -i, the i-th argument had an illegal value.
 *          > 0:  if INFO = 1, an singular value did not converge  *          > 0:  if INFO = 1, a singular value did not converge
 *  *
 *  Further Details  *  Further Details
 *  ===============  *  ===============
Line 282 Line 282
       CALL DLASD8( ICOMPQ, K, D, Z, VF, VL, DIFL, DIFR, LDGNUM,        CALL DLASD8( ICOMPQ, K, D, Z, VF, VL, DIFL, DIFR, LDGNUM,
      $             WORK( ISIGMA ), WORK( IW ), INFO )       $             WORK( ISIGMA ), WORK( IW ), INFO )
 *  *
   *     Handle error returned
   *
         IF( INFO.NE.0 ) THEN
            CALL XERBLA( 'DLASD8', -INFO )
            RETURN
         END IF
   *
 *     Save the poles if ICOMPQ = 1.  *     Save the poles if ICOMPQ = 1.
 *  *
       IF( ICOMPQ.EQ.1 ) THEN        IF( ICOMPQ.EQ.1 ) THEN
Removed from v.1.1.1.1  
changed lines
  Added in v.1.10

CVSweb interface <joel.bertrand@systella.fr>