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

version 1.7, 2010/12/21 13:53:34 version 1.8, 2011/07/22 07:38:08
Line 1 Line 1
       SUBROUTINE DLAUU2( UPLO, N, A, LDA, INFO )        SUBROUTINE DLAUU2( UPLO, N, A, LDA, INFO )
 *  *
 *  -- LAPACK auxiliary routine (version 3.2) --  *  -- LAPACK auxiliary 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 16 Line 16
 *  Purpose  *  Purpose
 *  =======  *  =======
 *  *
 *  DLAUU2 computes the product U * U' or L' * L, where the triangular  *  DLAUU2 computes the product U * U**T or L**T * L, where the triangular
 *  factor U or L is stored in the upper or lower triangular part of  *  factor U or L is stored in the upper or lower triangular part of
 *  the array A.  *  the array A.
 *  *
Line 42 Line 42
 *  A       (input/output) DOUBLE PRECISION array, dimension (LDA,N)  *  A       (input/output) DOUBLE PRECISION array, dimension (LDA,N)
 *          On entry, the triangular factor U or L.  *          On entry, the triangular factor U or L.
 *          On exit, if UPLO = 'U', the upper triangle of A is  *          On exit, if UPLO = 'U', the upper triangle of A is
 *          overwritten with the upper triangle of the product U * U';  *          overwritten with the upper triangle of the product U * U**T;
 *          if UPLO = 'L', the lower triangle of A is overwritten with  *          if UPLO = 'L', the lower triangle of A is overwritten with
 *          the lower triangle of the product L' * L.  *          the lower triangle of the product L**T * L.
 *  *
 *  LDA     (input) INTEGER  *  LDA     (input) INTEGER
 *          The leading dimension of the array A.  LDA >= max(1,N).  *          The leading dimension of the array A.  LDA >= max(1,N).
Line 100 Line 100
 *  *
       IF( UPPER ) THEN        IF( UPPER ) THEN
 *  *
 *        Compute the product U * U'.  *        Compute the product U * U**T.
 *  *
          DO 10 I = 1, N           DO 10 I = 1, N
             AII = A( I, I )              AII = A( I, I )
Line 115 Line 115
 *  *
       ELSE        ELSE
 *  *
 *        Compute the product L' * L.  *        Compute the product L**T * L.
 *  *
          DO 20 I = 1, N           DO 20 I = 1, N
             AII = A( I, I )              AII = A( I, I )

Removed from v.1.7  
changed lines
  Added in v.1.8


CVSweb interface <joel.bertrand@systella.fr>