-
Notifications
You must be signed in to change notification settings - Fork 206
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
[CUSOLVER] Add ordering functions #1198
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.
Thanks! Couple of comments (that apply to all of these functions).
lib/cusolver/ordering.jl
Outdated
rowsA = A.rowPtr |> Vector{Cint} | ||
colsA = A.colVal |> Vector{Cint} | ||
p = zeros(Cint, n) | ||
cusolverSpXcsrsymrcmHost(sparse_handle(), n, nnz(A), descA, rowsA, colsA, p) |
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 do you use the host APIs here?
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.
symrcm
,symmdq
,symamd
andmetisnd
are used insidecusolverRF
(Add a wrapper to cusolverRF in CUSOLVER #828) but no low-level wrappers are implemented.zfd
should be combined in practice withilu0
routine (incomplete LU factorization) because the factorization is done without pivoting and if one diagonal element is zero the factorization failed.
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.
Ah, there aren't any device equivalents for these APIs...
lib/cusolver/ordering.jl
Outdated
n, m = A.dims | ||
(m ≠ m) && throw(DimensionMismatch("A must be square, but has dimensions ($n,$m)!")) | ||
descA = CUSPARSE.CuMatrixDescriptor('G', 'L', 'N', index) | ||
rowsA = A.rowPtr |> Vector{Cint} |
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.
We don't want to do this.
@maleadt I take into account all your comments. |
Codecov Report
@@ Coverage Diff @@
## master #1198 +/- ##
===========================================
+ Coverage 0.49% 80.45% +79.96%
===========================================
Files 118 119 +1
Lines 7953 8374 +421
===========================================
+ Hits 39 6737 +6698
+ Misses 7914 1637 -6277
Continue to review full report at Codecov.
|
Thanks, looks good! |
No description provided.