diff options
| author | julie <julielangou@users.noreply.github.com> | 2012-08-03 18:17:41 +0000 |
|---|---|---|
| committer | julie <julielangou@users.noreply.github.com> | 2012-08-03 18:17:41 +0000 |
| commit | d797a2b9defae2991b28ca622acce3b72ae68494 (patch) | |
| tree | 40cd4a012c75840d7bdeef63d9445de59bc744c7 /SRC | |
| parent | 199c646448f65fb375a8aaf7383c705dab58f550 (diff) | |
| download | lapack-d797a2b9defae2991b28ca622acce3b72ae68494.tar.gz lapack-d797a2b9defae2991b28ca622acce3b72ae68494.tar.bz2 lapack-d797a2b9defae2991b28ca622acce3b72ae68494.zip | |
error in comment detected by Mathieu (PLASMA)
Routine affected: [CZ]HEGST et [CZ]HEGS2,
B is IN in comment, but in the code, B is actually changed, but after restored.
Here is the code portion that is to blame. [CZ]HEGS2
00203 IF( K.LT.N ) THEN
00204 CALL ZDSCAL( N-K, ONE / BKK, A( K, K+1 ), LDA )
00205 CT = -HALF*AKK
00206 CALL ZLACGV( N-K, A( K, K+1 ), LDA )
00207 CALL ZLACGV( N-K, B( K, K+1 ), LDB )
00208 CALL ZAXPY( N-K, CT, B( K, K+1 ), LDB, A( K, K+1 ),
00209 $ LDA )
00210 CALL ZHER2( UPLO, N-K, -CONE, A( K, K+1 ), LDA,
00211 $ B( K, K+1 ), LDB, A( K+1, K+1 ), LDA )
00212 CALL ZAXPY( N-K, CT, B( K, K+1 ), LDB, A( K, K+1 ),
00213 $ LDA )
00214 CALL ZLACGV( N-K, B( K, K+1 ), LDB )
00215 CALL ZTRSV( UPLO, 'Conjugate transpose', 'Non-unit',
00216 $ N-K, B( K+1, K+1 ), LDB, A( K, K+1 ),
00217 $ LDA )
00218 CALL ZLACGV( N-K, A( K, K+1 ), LDA )
Diffstat (limited to 'SRC')
| -rw-r--r-- | SRC/chegs2.f | 2 | ||||
| -rw-r--r-- | SRC/chegst.f | 2 | ||||
| -rw-r--r-- | SRC/zhegs2.f | 2 | ||||
| -rw-r--r-- | SRC/zhegst.f | 2 |
4 files changed, 4 insertions, 4 deletions
diff --git a/SRC/chegs2.f b/SRC/chegs2.f index 9c59ecce..421e21d1 100644 --- a/SRC/chegs2.f +++ b/SRC/chegs2.f @@ -92,7 +92,7 @@ *> The leading dimension of the array A. LDA >= max(1,N). *> \endverbatim *> -*> \param[in] B +*> \param[in,out] B *> \verbatim *> B is COMPLEX array, dimension (LDB,N) *> The triangular factor from the Cholesky factorization of B, diff --git a/SRC/chegst.f b/SRC/chegst.f index 93ab7ba5..a591dcdd 100644 --- a/SRC/chegst.f +++ b/SRC/chegst.f @@ -92,7 +92,7 @@ *> The leading dimension of the array A. LDA >= max(1,N). *> \endverbatim *> -*> \param[in] B +*> \param[in,out] B *> \verbatim *> B is COMPLEX array, dimension (LDB,N) *> The triangular factor from the Cholesky factorization of B, diff --git a/SRC/zhegs2.f b/SRC/zhegs2.f index 95b1e075..f131c3aa 100644 --- a/SRC/zhegs2.f +++ b/SRC/zhegs2.f @@ -92,7 +92,7 @@ *> The leading dimension of the array A. LDA >= max(1,N). *> \endverbatim *> -*> \param[in] B +*> \param[in,out] B *> \verbatim *> B is COMPLEX*16 array, dimension (LDB,N) *> The triangular factor from the Cholesky factorization of B, diff --git a/SRC/zhegst.f b/SRC/zhegst.f index 93dc9485..314fbec9 100644 --- a/SRC/zhegst.f +++ b/SRC/zhegst.f @@ -92,7 +92,7 @@ *> The leading dimension of the array A. LDA >= max(1,N). *> \endverbatim *> -*> \param[in] B +*> \param[in,out] B *> \verbatim *> B is COMPLEX*16 array, dimension (LDB,N) *> The triangular factor from the Cholesky factorization of B, |