--- rpl/lapack/lapack/ztrevc3.f 2016/08/27 15:31:05 1.1 +++ rpl/lapack/lapack/ztrevc3.f 2023/08/07 08:39:42 1.7 @@ -2,24 +2,24 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * *> \htmlonly -*> Download ZTREVC3 + dependencies -*> -*> [TGZ] -*> -*> [ZIP] -*> +*> Download ZTREVC3 + dependencies +*> +*> [TGZ] +*> +*> [ZIP] +*> *> [TXT] -*> \endhtmlonly +*> \endhtmlonly * * Definition: * =========== * -* SUBROUTINE ZTREVC3( SIDE, HOWMNY, SELECT, N, T, LDT, VL, LDVL, -* VR, LDVR, MM, M, WORK, LWORK, RWORK, INFO ) +* SUBROUTINE ZTREVC3( SIDE, HOWMNY, SELECT, N, T, LDT, VL, LDVL, VR, +* $ LDVR, MM, M, WORK, LWORK, RWORK, LRWORK, INFO) * * .. Scalar Arguments .. * CHARACTER HOWMNY, SIDE @@ -222,10 +222,6 @@ *> \author Univ. of Colorado Denver *> \author NAG Ltd. * -*> \date November 2011 -* -* @precisions fortran z -> c -* *> \ingroup complex16OTHERcomputational * *> \par Further Details: @@ -247,10 +243,9 @@ $ LDVR, MM, M, WORK, LWORK, RWORK, LRWORK, INFO) IMPLICIT NONE * -* -- LAPACK computational routine (version 3.4.0) -- +* -- LAPACK computational routine -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* November 2011 * * .. Scalar Arguments .. CHARACTER HOWMNY, SIDE @@ -287,16 +282,17 @@ EXTERNAL LSAME, ILAENV, IZAMAX, DLAMCH, DZASUM * .. * .. External Subroutines .. - EXTERNAL XERBLA, ZCOPY, ZDSCAL, ZGEMV, ZLATRS + EXTERNAL XERBLA, ZCOPY, ZDSCAL, ZGEMV, ZLATRS, + $ ZGEMM, DLABAD, ZLASET, ZLACPY * .. * .. Intrinsic Functions .. - INTRINSIC ABS, DBLE, DCMPLX, CONJG, AIMAG, MAX + INTRINSIC ABS, DBLE, DCMPLX, CONJG, DIMAG, MAX * .. * .. Statement Functions .. DOUBLE PRECISION CABS1 * .. * .. Statement Function definitions .. - CABS1( CDUM ) = ABS( DBLE( CDUM ) ) + ABS( AIMAG( CDUM ) ) + CABS1( CDUM ) = ABS( DBLE( CDUM ) ) + ABS( DIMAG( CDUM ) ) * .. * .. Executable Statements .. * @@ -593,7 +589,7 @@ * When the number of vectors stored reaches NB, * or if this was last vector, do the GEMM IF( (IV.EQ.NB) .OR. (KI.EQ.N) ) THEN - CALL ZGEMM( 'N', 'N', N, IV, N-KI+IV, ONE, + CALL ZGEMM( 'N', 'N', N, IV, N-KI+IV, CONE, $ VL( 1, KI-IV+1 ), LDVL, $ WORK( KI-IV+1 + (1)*N ), N, $ CZERO,