[julia-users] Re: Eigenvalues of symmetric dense n=10^6 matrix (ScaLAPACK.jl?)

2016-03-08 Thread Jack Poulson
ue would likely be preferred. Jack On Tuesday, March 8, 2016 at 8:20:10 PM UTC-8, Jack Poulson wrote: > > There were several recent extensions to Elemental's SVD support, as > detailed here: > https://github.com/elemental/Elemental/issues/125 > > In particular, FULL_SVD,

[julia-users] Re: Eigenvalues of symmetric dense n=10^6 matrix (ScaLAPACK.jl?)

2016-03-08 Thread Jack Poulson
There were several recent extensions to Elemental's SVD support, as detailed here: https://github.com/elemental/Elemental/issues/125 In particular, FULL_SVD, THIN_SVD, COMPACT_SVD, and PRODUCT_SVD are now all supported, where the latter should used if you do not need small triplets as it uses a

[julia-users] Re: higher rank sparse matrices

2015-09-15 Thread Jack Poulson
I believe that Tony is suggesting manually applying the sparse operator rather than explicitly constructing it and then applying it. This is a common (and significant) performance optimization when a sparse operator is only used once or twice, as the construction of the sparse matrix is often *

Re: [julia-users] Can't figure out why Julia code takes 2x time of identical Python

2014-03-09 Thread Jack Poulson
On Sunday, March 9, 2014 3:08:53 PM UTC-4, Milan Bouchet-Valat wrote: > > Le dimanche 09 mars 2014 à 11:39 -0700, Jack Poulson a écrit : > > On Tuesday, December 24, 2013 7:55:23 PM UTC-5, Ivar Nesje wrote: > > There is a SubArray implementation, sub(

Re: [julia-users] Can't figure out why Julia code takes 2x time of identical Python

2014-03-09 Thread Jack Poulson
On Tuesday, December 24, 2013 7:55:23 PM UTC-5, Ivar Nesje wrote: > > There is a SubArray implementation, sub(), but it has some performance > issues related to indexing that should be fixed before it becomes default > for slicing. > > For small arrays I would guess that hardcoded multiplication