Skip to content

Commit bce48d4

Browse files
martin-frbgmartin-frbg
authored
Fix typos and sytrd boundary workspace (Reference-LAPACK PR 1030)
1 parent 14a8a9a commit bce48d4

File tree

7 files changed

+27
-20
lines changed

7 files changed

+27
-20
lines changed

‎lapack-netlib/SRC/cgehrd.f‎

Lines changed: 5 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -163,7 +163,8 @@
163163
*> \endverbatim
164164
*>
165165
* =====================================================================
166-
SUBROUTINE CGEHRD( N, ILO, IHI, A, LDA, TAU, WORK, LWORK, INFO )
166+
SUBROUTINE CGEHRD( N, ILO, IHI, A, LDA, TAU, WORK, LWORK,
167+
$ INFO )
167168
*
168169
* -- LAPACK computational routine --
169170
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -193,7 +194,8 @@ SUBROUTINE CGEHRD( N, ILO, IHI, A, LDA, TAU, WORK, LWORK, INFO )
193194
COMPLEX EI
194195
* ..
195196
* .. External Subroutines ..
196-
EXTERNAL CAXPY, CGEHD2, CGEMM, CLAHR2, CLARFB, CTRMM,
197+
EXTERNAL CAXPY, CGEHD2, CGEMM, CLAHR2, CLARFB,
198+
$ CTRMM,
197199
$ XERBLA
198200
* ..
199201
* .. Intrinsic Functions ..
@@ -230,7 +232,7 @@ SUBROUTINE CGEHRD( N, ILO, IHI, A, LDA, TAU, WORK, LWORK, INFO )
230232
IF( NH.LE.1 ) THEN
231233
LWKOPT = 1
232234
ELSE
233-
NB = MIN( NBMAX, ILAENV( 1, 'DGEHRD', ' ', N, ILO, IHI,
235+
NB = MIN( NBMAX, ILAENV( 1, 'CGEHRD', ' ', N, ILO, IHI,
234236
$ -1 ) )
235237
LWKOPT = N*NB + TSIZE
236238
END IF

‎lapack-netlib/SRC/chetrd.f‎

Lines changed: 8 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -139,7 +139,7 @@
139139
*> \author Univ. of Colorado Denver
140140
*> \author NAG Ltd.
141141
*
142-
*> \ingroup complexHEcomputational
142+
*> \ingroup hetrd
143143
*
144144
*> \par Further Details:
145145
* =====================
@@ -188,7 +188,8 @@
188188
*> \endverbatim
189189
*>
190190
* =====================================================================
191-
SUBROUTINE CHETRD( UPLO, N, A, LDA, D, E, TAU, WORK, LWORK, INFO )
191+
SUBROUTINE CHETRD( UPLO, N, A, LDA, D, E, TAU, WORK, LWORK,
192+
$ INFO )
192193
*
193194
* -- LAPACK computational routine --
194195
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -225,7 +226,8 @@ SUBROUTINE CHETRD( UPLO, N, A, LDA, D, E, TAU, WORK, LWORK, INFO )
225226
* .. External Functions ..
226227
LOGICAL LSAME
227228
INTEGER ILAENV
228-
EXTERNAL LSAME, ILAENV
229+
REAL SROUNDUP_LWORK
230+
EXTERNAL LSAME, ILAENV, SROUNDUP_LWORK
229231
* ..
230232
* .. Executable Statements ..
231233
*
@@ -249,8 +251,8 @@ SUBROUTINE CHETRD( UPLO, N, A, LDA, D, E, TAU, WORK, LWORK, INFO )
249251
* Determine the block size.
250252
*
251253
NB = ILAENV( 1, 'CHETRD', UPLO, N, -1, -1, -1 )
252-
LWKOPT = N*NB
253-
WORK( 1 ) = LWKOPT
254+
LWKOPT = MAX( 1, N*NB )
255+
WORK( 1 ) = SROUNDUP_LWORK(LWKOPT)
254256
END IF
255257
*
256258
IF( INFO.NE.0 ) THEN
@@ -367,7 +369,7 @@ SUBROUTINE CHETRD( UPLO, N, A, LDA, D, E, TAU, WORK, LWORK, INFO )
367369
$ TAU( I ), IINFO )
368370
END IF
369371
*
370-
WORK( 1 ) = LWKOPT
372+
WORK( 1 ) = SROUNDUP_LWORK(LWKOPT)
371373
RETURN
372374
*
373375
* End of CHETRD

‎lapack-netlib/SRC/sgelqt.f‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -18,7 +18,7 @@
1818
*>
1919
*> \verbatim
2020
*>
21-
*> DGELQT computes a blocked LQ factorization of a real M-by-N matrix A
21+
*> SGELQT computes a blocked LQ factorization of a real M-by-N matrix A
2222
*> using the compact WY representation of Q.
2323
*> \endverbatim
2424
*
@@ -93,7 +93,7 @@
9393
*> \author Univ. of Colorado Denver
9494
*> \author NAG Ltd.
9595
*
96-
*> \ingroup doubleGEcomputational
96+
*> \ingroup gelqt
9797
*
9898
*> \par Further Details:
9999
* =====================

‎lapack-netlib/SRC/sgemlq.f‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -74,7 +74,7 @@
7474
*> A is REAL array, dimension
7575
*> (LDA,M) if SIDE = 'L',
7676
*> (LDA,N) if SIDE = 'R'
77-
*> Part of the data structure to represent Q as returned by DGELQ.
77+
*> Part of the data structure to represent Q as returned by SGELQ.
7878
*> \endverbatim
7979
*>
8080
*> \param[in] LDA

‎lapack-netlib/SRC/sgemlqt.f‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -20,7 +20,7 @@
2020
*>
2121
*> \verbatim
2222
*>
23-
*> DGEMLQT overwrites the general real M-by-N matrix C with
23+
*> SGEMLQT overwrites the general real M-by-N matrix C with
2424
*>
2525
*> SIDE = 'L' SIDE = 'R'
2626
*> TRANS = 'N': Q C C Q
@@ -145,7 +145,7 @@
145145
*> \author Univ. of Colorado Denver
146146
*> \author NAG Ltd.
147147
*
148-
*> \ingroup doubleGEcomputational
148+
*> \ingroup gemlqt
149149
*
150150
* =====================================================================
151151
SUBROUTINE SGEMLQT( SIDE, TRANS, M, N, K, MB, V, LDV, T, LDT,

‎lapack-netlib/SRC/ssytrd.f‎

Lines changed: 5 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -188,7 +188,8 @@
188188
*> \endverbatim
189189
*>
190190
* =====================================================================
191-
SUBROUTINE SSYTRD( UPLO, N, A, LDA, D, E, TAU, WORK, LWORK, INFO )
191+
SUBROUTINE SSYTRD( UPLO, N, A, LDA, D, E, TAU, WORK, LWORK,
192+
$ INFO )
192193
*
193194
* -- LAPACK computational routine --
194195
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -248,7 +249,7 @@ SUBROUTINE SSYTRD( UPLO, N, A, LDA, D, E, TAU, WORK, LWORK, INFO )
248249
* Determine the block size.
249250
*
250251
NB = ILAENV( 1, 'SSYTRD', UPLO, N, -1, -1, -1 )
251-
LWKOPT = N*NB
252+
LWKOPT = MAX( 1, N*NB )
252253
WORK( 1 ) = SROUNDUP_LWORK(LWKOPT)
253254
END IF
254255
*
@@ -316,7 +317,8 @@ SUBROUTINE SSYTRD( UPLO, N, A, LDA, D, E, TAU, WORK, LWORK, INFO )
316317
* Update the unreduced submatrix A(1:i-1,1:i-1), using an
317318
* update of the form: A := A - V*W**T - W*V**T
318319
*
319-
CALL SSYR2K( UPLO, 'No transpose', I-1, NB, -ONE, A( 1, I ),
320+
CALL SSYR2K( UPLO, 'No transpose', I-1, NB, -ONE, A( 1,
321+
$ I ),
320322
$ LDA, WORK, LDWORK, ONE, A, LDA )
321323
*
322324
* Copy superdiagonal elements back into A, and diagonal

‎lapack-netlib/SRC/zhetrd.f‎

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -139,7 +139,7 @@
139139
*> \author Univ. of Colorado Denver
140140
*> \author NAG Ltd.
141141
*
142-
*> \ingroup complex16HEcomputational
142+
*> \ingroup hetrd
143143
*
144144
*> \par Further Details:
145145
* =====================
@@ -188,7 +188,8 @@
188188
*> \endverbatim
189189
*>
190190
* =====================================================================
191-
SUBROUTINE ZHETRD( UPLO, N, A, LDA, D, E, TAU, WORK, LWORK, INFO )
191+
SUBROUTINE ZHETRD( UPLO, N, A, LDA, D, E, TAU, WORK, LWORK,
192+
$ INFO )
192193
*
193194
* -- LAPACK computational routine --
194195
* -- LAPACK is a software package provided by Univ. of Tennessee, --
@@ -249,7 +250,7 @@ SUBROUTINE ZHETRD( UPLO, N, A, LDA, D, E, TAU, WORK, LWORK, INFO )
249250
* Determine the block size.
250251
*
251252
NB = ILAENV( 1, 'ZHETRD', UPLO, N, -1, -1, -1 )
252-
LWKOPT = N*NB
253+
LWKOPT = MAX( 1, N*NB )
253254
WORK( 1 ) = LWKOPT
254255
END IF
255256
*

0 commit comments

Comments
 (0)