-
Notifications
You must be signed in to change notification settings - Fork 452
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.
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
Fix error on LAPACKE_*tpmqrt_work for ROW MAJOR matrices #540
Fix error on LAPACKE_*tpmqrt_work for ROW MAJOR matrices #540
Conversation
Codecov ReportBase: 82.37% // Head: 82.37% // No change to project coverage 👍
Additional details and impacted files@@ Coverage Diff @@
## master #540 +/- ##
=======================================
Coverage 82.37% 82.37%
=======================================
Files 1894 1894
Lines 190681 190681
=======================================
Hits 157067 157067
Misses 33614 33614 Help us with your feedback. Take ten seconds to tell us how you rate us. Have a feature suggestion? Share it here. ☔ View full report at Codecov. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The fixes of the leading dimensions look good. I like your solution to abstract from side
by setting nrowsA
, ncolsA
and nrowsV
.
Thanks a lot for reviewing this PR, @angsch! I applied all your fixes in this last commit. I will now apply the changes to the other |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good to me 👍
The fix likely applies toLAPACKE_?tprfb_work
, too. The dimensions of A
depend on side
. When side = 'Left'
, A
has minimum dimensions (K, N)
. When side = 'Right'
, A
is (M, K)
. The row major interface seems to only support side = 'Left'
:
lapack/LAPACKE/src/lapacke_ctprfb_work.c
Lines 52 to 53 in 28f7e83
} else if( matrix_layout == LAPACK_ROW_MAJOR ) { | |
lapack_int lda_t = MAX(1,k); |
Perhaps it make sense to move this to a separate issue.
Fixes #404
Description
There are many bugs in
LAPACKE_*tpmqrt_work
routines for ROW MAJOR matrices, as pointed out by @haldaas in #404.Checklist
LAPACKE_dtpmqrt_work
LAPACKE_*tpmqrt_work
whenever my solution is validated for the double algorithm