search

achirkin / easytensor

Many-dimensional type-safe numeric ops
https://hackage.haskell.org/package/easytensor
BSD 3-Clause "New" or "Revised" License
46 stars 2 forks source link

Stability of LU factorization #15

Closed achirkin closed 5 years ago

achirkin commented 5 years ago

Sometimes LU factorization is very unstable, especially for Float matrices. We should fix this some day.

https://github.com/achirkin/easytensor/blob/4cc9a381ad09c7a013a121cc02231f0ac2de649d/easytensor/test/Numeric/MatrixFloatTest.hs#L80-L90

achirkin commented 5 years ago

Solved by 3b39e4a644c3c40d03ea1fad04eeeb22400f5d48 There must have been a bug before, but the new version seems to be far more robust.