morning all, I've just released a new module PDL::SVDSLEPc <https://metacpan.org/release/PDL-SVDSLEPc> to CPAN: it's basically just a thin wrapper around the (reduced) SVD solver routines from the SLEPc c library ("Scalable Library for Eigenvalue Problem Computations", http://slepc.upv.es/). Preliminary tests with the default configuration run about 4-5x faster than PDL::SVDLIBC <https://metacpan.org/release/PDL-SVDLIBC> on my machine here (I believe it's using ARPACK under the hood, which was mentioned on this list a few years back, but my searches have failed to turn up any PDL-ARPACK interface code).
It would be nice to have a more general interface to all the SLEPc functions, but doing that the "right" way would seem to involve a low-level PDL interface to the underlying PETSc matrix representation classes ( http://www.mcs.anl.gov/petsc/ ; these guys are using OpenMPI for parallel computation), and I'm running low on tuits at the moment. Just something on my wish list for future PDL incarnations :-) marmosets, Bryan -- Bryan Jurish "There is *always* one more bug." [email protected] -Lubarsky's Law of Cybernetic Entomology
------------------------------------------------------------------------------
_______________________________________________ pdl-general mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/pdl-general
