Diff for /rpl/lapack/lapack/zlatrd.f between versions 1.17 and 1.19

version 1.17, 2017/06/17 11:06:58 version 1.19, 2023/08/07 08:39:32
Line 140 Line 140
 *> \author Univ. of Colorado Denver  *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.  *> \author NAG Ltd.
 *  *
 *> \date December 2016  
 *  
 *> \ingroup complex16OTHERauxiliary  *> \ingroup complex16OTHERauxiliary
 *  *
 *> \par Further Details:  *> \par Further Details:
Line 199 Line 197
 *  =====================================================================  *  =====================================================================
       SUBROUTINE ZLATRD( UPLO, N, NB, A, LDA, E, TAU, W, LDW )        SUBROUTINE ZLATRD( UPLO, N, NB, A, LDA, E, TAU, W, LDW )
 *  *
 *  -- LAPACK auxiliary routine (version 3.7.0) --  *  -- LAPACK auxiliary routine --
 *  -- 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..--
 *     December 2016  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          UPLO        CHARACTER          UPLO
Line 271 Line 268
 *  *
                ALPHA = A( I-1, I )                 ALPHA = A( I-1, I )
                CALL ZLARFG( I-1, ALPHA, A( 1, I ), 1, TAU( I-1 ) )                 CALL ZLARFG( I-1, ALPHA, A( 1, I ), 1, TAU( I-1 ) )
                E( I-1 ) = ALPHA                 E( I-1 ) = DBLE( ALPHA )
                A( I-1, I ) = ONE                 A( I-1, I ) = ONE
 *  *
 *              Compute W(1:i-1,i)  *              Compute W(1:i-1,i)
Line 325 Line 322
                ALPHA = A( I+1, I )                 ALPHA = A( I+1, I )
                CALL ZLARFG( N-I, ALPHA, A( MIN( I+2, N ), I ), 1,                 CALL ZLARFG( N-I, ALPHA, A( MIN( I+2, N ), I ), 1,
      $                      TAU( I ) )       $                      TAU( I ) )
                E( I ) = ALPHA                 E( I ) = DBLE( ALPHA )
                A( I+1, I ) = ONE                 A( I+1, I ) = ONE
 *  *
 *              Compute W(i+1:n,i)  *              Compute W(i+1:n,i)

Removed from v.1.17  
changed lines
  Added in v.1.19


CVSweb interface <joel.bertrand@systella.fr>