IncrementalModel_step6 - gama-platform/gama GitHub Wiki
This step illustrates how to define a multi-level model.
We propose to let the buildings manage what happens when the people are inside buildings. In this context, we will use the multi-level properties of GAMA: when a people agent will be inside a building, it will be captured by it and its species will be modified. It will be not anymore the people agent that will decide when to leave the building, but the building itself that will release it.
We will need to:
- Define a micro-species of people inside the building species (
people_in_building
). - Define two new behaviors for building:
let_people_leave
andlet_people_enter
. - Modify the aspect of the building.
- Modify some global variables for counting the number of infected people.
First, we define a new species called people_in_building
inside the building
species. Thus, a building could have agents of this species as members
and control them. The people_in_building
species has for parent the people
species, which means that a people_in_building
agent has all the attributes, aspects and behaviors of a people
agent.
In our case, we want a people agent inside a building does not do anything. Thus, we use the schedules
facet of the species to remove the people_in_building
from the scheduler.
species building {
...
species people_in_building parent: people schedules: [] {
}
...
}
We define a first reflex
for the buildings that will be activated at each simulation step and that will allow the building to capture all the people that are inside its geometry and that are not moving (target = nil
). Capturing agents means putting them inside its members
list and changing their species: here the people
agents become people_in_building
agents.
species building {
...
reflex let_people_enter {
capture (people inside self where (each.target = nil)) as: people_in_building;
}
....
}
We define a second reflex for the buildings that will be activated at each simulation step and that will allow the building to release some of the people_in_building
agents. First, it increments the staying counter of all the people_in_building
agents. Then it builds the list of leaving people by testing the same probability as before for all the people_in_building
agents. Finally, if this list is not empty, it releases them as people
agents (and gives them a new target point).
species building {
...
reflex let_people_leave {
ask people_in_building {
staying_counter <- staying_counter + 1;
}
release people_in_building where (flip(each.staying_counter / staying_coeff)) as: people in: world {
target <- any_location_in(one_of(building));
}
}
....
}
At last, we refine the aspect of the buildings: if there are no people inside the building, we draw it with gray color. If the number of people_in_building
infected is higher than the number of people_in_building
not infected, we draw it in red; otherwise in green. The number of infected people_in_building
and its total number will be computed once a step (through the update
facet of building
attribute).
species building {
int nb_infected <- 0 update: self.people_in_building count each.is_infected;
int nb_total <- 0 update: length(self.people_in_building);
aspect default {
draw shape color: nb_total = 0 ? #gray : (float(nb_infected) / nb_total > 0.5 ? #red : #green) border: #black depth: height;
}
}
In order to take into account the people that are inside the buildings for the computation of nb_people_infected
, we first build the list of people_in_building
. As people_in_building
is microspecies of building
, we cannot compute it directly like for the other species, we then aggregate all the list people_in_building
of all building in a single list (list_people_in_buildings
). Then, we compute the number of infected people as the number of people infected outside the building + the number of people infected inside them.
global {
...
list<people_in_building> list_people_in_buildings update: (building accumulate each.people_in_building);
int nb_people_infected <- nb_infected_init update: (people + list_people_in_buildings) count (each.is_infected);
...
}
https://github.com/gama-platform/gama.old/blob/GAMA_1.9.2/msi.gama.models/models/Tutorials/Incremental%20Model/models/Incremental%20Model%206.gaml