-
Notifications
You must be signed in to change notification settings - Fork 26.3k
Test torch.svd using complex float and double numbers (take 2) #45795
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
Conversation
Array allocated by ALLOCATE_ARRAY will be freed when it goes out of scope
…float and double numbers" This reverts commit 5a47a21.
facebook-github-bot
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.
@malfet has imported this pull request. If you are a Facebook employee, you can view this diff on Phabricator.
facebook-github-bot
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.
@malfet has imported this pull request. If you are a Facebook employee, you can view this diff on Phabricator.
IvanYashchuk
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.
Hi! I left a few suggestions regarding lrwork that come from the MAGMA documentation.
| int64_t lrwork; // These settings are valid for on LAPACK 3.6+ | ||
| if (jobz == MagmaNoVec){ | ||
| lrwork = 7 * mn; | ||
| } else if (mx > 10 * mn){ |
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.
mx > 10 * n was an arbitrary decision, since LAPACK seems to not document this number #27653 (comment).
But MAGMA documents this condition. The threshold for mx >> mn is currently mx >= mn*17/9.
In code, this condition is calculated with (magma_int_t) (mn * 17. / 9.), maybe we can make it mx >= mn * 2?
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.
Do you know how 17/9 ration is computed? As in the documentation is simply says if mx >> mn
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 don't know.
Adds support for magmaSvd for complex numbers
Fixes use-after-free error in
apply_symeig