Merge remote-tracking branch 'origin/master'

Author: khvilaboa

epidemic_roman/README.md

@@ -0,0 +1,12 @@
+# Susceptible-Infected-Recovered Disease Spread Model
+
+Yet another Cell-DEVS implementation of the cellular automata model presented by [Hoya White et al.](https://www.sciencedirect.com/science/article/pii/S0096300306009295).
+However, rules have been re-interpreted by [Román Cárdenas](mailto:[email protected]).
+
+
+## How to run the model:
+
+In the linux terminal, type:  
+```shell
+> ./../cd++ -msir_roman.ma -lresults/sir_roman.log
+```

epidemic_roman/sir_roman.inc

@@ -0,0 +1,19 @@
+#BeginMacro(recovery)
+0.4
+#EndMacro
+
+
+#BeginMacro(i_effect)
+min
+(
+  (0,0)~sus,
+  (0,0)~sus * $virulence *
+  (
+    (0,0)~infec + $connection * $movement / (0,0)~pop *
+    (
+      (1,0)~pop * (1,0)~infec + (-1,0)~pop * (-1,0)~infec +
+      (0,1)~pop * (0,1)~infec + (0,-1)~pop * (0,-1)~infec
+    )
+  )
+)
+#EndMacro

epidemic_roman/sir_roman.ma

@@ -0,0 +1,35 @@
+#include(sir_roman.inc)
+[top]
+components : area
+
+[area]
+type : cell
+width : 50
+height : 50
+delay : transport
+defaultDelayTime : 1
+border : wrapped
+neighbors : 			area(-1,0)
+neighbors : area(0,-1)  area(0,0)  area(0,1)
+neighbors : 			area(1,0)
+initialvalue : -1
+
+
+localtransition : sir-roman-rule
+
+statevariables: population virulence connection movement i_sus i_infec i_rec
+statevalues: 100 0.6 1.0 0.5 1.0 0.0 0.0
+initialvariablesvalue: sir_roman.var
+
+neighborports : initial infec rec pop sus
+
+[sir-roman-rule]
+rule : {~initial := 0; ~pop := $population; ~infec := $i_infec; ~sus := $i_sus; ~rec := $i_rec; }
+ 1 {(0,0)~initial = -1}
+
+rule : {~pop := $population; ~infec := $i_infec; ~sus := $i_sus; ~rec := $i_rec; }
+{
+$i_sus := round(((0,0)~sus - #macro(i_effect))*100)/100;
+$i_infec:= round(((1 - #macro(recovery)) * (0,0)~infec + (0,0)~sus - $i_sus)*100)/100;
+$i_rec := 1 - $i_sus - $i_infec;}
+1 { (0,0)~initial != -1 }

epidemic_roman/sir_roman.pal

@@ -0,0 +1,7 @@
+[-100; -0.001] 255 0 0
+[0; 0.200] 0 0 0
+[0.201; 0.400] 50 50 50
+[0.401; 0.600] 100 100 100
+[0.601; 0.800] 150 150 150
+[0.801; 1.00] 200 200 200
+[1.001;100] 255 0 0

epidemic_roman/sir_roman.var

@@ -0,0 +1 @@
+(25,25) = 10 0.6 1 1 0.3 0.7 0.0