module Simulation.Aivika.Process.Random
(randomUniformProcess,
randomUniformProcess_,
randomUniformIntProcess,
randomUniformIntProcess_,
randomTriangularProcess,
randomTriangularProcess_,
randomNormalProcess,
randomNormalProcess_,
randomLogNormalProcess,
randomLogNormalProcess_,
randomExponentialProcess,
randomExponentialProcess_,
randomErlangProcess,
randomErlangProcess_,
randomPoissonProcess,
randomPoissonProcess_,
randomBinomialProcess,
randomBinomialProcess_,
randomGammaProcess,
randomGammaProcess_,
randomBetaProcess,
randomBetaProcess_,
randomWeibullProcess,
randomWeibullProcess_,
randomDiscreteProcess,
randomDiscreteProcess_) where
import Control.Monad
import Control.Monad.Trans
import Simulation.Aivika.Generator
import Simulation.Aivika.Parameter
import Simulation.Aivika.Parameter.Random
import Simulation.Aivika.Process
randomUniformProcess :: Double
-> Double
-> Process Double
randomUniformProcess :: Double -> Double -> Process Double
randomUniformProcess Double
min Double
max =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Parameter Double
randomUniform Double
min Double
max
Double -> Process ()
holdProcess Double
t
Double -> Process Double
forall a. a -> Process a
forall (m :: * -> *) a. Monad m => a -> m a
return Double
t
randomUniformProcess_ :: Double
-> Double
-> Process ()
randomUniformProcess_ :: Double -> Double -> Process ()
randomUniformProcess_ Double
min Double
max =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Parameter Double
randomUniform Double
min Double
max
Double -> Process ()
holdProcess Double
t
randomUniformIntProcess :: Int
-> Int
-> Process Int
randomUniformIntProcess :: Int -> Int -> Process Int
randomUniformIntProcess Int
min Int
max =
do Int
t <- Parameter Int -> Process Int
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Int -> Process Int) -> Parameter Int -> Process Int
forall a b. (a -> b) -> a -> b
$ Int -> Int -> Parameter Int
randomUniformInt Int
min Int
max
Double -> Process ()
holdProcess (Double -> Process ()) -> Double -> Process ()
forall a b. (a -> b) -> a -> b
$ Int -> Double
forall a b. (Integral a, Num b) => a -> b
fromIntegral Int
t
Int -> Process Int
forall a. a -> Process a
forall (m :: * -> *) a. Monad m => a -> m a
return Int
t
randomUniformIntProcess_ :: Int
-> Int
-> Process ()
randomUniformIntProcess_ :: Int -> Int -> Process ()
randomUniformIntProcess_ Int
min Int
max =
do Int
t <- Parameter Int -> Process Int
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Int -> Process Int) -> Parameter Int -> Process Int
forall a b. (a -> b) -> a -> b
$ Int -> Int -> Parameter Int
randomUniformInt Int
min Int
max
Double -> Process ()
holdProcess (Double -> Process ()) -> Double -> Process ()
forall a b. (a -> b) -> a -> b
$ Int -> Double
forall a b. (Integral a, Num b) => a -> b
fromIntegral Int
t
randomTriangularProcess :: Double
-> Double
-> Double
-> Process Double
randomTriangularProcess :: Double -> Double -> Double -> Process Double
randomTriangularProcess Double
min Double
median Double
max =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Double -> Parameter Double
randomTriangular Double
min Double
median Double
max
Double -> Process ()
holdProcess Double
t
Double -> Process Double
forall a. a -> Process a
forall (m :: * -> *) a. Monad m => a -> m a
return Double
t
randomTriangularProcess_ :: Double
-> Double
-> Double
-> Process ()
randomTriangularProcess_ :: Double -> Double -> Double -> Process ()
randomTriangularProcess_ Double
min Double
median Double
max =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Double -> Parameter Double
randomTriangular Double
min Double
median Double
max
Double -> Process ()
holdProcess Double
t
randomNormalProcess :: Double
-> Double
-> Process Double
randomNormalProcess :: Double -> Double -> Process Double
randomNormalProcess Double
mu Double
nu =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Parameter Double
randomNormal Double
mu Double
nu
Bool -> Process () -> Process ()
forall (f :: * -> *). Applicative f => Bool -> f () -> f ()
when (Double
t Double -> Double -> Bool
forall a. Ord a => a -> a -> Bool
> Double
0) (Process () -> Process ()) -> Process () -> Process ()
forall a b. (a -> b) -> a -> b
$
Double -> Process ()
holdProcess Double
t
Double -> Process Double
forall a. a -> Process a
forall (m :: * -> *) a. Monad m => a -> m a
return Double
t
randomNormalProcess_ :: Double
-> Double
-> Process ()
randomNormalProcess_ :: Double -> Double -> Process ()
randomNormalProcess_ Double
mu Double
nu =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Parameter Double
randomNormal Double
mu Double
nu
Bool -> Process () -> Process ()
forall (f :: * -> *). Applicative f => Bool -> f () -> f ()
when (Double
t Double -> Double -> Bool
forall a. Ord a => a -> a -> Bool
> Double
0) (Process () -> Process ()) -> Process () -> Process ()
forall a b. (a -> b) -> a -> b
$
Double -> Process ()
holdProcess Double
t
randomLogNormalProcess :: Double
-> Double
-> Process Double
randomLogNormalProcess :: Double -> Double -> Process Double
randomLogNormalProcess Double
mu Double
nu =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Parameter Double
randomLogNormal Double
mu Double
nu
Double -> Process ()
holdProcess Double
t
Double -> Process Double
forall a. a -> Process a
forall (m :: * -> *) a. Monad m => a -> m a
return Double
t
randomLogNormalProcess_ :: Double
-> Double
-> Process ()
randomLogNormalProcess_ :: Double -> Double -> Process ()
randomLogNormalProcess_ Double
mu Double
nu =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Parameter Double
randomLogNormal Double
mu Double
nu
Double -> Process ()
holdProcess Double
t
randomExponentialProcess :: Double
-> Process Double
randomExponentialProcess :: Double -> Process Double
randomExponentialProcess Double
mu =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Parameter Double
randomExponential Double
mu
Double -> Process ()
holdProcess Double
t
Double -> Process Double
forall a. a -> Process a
forall (m :: * -> *) a. Monad m => a -> m a
return Double
t
randomExponentialProcess_ :: Double
-> Process ()
randomExponentialProcess_ :: Double -> Process ()
randomExponentialProcess_ Double
mu =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Parameter Double
randomExponential Double
mu
Double -> Process ()
holdProcess Double
t
randomErlangProcess :: Double
-> Int
-> Process Double
randomErlangProcess :: Double -> Int -> Process Double
randomErlangProcess Double
beta Int
m =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Int -> Parameter Double
randomErlang Double
beta Int
m
Double -> Process ()
holdProcess Double
t
Double -> Process Double
forall a. a -> Process a
forall (m :: * -> *) a. Monad m => a -> m a
return Double
t
randomErlangProcess_ :: Double
-> Int
-> Process ()
randomErlangProcess_ :: Double -> Int -> Process ()
randomErlangProcess_ Double
beta Int
m =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Int -> Parameter Double
randomErlang Double
beta Int
m
Double -> Process ()
holdProcess Double
t
randomPoissonProcess :: Double
-> Process Int
randomPoissonProcess :: Double -> Process Int
randomPoissonProcess Double
mu =
do Int
t <- Parameter Int -> Process Int
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Int -> Process Int) -> Parameter Int -> Process Int
forall a b. (a -> b) -> a -> b
$ Double -> Parameter Int
randomPoisson Double
mu
Double -> Process ()
holdProcess (Double -> Process ()) -> Double -> Process ()
forall a b. (a -> b) -> a -> b
$ Int -> Double
forall a b. (Integral a, Num b) => a -> b
fromIntegral Int
t
Int -> Process Int
forall a. a -> Process a
forall (m :: * -> *) a. Monad m => a -> m a
return Int
t
randomPoissonProcess_ :: Double
-> Process ()
randomPoissonProcess_ :: Double -> Process ()
randomPoissonProcess_ Double
mu =
do Int
t <- Parameter Int -> Process Int
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Int -> Process Int) -> Parameter Int -> Process Int
forall a b. (a -> b) -> a -> b
$ Double -> Parameter Int
randomPoisson Double
mu
Double -> Process ()
holdProcess (Double -> Process ()) -> Double -> Process ()
forall a b. (a -> b) -> a -> b
$ Int -> Double
forall a b. (Integral a, Num b) => a -> b
fromIntegral Int
t
randomBinomialProcess :: Double
-> Int
-> Process Int
randomBinomialProcess :: Double -> Int -> Process Int
randomBinomialProcess Double
prob Int
trials =
do Int
t <- Parameter Int -> Process Int
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Int -> Process Int) -> Parameter Int -> Process Int
forall a b. (a -> b) -> a -> b
$ Double -> Int -> Parameter Int
randomBinomial Double
prob Int
trials
Double -> Process ()
holdProcess (Double -> Process ()) -> Double -> Process ()
forall a b. (a -> b) -> a -> b
$ Int -> Double
forall a b. (Integral a, Num b) => a -> b
fromIntegral Int
t
Int -> Process Int
forall a. a -> Process a
forall (m :: * -> *) a. Monad m => a -> m a
return Int
t
randomBinomialProcess_ :: Double
-> Int
-> Process ()
randomBinomialProcess_ :: Double -> Int -> Process ()
randomBinomialProcess_ Double
prob Int
trials =
do Int
t <- Parameter Int -> Process Int
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Int -> Process Int) -> Parameter Int -> Process Int
forall a b. (a -> b) -> a -> b
$ Double -> Int -> Parameter Int
randomBinomial Double
prob Int
trials
Double -> Process ()
holdProcess (Double -> Process ()) -> Double -> Process ()
forall a b. (a -> b) -> a -> b
$ Int -> Double
forall a b. (Integral a, Num b) => a -> b
fromIntegral Int
t
randomGammaProcess :: Double
-> Double
-> Process Double
randomGammaProcess :: Double -> Double -> Process Double
randomGammaProcess Double
kappa Double
theta =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Parameter Double
randomGamma Double
kappa Double
theta
Double -> Process ()
holdProcess Double
t
Double -> Process Double
forall a. a -> Process a
forall (m :: * -> *) a. Monad m => a -> m a
return Double
t
randomGammaProcess_ :: Double
-> Double
-> Process ()
randomGammaProcess_ :: Double -> Double -> Process ()
randomGammaProcess_ Double
kappa Double
theta =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Parameter Double
randomGamma Double
kappa Double
theta
Double -> Process ()
holdProcess Double
t
randomBetaProcess :: Double
-> Double
-> Process Double
randomBetaProcess :: Double -> Double -> Process Double
randomBetaProcess Double
alpha Double
beta =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Parameter Double
randomBeta Double
alpha Double
beta
Double -> Process ()
holdProcess Double
t
Double -> Process Double
forall a. a -> Process a
forall (m :: * -> *) a. Monad m => a -> m a
return Double
t
randomBetaProcess_ :: Double
-> Double
-> Process ()
randomBetaProcess_ :: Double -> Double -> Process ()
randomBetaProcess_ Double
alpha Double
beta =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Parameter Double
randomBeta Double
alpha Double
beta
Double -> Process ()
holdProcess Double
t
randomWeibullProcess :: Double
-> Double
-> Process Double
randomWeibullProcess :: Double -> Double -> Process Double
randomWeibullProcess Double
alpha Double
beta =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Parameter Double
randomWeibull Double
alpha Double
beta
Double -> Process ()
holdProcess Double
t
Double -> Process Double
forall a. a -> Process a
forall (m :: * -> *) a. Monad m => a -> m a
return Double
t
randomWeibullProcess_ :: Double
-> Double
-> Process ()
randomWeibullProcess_ :: Double -> Double -> Process ()
randomWeibullProcess_ Double
alpha Double
beta =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ Double -> Double -> Parameter Double
randomWeibull Double
alpha Double
beta
Double -> Process ()
holdProcess Double
t
randomDiscreteProcess :: DiscretePDF Double
-> Process Double
randomDiscreteProcess :: DiscretePDF Double -> Process Double
randomDiscreteProcess DiscretePDF Double
dpdf =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ DiscretePDF Double -> Parameter Double
forall a. DiscretePDF a -> Parameter a
randomDiscrete DiscretePDF Double
dpdf
Double -> Process ()
holdProcess Double
t
Double -> Process Double
forall a. a -> Process a
forall (m :: * -> *) a. Monad m => a -> m a
return Double
t
randomDiscreteProcess_ :: DiscretePDF Double
-> Process ()
randomDiscreteProcess_ :: DiscretePDF Double -> Process ()
randomDiscreteProcess_ DiscretePDF Double
dpdf =
do Double
t <- Parameter Double -> Process Double
forall a. Parameter a -> Process a
forall (m :: * -> *) a. ParameterLift m => Parameter a -> m a
liftParameter (Parameter Double -> Process Double)
-> Parameter Double -> Process Double
forall a b. (a -> b) -> a -> b
$ DiscretePDF Double -> Parameter Double
forall a. DiscretePDF a -> Parameter a
randomDiscrete DiscretePDF Double
dpdf
Double -> Process ()
holdProcess Double
t