[Haskell] Linear Equation Solver Using Arrays

Xiao-Yong Jin xj2106 at columbia.edu
Sun Sep 16 01:22:03 EDT 2007

Dear Haskellers,

My thanks to people on the irc channel, especially `int-e'.
With their help, I managed to write a linear equation solver
using STUArrays.

-------------- next part --------------
A non-text attachment was scrubbed...
Name: MatrixOp.hs.gz
Type: application/octet-stream
Size: 1379 bytes
Desc: MatrixOp, only contains a linear solver, yet
Url : http://www.haskell.org/pipermail/haskell/attachments/20070916/fc035081/MatrixOp.hs.obj
-------------- next part --------------

It is basically a direct rewrite of a C version, which
adopts Crout's LU decomposition method with implicit
pivoting, in chapter 2.3 of the book Numerical Recipes.

I must admit that the code is in a really bad style,
somewhat ugly, because of the limit of my Haskell knowledge.
The code is attached in this email for your amusement.  I
would be really thankful if you can give me any kind of
comments, idea, improvement, criticism, because I sincerely
want to make the code better.

Best regards,
    c/*    __o/*
    <\     * (__
    */\      <

More information about the Haskell mailing list