--- rpl/lapack/lapack/zgesvj.f 2017/06/17 11:06:44 1.5 +++ rpl/lapack/lapack/zgesvj.f 2018/05/29 06:55:22 1.6 @@ -52,7 +52,7 @@ * *> \param[in] JOBA *> \verbatim -*> JOBA is CHARACTER* 1 +*> JOBA is CHARACTER*1 *> Specifies the structure of A. *> = 'L': The input matrix A is lower triangular; *> = 'U': The input matrix A is upper triangular; @@ -206,7 +206,7 @@ *> *> \param[in,out] CWORK *> \verbatim -*> CWORK is COMPLEX*16 array, dimension max(1,LWORK). +*> CWORK is COMPLEX*16 array, dimension (max(1,LWORK)) *> Used as workspace. *> If on entry LWORK .EQ. -1, then a workspace query is assumed and *> no computation is done; CWORK(1) is set to the minial (and optimal) @@ -221,7 +221,7 @@ *> *> \param[in,out] RWORK *> \verbatim -*> RWORK is DOUBLE PRECISION array, dimension max(6,LRWORK). +*> RWORK is DOUBLE PRECISION array, dimension (max(6,LRWORK)) *> On entry, *> If JOBU .EQ. 'C' : *> RWORK(1) = CTOL, where CTOL defines the threshold for convergence. @@ -318,6 +318,8 @@ *> \par References: * ================ *> +*> \verbatim +*> *> [1] P. P. M. De Rijk: A one-sided Jacobi algorithm for computing the *> singular value decomposition on a vector computer. *> SIAM J. Sci. Stat. Comp., Vol. 10 (1998), pp. 359-371. @@ -349,7 +351,7 @@ SUBROUTINE ZGESVJ( JOBA, JOBU, JOBV, M, N, A, LDA, SVA, MV, V, $ LDV, CWORK, LWORK, RWORK, LRWORK, INFO ) * -* -- LAPACK computational routine (version 3.7.0) -- +* -- LAPACK computational routine (version 3.8.0) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- * June 2016 @@ -407,7 +409,7 @@ * .. External Subroutines .. * .. * from BLAS - EXTERNAL ZCOPY, ZROT, ZDSCAL, ZSWAP + EXTERNAL ZCOPY, ZROT, ZDSCAL, ZSWAP, ZAXPY * from LAPACK EXTERNAL DLASCL, ZLASCL, ZLASET, ZLASSQ, XERBLA EXTERNAL ZGSVJ0, ZGSVJ1