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

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

version 1.7, 2010/12/21 13:53:39 version 1.8, 2011/07/22 07:38:11
Line 1 Line 1
       SUBROUTINE DSYTRD( UPLO, N, A, LDA, D, E, TAU, WORK, LWORK, INFO )        SUBROUTINE DSYTRD( UPLO, N, A, LDA, D, E, TAU, 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                                                      --
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          UPLO        CHARACTER          UPLO
Line 92 Line 92
 *  *
 *  Each H(i) has the form  *  Each H(i) has the form
 *  *
 *     H(i) = I - tau * v * v'  *     H(i) = I - tau * v * v**T
 *  *
 *  where tau is a real scalar, and v is a real vector with  *  where tau is a real scalar, and v is a real vector with
 *  v(i+1:n) = 0 and v(i) = 1; v(1:i-1) is stored on exit in  *  v(i+1:n) = 0 and v(i) = 1; v(1:i-1) is stored on exit in
Line 105 Line 105
 *  *
 *  Each H(i) has the form  *  Each H(i) has the form
 *  *
 *     H(i) = I - tau * v * v'  *     H(i) = I - tau * v * v**T
 *  *
 *  where tau is a real scalar, and v is a real vector with  *  where tau is a real scalar, and v is a real vector with
 *  v(1:i) = 0 and v(i+1) = 1; v(i+2:n) is stored on exit in A(i+2:n,i),  *  v(1:i) = 0 and v(i+1) = 1; v(i+2:n) is stored on exit in A(i+2:n,i),
Line 235 Line 235
      $                   LDWORK )       $                   LDWORK )
 *  *
 *           Update the unreduced submatrix A(1:i-1,1:i-1), using an  *           Update the unreduced submatrix A(1:i-1,1:i-1), using an
 *           update of the form:  A := A - V*W' - W*V'  *           update of the form:  A := A - V*W**T - W*V**T
 *  *
             CALL DSYR2K( UPLO, 'No transpose', I-1, NB, -ONE, A( 1, I ),              CALL DSYR2K( UPLO, 'No transpose', I-1, NB, -ONE, A( 1, I ),
      $                   LDA, WORK, LDWORK, ONE, A, LDA )       $                   LDA, WORK, LDWORK, ONE, A, LDA )
Line 266 Line 266
      $                   TAU( I ), WORK, LDWORK )       $                   TAU( I ), WORK, LDWORK )
 *  *
 *           Update the unreduced submatrix A(i+ib:n,i+ib:n), using  *           Update the unreduced submatrix A(i+ib:n,i+ib:n), using
 *           an update of the form:  A := A - V*W' - W*V'  *           an update of the form:  A := A - V*W**T - W*V**T
 *  *
             CALL DSYR2K( UPLO, 'No transpose', N-I-NB+1, NB, -ONE,              CALL DSYR2K( UPLO, 'No transpose', N-I-NB+1, NB, -ONE,
      $                   A( I+NB, I ), LDA, WORK( NB+1 ), LDWORK, ONE,       $                   A( I+NB, I ), LDA, WORK( NB+1 ), LDWORK, ONE,
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>