diff options
| author | Bas van Dijk <v.dijk.bas@gmail.com> | 2011-12-24 20:15:43 +0100 |
|---|---|---|
| committer | Bas van Dijk <v.dijk.bas@gmail.com> | 2011-12-24 20:15:43 +0100 |
| commit | a4f38eb196209436f72b938f6355f6e28474bef3 (patch) | |
| tree | 48bd5f980bc248c30e4b2dba8a7546b04ea7ba18 /lib/Numeric/LinearAlgebra | |
| parent | c3bda2d38c432fb53ce456cba295b097fd4d6ad1 (diff) | |
Fix build on GHC-7.4.1-rc1
Diffstat (limited to 'lib/Numeric/LinearAlgebra')
| -rw-r--r-- | lib/Numeric/LinearAlgebra/Algorithms.hs | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/lib/Numeric/LinearAlgebra/Algorithms.hs b/lib/Numeric/LinearAlgebra/Algorithms.hs index bea33ea..a6b4dae 100644 --- a/lib/Numeric/LinearAlgebra/Algorithms.hs +++ b/lib/Numeric/LinearAlgebra/Algorithms.hs | |||
| @@ -567,6 +567,11 @@ epslist = [ (fromIntegral k, golubeps k k) | k <- [1..]] | |||
| 567 | 567 | ||
| 568 | geps delta = head [ k | (k,g) <- epslist, g<delta] | 568 | geps delta = head [ k | (k,g) <- epslist, g<delta] |
| 569 | 569 | ||
| 570 | expGolub :: ( Fractional t, Element t, Field t | ||
| 571 | , Normed Matrix t | ||
| 572 | , RealFrac (RealOf t) | ||
| 573 | , Floating (RealOf t) | ||
| 574 | ) => Matrix t -> Matrix t | ||
| 570 | expGolub m = iterate msq f !! j | 575 | expGolub m = iterate msq f !! j |
| 571 | where j = max 0 $ floor $ logBase 2 $ pnorm Infinity m | 576 | where j = max 0 $ floor $ logBase 2 $ pnorm Infinity m |
| 572 | a = m */ fromIntegral ((2::Int)^j) | 577 | a = m */ fromIntegral ((2::Int)^j) |