diff options
author | Julien Langou <julien.langou@ucdenver.edu> | 2016-05-23 10:41:12 -0600 |
---|---|---|
committer | Julien Langou <julien.langou@ucdenver.edu> | 2016-05-23 10:41:12 -0600 |
commit | 765ec982cc6bd6940744b02e08b81bae71065f40 (patch) | |
tree | 2aa0c917b2195b2fff7f405b8e85b224903e1272 /SRC/zgelss.f | |
parent | 6cfa4785be8368938ab4da6b701ebf784910c156 (diff) | |
download | lapack-765ec982cc6bd6940744b02e08b81bae71065f40.tar.gz lapack-765ec982cc6bd6940744b02e08b81bae71065f40.tar.bz2 lapack-765ec982cc6bd6940744b02e08b81bae71065f40.zip |
apply 1742, 1743 and 1744
------------------------------------------------------------------------
r1744 | langou | 2016-04-29 15:39:17 -0600 (Fri, 29 Apr 2016) | 15 lines
See post of Nathan Whitehead on the forum (4/29/2016)
http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=4958
Addressing second comment (out of two):
"For GELSS, there is a DUM dummy argument but it is complex, [CZ]GEBRD expects
a real array as 6th argument. I added a "REAL DUMMY(1)" in local arrays, then
passed DUMMY for 6th argument of [CZ]GEBRD."
(Use the array S instead of creating a REAL DUMMY(1) array.)
(Update the thank-you file.)
------------------------------------------------------------------------
r1743 | langou | 2016-04-29 15:12:04 -0600 (Fri, 29 Apr 2016) | 4 lines
( minor: re-indent the file zgges3.f )
------------------------------------------------------------------------
r1742 | langou | 2016-04-29 15:11:31 -0600 (Fri, 29 Apr 2016) | 14 lines
See post of Nathan Whitehead on the forum (4/29/2016)
http://icl.cs.utk.edu/lapack-forum/viewtopic.php?f=13&t=4958
Addressing first comment (out of two):
"For CGGES3, there is an argument that is complex where a real array is expected.
SRC/cgges3.f:397
I changed WORK to RWORK in the second to last argument to CHGEQZ."
(Update the thank-you file.)
Diffstat (limited to 'SRC/zgelss.f')
-rw-r--r-- | SRC/zgelss.f | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/SRC/zgelss.f b/SRC/zgelss.f index 56e58ddf..a6a522b6 100644 --- a/SRC/zgelss.f +++ b/SRC/zgelss.f @@ -285,8 +285,8 @@ * Path 1 - overdetermined or exactly determined * * Compute space needed for ZGEBRD - CALL ZGEBRD( MM, N, A, LDA, S, DUM(1), DUM(1), - $ DUM(1), DUM(1), -1, INFO ) + CALL ZGEBRD( MM, N, A, LDA, S, S, DUM(1), DUM(1), DUM(1), + $ -1, INFO ) LWORK_ZGEBRD=DUM(1) * Compute space needed for ZUNMBR CALL ZUNMBR( 'Q', 'L', 'C', MM, NRHS, N, A, LDA, DUM(1), @@ -315,8 +315,8 @@ $ -1, INFO ) LWORK_ZGELQF=DUM(1) * Compute space needed for ZGEBRD - CALL ZGEBRD( M, M, A, LDA, S, DUM(1), DUM(1), - $ DUM(1), DUM(1), -1, INFO ) + CALL ZGEBRD( M, M, A, LDA, S, S, DUM(1), DUM(1), + $ DUM(1), -1, INFO ) LWORK_ZGEBRD=DUM(1) * Compute space needed for ZUNMBR CALL ZUNMBR( 'Q', 'L', 'C', M, NRHS, N, A, LDA, @@ -346,8 +346,8 @@ * Path 2 - underdetermined * * Compute space needed for ZGEBRD - CALL ZGEBRD( M, N, A, LDA, S, DUM(1), DUM(1), - $ DUM(1), DUM(1), -1, INFO ) + CALL ZGEBRD( M, N, A, LDA, S, S, DUM(1), DUM(1), + $ DUM(1), -1, INFO ) LWORK_ZGEBRD=DUM(1) * Compute space needed for ZUNMBR CALL ZUNMBR( 'Q', 'L', 'C', M, NRHS, M, A, LDA, |