
4 Feb
2009
4 Feb
'09
2:09 p.m.
2009/2/3 Rafael Gustavo da Cunha Pereira Pinto
After a discussion on whether is possible to compile hmatrix in Windows, I decided to go crazy and do a LU decomposition entirely in Haskell... import Data.Array.IArray
...
e_an i j=a!(i,j)-(lik i)*a!(k,j)
There are three different representations of arrays in this code: arrays, lists and a functional one where f represents an array with elements (f i j). Looks like a candidate for a stream fusion type thing so the user only writes code using one representation and some RULES switch back and forth between representations behind the scenes. -- Dan