- Haskell-Cafe - Haskell.org

Re[2]: [Haskell] GHC inferred type different for pointed/point-free defs?
by Bulat Ziganshin 04 Nov '06

04 Nov '06

Hello Dan, Saturday, November 4, 2006, 3:49:13 AM, you wrote: > You might want to rename your (&&&) and (|||), though, since they clash > with those defined in Control.Arrow. > I found that by defining (???) then way I did, I could use > (Control.Arrow.|||) as-is for the second part. i just not use arrows because i don't understand them. my definitions of "&&&" and "|||" is like those in C/perl/ruby and most times i use them alone: debugging &&& debug_info, first_time ||| putStrLn "" in the far future, i will be glad to see && and || redefined in this way (or some like this; it should just provide ruby's look&feel plus extensibility). although for my own program this may be solved by just not importing Prelude. hmmm, i should think about it.. but && and || functions is not exactly the same as ?: operator (nor in C, nor in standard Haskell, nor in my implementation). with these operations condition is tested _twice_ and if second argument has "false" value, third one will be returned anyway: True && False || True => True True ? False : True => False btw, it may be extended slightly more by adding application version: a &&. f = a &&& (f a) example: str &&. ("debug info: "++) now such function named unlessNull that is not so neat -- Best regards, Bulat mailto:Bulat.Ziganshin@gmail.com

1 0

Re: [Haskell] GHC inferred type different for pointed/point-free defs?
by Neil Mitchell 03 Nov '06

03 Nov '06

Hi, You probably want to read up on: Monomorphism restriction (I don't think this applies here, but I'm never too sure!): http://www.haskell.org/hawiki/MonomorphismRestriction Defaulting: http://www.haskell.org/onlinereport/decls.html#sect4.3.4 I don't understand either to any great degree, but when points free starts changing things, its usually one of those two kicking in. And as a side note, haskell@ tends to be used more for announcements, and haskell-cafe@ is used more for asking questions etc. Thanks Neil On 11/3/06, Dan Weston <westondan(a)imageworks.com> wrote: > Help! One of two things is going on: > > 1) I don't understand what I'm doing > 2) GHC is inferring different types for pointed and > point-free function definition. > > I wanted to define Haskell equivalents to the C ternary operator. > Basically, predicate ??? doIfTrue ||| doIfFalse > > For some reason, though, in GHC 6.4 and 6.7, the types > inferred for abs1 and ab2 (defined below) are different: > > *Main> :t abs1 > abs1 :: (Num a, Ord a) => a -> a > > *Main> :t abs2 > abs2 :: Integer -> Integer > > How is abs2 less general just because I use point-free > notation instead of pointed notation? They should have the same > type, no? I'm stumped... > > > import Control.Arrow > > import Data.Either > > > > infix 1 ?, ??, ???, ???? > > > > (?) True = Left > > (?) False = Right > > > > (??) True = Left . fst > > (??) False = Right . snd > > > > p ??? q = (p &&& arr id) >>> uncurry (?) >>> q > > p ???? q = (p &&& arr id) >>> uncurry (??) >>> q > > > > abs1 x = (< 0) ??? negate ||| id $ x > > abs2 = (< 0) ??? negate ||| id > > > > checks :: [Bool] > > checks = > > [ > > (True ? 3 ) == (Left 3), > > (False ? 4 ) == (Right 4), > > (True ?? (3,4)) == (Left 3), > > (False ?? (3,4)) == (Right 4), > > (even ??? (`div` 2) ||| (+1).(*3) $ 3 ) == ( 10), > > (even ??? (`div` 2) ||| (+1).(*3) $ 4 ) == ( 2), > > (uncurry (==) ???? (+1) ||| (*3) $ (3,3)) == ( 4), > > (uncurry (==) ???? (+1) ||| (*3) $ (3,4)) == ( 12), > > (uncurry (==) ???? (+1) +++ (*3) $ (3,3)) == (Left 4), > > (uncurry (==) ???? (+1) +++ (*3) $ (3,4)) == (Right 12), > > (abs1 5 ) == ( 5), > > (abs1 (-5) ) == ( 5), > > (abs2 5 ) == ( 5), > > (abs2 (-5) ) == ( 5) > > ] > > > > main = print (and checks) > > _______________________________________________ > Haskell mailing list > Haskell(a)haskell.org > http://www.haskell.org/mailman/listinfo/haskell >

1 0

Finding Memory Leaks
by Matthew Bromberg 03 Nov '06

03 Nov '06

I have a rather complex Haskell program that performs an engineering simulation using a lot of large complex valued matrices. The imperative components of the simulation are written in C and interfaced into Haskell which handles the higher level logic. This has served me fairly well, and I've been pleased with both the development style and the performance of my code. However I've run into a fairly major issue with an unresolvable memory leak in the Haskell code. The memory leak is apparent since both private allocated bytes and page file bytes associated with the process increase linearly over time while the process is running. Eventually it causes the process to crash. I believe it's in the Haskell code because I have replaced malloc, free, realloc etc. in the C code with a memory leak detector and the C code appears to free all allocated resources. I also know that the memory leak is due to unmatched free's from a call to malloc from one of the windows XP system libraries. Unfortunately GHC's heap profiler does not reveal the memory leak. All my code is compiled using the -prof -auto flag and then run using +RTS -hc -RTS. The resulting plots do not show a linear increase in heap usage, although the Windows XP operating system does report such an increase. This leaves me in a difficult position in terms of finding the offending code. I am currently resorting to commenting out code line by line, but it's very tedious since I have to preserve the state functionality of the code to get it to run. Is there a way to compile Haskell into C using GHC and then forcing the C code to hook up to the memory leak detecting versions of malloc? I know that GHC can produce intermediary C code, but I have no idea how to access it and how to setup all the run time libraries etc, so that I can compile the C code separately. Any advice for me in this situation would be appreciated.

4 5

[Haskell-cafe]Prime Generator time limit exceeded
by alaiyeshi 03 Nov '06

03 Nov '06

Hi I'm new to Haskell. I found this site on the Haskell wiki https://www.spoj.pl. But I got some trouble on trying to solve the problem titled "Prime Generator" https://www.spoj.pl/problems/PRIME1. The online-judge system tells me "time limit excedded" Would you be so kind to tell me how to make it more faster? And any other suggestion is welcome. Thanks in advance. --------------------------------------Code begin------------------------------------------------------------ module Main where import IO import List main = do input_size<-getLine content<-get_contents (read input_size) mapM_ (\r-> do mapM_ (print) (primeGenerator (parse r)); putStrLn "") content get_contents n | n == 0 = return [] | otherwise = do content<-getLine rests<-get_contents (n-1) return ([content]++rests) primeGenerator [start,end] = [x | x<-[start..end], all (== 1) (map (gcd x) [2..(x-1)]), x/=1] parse s = unfoldr (\x-> case x of [] -> Nothing _ -> Just (head (reads x))) s -------------------------------Code ends-------------------------------------------------------------------------------- (BTW: I'm new to this mailling list also, forgive my rudeness if I am, and forgive my poor English)

5 11

Cost of seq (sorry about the old subject)
by Chad Scherrer 03 Nov '06

03 Nov '06

On 11/2/06, Bulat Ziganshin <bulat.ziganshin(a)gmail.com> wrote: > Hello Chad, > > Thursday, November 2, 2006, 8:49:17 PM, you [Chad] wrote: > > Would it be possible for the type system to distinguish at compile > > time whether something would need to be evaluated, and optimize away > > redundant `seq`s? Maybe this is what the strictness analyzer does > > already. > > afaik, seq cost is zero (except for evaluating expressions, of course) > So if x has already been evaluated, does x `seq` y evaluate just as quickly as y alone, or does it require extra cycles to make sure x has been evaluated? -- Chad Scherrer "Time flies like an arrow; fruit flies like a banana" -- Groucho Marx

2 1

Re: [Haskell-cafe] Class
by Slavomir Kaslev 03 Nov '06

03 Nov '06

On 11/2/06, Sebastian Sylvan <sylvan(a)student.chalmers.se> wrote: > On 11/2/06, Slavomir Kaslev <slavomir.kaslev(a)gmail.com> wrote: > > Little by little, I think I am getting closer. > > > > class Show a => Visible a where > > toString :: a -> String > > toString = show > > size :: a -> Int > > size = length . show > > > > Is just declaration, not definition. It doesn't define anything, even > > though it has default implementations for toString and size one still > > needs to define instance of it. Right? > > > > As Sebastian Sylvan proposed, I probably need something like this: > > > > class Visible a where > > toString :: a -> String > > size :: a -> Int > > > > instance Show a => Visible a where > > toString = show > > size = length . show > > > > But it seems ghc doesn't like instance definitions like 'instance Show > > a => Visible a where ...'. Why? > > This is not Haskell98 since instance declarations of this form could > cause the type checker to go into an infinite loop. This particular > one is okay, though, but you have to start ghc with > -fallow-undecidable-instances and -fglasgow-exts I'm afraid. > > > -- > Sebastian Sylvan > +46(0)736-818655 > UIN: 44640862 > Thanks, Sebastian. That was helpful. Are there any papers on the subject? -- Slavomir Kaslev

3 3

How to improve speed? (MersenneTwister is several times slower than C version)
by isto 03 Nov '06

03 Nov '06

Hi all, On HaWiki was an announcement of MersenneTwister made by Lennart Augustsson. On a typical run to find out 10000000th rnd num the output is (code shown below): $ time ./testMTla Testing Mersenne Twister. Result is [3063349438] real 0m4.925s user 0m4.856s I was exercising with the very same algorithm and tried to make it efficient (by using IOUArray): now a typical run looks like (code shown below): $ time ./testMT Testing Mersenne Twister. 3063349438 real 0m3.032s user 0m3.004s The original C-version (modified so that only the last number is shown) gives typically $ time ./mt19937ar outputs of genrand_int32() 3063349438 real 0m0.624s user 0m0.616s Results are similar with 64 bit IOUArray against 64 bit C variant. C seems to work about 5 to 10 times faster in this case. I have tried to do different things but now I'm stuck. unsafeRead and unsafeWrite improved a bit the lazy (STUArray-version) and IOUArray-versions but not very much. I took a look of Core file but then, I'm not sure where the boxed values are ok. E.g. should IOUArray Int Word64 be replaced with something else? Any hints and comments on how to improve the efficiency and make everything better will be appreciated a lot! br, Isto ----------------------------- testMTla.hs (MersenneTwister, see HaWiki) module Main where -- ghc -O3 -optc-O3 -optc-ffast-math -fexcess-precision --make testMTla import MersenneTwister main = do putStrLn "Testing Mersenne Twister." let mt = mersenneTwister 100 w = take 1 (drop 9999999 mt) -- w = take 1 (drop 99 mt) putStrLn $ "Result is " ++ (show w) ----------------------------- ----------------------------- testMT.hs module Main where -- Compile eg with -- ghc -O3 -optc-O3 -optc-ffast-math -fexcess-precision --make testMT import Mersenne genRNums32 :: MT32 -> Int -> IO (MT32) genRNums32 mt nCnt = gRN mt nCnt where gRN :: MT32 -> Int -> IO (MT32) gRN mt nCnt | mt `seq` nCnt `seq` False = undefined gRN mt 1 = do (r,mt') <- next32 mt putStrLn $ (show r) return mt' gRN mt nCnt = do (r,mt') <- next32 mt gRN mt' $! (nCnt-1) main = do putStrLn "Testing Mersenne Twister." mt32 <- initialiseGenerator32 100 genRNums32 mt32 10000000 ----------------------------- ----------------------------- Mersenne.hs (sorry for linewraps) module Mersenne where import Data.Bits import Data.Word import Data.Array.Base import Data.Array.MArray import Data.Array.IO -- import System.Random data MT32 = MT32 (IOUArray Int Word32) Int data MT64 = MT64 (IOUArray Int Word64) Int last32bitsof :: Word32 -> Word32 last32bitsof a = a .&. 0xffffffff -- == (2^32-1) lm32 = 0x7fffffff :: Word32 um32 = 0x80000000 :: Word32 mA32 = 0x9908b0df :: Word32 -- == 2567483615 -- Array of length 624. initialiseGenerator32 :: Int -> IO MT32 initialiseGenerator32 seed = do let s = last32bitsof (fromIntegral seed)::Word32 mt <- newArray (0,623) (0::Word32) unsafeWrite mt 0 s iG mt s 1 mt' <- generateNumbers32 mt return (MT32 mt' 0) where iG :: (IOUArray Int Word32) -> Word32 -> Int -> IO (IOUArray Int Word32) iG mt lastNro n | n == 624 = return mt | otherwise = do let n1 = lastNro `xor` (shiftR lastNro 30) new = (1812433253 * n1 + (fromIntegral n)::Word32) unsafeWrite mt n new iG mt new (n+1) generateNumbers32 :: (IOUArray Int Word32) -> IO (IOUArray Int Word32) generateNumbers32 mt = gLoop 0 mt where gLoop :: Int -> (IOUArray Int Word32) -> IO (IOUArray Int Word32) gLoop i mt | i==623 = do wL <- unsafeRead mt 623 w0 <- unsafeRead mt 0 w396 <- unsafeRead mt 396 let y = (wL .&. um32) .|. (w0 .&. lm32) :: Word32 if even y then unsafeWrite mt 623 (w396 `xor` (shiftR y 1)) else unsafeWrite mt 623 (w396 `xor` (shiftR y 1) `xor` mA32) return mt | otherwise = do wi <- unsafeRead mt i wi1 <- unsafeRead mt (i+1) w3 <- unsafeRead mt ((i+397) `mod` 624) let y = (wi .&. um32) .|. (wi1 .&. lm32) if even y then unsafeWrite mt i (w3 `xor` (shiftR y 1)) else unsafeWrite mt i (w3 `xor` (shiftR y 1) `xor` mA32) gLoop (i+1) mt next32 :: MT32 -> IO (Word32, MT32) next32 (MT32 mt i) | i >= 624 = do mt' <- generateNumbers32 mt let m = MT32 mt' (i `mod` 624) (w,m') <- next32 m return (w,m') | otherwise = do y <- unsafeRead mt i let y1 = y `xor` (shiftR y 11) y2 = y1 `xor` ((shiftL y1 7 ) .&. 0x9d2c5680) -- == 2636928640 y3 = y2 `xor` ((shiftL y2 15) .&. 0xefc60000) -- == 4022730752 y4 = y3 `xor` (shiftR y3 18) return $ (y4, MT32 mt (i+1)) mA64 = 0xB5026F5AA96619E9 :: Word64 um64 = 0xFFFFFFFF80000000 :: Word64 lm64 = 0x7FFFFFFF :: Word64 initialiseGenerator64 :: Int -> IO (MT64) initialiseGenerator64 seed = do let s = (fromIntegral seed)::Word64 mt <- newArray (0,311) (0::Word64) unsafeWrite mt 0 s iG mt s 1 generateNumbers64 mt return (MT64 mt 0) where iG :: (IOUArray Int Word64) -> Word64 -> Int -> IO (IOUArray Int Word64) iG mt lN i | mt `seq` lN `seq` i `seq` False = undefined iG mt lastNro 312 = return mt iG mt lastNro n = do let n1 = lastNro `xor` (shiftR lastNro 62) new = (6364136223846793005 * n1 + (fromIntegral n)::Word64) unsafeWrite mt n new iG mt new $! (n+1) generateNumbers64 :: (IOUArray Int Word64) -> IO () generateNumbers64 mt = gLoop 0 where gLoop :: Int -> IO () gLoop i | i `seq` False = undefined gLoop 311 = do wL <- unsafeRead mt 311 w0 <- unsafeRead mt 0 w155 <- unsafeRead mt 155 let y = (wL .&. um64) .|. (w0 .&. lm64) :: Word64 if even y then unsafeWrite mt 311 (w155 `xor` (shiftR y 1)) else unsafeWrite mt 311 (w155 `xor` (shiftR y 1) `xor` mA64) return () gLoop i = do wi <- unsafeRead mt i wi1 <- unsafeRead mt (i+1) w3 <- unsafeRead mt ((i+156) `mod` 312) let y = (wi .&. um64) .|. (wi1 .&. lm64) if even y then unsafeWrite mt i (w3 `xor` (shiftR y 1)) else unsafeWrite mt i (w3 `xor` (shiftR y 1) `xor` mA64) gLoop $! (i+1) next64 :: MT64 -> IO (Word64, MT64) next64 (MT64 mt 312) = do generateNumbers64 mt let m = MT64 mt 0 (w,m') <- next64 m return (w,m') next64 (MT64 mt i) = do y <- unsafeRead mt i let y1 = y `xor` ((shiftR y 29) .&. 0x5555555555555555) y2 = y1 `xor` ((shiftL y1 17) .&. 0x71D67FFFEDA60000) y3 = y2 `xor` ((shiftL y2 37) .&. 0xFFF7EEE000000000) y4 = y3 `xor` (shiftR y3 43) return $! (y4, MT64 mt (i+1))

5 14

Non-integer part of a Double
by Maurício 03 Nov '06

03 Nov '06

Hi, Is there a function to get that? I'm using \x -> x - fromIntegral(floor x) since I was not able to find something better, but I guess I have missed something in the standard library since there are functions with similar funcionality. Thanks, Maurício

3 3

Fwd: [Haskell-cafe] Class
by Slavomir Kaslev 03 Nov '06

03 Nov '06

Little by little, I think I am getting closer. class Show a => Visible a where toString :: a -> String toString = show size :: a -> Int size = length . show Is just declaration, not definition. It doesn't define anything, even though it has default implementations for toString and size one still needs to define instance of it. Right? As Sebastian Sylvan proposed, I probably need something like this: class Visible a where toString :: a -> String size :: a -> Int instance Show a => Visible a where toString = show size = length . show But it seems ghc doesn't like instance definitions like 'instance Show a => Visible a where ...'. Why? -- Slavomir Kaslev

2 1

Re: [Haskell-cafe] Class
by Slavomir Kaslev 03 Nov '06

03 Nov '06

On 11/2/06, Sebastian Sylvan <sylvan(a)student.chalmers.se> wrote: > On 11/2/06, Slavomir Kaslev <slavomir.kaslev(a)gmail.com> wrote: > > > > > Visible 'subclasses' Show. Doesn't this mean that Visible should be > > defined for all types with Show instances? > > No, it just means that any type that you want to instantiate in > Visible needs to also have a Show instance (is this a restriction you > really want? See below). Well, yes that's what I meant. I want class Visible for objects that have Show instance and I want to be able to use this Show instance in the definition of Visible. Again this is not some practical example, I am just trying to grasp Haskell classes and instances. > > You could do something like > > instance Show a => Visible a > > (the default implementations that you provided in the class should be > enough, you could possibly move the implementations to the "Show > instance" and remove the Show super class from the Visible class - > then you could make even non-Show types instances of Visible without > having to define their instance for Show as well). > > /S > > > -- > Sebastian Sylvan > +46(0)736-818655 > UIN: 44640862 > This doesn't work either: import List class Visible a where toString :: a -> String size :: a -> Int instance Show a => Visible a where toString = show size = length . show instance Visible a => Visible [a] where toString = concat . map toString size = foldl (+) 0 . map size vSort :: (Visible a, Ord a) => [a] -> String vSort = toString . List.sort s = vSort [1..3] gives: Illegal instance declaration for `Visible a' (The instance type must be of form (T a b c) where T is not a synonym, and a,b,c are distinct type variables) In the instance declaration for `Visible a' Failed, modules loaded: none. at instance Show a => Visible a where. Probably I should reconsider my expectations? How should something like this designed? -- Slavomir Kaslev

2 3