Skip to content

Fix xORBDB6 performs numerically unadvisable operations #647

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 10 commits into from
Jul 15, 2022
Merged

Fix xORBDB6 performs numerically unadvisable operations #647

Changes from 1 commit
Commits
Show all changes
10 commits
Select commit Hold shift + click to select a range
d245b4f
SORCSD2BY1: remove dead code
christoph-conrads Oct 25, 2021
bdcd890
SORCSD: fix documentation on matrix dimensions
christoph-conrads Oct 25, 2021
643e7e3
SORBDB6: fix indexing, set vectors to zero
christoph-conrads Oct 26, 2021
774e5a9
SORBDB6: document algorithm, when output is zero
christoph-conrads Oct 26, 2021
deed8f3
SORBDB6: improve numerical stability
christoph-conrads Nov 11, 2021
a015b21
SORBDB6: replace integer literal with constant
christoph-conrads May 29, 2022
94419e8
DORBDB6: improve numerical stability, documentation
christoph-conrads Nov 11, 2021
6479e0f
CUNBDB6: improve numerical stability, documentation
christoph-conrads Nov 11, 2021
54b3964
ZUNBDB6: improve numerical stability, documentation
christoph-conrads Nov 11, 2021
145b1c3
xORBDB6/xUNBDB6: fix read of uninitialized variable
christoph-conrads Jul 14, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Next Next commit
SORCSD2BY1: remove dead code 
The zeros are overwritten by the following SORGQR call.
  • Loading branch information
@christoph-conrads
christoph-conrads committed Jul 10, 2022
commit d245b4f6ef5ed18cff4ef53d75a96b49f259bc3a
3 changes: 0 additions & 3 deletions SRC/sorcsd2by1.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -676,9 +676,6 @@ SUBROUTINE SORCSD2BY1( JOBU1, JOBU2, JOBV1T, M, P, Q, X11, LDX11,
END IF
IF( WANTU1 .AND. P .GT. 0 ) THEN
CALL SCOPY( P, WORK(IORBDB), 1, U1, 1 )
DO J = 2, P
U1(1,J) = ZERO
END DO
CALL SLACPY( 'L', P-1, M-Q-1, X11(2,1), LDX11, U1(2,2),
$ LDU1 )
CALL SORGQR( P, P, M-Q, U1, LDU1, WORK(ITAUP1),
Expand Down