Hi all, hi Manuel,
thanks for the help. Yes, i read the user docs for "earlier, faster,
thriftier, smaller ..." or a permutation of this,
but i do not _really_ remember from my university time what strict
functions are. Is by strict meant as in Formal OO that
a function is trict iff its result is _|_ iff any of its arguments is _|_ ?
I will try to use
data Genauigkeit = Absolut {wert, intervall :: !Double}
| Intervall !Double
is that what you meant? I will read the strictness flag documentation.
Thanks for the hints.
Regards,
Andreas Schroeder
|---------+----------------------------------------->
| | Manuel M T Chakravarty |
| |
-------------------------------------------------------------------------------------------------------------------| | | | An: Andreas.Schroeder@gillardon.de | | Kopie: glasgow-haskell-users@haskell.org | | Thema: Re: GHC with MS .Net 2003 C compiler | -------------------------------------------------------------------------------------------------------------------|
Andreas.Schroeder@gillardon.de wrote, Generally, did you look at http://haskell.org/ghc/docs/latest/html/users_guide/faster.html
module AIBD(aibd, rechneZins, rechneRate, rechneKapital) where import Foreign(unsafePerformIO)
type Funktion = Double -> Double
data Genauigkeit = Absolut {wert, intervall :: Double} | Intervall Double
See the paragraph under "Use strictness annotations" on the above web page on how to use strictness annotations to speed this data type up.
istOk :: Genauigkeit -> Double -> Double -> Bool istOk (Absolut w i) x _ = abs (x-w) <= i istOk (Intervall i) _ dx = abs dx <= i
sekantenVerfahren :: Funktion -> Genauigkeit -> Genauigkeit -> Int -> Double -> Double sekantenVerfahren f gx gy tiefe start = sekantenIter f tiefe gx gy x1 x2 y1 y2 where x1 = start x2 = start + start * 0.1 y1 = f x1 y2 = f x2 sekantenIter _ 0 _ _ _ x2 _ _ = x2 sekantenIter f tiefe gx gy x1 x2 y1 y2 = let x3 = x2 - y2 * (x2 - x1) / (y2 - y1) y3 = f x3 dy = y3 - y2 dx = x3 - x2 in if (istOk gx x3 dx) && (istOk gy y3 dy) then x3 else sekantenIter f (tiefe-1) gx gy x2 x3 y2 y3
Code like this may benefit from using the option "-O2" and possibly also "-fliberate-case-threshold100" http://haskell.org/ghc/docs/latest/html/users_guide/flag-reference.html#AEN5...
foreign import ccall "HAIBDUtil.h caibd" caibd :: Double -> Double -> Double -> Int -> IO Double
aibd:: Double -> Double -> Double -> Int -> Double aibd kapital zins rate jahre = unsafePerformIO (caibd kapital zins rate jahre) {- does the same than aibd:: Double -> Double -> Double -> Int -> Double aibd kapital _ _ 0 = kapital aibd kapital zins rate jahre = aibd (kapital + kapital * zins - rate) zins rate (jahre-1) -}
Two comments with -O2 -fliberate-case-threshold100 the Haskell version may be as fast as the C version. If not, better foreign import as a pure functions foreign import ccall "HAIBDUtil.h caibd" caibd :: Double -> Double -> Double -> Int -> Double and omit the unsafePerformIO. Cheers, Manuel _______________________________________________ Glasgow-haskell-users mailing list Glasgow-haskell-users@haskell.org http://www.haskell.org/mailman/listinfo/glasgow-haskell-users