loom

1
// Encapsulating the model file into the LorenzSystem type.
2
// This way you can create many objects from one model.
3
import "diff/LorenzSystem.diff" type LorenzSystem
4

5
// Model object initialization
6
LorenzSystem ls = {
7
    sigma : 10.0,
8
    b     : 8.0/3,
9
    R     : 28.0,
10
    x : 10.0,
11
    y : 1.0,
12
    z : 1.0,
13
}
14

15
// Charting module
16
import "Chart" chart = {title : "Lorenz System"}
17

18
// Delayed execution
19
finally chart.show()
20

21
// Simulation execution module
22
import "Scene" scene = {
23
    // Initializing simulation data
24
    t  : 0.0,
25
    tk : 50.0,
26
    dt : 0.01,
27
    // Lambda with a closure as a callback function
28
    callback : func [chart : var, ls] ()
29
    {
30
        call chart.addPoint("x-y projection", ls.x, ls.y)
31
        call chart.addPoint("x-z projection", ls.x, ls.z)
32
    },
33
    callback_interval : 100,
34
    // Populating an array of scene objects
35
    actors : [ls],
36
}
37

38
// Run simulation
39
run scene
40