Skip to content
This repository was archived by the owner on Mar 12, 2021. It is now read-only.

Fix sparse mul!#692

Merged
bors[bot] merged 3 commits into
JuliaGPU:masterfrom
amontoison:adjoint-sparse-mul-fix
Apr 27, 2020
Merged

Fix sparse mul!#692
bors[bot] merged 3 commits into
JuliaGPU:masterfrom
amontoison:adjoint-sparse-mul-fix

Conversation

@amontoison

@amontoison amontoison commented Apr 22, 2020

Copy link
Copy Markdown
Member

@maleadt

maleadt commented Apr 23, 2020

Copy link
Copy Markdown
Member

Thanks!
bors r+

bors Bot added a commit that referenced this pull request Apr 23, 2020
@bors @amontoison
Merge #692
991ecf8 
692: Fix sparse mul! r=maleadt a=amontoison

close #629 
#630 , #637 
@haampie 

Co-authored-by: Alexis Montoison <alexis.montoison@polymtl.ca>
@bors

bors Bot commented Apr 23, 2020

Copy link
Copy Markdown
Contributor

Build failed:

@amontoison

amontoison commented Apr 23, 2020

Copy link
Copy Markdown
Member Author

Tests failed with A' * x when A is a CuSparseMatrixCSC{ComplexF32} or a CuSparseMatrixCSC{ComplexF64}. All tests passed with CuSparseMatrixCSR{...}.

@amontoison amontoison force-pushed the adjoint-sparse-mul-fix branch 2 times, most recently from faca527 to dd44454 Compare April 23, 2020 23:38
@amontoison

amontoison commented Apr 23, 2020

Copy link
Copy Markdown
Member Author

The culprit should be mv! : https://github.com/JuliaGPU/CuArrays.jl/blob/master/src/sparse/wrappers.jl#L554-L584

@amontoison

amontoison commented Apr 23, 2020 •
edited
Loading

Copy link
Copy Markdown
Member Author

The conjugate-transpose of a transposed complex matrix A is not A but conj(A), it's why tests failed with complex matrices.

@amontoison amontoison force-pushed the adjoint-sparse-mul-fix branch 2 times, most recently from 97c71dd to 713422b Compare April 26, 2020 20:06
@amontoison

amontoison commented Apr 26, 2020

Copy link
Copy Markdown
Member Author

I read the documentation of csrmv function and it seems that they forgot this corner case.
https://docs.nvidia.com/cuda/cusparse/index.html#csrmv

@maleadt, it's ready to be tested.

amontoison
amontoison commented Apr 26, 2020
View reviewed changes
Comment thread src/sparse/wrappers.jl
ctransa = 'T'
end
cutransa = cusparseop(ctransa)
cuind = cusparseindex(index)

@amontoison amontoison Apr 26, 2020

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

cuind variable is not used. I removed it.

@amontoison amontoison force-pushed the adjoint-sparse-mul-fix branch from 713422b to 2ec414f Compare April 26, 2020 20:13
@maleadt

maleadt commented Apr 27, 2020

Copy link
Copy Markdown
Member

bors d+

@bors

bors Bot commented Apr 27, 2020

Copy link
Copy Markdown
Contributor

✌️ amontoison can now approve this pull request. To approve and merge a pull request, simply reply with bors r+. More detailed instructions are available here.

@maleadt

maleadt commented Apr 27, 2020

Copy link
Copy Markdown
Member

bors try

bors Bot added a commit that referenced this pull request Apr 27, 2020
@bors
Try #692:
602e5fc
@maleadt

maleadt commented Apr 27, 2020

Copy link
Copy Markdown
Member

What about:

[[DEPRECATED]] use cusparseSpMV() instead. The routine will be removed in the next major release

@bors

bors Bot commented Apr 27, 2020

Copy link
Copy Markdown
Contributor

try

Build failed:

@maleadt

maleadt commented Apr 27, 2020

Copy link
Copy Markdown
Member

bors try

bors Bot added a commit that referenced this pull request Apr 27, 2020
@bors
Try #692:
bd4aa17
@bors

bors Bot commented Apr 27, 2020

Copy link
Copy Markdown
Contributor

try

Build succeeded:

@amontoison

amontoison commented Apr 27, 2020

Copy link
Copy Markdown
Member Author

About the DEPRECATED csrmv, it concerned numerous functions in CUSPARSE module. 😟

@amontoison

amontoison commented Apr 27, 2020

Copy link
Copy Markdown
Member Author

bors r+

@bors

bors Bot commented Apr 27, 2020

Copy link
Copy Markdown
Contributor

Build succeeded:

@bors bors Bot merged commit 16c0080 into JuliaGPU:master Apr 27, 2020
Sign up for free to subscribe to this conversation on GitHub. Already have an account? Sign in.

Reviewers

No reviews

Assignees

No one assigned

Labels

None yet

Projects

None yet

Milestone

No milestone

Development

Successfully merging this pull request may close these issues.

mul! failing for adjoint sparse matrices

2 participants

@amontoison @maleadt