<?xml version="1.0" encoding="UTF-8"?> <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"> <html> <head> <!-- Generated by HsColour, http://www.cs.york.ac.uk/fp/darcs/hscolour/ --> <title>src/ForSyDe/Shallow/StochasticLib.hs</title> <link type='text/css' rel='stylesheet' href='hscolour.css' /> </head> <body> <pre><a name="line-1"></a><span class='hs-comment'>-----------------------------------------------------------------------------</span> <a name="line-2"></a><span class='hs-comment'>-- |</span> <a name="line-3"></a><span class='hs-comment'>-- Module : ForSyDe.Shallow.StochasticLib</span> <a name="line-4"></a><span class='hs-comment'>-- Copyright : (c) SAM Group, KTH/ICT/ECS 2007-2008</span> <a name="line-5"></a><span class='hs-comment'>-- License : BSD-style (see the file LICENSE)</span> <a name="line-6"></a><span class='hs-comment'>-- </span> <a name="line-7"></a><span class='hs-comment'>-- Maintainer : forsyde-dev@ict.kth.se</span> <a name="line-8"></a><span class='hs-comment'>-- Stability : experimental</span> <a name="line-9"></a><span class='hs-comment'>-- Portability : portable</span> <a name="line-10"></a><span class='hs-comment'>--</span> <a name="line-11"></a><span class='hs-comment'>-- The stochastic library provides a few stochastic skeletons, which are</span> <a name="line-12"></a><span class='hs-comment'>-- relatives to the skeletons of the synchronous library. These skeletons are</span> <a name="line-13"></a><span class='hs-comment'>-- based on two elementary functions, 'sigmaUn' and 'sigmaGe'</span> <a name="line-14"></a><span class='hs-comment'>-- which provide stochastic signals. The background and motivation for this</span> <a name="line-15"></a><span class='hs-comment'>-- approach is described in the paper </span> <a name="line-16"></a><span class='hs-comment'>--</span> <a name="line-17"></a><span class='hs-comment'>-- Axel Jantsch, Ingo Sander, and Wenbiao Wu,</span> <a name="line-18"></a><span class='hs-comment'>-- \"The usage of stochastic processes in embedded system specifications\",</span> <a name="line-19"></a><span class='hs-comment'>-- In /Proceedings of the Ninth International Symposium on Hardware and Software Codesign/, </span> <a name="line-20"></a><span class='hs-comment'>-- April 2001 (<<a href="http://web.it.kth.se/~axel/papers/2001/codes-2001.pdf">http://web.it.kth.se/~axel/papers/2001/codes-2001.pdf</a>>). </span> <a name="line-21"></a><span class='hs-comment'>--</span> <a name="line-22"></a><span class='hs-comment'>-- Unfortunately, not all of the suggested skeletons are</span> <a name="line-23"></a><span class='hs-comment'>-- implemented. In particular, consolidation-based process</span> <a name="line-24"></a><span class='hs-comment'>-- constructors and all constructors for the untimed and the discrete</span> <a name="line-25"></a><span class='hs-comment'>-- timed MoCs are missing.</span> <a name="line-26"></a><span class='hs-comment'>-----------------------------------------------------------------------------</span> <a name="line-27"></a> <a name="line-28"></a><span class='hs-keyword'>module</span> <span class='hs-conid'>ForSyDe</span><span class='hs-varop'>.</span><span class='hs-conid'>Shallow</span><span class='hs-varop'>.</span><span class='hs-conid'>StochasticLib</span><span class='hs-layout'>(</span> <a name="line-29"></a> <span class='hs-comment'>-- * Select based synchronous process constructors</span> <a name="line-30"></a> <span class='hs-varid'>selMapSY</span><span class='hs-layout'>,</span> <span class='hs-varid'>selScanlSY</span><span class='hs-layout'>,</span> <span class='hs-varid'>selMealySY</span><span class='hs-layout'>,</span> <span class='hs-varid'>selMooreSY</span> <a name="line-31"></a> <span class='hs-comment'>-- * Elementary stochastic processes</span> <a name="line-32"></a> <span class='hs-layout'>,</span> <span class='hs-varid'>sigmaUn</span><span class='hs-layout'>,</span> <span class='hs-varid'>sigmaGe</span><span class='hs-layout'>)</span> <span class='hs-keyword'>where</span> <a name="line-33"></a> <a name="line-34"></a><span class='hs-keyword'>import</span> <span class='hs-conid'>ForSyDe</span><span class='hs-varop'>.</span><span class='hs-conid'>Shallow</span><span class='hs-varop'>.</span><span class='hs-conid'>SynchronousLib</span> <a name="line-35"></a><span class='hs-keyword'>import</span> <span class='hs-conid'>ForSyDe</span><span class='hs-varop'>.</span><span class='hs-conid'>Shallow</span><span class='hs-varop'>.</span><span class='hs-conid'>Signal</span> <a name="line-36"></a><span class='hs-keyword'>import</span> <span class='hs-conid'>System</span><span class='hs-varop'>.</span><span class='hs-conid'>Random</span> <a name="line-37"></a> <a name="line-38"></a><a name="selMapSY"></a><span class='hs-comment'>-- | The skeleton 'selMapSY' is a stochastic variant of</span> <a name="line-39"></a><span class='hs-comment'>-- 'mapSY'. It has an internal stochastic process and selects one</span> <a name="line-40"></a><span class='hs-comment'>-- out of two combinatorial functions depending on the output of the</span> <a name="line-41"></a><span class='hs-comment'>-- stochastic process.</span> <a name="line-42"></a><span class='hs-definition'>selMapSY</span> <span class='hs-keyglyph'>::</span> <span class='hs-conid'>Int</span> <span class='hs-comment'>-- ^The seed for the stochastic process</span> <a name="line-43"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^The first alternative function</span> <a name="line-44"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^The second alternative function</span> <a name="line-45"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-varid'>a</span> <span class='hs-comment'>-- ^The input signal</span> <a name="line-46"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-varid'>b</span> <span class='hs-comment'>-- ^The output signal of the process</span> <a name="line-47"></a><span class='hs-definition'>selMapSY</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-conid'>NullS</span> <span class='hs-keyglyph'>=</span> <span class='hs-conid'>NullS</span> <a name="line-48"></a><span class='hs-definition'>selMapSY</span> <span class='hs-varid'>seed</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>xs</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>selmap1</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>f1</span> <span class='hs-layout'>(</span><span class='hs-varid'>sigmaUn</span> <span class='hs-varid'>seed</span> <span class='hs-layout'>(</span><span class='hs-num'>0</span><span class='hs-layout'>,</span><span class='hs-num'>1</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span> <span class='hs-varid'>xs</span> <a name="line-49"></a> <span class='hs-keyword'>where</span> <span class='hs-varid'>selmap1</span> <span class='hs-keyglyph'>::</span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span><span class='hs-keyglyph'>-></span><span class='hs-varid'>b</span><span class='hs-layout'>)</span><span class='hs-keyglyph'>-></span><span class='hs-layout'>(</span><span class='hs-varid'>a</span><span class='hs-keyglyph'>-></span><span class='hs-varid'>b</span><span class='hs-layout'>)</span><span class='hs-keyglyph'>-></span><span class='hs-layout'>(</span><span class='hs-conid'>Signal</span> <span class='hs-conid'>Int</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-varid'>b</span> <a name="line-50"></a> <span class='hs-varid'>selmap1</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-conid'>NullS</span> <span class='hs-keyglyph'>=</span> <span class='hs-conid'>NullS</span> <a name="line-51"></a> <span class='hs-varid'>selmap1</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>f1</span> <span class='hs-layout'>(</span><span class='hs-varid'>s</span><span class='hs-conop'>:-</span><span class='hs-varid'>ss</span><span class='hs-layout'>)</span> <span class='hs-layout'>(</span><span class='hs-varid'>x</span><span class='hs-conop'>:-</span><span class='hs-varid'>xs</span><span class='hs-layout'>)</span> <a name="line-52"></a> <span class='hs-keyglyph'>=</span> <span class='hs-layout'>(</span><span class='hs-varid'>select1</span> <span class='hs-varid'>s</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>x</span><span class='hs-layout'>)</span> <span class='hs-conop'>:-</span> <span class='hs-layout'>(</span><span class='hs-varid'>selmap1</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>ss</span> <span class='hs-varid'>xs</span><span class='hs-layout'>)</span> <a name="line-53"></a> <span class='hs-varid'>selmap1</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-conid'>NullS</span> <span class='hs-keyword'>_</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>error</span> <span class='hs-str'>"selMapSY: empty seed signal."</span> <a name="line-54"></a> <a name="line-55"></a><a name="selScanlSY"></a><span class='hs-comment'>-- | The skeleton 'selScanlSY' is a stochastic variant of 'scanlSY'.</span> <a name="line-56"></a><span class='hs-definition'>selScanlSY</span> <span class='hs-keyglyph'>::</span> <span class='hs-conid'>Int</span> <span class='hs-comment'>-- ^The seed</span> <a name="line-57"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^The first alternative next-state function</span> <a name="line-58"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^The second alternative function</span> <a name="line-59"></a> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span> <span class='hs-comment'>-- ^The initial state</span> <a name="line-60"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-varid'>b</span> <span class='hs-comment'>-- ^The input signal</span> <a name="line-61"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-varid'>a</span> <span class='hs-comment'>-- ^The output signal</span> <a name="line-62"></a><span class='hs-definition'>selScanlSY</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-conid'>NullS</span> <span class='hs-keyglyph'>=</span> <span class='hs-conid'>NullS</span> <a name="line-63"></a><span class='hs-definition'>selScanlSY</span> <span class='hs-varid'>seed</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>mem</span> <span class='hs-varid'>xs</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>selscan1</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>mem</span> <span class='hs-layout'>(</span><span class='hs-varid'>sigmaUn</span> <span class='hs-varid'>seed</span> <span class='hs-layout'>(</span><span class='hs-num'>0</span><span class='hs-layout'>,</span><span class='hs-num'>1</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span> <span class='hs-varid'>xs</span> <a name="line-64"></a> <span class='hs-keyword'>where</span> <span class='hs-varid'>selscan1</span> <span class='hs-keyglyph'>::</span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span> <a name="line-65"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-conid'>Int</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-varid'>b</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-varid'>a</span> <a name="line-66"></a> <span class='hs-varid'>selscan1</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-conid'>NullS</span> <span class='hs-keyglyph'>=</span> <span class='hs-conid'>NullS</span> <a name="line-67"></a> <span class='hs-varid'>selscan1</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>mem</span> <span class='hs-layout'>(</span><span class='hs-varid'>s</span><span class='hs-conop'>:-</span><span class='hs-varid'>ss</span><span class='hs-layout'>)</span> <span class='hs-layout'>(</span><span class='hs-varid'>x</span><span class='hs-conop'>:-</span><span class='hs-varid'>xs</span><span class='hs-layout'>)</span> <a name="line-68"></a> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>newmem</span> <span class='hs-conop'>:-</span> <span class='hs-layout'>(</span><span class='hs-varid'>selscan1</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>newmem</span> <span class='hs-varid'>ss</span> <span class='hs-varid'>xs</span><span class='hs-layout'>)</span> <a name="line-69"></a> <span class='hs-keyword'>where</span> <span class='hs-varid'>newmem</span> <span class='hs-keyglyph'>=</span> <span class='hs-layout'>(</span><span class='hs-varid'>select2</span> <span class='hs-varid'>s</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>mem</span> <span class='hs-varid'>x</span><span class='hs-layout'>)</span> <a name="line-70"></a> <span class='hs-varid'>selscan1</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-conid'>NullS</span> <span class='hs-keyword'>_</span> <a name="line-71"></a> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>error</span> <span class='hs-str'>"selScanlSY: empty seed signal"</span> <a name="line-72"></a> <a name="line-73"></a><a name="select1"></a><span class='hs-definition'>select1</span> <span class='hs-keyglyph'>::</span> <span class='hs-conid'>Int</span> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span><span class='hs-keyglyph'>-></span><span class='hs-varid'>b</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span> <a name="line-74"></a><span class='hs-definition'>select1</span> <span class='hs-num'>0</span> <span class='hs-varid'>f0</span> <span class='hs-keyword'>_</span> <span class='hs-varid'>x</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>x</span> <a name="line-75"></a><span class='hs-definition'>select1</span> <span class='hs-num'>1</span> <span class='hs-keyword'>_</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>x</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>x</span> <a name="line-76"></a><span class='hs-definition'>select1</span> <span class='hs-varid'>s</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>error</span> <span class='hs-layout'>(</span><span class='hs-str'>"select1: seed value neither 0 nor 1: "</span> <a name="line-77"></a> <span class='hs-varop'>++</span> <span class='hs-layout'>(</span><span class='hs-varid'>show</span> <span class='hs-varid'>s</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span> <a name="line-78"></a> <a name="line-79"></a><a name="select2"></a><span class='hs-definition'>select2</span> <span class='hs-keyglyph'>::</span> <span class='hs-conid'>Int</span> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>c</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span><span class='hs-keyglyph'>-></span><span class='hs-varid'>b</span><span class='hs-keyglyph'>-></span><span class='hs-varid'>c</span><span class='hs-layout'>)</span> <a name="line-80"></a> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>c</span> <a name="line-81"></a><span class='hs-definition'>select2</span> <span class='hs-num'>0</span> <span class='hs-varid'>f0</span> <span class='hs-keyword'>_</span> <span class='hs-varid'>x</span> <span class='hs-varid'>y</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>x</span> <span class='hs-varid'>y</span> <a name="line-82"></a><span class='hs-definition'>select2</span> <span class='hs-num'>1</span> <span class='hs-keyword'>_</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>x</span> <span class='hs-varid'>y</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>x</span> <span class='hs-varid'>y</span> <a name="line-83"></a><span class='hs-definition'>select2</span> <span class='hs-varid'>s</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>error</span> <span class='hs-layout'>(</span><span class='hs-str'>"select2: seed value neither 0 nor 1: "</span> <a name="line-84"></a> <span class='hs-varop'>++</span> <span class='hs-layout'>(</span><span class='hs-varid'>show</span> <span class='hs-varid'>s</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span> <a name="line-85"></a> <a name="line-86"></a><a name="selMooreSY"></a><span class='hs-comment'>-- | 'selMooreSY' is the stochastic variant of mooreSY. Both the </span> <a name="line-87"></a><span class='hs-comment'>-- next-state and the output function is randomly selected based on a </span> <a name="line-88"></a><span class='hs-comment'>-- uniform distribution.</span> <a name="line-89"></a><span class='hs-definition'>selMooreSY</span> <span class='hs-keyglyph'>::</span> <span class='hs-conid'>Int</span> <span class='hs-comment'>-- ^The seed for the next-state function</span> <a name="line-90"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Int</span> <span class='hs-comment'>-- ^The seed for the output function</span> <a name="line-91"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^First alternative for the next-state function</span> <a name="line-92"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^Second alternative for the next-state function</span> <a name="line-93"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>c</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^First alternative for the output function</span> <a name="line-94"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>c</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^Second alternative for the output function</span> <a name="line-95"></a> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span> <span class='hs-comment'>-- ^The initial state</span> <a name="line-96"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-varid'>b</span> <span class='hs-comment'>-- ^The input signal</span> <a name="line-97"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-varid'>c</span> <span class='hs-comment'>-- ^The output signal</span> <a name="line-98"></a><span class='hs-definition'>selMooreSY</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-conid'>NullS</span> <span class='hs-keyglyph'>=</span> <span class='hs-conid'>NullS</span> <a name="line-99"></a><span class='hs-definition'>selMooreSY</span> <span class='hs-varid'>seedg</span> <span class='hs-varid'>seedf</span> <span class='hs-varid'>g0</span> <span class='hs-varid'>g1</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>w0</span> <span class='hs-varid'>s</span> <a name="line-100"></a> <span class='hs-keyglyph'>=</span> <span class='hs-layout'>(</span><span class='hs-layout'>(</span><span class='hs-varid'>selMapSY</span> <span class='hs-varid'>seedf</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>f1</span> <span class='hs-layout'>)</span> <span class='hs-varop'>.</span> <span class='hs-layout'>(</span><span class='hs-varid'>selScanlSY</span> <span class='hs-varid'>seedg</span> <span class='hs-varid'>g0</span> <span class='hs-varid'>g1</span> <span class='hs-varid'>w0</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span> <span class='hs-varid'>s</span> <a name="line-101"></a> <a name="line-102"></a><a name="selMealySY"></a><span class='hs-comment'>-- | 'selMealySY' is the stochastic variant of mealySY. Both the </span> <a name="line-103"></a><span class='hs-comment'>-- next-state and the output function is randomly selected based on a </span> <a name="line-104"></a><span class='hs-comment'>-- uniform distribution.</span> <a name="line-105"></a><span class='hs-definition'>selMealySY</span> <span class='hs-keyglyph'>::</span> <span class='hs-conid'>Int</span> <span class='hs-comment'>-- ^The seed for the next-state function</span> <a name="line-106"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Int</span> <span class='hs-comment'>-- ^The seed for the output function</span> <a name="line-107"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^First alternative for the next-state function</span> <a name="line-108"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^Second alternative for the next-state function</span> <a name="line-109"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>c</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^First alternative for the output function</span> <a name="line-110"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-varid'>a</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>b</span> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>c</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^Second alternative for the output function</span> <a name="line-111"></a> <span class='hs-keyglyph'>-></span> <span class='hs-varid'>a</span> <span class='hs-comment'>-- ^The initial state</span> <a name="line-112"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-varid'>b</span> <span class='hs-comment'>-- ^The input signal</span> <a name="line-113"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-varid'>c</span> <span class='hs-comment'>-- ^The output signal</span> <a name="line-114"></a><span class='hs-definition'>selMealySY</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-keyword'>_</span> <span class='hs-conid'>NullS</span> <span class='hs-keyglyph'>=</span> <span class='hs-conid'>NullS</span> <a name="line-115"></a><span class='hs-definition'>selMealySY</span> <span class='hs-varid'>seedg</span> <span class='hs-varid'>seedf</span> <span class='hs-varid'>g0</span> <span class='hs-varid'>g1</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>w0</span> <span class='hs-varid'>s</span> <a name="line-116"></a> <span class='hs-keyglyph'>=</span> <span class='hs-layout'>(</span><span class='hs-layout'>(</span><span class='hs-varid'>selMapSY</span> <span class='hs-varid'>seedf</span> <span class='hs-varid'>f0'</span> <span class='hs-varid'>f1'</span> <span class='hs-layout'>)</span> <span class='hs-varop'>.</span> <span class='hs-layout'>(</span><span class='hs-varid'>zipSY</span> <span class='hs-varid'>s</span><span class='hs-layout'>)</span> <span class='hs-varop'>.</span> <span class='hs-layout'>(</span><span class='hs-varid'>selScanlSY</span> <span class='hs-varid'>seedg</span> <span class='hs-varid'>g0</span> <span class='hs-varid'>g1</span> <span class='hs-varid'>w0</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span> <span class='hs-varid'>s</span> <a name="line-117"></a> <span class='hs-keyword'>where</span> <span class='hs-varid'>f0'</span> <span class='hs-layout'>(</span><span class='hs-varid'>b</span><span class='hs-layout'>,</span> <span class='hs-varid'>a</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>f0</span> <span class='hs-varid'>a</span> <span class='hs-varid'>b</span> <a name="line-118"></a> <span class='hs-varid'>f1'</span> <span class='hs-layout'>(</span><span class='hs-varid'>b</span><span class='hs-layout'>,</span> <span class='hs-varid'>a</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>f1</span> <span class='hs-varid'>a</span> <span class='hs-varid'>b</span> <a name="line-119"></a><a name="sigmaUn"></a><span class='hs-comment'>-- |'sigmaUn' generates a signal list of uniformly distributed Int within</span> <a name="line-120"></a><span class='hs-comment'>-- the given range and with a given seed. </span> <a name="line-121"></a><span class='hs-definition'>sigmaUn</span> <span class='hs-keyglyph'>::</span> <span class='hs-conid'>Int</span> <span class='hs-comment'>-- ^The seed</span> <a name="line-122"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-conid'>Int</span><span class='hs-layout'>,</span> <span class='hs-conid'>Int</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^The interval from which the stochastic values are </span> <a name="line-123"></a> <span class='hs-comment'>-- taken</span> <a name="line-124"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-conid'>Int</span> <span class='hs-comment'>-- ^The sequence of stochastic values</span> <a name="line-125"></a><span class='hs-definition'>sigmaUn</span> <span class='hs-varid'>seed</span> <span class='hs-varid'>range</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>signal</span> <span class='hs-layout'>(</span><span class='hs-varid'>stoch</span> <span class='hs-varid'>range</span> <span class='hs-layout'>(</span><span class='hs-varid'>mkStdGen</span> <span class='hs-varid'>seed</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span> <a name="line-126"></a> <span class='hs-keyword'>where</span> <span class='hs-varid'>stoch</span> <span class='hs-keyglyph'>::</span> <span class='hs-layout'>(</span><span class='hs-conid'>Int</span><span class='hs-layout'>,</span> <span class='hs-conid'>Int</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>StdGen</span> <span class='hs-keyglyph'>-></span> <span class='hs-keyglyph'>[</span><span class='hs-conid'>Int</span><span class='hs-keyglyph'>]</span> <a name="line-127"></a> <span class='hs-varid'>stoch</span> <span class='hs-varid'>range</span> <span class='hs-varid'>g</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>newNum</span> <span class='hs-varop'>`seq`</span> <a name="line-128"></a> <span class='hs-layout'>(</span><span class='hs-varid'>newNum</span> <span class='hs-conop'>:</span> <span class='hs-layout'>(</span><span class='hs-varid'>stoch</span> <span class='hs-varid'>range</span> <span class='hs-varid'>newGen</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span> <a name="line-129"></a> <span class='hs-keyword'>where</span> <span class='hs-varid'>newNum</span> <span class='hs-keyglyph'>=</span> <span class='hs-layout'>(</span><span class='hs-varid'>fst</span> <span class='hs-layout'>(</span><span class='hs-varid'>randomR</span> <span class='hs-varid'>range</span> <span class='hs-varid'>g</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span> <a name="line-130"></a> <span class='hs-varid'>newGen</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>snd</span> <span class='hs-layout'>(</span><span class='hs-varid'>next</span> <span class='hs-varid'>g</span><span class='hs-layout'>)</span> <a name="line-131"></a> <a name="line-132"></a><a name="sigmaGe"></a><span class='hs-comment'>-- |'sigmaGe' is a more general stochastic process. The first argument is a</span> <a name="line-133"></a><span class='hs-comment'>-- function f which describes the distribution. For each value v in the</span> <a name="line-134"></a><span class='hs-comment'>-- given range (r1,r2), f(v) is the probability that v is generated. </span> <a name="line-135"></a><span class='hs-comment'>--</span> <a name="line-136"></a><span class='hs-comment'>-- Note, that the user has to make sure that sum(f(v))=1 for v in (r1,r2).</span> <a name="line-137"></a><span class='hs-comment'>--</span> <a name="line-138"></a><span class='hs-comment'>-- For illustration consider the following example.</span> <a name="line-139"></a><span class='hs-comment'>--</span> <a name="line-140"></a><span class='hs-comment'>-- > pdist :: Float -> Float</span> <a name="line-141"></a><span class='hs-comment'>-- > pdist d = 1\/\(2**d\)</span> <a name="line-142"></a><span class='hs-comment'>-- > pdistsum 1 = pdist 1</span> <a name="line-143"></a><span class='hs-comment'>-- > pdistsum d = \(pdist d\) + \(pdistsum \(d-1\)\)</span> <a name="line-144"></a><span class='hs-comment'>--</span> <a name="line-145"></a><span class='hs-comment'>-- > pdistnorm :: Float -> Float -> Float</span> <a name="line-146"></a><span class='hs-comment'>-- > pdistnorm dmax d = 1\/((pdistsum dmax) * (2**d))</span> <a name="line-147"></a><span class='hs-comment'>--</span> <a name="line-148"></a><span class='hs-comment'>-- @pdistnorm dmax d@ gives the probability of a value <= d;</span> <a name="line-149"></a><span class='hs-comment'>--</span> <a name="line-150"></a><span class='hs-comment'>-- @pdistnorm dmax dmax@ is always 1.0</span> <a name="line-151"></a><span class='hs-comment'>--</span> <a name="line-152"></a><span class='hs-comment'>-- Hence, using pdistnorm as a function in 'sigmaGe' gives an exponantial</span> <a name="line-153"></a><span class='hs-comment'>-- distribution for values in the range \[0, dmax\].</span> <a name="line-154"></a><span class='hs-definition'>sigmaGe</span> <span class='hs-keyglyph'>::</span> <span class='hs-layout'>(</span><span class='hs-conid'>Float</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Float</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^The stochastic distribution</span> <a name="line-155"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Int</span> <span class='hs-comment'>-- ^The seed</span> <a name="line-156"></a> <span class='hs-keyglyph'>-></span> <span class='hs-layout'>(</span><span class='hs-conid'>Int</span><span class='hs-layout'>,</span> <span class='hs-conid'>Int</span><span class='hs-layout'>)</span> <span class='hs-comment'>-- ^The range</span> <a name="line-157"></a> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Signal</span> <span class='hs-conid'>Int</span> <span class='hs-comment'>-- ^The sequence of stochastic values</span> <a name="line-158"></a><span class='hs-definition'>sigmaGe</span> <span class='hs-varid'>f</span> <span class='hs-varid'>seed</span> <span class='hs-layout'>(</span><span class='hs-varid'>r1</span><span class='hs-layout'>,</span><span class='hs-varid'>r2</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>sigma2</span> <span class='hs-layout'>(</span><span class='hs-varid'>checkSum</span> <span class='hs-varid'>f</span> <span class='hs-layout'>(</span><span class='hs-varid'>fromIntegral</span> <span class='hs-varid'>r1</span><span class='hs-layout'>)</span> <a name="line-159"></a> <span class='hs-layout'>(</span><span class='hs-varid'>fromIntegral</span> <span class='hs-varid'>r2</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span> <span class='hs-varid'>f</span> <span class='hs-varid'>seed</span> <span class='hs-layout'>(</span><span class='hs-varid'>r1</span><span class='hs-layout'>,</span><span class='hs-varid'>r2</span><span class='hs-layout'>)</span> <a name="line-160"></a> <span class='hs-keyword'>where</span> <span class='hs-varid'>sigma2</span> <span class='hs-varid'>s</span> <span class='hs-varid'>f</span> <span class='hs-varid'>seed</span> <span class='hs-layout'>(</span><span class='hs-varid'>r1</span><span class='hs-layout'>,</span><span class='hs-varid'>r2</span><span class='hs-layout'>)</span> <a name="line-161"></a> <span class='hs-keyglyph'>|</span> <span class='hs-varid'>s</span> <span class='hs-varop'>></span> <span class='hs-num'>0.999</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>signal</span> <span class='hs-layout'>(</span><span class='hs-varid'>sigma1</span> <span class='hs-layout'>(</span><span class='hs-varid'>mkStdGen</span> <span class='hs-varid'>seed</span><span class='hs-layout'>)</span> <a name="line-162"></a> <span class='hs-layout'>(</span><span class='hs-varid'>mkdlist</span> <span class='hs-varid'>f</span> <span class='hs-layout'>(</span><span class='hs-varid'>fromIntegral</span> <span class='hs-layout'>(</span><span class='hs-varid'>r2</span><span class='hs-comment'>-</span><span class='hs-varid'>r1</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span> <a name="line-163"></a> <span class='hs-keyglyph'>|</span> <span class='hs-varid'>otherwise</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>error</span> <a name="line-164"></a> <span class='hs-layout'>(</span><span class='hs-str'>"sigmaGe: sum of probabilitites is "</span> <a name="line-165"></a> <span class='hs-varop'>++</span> <span class='hs-layout'>(</span><span class='hs-varid'>show</span> <span class='hs-varid'>s</span><span class='hs-layout'>)</span> <span class='hs-varop'>++</span> <span class='hs-str'>". It must be 1."</span><span class='hs-layout'>)</span> <a name="line-166"></a> <span class='hs-varid'>checkSum</span> <span class='hs-keyglyph'>::</span> <span class='hs-layout'>(</span><span class='hs-conid'>Float</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Float</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Float</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Float</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Float</span> <a name="line-167"></a> <span class='hs-varid'>checkSum</span> <span class='hs-varid'>f</span> <span class='hs-varid'>c</span> <span class='hs-varid'>max</span> <span class='hs-keyglyph'>|</span> <span class='hs-varid'>c</span> <span class='hs-varop'>==</span> <span class='hs-varid'>max</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>f</span> <span class='hs-varid'>c</span> <a name="line-168"></a> <span class='hs-keyglyph'>|</span> <span class='hs-varid'>otherwise</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>f</span><span class='hs-layout'>(</span><span class='hs-varid'>c</span><span class='hs-layout'>)</span> <span class='hs-varop'>+</span> <span class='hs-layout'>(</span><span class='hs-varid'>checkSum</span> <span class='hs-varid'>f</span> <span class='hs-layout'>(</span><span class='hs-varid'>c</span><span class='hs-varop'>+</span><span class='hs-num'>1</span><span class='hs-layout'>)</span> <span class='hs-varid'>max</span><span class='hs-layout'>)</span> <a name="line-169"></a> <a name="line-170"></a> <span class='hs-varid'>sigma1</span> <span class='hs-keyglyph'>::</span> <span class='hs-conid'>StdGen</span> <span class='hs-keyglyph'>-></span> <span class='hs-keyglyph'>[</span><span class='hs-conid'>Float</span><span class='hs-keyglyph'>]</span> <span class='hs-keyglyph'>-></span> <span class='hs-keyglyph'>[</span><span class='hs-conid'>Int</span><span class='hs-keyglyph'>]</span> <a name="line-171"></a> <span class='hs-varid'>sigma1</span> <span class='hs-varid'>g</span> <span class='hs-varid'>fl</span> <span class='hs-keyglyph'>=</span> <span class='hs-layout'>(</span><span class='hs-varid'>findk</span> <span class='hs-layout'>(</span><span class='hs-varid'>fst</span> <span class='hs-layout'>(</span><span class='hs-varid'>randomR</span> <span class='hs-layout'>(</span><span class='hs-num'>0.0</span><span class='hs-layout'>,</span><span class='hs-num'>1.0</span><span class='hs-layout'>)</span> <span class='hs-varid'>g</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span> <span class='hs-varid'>fl</span><span class='hs-layout'>)</span> <a name="line-172"></a> <span class='hs-conop'>:</span> <span class='hs-layout'>(</span><span class='hs-varid'>sigma1</span> <span class='hs-layout'>(</span><span class='hs-varid'>snd</span> <span class='hs-layout'>(</span><span class='hs-varid'>next</span> <span class='hs-varid'>g</span><span class='hs-layout'>)</span><span class='hs-layout'>)</span> <span class='hs-varid'>fl</span><span class='hs-layout'>)</span> <a name="line-173"></a> <a name="line-174"></a> <span class='hs-varid'>findk</span> <span class='hs-keyglyph'>::</span> <span class='hs-conid'>Float</span> <span class='hs-keyglyph'>-></span> <span class='hs-keyglyph'>[</span><span class='hs-conid'>Float</span><span class='hs-keyglyph'>]</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Int</span> <a name="line-175"></a> <span class='hs-varid'>findk</span> <span class='hs-varid'>r</span> <span class='hs-varid'>fs</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>findk1</span> <span class='hs-num'>0</span> <span class='hs-varid'>r</span> <span class='hs-varid'>fs</span> <a name="line-176"></a> <a name="line-177"></a> <span class='hs-varid'>findk1</span> <span class='hs-varid'>k</span> <span class='hs-varid'>r</span> <span class='hs-layout'>(</span><span class='hs-varid'>f</span><span class='hs-conop'>:</span><span class='hs-varid'>fs</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>|</span> <span class='hs-varid'>r</span> <span class='hs-varop'><</span> <span class='hs-varid'>f</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>k</span> <a name="line-178"></a> <span class='hs-keyglyph'>|</span> <span class='hs-varid'>otherwise</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>findk1</span> <span class='hs-layout'>(</span><span class='hs-varid'>k</span><span class='hs-varop'>+</span><span class='hs-num'>1</span><span class='hs-layout'>)</span> <span class='hs-varid'>r</span> <span class='hs-varid'>fs</span> <a name="line-179"></a> <span class='hs-varid'>findk1</span> <span class='hs-varid'>k</span> <span class='hs-keyword'>_</span> <span class='hs-conid'>[]</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>k</span> <a name="line-180"></a> <a name="line-181"></a> <span class='hs-varid'>mkdlist</span> <span class='hs-keyglyph'>::</span> <span class='hs-layout'>(</span><span class='hs-conid'>Float</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Float</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Float</span> <span class='hs-keyglyph'>-></span> <span class='hs-keyglyph'>[</span><span class='hs-conid'>Float</span><span class='hs-keyglyph'>]</span> <a name="line-182"></a> <span class='hs-varid'>mkdlist</span> <span class='hs-varid'>f</span> <span class='hs-varid'>d</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>scanl</span> <span class='hs-layout'>(</span><span class='hs-varid'>sumf</span> <span class='hs-varid'>f</span><span class='hs-layout'>)</span> <span class='hs-num'>0.0</span> <span class='hs-keyglyph'>[</span><span class='hs-num'>1</span><span class='hs-keyglyph'>..</span><span class='hs-varid'>d</span><span class='hs-keyglyph'>]</span> <a name="line-183"></a> <a name="line-184"></a> <span class='hs-varid'>sumf</span> <span class='hs-keyglyph'>::</span> <span class='hs-layout'>(</span><span class='hs-conid'>Float</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Float</span><span class='hs-layout'>)</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Float</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Float</span> <span class='hs-keyglyph'>-></span> <span class='hs-conid'>Float</span> <a name="line-185"></a> <span class='hs-varid'>sumf</span> <span class='hs-varid'>g</span> <span class='hs-varid'>x</span> <span class='hs-varid'>y</span> <span class='hs-keyglyph'>=</span> <span class='hs-varid'>x</span> <span class='hs-varop'>+</span> <span class='hs-layout'>(</span><span class='hs-varid'>g</span> <span class='hs-varid'>y</span><span class='hs-layout'>)</span> <a name="line-186"></a> <a name="line-187"></a><span class='hs-comment'>--pdist :: Float -> Float</span> <a name="line-188"></a><span class='hs-comment'>--pdist d = 1/(2**d)</span> <a name="line-189"></a> <a name="line-190"></a><span class='hs-comment'>--pdistsum 1 = pdist 1</span> <a name="line-191"></a><span class='hs-comment'>--pdistsum d = (pdist d) + (pdistsum (d-1))</span> <a name="line-192"></a> <a name="line-193"></a><span class='hs-comment'>-- The function pdistnorm can be used as a function in sigmaGe for an</span> <a name="line-194"></a><span class='hs-comment'>-- exponantial distribution of values in the range [0, dmax]:</span> <a name="line-195"></a><span class='hs-comment'>--pdistnorm :: Float -> Float -> Float</span> <a name="line-196"></a><span class='hs-comment'>--pdistnorm dmax d = 1/((pdistsum dmax) * (2**d))</span> <a name="line-197"></a> <a name="line-198"></a><span class='hs-comment'>--pdnormsum dmax 1 = pdistnorm dmax 1</span> <a name="line-199"></a><span class='hs-comment'>--pdnormsum dmax d = (pdistnorm dmax d) + (pdnormsum dmax (d-1))</span> <a name="line-200"></a> <a name="line-201"></a> <a name="line-202"></a><span class='hs-comment'>-----------------------------------------------------------------------------</span> <a name="line-203"></a><span class='hs-comment'>-- Test section:</span> <a name="line-204"></a><span class='hs-comment'>--</span> <a name="line-205"></a><span class='hs-comment'>-- These tests are commented to avoid warnings about not-used functions.</span> <a name="line-206"></a><span class='hs-comment'>-- But the test functions work and are useful.</span> <a name="line-207"></a><span class='hs-comment'>-- testAll = "test selMapSY: " ++ testSelMap </span> <a name="line-208"></a><span class='hs-comment'>-- ++ ", test selMooreSY: " ++ testSelMoore</span> <a name="line-209"></a><span class='hs-comment'>-- ++ ", test selMealySY: " ++ testSelMealy</span> <a name="line-210"></a> <a name="line-211"></a><span class='hs-comment'>-- testSelMap = show so ++ " -> " ++ (cmpSig so (signal [0,3,4,5,4,5,8,9,8,11]))</span> <a name="line-212"></a><span class='hs-comment'>-- where f0 x = x + 1</span> <a name="line-213"></a><span class='hs-comment'>-- f1 x = x - 1</span> <a name="line-214"></a><span class='hs-comment'>-- so = takeS 10 (selMapSY 876876 f0 f1 (signal [1,2..]))</span> <a name="line-215"></a> <a name="line-216"></a><span class='hs-comment'>-- testSelMoore = show so ++ " -> " </span> <a name="line-217"></a><span class='hs-comment'>-- ++ (cmpSig so (signal [10,2,3,-40,-5,0,7,-80,0,-100]))</span> <a name="line-218"></a><span class='hs-comment'>-- where so = takeS 10 (selMooreSY 7667567 123234 g0 g1 f0 f1 w0 </span> <a name="line-219"></a><span class='hs-comment'>-- (signal [1,2..]))</span> <a name="line-220"></a><span class='hs-comment'>-- g0 (0,y) x | even x = (0,x)</span> <a name="line-221"></a><span class='hs-comment'>-- | otherwise = (1,x)</span> <a name="line-222"></a><span class='hs-comment'>-- g0 (1,y) x | x `mod` 3 == 0 = (0,x)</span> <a name="line-223"></a><span class='hs-comment'>-- | otherwise = (1,x)</span> <a name="line-224"></a><span class='hs-comment'>-- g1 (0,y) x | even x = (1,x)</span> <a name="line-225"></a><span class='hs-comment'>-- | otherwise = (0, x)</span> <a name="line-226"></a><span class='hs-comment'>-- g1 (1,y) x | x `mod` 3 == 0 = (0,0)</span> <a name="line-227"></a><span class='hs-comment'>-- | otherwise = (1,x)</span> <a name="line-228"></a><span class='hs-comment'>-- f0 (0,y) = y</span> <a name="line-229"></a><span class='hs-comment'>-- f0 (1,y) = -1 * y</span> <a name="line-230"></a><span class='hs-comment'>-- f1 (0,y) = 10 * y</span> <a name="line-231"></a><span class='hs-comment'>-- f1 (1,y) = -10 * y</span> <a name="line-232"></a><span class='hs-comment'>-- w0 = (0,0)</span> <a name="line-233"></a> <a name="line-234"></a><span class='hs-comment'>-- testSelMealy = show so ++ " -> " </span> <a name="line-235"></a><span class='hs-comment'>-- ++ (cmpSig so (signal [10,2,3,-40,-5,0,7,-80,0,-100]))</span> <a name="line-236"></a><span class='hs-comment'>-- where so = takeS 10 (selMealySY 7667567 123234 g0 g1 f0 f1 w0 </span> <a name="line-237"></a><span class='hs-comment'>-- (signal [1,2..]))</span> <a name="line-238"></a><span class='hs-comment'>-- g0 (0,y) x | even x = (0,x)</span> <a name="line-239"></a><span class='hs-comment'>-- | otherwise = (1,x)</span> <a name="line-240"></a><span class='hs-comment'>-- g0 (1,y) x | x `mod` 3 == 0 = (0,x)</span> <a name="line-241"></a><span class='hs-comment'>-- | otherwise = (1,x)</span> <a name="line-242"></a><span class='hs-comment'>-- g1 (0,y) x | even x = (1,x)</span> <a name="line-243"></a><span class='hs-comment'>-- | otherwise = (0, x)</span> <a name="line-244"></a><span class='hs-comment'>-- g1 (1,y) x | x `mod` 3 == 0 = (0,0)</span> <a name="line-245"></a><span class='hs-comment'>-- | otherwise = (1,x)</span> <a name="line-246"></a><span class='hs-comment'>-- f0 (0,y) x = y</span> <a name="line-247"></a><span class='hs-comment'>-- f0 (1,y) x = -1 * y</span> <a name="line-248"></a><span class='hs-comment'>-- f1 (0,y) x = 10 * y</span> <a name="line-249"></a><span class='hs-comment'>-- f1 (1,y) x = -10 * y</span> <a name="line-250"></a><span class='hs-comment'>-- w0 = (0,0)</span> <a name="line-251"></a><span class='hs-comment'>--</span> <a name="line-252"></a><span class='hs-comment'>--</span> <a name="line-253"></a><span class='hs-comment'>-- cmpSig :: Eq a => Signal a -> Signal a -> String</span> <a name="line-254"></a><span class='hs-comment'>-- cmpSig s1 s2 | s1 == s2 = "OK"</span> <a name="line-255"></a><span class='hs-comment'>-- | otherwise = "Not OK"</span> </pre></body> </html>