diff options
| author | julie <julielangou@users.noreply.github.com> | 2011-04-02 11:08:56 +0000 |
|---|---|---|
| committer | julie <julielangou@users.noreply.github.com> | 2011-04-02 11:08:56 +0000 |
| commit | f2953573ede24d7f8c01fdb18de48f65f00a9943 (patch) | |
| tree | 53172aa9083b9aa1abe2d6c130f7c173d8d8725b /SRC/dgbbrd.f | |
| parent | 53b71f5605f83d116ab6bcf477bfb6d2ca757de1 (diff) | |
| download | lapack-f2953573ede24d7f8c01fdb18de48f65f00a9943.tar.gz lapack-f2953573ede24d7f8c01fdb18de48f65f00a9943.tar.bz2 lapack-f2953573ede24d7f8c01fdb18de48f65f00a9943.zip | |
First pass to homgenize notation for transpose (**T) and conjugate transpose (**H)
Corresponds to bug0024
Diffstat (limited to 'SRC/dgbbrd.f')
| -rw-r--r-- | SRC/dgbbrd.f | 18 |
1 files changed, 9 insertions, 9 deletions
diff --git a/SRC/dgbbrd.f b/SRC/dgbbrd.f index 5cefa7f1..12f0d8f8 100644 --- a/SRC/dgbbrd.f +++ b/SRC/dgbbrd.f @@ -19,20 +19,20 @@ * ======= * * DGBBRD reduces a real general m-by-n band matrix A to upper -* bidiagonal form B by an orthogonal transformation: Q' * A * P = B. +* bidiagonal form B by an orthogonal transformation: Q**T * A * P = B. * -* The routine computes B, and optionally forms Q or P', or computes -* Q'*C for a given matrix C. +* The routine computes B, and optionally forms Q or P**T, or computes +* Q**T*C for a given matrix C. * * Arguments * ========= * * VECT (input) CHARACTER*1 -* Specifies whether or not the matrices Q and P' are to be +* Specifies whether or not the matrices Q and P**T are to be * formed. -* = 'N': do not form Q or P'; +* = 'N': do not form Q or P**T; * = 'Q': form Q only; -* = 'P': form P' only; +* = 'P': form P**T only; * = 'B': form both. * * M (input) INTEGER @@ -85,7 +85,7 @@ * * C (input/output) DOUBLE PRECISION array, dimension (LDC,NCC) * On entry, an m-by-ncc matrix C. -* On exit, C is overwritten by Q'*C. +* On exit, C is overwritten by Q**T*C. * C is not referenced if NCC = 0. * * LDC (input) INTEGER @@ -157,7 +157,7 @@ RETURN END IF * -* Initialize Q and P' to the unit matrix, if needed +* Initialize Q and P**T to the unit matrix, if needed * IF( WANTQ ) $ CALL DLASET( 'Full', M, M, ZERO, ONE, Q, LDQ ) @@ -334,7 +334,7 @@ * IF( WANTPT ) THEN * -* accumulate product of plane rotations in P' +* accumulate product of plane rotations in P**T * DO 60 J = J1, J2, KB1 CALL DROT( N, PT( J+KUN-1, 1 ), LDPT, |