--- rpl/lapack/lapack/zlarft.f 2012/07/31 11:06:39 1.11 +++ rpl/lapack/lapack/zlarft.f 2016/08/27 15:35:00 1.18 @@ -1,4 +1,4 @@ -*> \brief \b ZLARFT +*> \brief \b ZLARFT forms the triangular factor T of a block reflector H = I - vtvH * * =========== DOCUMENTATION =========== * @@ -130,7 +130,7 @@ *> \author Univ. of Colorado Denver *> \author NAG Ltd. * -*> \date April 2012 +*> \date June 2016 * *> \ingroup complex16OTHERauxiliary * @@ -163,10 +163,10 @@ * ===================================================================== SUBROUTINE ZLARFT( DIRECT, STOREV, N, K, V, LDV, TAU, T, LDT ) * -* -- LAPACK auxiliary routine (version 3.4.1) -- +* -- LAPACK auxiliary routine (version 3.6.1) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* April 2012 +* June 2016 * * .. Scalar Arguments .. CHARACTER DIRECT, STOREV @@ -187,7 +187,7 @@ INTEGER I, J, PREVLASTV, LASTV * .. * .. External Subroutines .. - EXTERNAL ZGEMV, ZLACGV, ZTRMV + EXTERNAL ZGEMV, ZLACGV, ZTRMV, ZGEMM * .. * .. External Functions .. LOGICAL LSAME @@ -303,7 +303,7 @@ * CALL ZGEMM( 'N', 'C', K-I, 1, N-K+I-J, -TAU( I ), $ V( I+1, J ), LDV, V( I, J ), LDV, - $ ONE, T( I+1, I ), LDT ) + $ ONE, T( I+1, I ), LDT ) END IF * * T(i+1:k,i) := T(i+1:k,i+1:k) * T(i+1:k,i)