regress/sys/uvm/pdsim/genpat5.hs

   1 {-
   2 /*      $NetBSD: genpat5.hs,v 1.3 2006/09/30 16:14:00 yamt Exp $        */
   3
   4 /*-
   5  * Copyright (c)2005 YAMAMOTO Takashi,
   6  * All rights reserved.
   7  *
   8  * Redistribution and use in source and binary forms, with or without
   9  * modification, are permitted provided that the following conditions
  10  * are met:
  11  * 1. Redistributions of source code must retain the above copyright
  12  *    notice, this list of conditions and the following disclaimer.
  13  * 2. Redistributions in binary form must reproduce the above copyright
  14  *    notice, this list of conditions and the following disclaimer in the
  15  *    documentation and/or other materials provided with the distribution.
  16  *
  17  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
  18  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  19  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  20  * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
  21  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  22  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  23  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  24  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  25  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  26  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  27  * SUCH DAMAGE.
  28  */
  29 -}
  30
  31 import System.Environment
  32 import List
  33 import Numeric
  34 import Random
  35
  36 byturns (x:xs) (y:ys) = x:y:byturns xs ys
  37
  38 acc0 _ r1 = cycle [start..(start+n-1)] where
  39         (start, n) = r1
  40
  41 acc1 seed r1 = map f rands where
  42         (start, n) = r1
  43         f x = start + x `mod` n
  44         rands = unfoldr (Just . next) $ mkStdGen seed
  45
  46 acc2 seed r1 = map f rands where
  47         (start, n) = r1
  48         f x = start + intsqrt (x `mod` n ^ 2)
  49         rands = unfoldr (Just . next) $ mkStdGen seed where
  50         intsqrt = truncate . (sqrt . ((flip encodeFloat) 0)) . toInteger
  51
  52 acc3 seed r1 = map f rands where
  53         (start, n) = r1
  54         f x = start + g (mod x 1024) n
  55         g x n = x * x * n `div` (1024 * 1024)
  56         rands = unfoldr (Just . next) $ mkStdGen seed where
  57
  58 acc4 seed r1 = map g $ zip [1..] $ map f rands where
  59         f x = apply r1 x where d = x `mod` 2
  60         g (x,y) = x+y
  61         apply r x = fst r + intsqrt (x `mod` (snd r) ^ 2)
  62         intsqrt = truncate . (sqrt . ((flip encodeFloat) 0)) . toInteger
  63         rands = unfoldr (Just . next) $ mkStdGen seed
  64
  65 patterns = [acc0, acc1, acc2, acc3, acc4]
  66
  67 mklist t a b = take t $ byturns a b
  68
  69 main = getArgs >>= \ args ->
  70         let
  71                 pat = read (args !! 0)
  72                 t = read (args !! 1)
  73                 start = read (args !! 2)
  74                 n = read (args !! 3)
  75                 seed = 0
  76         in
  77         mapM_ print $ take t $ (patterns !! pat) seed (start, n)