aivika-transformers-4.3.2: Transformers for the Aivika simulation library

CopyrightCopyright (c) 2009-2015, David Sorokin <david.sorokin@gmail.com>
LicenseBSD3
MaintainerDavid Sorokin <david.sorokin@gmail.com>
Stabilityexperimental
Safe HaskellNone
LanguageHaskell2010

Simulation.Aivika.Trans.Circuit

Contents

Description

Tested with: GHC 7.10.1

It represents a circuit synchronized with the event queue. Also it allows creating the recursive links with help of the proc-notation.

The implementation is based on the Arrow Tutorial.

Synopsis

The Circuit Arrow

newtype Circuit m a b Source

Represents a circuit synchronized with the event queue. Besides, it allows creating the recursive links with help of the proc-notation.

Constructors

Circuit 

Fields

runCircuit :: a -> Event m (b, Circuit m a b)

Run the circuit.

iterateCircuitInIntegTimes :: MonadDES m => Circuit m a a -> a -> Event m (Task m a) Source

Iterate the circuit in the integration time points returning a task which completes after the final output of the circuit is received.

iterateCircuitInIntegTimes_ :: MonadDES m => Circuit m a a -> a -> Event m () Source

Iterate the circuit in the integration time points.

iterateCircuitInIntegTimesMaybe :: MonadDES m => Circuit m a (Maybe a) -> a -> Event m () Source

Iterate the circuit in the integration time points, interrupting the iteration immediately if Nothing is returned within the Circuit computation.

iterateCircuitInIntegTimesEither :: MonadDES m => Circuit m a (Either b a) -> a -> Event m (Task m (Either b a)) Source

Iterate the circuit in the integration time points returning a task that computes the final output of the circuit either after all points are exhausted, or after the Left result of type b is received, which interrupts the computation immediately.

iterateCircuitInTimes :: MonadDES m => [Double] -> Circuit m a a -> a -> Event m (Task m a) Source

Iterate the circuit in the specified time points returning a task which completes after the final output of the circuit is received.

iterateCircuitInTimes_ :: MonadDES m => [Double] -> Circuit m a a -> a -> Event m () Source

Iterate the circuit in the specified time points.

iterateCircuitInTimesMaybe :: MonadDES m => [Double] -> Circuit m a (Maybe a) -> a -> Event m () Source

Iterate the circuit in the specified time points, interrupting the iteration immediately if Nothing is returned within the Circuit computation.

iterateCircuitInTimesEither :: MonadDES m => [Double] -> Circuit m a (Either b a) -> a -> Event m (Task m (Either b a)) Source

Iterate the circuit in the specified time points returning a task that computes the final output of the circuit either after all points are exhausted, or after the Left result of type b is received, which interrupts the computation immediately.

Circuit Primitives

arrCircuit :: MonadDES m => (a -> Event m b) -> Circuit m a b Source

Create a simple circuit by the specified handling function that runs the computation for each input value to get an output.

accumCircuit :: MonadDES m => (acc -> a -> Event m (acc, b)) -> acc -> Circuit m a b Source

Accumulator that outputs a value determined by the supplied function.

The Arrival Circuit

arrivalCircuit :: MonadDES m => Circuit m a (Arrival a) Source

A circuit that adds the information about the time points at which the values were received.

Delaying the Circuit

delayCircuit :: MonadDES m => a -> Circuit m a a Source

Delay the input by one step using the specified initial value.

The Time Circuit

timeCircuit :: MonadDES m => Circuit m a Double Source

A circuit that returns the current modeling time.

Conditional Computation

(<?<) Source

Arguments

:: MonadDES m 
=> Circuit m b c

process the event if it presents

-> Circuit m a (Maybe b)

whether there is an event

-> Circuit m a (Maybe c)

the resulting circuit that processes only the represented events

Like <<< but processes only the represented events.

(>?>) Source

Arguments

:: MonadDES m 
=> Circuit m a (Maybe b)

whether there is an event

-> Circuit m b c

process the event if it presents

-> Circuit m a (Maybe c)

the resulting circuit that processes only the represented events

Like >>> but processes only the represented events.

filterCircuit :: MonadDES m => (a -> Bool) -> Circuit m a b -> Circuit m a (Maybe b) Source

Filter the circuit, calculating only those parts of the circuit that satisfy the specified predicate.

filterCircuitM :: MonadDES m => (a -> Event m Bool) -> Circuit m a b -> Circuit m a (Maybe b) Source

Filter the circuit within the Event computation, calculating only those parts of the circuit that satisfy the specified predicate.

neverCircuit :: MonadDES m => Circuit m a (Maybe b) Source

The source of events that never occur.

Converting to Signals and Processors

circuitSignaling :: MonadDES m => Circuit m a b -> Signal m a -> Signal m b Source

Get a signal transform by the specified circuit.

circuitProcessor :: MonadDES m => Circuit m a b -> Processor m a b Source

Transform the circuit to a processor.

Integrals and Difference Equations

integCircuit Source

Arguments

:: MonadDES m 
=> Double

the initial value

-> Circuit m Double Double

map the derivative to an integral

An approximation of the integral using Euler's method.

This function can be rather inaccurate as it depends on the time points at wich the Circuit computation is actuated. Also Euler's method per se is not most accurate, although simple enough for implementation.

Consider using the integ function whenever possible. That function can integrate with help of the Runge-Kutta method by the specified integration time points that are passed in the simulation specs to every Simulation, when running the model.

At the same time, the integCircuit function has no mutable state unlike the former. The latter consumes less memory but at the cost of inaccuracy and relatively more slow simulation, had we requested the integral in the same time points.

Regarding the recursive equations, the both functions allow defining them but whithin different computations (either with help of the recursive do-notation or the proc-notation).

integCircuitEither Source

Arguments

:: MonadDES m 
=> Double

the initial value

-> Circuit m (Either Double Double) Double

map either a new Left value or the Right derivative to an integral

Like integCircuit but allows either setting a new Left integral value, or using the Right derivative when integrating by Euler's method.

sumCircuit Source

Arguments

:: (MonadDES m, Num a) 
=> a

the initial value

-> Circuit m a a

map the difference to a sum

A sum of differences starting from the specified initial value.

Consider using the more accurate diffsum function whener possible as it is calculated in every integration time point specified by specs passed in to every Simulation, when running the model.

At the same time, the sumCircuit function has no mutable state and it consumes less memory than the former.

Regarding the recursive equations, the both functions allow defining them but whithin different computations (either with help of the recursive do-notation or the proc-notation).

sumCircuitEither Source

Arguments

:: (MonadDES m, Num a) 
=> a

the initial value

-> Circuit m (Either a a) a

map either a new Left value or the Right difference to a sum

Like sumCircuit but allows either setting a new Left value for the sum, or updating it by specifying the Right difference.

The Circuit Transform

circuitTransform :: (MonadSD m, MonadDES m) => Circuit m a b -> Transform m a b Source

Approximate the circuit as a transform of time varying function, calculating the values in the integration time points and then interpolating in all other time points. The resulting transform computation is synchronized with the event queue.

This procedure consumes memory as the underlying memoization allocates an array to store the calculated values.

Debugging

traceCircuit Source

Arguments

:: MonadDES m 
=> Maybe String

the request message

-> Maybe String

the response message

-> Circuit m a b

a circuit

-> Circuit m a b 

Show the debug messages with the current simulation time.