Diff for /rpl/lapack/lapack/zlauu2.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:33
Line 95 Line 95
 *> \author Univ. of Colorado Denver  *> \author Univ. of Colorado Denver
 *> \author NAG Ltd.  *> \author NAG Ltd.
 *  *
 *> \date December 2016  
 *  
 *> \ingroup complex16OTHERauxiliary  *> \ingroup complex16OTHERauxiliary
 *  *
 *  =====================================================================  *  =====================================================================
       SUBROUTINE ZLAUU2( UPLO, N, A, LDA, INFO )        SUBROUTINE ZLAUU2( UPLO, N, A, LDA, INFO )
 *  *
 *  -- 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 165 Line 162
 *        Compute the product U * U**H.  *        Compute the product U * U**H.
 *  *
          DO 10 I = 1, N           DO 10 I = 1, N
             AII = A( I, I )              AII = DBLE( A( I, I ) )
             IF( I.LT.N ) THEN              IF( I.LT.N ) THEN
                A( I, I ) = AII*AII + DBLE( ZDOTC( N-I, A( I, I+1 ), LDA,                 A( I, I ) = AII*AII + DBLE( ZDOTC( N-I, A( I, I+1 ), LDA,
      $                     A( I, I+1 ), LDA ) )       $                     A( I, I+1 ), LDA ) )
Line 184 Line 181
 *        Compute the product L**H * L.  *        Compute the product L**H * L.
 *  *
          DO 20 I = 1, N           DO 20 I = 1, N
             AII = A( I, I )              AII = DBLE( A( I, I ) )
             IF( I.LT.N ) THEN              IF( I.LT.N ) THEN
                A( I, I ) = AII*AII + DBLE( ZDOTC( N-I, A( I+1, I ), 1,                 A( I, I ) = AII*AII + DBLE( ZDOTC( N-I, A( I+1, I ), 1,
      $                     A( I+1, I ), 1 ) )       $                     A( I+1, I ), 1 ) )

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


CVSweb interface <joel.bertrand@systella.fr>