Skip to content
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

Documentation for the package #66

Open
ViralBShah opened this issue May 28, 2020 · 5 comments
Open

Documentation for the package #66

ViralBShah opened this issue May 28, 2020 · 5 comments

Comments

@ViralBShah
Copy link
Member

Filing a help needed issue to put together Documenter based documentation for this package.

@ViralBShah
Copy link
Member Author

@andreasnoack Could you list in the README, the current set of capabilities available in the package, and a sort of quick roadmap of the next set of things to add? This will help people pick up the project.

@alanedelman ^

@alanedelman
Copy link

alanedelman commented Jun 5, 2020

i'd love for performance guidance right on top of when an svd from elemental might be beneficial --- example a statement that suggests that at least on one machine with 16 processors, we found a performance increase over lapack with mkl and open blas.....

and if there is a way others can add performance anecdotes in a useful way...

@ViralBShah
Copy link
Member Author

ViralBShah commented Jun 5, 2020

Julia with MKL for SVD:

julia> a = randn(10^4,10^4);
julia> using LinearAlgebra
julia> @time u,s,v = svd(a)
136.560808 seconds (258.16 k allocations: 4.485 GiB, 0.16% gc time)

With Elemental using 20 processors on a single box using the example in the README:

julia> @mpi_do man Elemental.gaussian!(A, 10000, 10000);
julia> @time @mpi_do man U, s, V = svd(A);
 80.731189 seconds (2.68 k allocations: 141.296 KiB)

@andreasnoack
Copy link
Member

I'm a little surprised by this. Would it be easy to try MKL's LAPACK.gesvd!. Our svd calls LAPACK.gesdd!.

@ViralBShah
Copy link
Member Author

ViralBShah commented Jun 5, 2020

LAPACK+OpenBLAS (the julia default) is 170s.

Should be easy to try gesvd - just have to go through LAPACK docs and the wrappers.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants