-
Notifications
You must be signed in to change notification settings - Fork 111
Lapack geqrf #2858
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
base: develop
Are you sure you want to change the base?
Lapack geqrf #2858
Conversation
The wrappers will allow use of LAPACK, cusolver, rocsolver and magma Signed-off-by: Luc Berger-Vergiat <[email protected]>
Signed-off-by: Luc Berger-Vergiat <[email protected]>
Signed-off-by: Luc Berger-Vergiat <[email protected]>
jgfouca
left a comment
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. Just a few minor things.
| static_assert(std::is_same_v<typename InfoArray::non_const_value_type, int>, | ||
| "KokkosLapack::geqrf: Info must be an array of integers."); | ||
|
|
||
| int64_t m = A.extent(0); |
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.
These could all be consts.
| for (int index(k + 1); index < m; ++index) { | ||
| h_v[index] = h_A(index, k); | ||
| } | ||
| #if 0 // def HAVE_KOKKOSKERNELS_DEBUG |
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.
Is it worthwhile to keep these #if 0 blocks around? Should we switch these to something like #if QEQRF_DEBUG so that the programmer can turn these all back on easily?
| /// Geqrf | ||
| /// | ||
|
|
||
| void F77_BLAS_MANGLE(sgeqrf, SGEQRF)(int*, int*, float*, int*, float*, float*, int*, int*); |
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.
I noticed none of these are consts. Is every argument an output or inout?
Signed-off-by: Luc Berger-Vergiat <[email protected]>
Signed-off-by: Luc Berger-Vergiat <[email protected]>
Signed-off-by: Luc Berger-Vergiat <[email protected]>
Signed-off-by: Luc Berger-Vergiat <[email protected]>
Taking over PR #2205 rebasing on the current state of develop and re-running tests, will also have a look at the unit-tests and might add more checks.