I look forward to see what comes of this!<div><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Fri, Nov 30, 2012 at 4:58 PM, Mark Flamer <span dir="ltr">&lt;<a href="mailto:mark@flamerassoc.com" target="_blank">mark@flamerassoc.com</a>&gt;</span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Thanks for all the replies,<br>
 It sounds like there is enough interest and even some potential<br>
collaborators out there. I have created a few data structures to represent<br>
sparse vectors and matrices. The vector was a simple binary tree and the<br>
matrix a quad tree. As I suspected a standard IntMap was around 3X as fast<br>
as my binary tree, so I have switched to the IntMap for now. I was hoping to<br>
hold on to my quad tree for the matrix rep. because I like the elegance of<br>
the recursive algorithms like Strassen’s for multiplication. In the end it<br>
will most likely be best to have a few different matrix representations<br>
anyway. For instance, I know that compressed row form is most efficient for<br>
certain algorithms. I have a Gauss iterative solver working currently and am<br>
thinking of moving on to a parallel Conjugate gradient or direct solver<br>
using LU decomposition. I’m no expert in Haskell or numeric methods but I do<br>
my best. I’ll clean up what I have and open up the project on Github or<br>
Bitbucket. Any other thoughts or ideas are welcome.<br>
I’m currently using the Matrix Market files for testing. It would be nice to<br>
benchmark this against an industry standard solver in C or Fortran, without<br>
having to do a lot of work to set it up. Does anyone know of an easy way to<br>
do this? I’m just trying to get results to compare in orders of magnitude<br>
for now. It would be motivating to see these comparisons.<br>
<br>
<br>
<br>
<br>
--<br>
View this message in context: <a href="http://haskell.1045720.n5.nabble.com/Is-anyone-working-on-a-sparse-matrix-library-in-Haskell-tp5721452p5721556.html" target="_blank">http://haskell.1045720.n5.nabble.com/Is-anyone-working-on-a-sparse-matrix-library-in-Haskell-tp5721452p5721556.html</a><br>


<div class="HOEnZb"><div class="h5">Sent from the Haskell - Haskell-Cafe mailing list archive at Nabble.com.<br>
<br>
_______________________________________________<br>
Haskell-Cafe mailing list<br>
<a href="mailto:Haskell-Cafe@haskell.org">Haskell-Cafe@haskell.org</a><br>
<a href="http://www.haskell.org/mailman/listinfo/haskell-cafe" target="_blank">http://www.haskell.org/mailman/listinfo/haskell-cafe</a><br>
</div></div></blockquote></div><br></div>