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

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

version 1.7, 2010/12/21 13:53:43 version 1.8, 2011/07/22 07:38:14
Line 1 Line 1
       SUBROUTINE ZGEHD2( N, ILO, IHI, A, LDA, TAU, WORK, INFO )        SUBROUTINE ZGEHD2( N, ILO, IHI, A, LDA, TAU, WORK, 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 ..
       INTEGER            IHI, ILO, INFO, LDA, N        INTEGER            IHI, ILO, INFO, LDA, N
Line 16 Line 16
 *  =======  *  =======
 *  *
 *  ZGEHD2 reduces a complex general matrix A to upper Hessenberg form H  *  ZGEHD2 reduces a complex general matrix A to upper Hessenberg form H
 *  by a unitary similarity transformation:  Q' * A * Q = H .  *  by a unitary similarity transformation:  Q**H * A * Q = H .
 *  *
 *  Arguments  *  Arguments
 *  =========  *  =========
Line 63 Line 63
 *  *
 *  Each H(i) has the form  *  Each H(i) has the form
 *  *
 *     H(i) = I - tau * v * v'  *     H(i) = I - tau * v * v**H
 *  *
 *  where tau is a complex scalar, and v is a complex vector with  *  where tau is a complex scalar, and v is a complex vector with
 *  v(1:i) = 0, v(i+1) = 1 and v(ihi+1:n) = 0; v(i+2:ihi) is stored on  *  v(1:i) = 0, v(i+1) = 1 and v(ihi+1:n) = 0; v(i+2:ihi) is stored on
Line 134 Line 134
          CALL ZLARF( 'Right', IHI, IHI-I, A( I+1, I ), 1, TAU( I ),           CALL ZLARF( 'Right', IHI, IHI-I, A( I+1, I ), 1, TAU( I ),
      $               A( 1, I+1 ), LDA, WORK )       $               A( 1, I+1 ), LDA, WORK )
 *  *
 *        Apply H(i)' to A(i+1:ihi,i+1:n) from the left  *        Apply H(i)**H to A(i+1:ihi,i+1:n) from the left
 *  *
          CALL ZLARF( 'Left', IHI-I, N-I, A( I+1, I ), 1,           CALL ZLARF( 'Left', IHI-I, N-I, A( I+1, I ), 1,
      $               DCONJG( TAU( I ) ), A( I+1, I+1 ), LDA, WORK )       $               DCONJG( TAU( I ) ), A( I+1, I+1 ), LDA, WORK )
Removed from v.1.7  
changed lines
  Added in v.1.8

CVSweb interface <joel.bertrand@systella.fr>