2-by-1 CS decomposition fixes#405
Merged
langou merged 3 commits intoReference-LAPACK:masterfrom Feb 17, 2021
Merged
Conversation
Codecov Report
@@ Coverage Diff @@
## master #405 +/- ##
==========================================
- Coverage 83.33% 83.33% -0.01%
==========================================
Files 1820 1820
Lines 170849 170857 +8
==========================================
+ Hits 142378 142384 +6
- Misses 28471 28473 +2
Continue to review full report at Codecov.
|
thijssteel
approved these changes
Jan 9, 2021
Contributor
Author
|
Hi, do you want me to rebase? |
Contributor
|
Hi Christoph, yes, if you can rebase, this would be great. Cheers, Julien. |
9f5b5a8 to
552baf1
Compare
Contributor
Author
done |
552baf1 to
1318ba1
Compare
Contributor
Author
* check LRWORK * pass actual LRWORK value to xBBCSD
The xORCSD2BY1/xUNCSD2BY1 output matrix U2 was clearly not orthogonal/unitary for certain input matrix dimensions m, p, and q (e.g., m = 260, p=130, q=131). The reason was an accidental overwrite of data by xORGQR()/xUNGQR() when the WORK array was apparently large enough to use blocking.
Previously, xUNCSD2BY1 only allowed workspace queries by passing LWORK=-1 (note the missing "R"). The new commit makes xUNCSD2BY1 behave, e.g., like xHEEVD.
1318ba1 to
0c32d4b
Compare
Contributor
|
Thanks @christoph-conrads, that's great. Julien |
christoph-conrads
pushed a commit
to christoph-conrads/lapack
that referenced
this pull request
May 23, 2021
…1-fixes-20200422 2-by-1 CS decomposition fixes
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
3 participants
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
This pull request improves
xUNCSD2BY1andxORCSD2BY1.LRWORKargument toxUNCSD2BY()is checked.LRWORKvalue is passed toxBBCSD()insidexUNCSD2BY()instead of passing the minimum allowedLRWORKvalue.U2may be obviously not orthogonal or unitary, respectively, for certain input matrix dimension becausexORGQR()andxUNGQR()may overwrite the scalar factors of the elementary reflectors when assemblingU1. On the author's computer, this happens, e.g., withm=260,q=131, andp=130(these values may be dependent on the blocking size).xUNCSD2BY1assume workspace queries whenLRWORK = -1. This is also the behavior exhibited by, e.g.,xHEEVD.The
U2orthogonality issue is caused byxORGQR()/xUNGQR()using more than the minimum workspace available. The issue can be either fixed byP) orThe patch choses the latter approach.