Diff for /rpl/lapack/lapack/zlabrd.f between versions 1.18 and 1.20

version 1.18, 2018/05/29 06:55:24 version 1.20, 2023/08/07 08:39:28
Line 157 Line 157
 *> \author Univ. of Colorado Denver  *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.  *> \author NAG Ltd.
 *  *
 *> \date June 2017  
 *  
 *> \ingroup complex16OTHERauxiliary  *> \ingroup complex16OTHERauxiliary
 *  *
 *> \par Further Details:  *> \par Further Details:
Line 212 Line 210
       SUBROUTINE ZLABRD( M, N, NB, A, LDA, D, E, TAUQ, TAUP, X, LDX, Y,        SUBROUTINE ZLABRD( M, N, NB, A, LDA, D, E, TAUQ, TAUP, X, LDX, Y,
      $                   LDY )       $                   LDY )
 *  *
 *  -- LAPACK auxiliary routine (version 3.7.1) --  *  -- 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..--
 *     June 2017  
 *  *
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       INTEGER            LDA, LDX, LDY, M, N, NB        INTEGER            LDA, LDX, LDY, M, N, NB
Line 270 Line 267
             ALPHA = A( I, I )              ALPHA = A( I, I )
             CALL ZLARFG( M-I+1, ALPHA, A( MIN( I+1, M ), I ), 1,              CALL ZLARFG( M-I+1, ALPHA, A( MIN( I+1, M ), I ), 1,
      $                   TAUQ( I ) )       $                   TAUQ( I ) )
             D( I ) = ALPHA              D( I ) = DBLE( ALPHA )
             IF( I.LT.N ) THEN              IF( I.LT.N ) THEN
                A( I, I ) = ONE                 A( I, I ) = ONE
 *  *
Line 310 Line 307
                ALPHA = A( I, I+1 )                 ALPHA = A( I, I+1 )
                CALL ZLARFG( N-I, ALPHA, A( I, MIN( I+2, N ) ), LDA,                 CALL ZLARFG( N-I, ALPHA, A( I, MIN( I+2, N ) ), LDA,
      $                      TAUP( I ) )       $                      TAUP( I ) )
                E( I ) = ALPHA                 E( I ) = DBLE( ALPHA )
                A( I, I+1 ) = ONE                 A( I, I+1 ) = ONE
 *  *
 *              Compute X(i+1:m,i)  *              Compute X(i+1:m,i)
Line 354 Line 351
             ALPHA = A( I, I )              ALPHA = A( I, I )
             CALL ZLARFG( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,              CALL ZLARFG( N-I+1, ALPHA, A( I, MIN( I+1, N ) ), LDA,
      $                   TAUP( I ) )       $                   TAUP( I ) )
             D( I ) = ALPHA              D( I ) = DBLE( ALPHA )
             IF( I.LT.M ) THEN              IF( I.LT.M ) THEN
                A( I, I ) = ONE                 A( I, I ) = ONE
 *  *
Line 388 Line 385
                ALPHA = A( I+1, I )                 ALPHA = A( I+1, I )
                CALL ZLARFG( M-I, ALPHA, A( MIN( I+2, M ), I ), 1,                 CALL ZLARFG( M-I, ALPHA, A( MIN( I+2, M ), I ), 1,
      $                      TAUQ( I ) )       $                      TAUQ( I ) )
                E( I ) = ALPHA                 E( I ) = DBLE( ALPHA )
                A( I+1, I ) = ONE                 A( I+1, I ) = ONE
 *  *
 *              Compute Y(i+1:n,i)  *              Compute Y(i+1:n,i)

Removed from v.1.18  
changed lines
  Added in v.1.20


CVSweb interface <joel.bertrand@systella.fr>