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/dgelsy.f | |
| parent | 53b71f5605f83d116ab6bcf477bfb6d2ca757de1 (diff) | |
First pass to homgenize notation for transpose (**T) and conjugate transpose (**H)
Corresponds to bug0024
Diffstat (limited to 'SRC/dgelsy.f')
| -rw-r--r-- | SRC/dgelsy.f | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/SRC/dgelsy.f b/SRC/dgelsy.f index 9b5a7844..4200e4d9 100644 --- a/SRC/dgelsy.f +++ b/SRC/dgelsy.f @@ -42,8 +42,8 @@ * A * P = Q * [ T11 0 ] * Z * [ 0 0 ] * The minimum-norm solution is then -* X = P * Z' [ inv(T11)*Q1'*B ] -* [ 0 ] +* X = P * Z**T [ inv(T11)*Q1**T*B ] +* [ 0 ] * where Q1 consists of the first RANK columns of Q. * * This routine is basically identical to the original xGELSX except @@ -325,7 +325,7 @@ * workspace: 2*MN. * Details of Householder rotations stored in WORK(MN+1:2*MN) * -* B(1:M,1:NRHS) := Q' * B(1:M,1:NRHS) +* B(1:M,1:NRHS) := Q**T * B(1:M,1:NRHS) * CALL DORMQR( 'Left', 'Transpose', M, NRHS, MN, A, LDA, WORK( 1 ), $ B, LDB, WORK( 2*MN+1 ), LWORK-2*MN, INFO ) @@ -344,7 +344,7 @@ 30 CONTINUE 40 CONTINUE * -* B(1:N,1:NRHS) := Y' * B(1:N,1:NRHS) +* B(1:N,1:NRHS) := Y**T * B(1:N,1:NRHS) * IF( RANK.LT.N ) THEN CALL DORMRZ( 'Left', 'Transpose', N, NRHS, RANK, N-RANK, A, |