-
Notifications
You must be signed in to change notification settings - Fork 25
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
Random matrix sampling on the GPU #74
Conversation
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.
Minor issues should be fixed, the larger curand()
interface issue might be fixed in future.
curandommatrices/src/circular.jl
Outdated
function _qr_fix!(z::CuMatrix) | ||
q, r = CuArrays.qr!(z) | ||
ph = diag(r) | ||
len = min(length(ph), 1024) |
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.
Why the magic number 1024
? Explain it or at least point to some information.
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.
This is just an ugly hack as warpsize()
results in a segfault.
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.
Then document it.
I think we should merge this and state clearly, in big bold red letters in the readme that this is highly experimental and we would appreciate any feedback.
As there are probably some mistakes I think people should have access to this iteration of GPU support and state their comments