I look forward to see what comes of this!
On Fri, Nov 30, 2012 at 4:58 PM, Mark Flamer <m...@flamerassoc.com> wrote: > Thanks for all the replies, > It sounds like there is enough interest and even some potential > collaborators out there. I have created a few data structures to represent > sparse vectors and matrices. The vector was a simple binary tree and the > matrix a quad tree. As I suspected a standard IntMap was around 3X as fast > as my binary tree, so I have switched to the IntMap for now. I was hoping > to > hold on to my quad tree for the matrix rep. because I like the elegance of > the recursive algorithms like Strassen’s for multiplication. In the end it > will most likely be best to have a few different matrix representations > anyway. For instance, I know that compressed row form is most efficient for > certain algorithms. I have a Gauss iterative solver working currently and > am > thinking of moving on to a parallel Conjugate gradient or direct solver > using LU decomposition. I’m no expert in Haskell or numeric methods but I > do > my best. I’ll clean up what I have and open up the project on Github or > Bitbucket. Any other thoughts or ideas are welcome. > I’m currently using the Matrix Market files for testing. It would be nice > to > benchmark this against an industry standard solver in C or Fortran, without > having to do a lot of work to set it up. Does anyone know of an easy way to > do this? I’m just trying to get results to compare in orders of magnitude > for now. It would be motivating to see these comparisons. > > > > > -- > View this message in context: > http://haskell.1045720.n5.nabble.com/Is-anyone-working-on-a-sparse-matrix-library-in-Haskell-tp5721452p5721556.html > Sent from the Haskell - Haskell-Cafe mailing list archive at Nabble.com. > > _______________________________________________ > Haskell-Cafe mailing list > Haskell-Cafe@haskell.org > http://www.haskell.org/mailman/listinfo/haskell-cafe >
_______________________________________________ Haskell-Cafe mailing list Haskell-Cafe@haskell.org http://www.haskell.org/mailman/listinfo/haskell-cafe