Producer Consumer - Revision history

Fei at 21:02, 15 November 2006

2006-11-15T21:02:51Z

← Older revision		Revision as of 21:02, 15 November 2006
Line 1:		Line 1:
	~~==Summary==~~


	~~==Screenshot==~~

	~~[[Image:Consumer.gif]]~~

	~~==Download Model==~~

	~~Download the .scn, .sel and .lse files by clicking on the following link: [[Image:ProducerConsumer.zip]]~~

	~~==Model Code Exploration==~~

	In the following sections we will examine all of the model files for this model. Note that instead of downloading the zip file above, you could just copy the text in the boxes below into a text editor and save it with the appropriate name (Section title). Opening the resulting .scn file in the SELES simulator would run this model.

	~~===Consume.scn===~~
	~~Seles Scenario~~
	~~Model.sel~~
	~~SimPriority Low Priority~~

	~~===Model.sel===~~
	~~Seles Model~~
	~~Time Units: Step kiloStep 1000 1000~~
	~~Landscape Events:~~
	~~Producer.lse DEBUG~~
	~~Consumer.lse DEBUG~~
	~~Update.lse DEBUG~~
	~~Spatial Variables:~~
	~~LSState[3]~~
	~~PrevLSState[3]~~
	~~Legends:~~
	~~LSStateTypes = {0:Uncolonised, 1:Producer, 2:Consumer}~~
	~~Global Variables:~~
	~~pInitialProducer = 0.01~~
	~~pInitialConsumer = 0.01~~
	~~pSpreadProducer = 0.1~~
	~~pSpreadConsumer = 0.2~~
	~~pDeathProducer = 0.01~~
	~~pDeathConsumer = 0.05~~
	~~// SpreadRateProducer = 1~~
	~~// SpreadRateConsumer = 1~~
	~~Output Frequency: 1~~

	~~===consumer.lse===~~
	~~// Process based definition for consumber behaviour in producer/consumer model~~
	~~LSEVENT: Consumder~~
	~~DEFINITIONS~~
	~~LAYER: LSState, PrevLSState~~
	~~GLOBAL CONSTANT: Uncolonised, Producer, Consumer~~
	~~GLOBAL VARIABLE: pInitialConsumer, pSpreadConsumer, pDeathConsumer //, SpreadRateConsumer~~
	~~ENDDEF~~
	~~// Just start once~~
	~~RETURNTIME = 0~~
	~~// Initialize in an empty cell~~
	~~PROBINIT = IF LSState EQ Uncolonised THEN pInitialConsumer ELSE 0~~
	~~TRANSITIONS~~
	~~// Decide if this cell survives~~
	~~Dies = UNIFORM(0,1) < pDeathConsumer~~
	~~// If the consumer is dying~~
	~~IF Dies~~
	~~LSState = Uncolonised~~
	~~ENDFN~~
	~~TRANSITIONS = !Dies~~
	~~LSState = Consumer~~
	~~ENDTR~~
	~~SPREADTIMESTEP = 1 //SpreadRateConsumer~~
	~~SPREADLOCATION~~
	~~REGION CENTRED(0,1, WRAPPED)~~
	~~ENDSL~~
	~~SPREADPROB~~
	~~// Always spread to the same cell (self), and spread probabilistical to neighbouring producer cells~~
	~~SPREADPROB = IF Location EQ (SOURCE Location) THEN 1~~
	~~ELSE IF PrevLSState EQ Producer THEN pSpreadConsumer~~
	~~ELSE 0~~
	~~ENDSP~~

	~~===producer.lse===~~
	~~// Process based definiton of producer behaviour in producer/consumer model~~
	~~LSEVENT: Producer~~
	~~DEFINITIONS~~
	~~LAYER: LSState, PrevLSState~~
	~~GLOBAL CONSTANT: Uncolonised, Producer, Consumer~~
	~~GLOBAL VARIABLE: pInitialProducer, pSpreadProducer, pDeathProducer //, SpreadRateProducer~~
	~~ENDDEF~~
	~~// Just start once~~
	~~RETURNTIME = 0~~
	~~// Initialize in an empty cell~~
	~~PROBINIT = IF LSState EQ Uncolonised THEN pInitialProducer ELSE 0~~
	~~TRANSITIONS~~
	~~// Decide if this cell survives~~
	~~Dies = UNIFORM(0,1) < pDeathProducer~~
	~~// If the producer is dying~~
	~~IF Dies~~
	~~LSState = Uncolonised~~
	~~ENDFN~~
	~~TRANSITIONS = (!Dies)~~
	~~LSState = Producer~~
	~~ENDTR~~
	~~SPREADTIMESTEP = 1 //SpreadRateProducer~~
	~~SPREADLOCATION~~
	~~notEaten = PrevLSState EQ Producer // check that the producer wasn't eaten during the last step~~
	~~REGION CENTRED(0, 1, WRAPPED)~~
	~~DECISION notEaten // don't spread if eaten (i.e. stop process)~~
	~~ENDSL~~
	~~SPREADPROB~~
	~~// Always spread to the same cell (self), and spread probabilistically to empty neighbouring cells~~
	~~SPREADPROB = IF Location EQ (SOURCE Location) THEN 1 // spread to same cell~~
	~~ELSE IF PrevLSState EQ Uncolonised THEN pSpreadProducer~~
	~~ELSE 0~~
	~~ENDSP~~

	~~===update.lse===~~
	~~// Event to update time step~~
	~~LSEVENT: Update~~
	~~DEFINITIONS~~
	~~LAYER: LSState, PrevLSState~~
	~~ENDDEF~~
	~~RETURNTIME~~
	~~RETURNTIME = 1~~
	~~PrevLSState = LSState~~
	~~ENDRT~~
	~~NUMCLUSTERS = 0~~

	~~==Suggested Experiments==~~

	~~To explore this cellular automata model further, try the following:~~

Fei: /* update.lse */

2006-11-08T23:19:20Z

update.lse

← Older revision		Revision as of 23:19, 8 November 2006
Line 111:		Line 111:

	===update.lse===		===update.lse===
			// Event to update time step
			LSEVENT: Update
			DEFINITIONS
			LAYER: LSState, PrevLSState
			ENDDEF
			RETURNTIME
			RETURNTIME = 1
			PrevLSState = LSState
			ENDRT
			NUMCLUSTERS = 0

	==Suggested Experiments==		==Suggested Experiments==

	To explore this cellular automata model further, try the following:		To explore this cellular automata model further, try the following:

Fei: /* producer.lse */

2006-11-08T23:18:49Z

producer.lse

← Older revision		Revision as of 23:18, 8 November 2006
Line 76:		Line 76:

	===producer.lse===		===producer.lse===
			// Process based definiton of producer behaviour in producer/consumer model
			LSEVENT: Producer
			DEFINITIONS
			LAYER: LSState, PrevLSState
			GLOBAL CONSTANT: Uncolonised, Producer, Consumer
			GLOBAL VARIABLE: pInitialProducer, pSpreadProducer, pDeathProducer //, SpreadRateProducer
			ENDDEF
			// Just start once
			RETURNTIME = 0
			// Initialize in an empty cell
			PROBINIT = IF LSState EQ Uncolonised THEN pInitialProducer ELSE 0
			TRANSITIONS
			// Decide if this cell survives
			Dies = UNIFORM(0,1) < pDeathProducer
			// If the producer is dying
			IF Dies
			LSState = Uncolonised
			ENDFN
			TRANSITIONS = (!Dies)
			LSState = Producer
			ENDTR
			SPREADTIMESTEP = 1 //SpreadRateProducer
			SPREADLOCATION
			notEaten = PrevLSState EQ Producer // check that the producer wasn't eaten during the last step
			REGION CENTRED(0, 1, WRAPPED)
			DECISION notEaten // don't spread if eaten (i.e. stop process)
			ENDSL
			SPREADPROB
			// Always spread to the same cell (self), and spread probabilistically to empty neighbouring cells
			SPREADPROB = IF Location EQ (SOURCE Location) THEN 1 // spread to same cell
			ELSE IF PrevLSState EQ Uncolonised THEN pSpreadProducer
			ELSE 0
			ENDSP

	===update.lse===		===update.lse===

Fei: /* consumer.lse */

2006-11-08T23:17:27Z

consumer.lse

← Older revision		Revision as of 23:17, 8 November 2006
Line 43:		Line 43:

	===consumer.lse===		===consumer.lse===
			// Process based definition for consumber behaviour in producer/consumer model
			LSEVENT: Consumder
			DEFINITIONS
			LAYER: LSState, PrevLSState
			GLOBAL CONSTANT: Uncolonised, Producer, Consumer
			GLOBAL VARIABLE: pInitialConsumer, pSpreadConsumer, pDeathConsumer //, SpreadRateConsumer
			ENDDEF
			// Just start once
			RETURNTIME = 0
			// Initialize in an empty cell
			PROBINIT = IF LSState EQ Uncolonised THEN pInitialConsumer ELSE 0
			TRANSITIONS
			// Decide if this cell survives
			Dies = UNIFORM(0,1) < pDeathConsumer
			// If the consumer is dying
			IF Dies
			LSState = Uncolonised
			ENDFN
			TRANSITIONS = !Dies
			LSState = Consumer
			ENDTR
			SPREADTIMESTEP = 1 //SpreadRateConsumer
			SPREADLOCATION
			REGION CENTRED(0,1, WRAPPED)
			ENDSL
			SPREADPROB
			// Always spread to the same cell (self), and spread probabilistical to neighbouring producer cells
			SPREADPROB = IF Location EQ (SOURCE Location) THEN 1
			ELSE IF PrevLSState EQ Producer THEN pSpreadConsumer
			ELSE 0
			ENDSP

	===producer.lse===		===producer.lse===

Fei at 23:15, 8 November 2006

2006-11-08T23:15:55Z

@@ Line 4: / Line 4: @@
 ==Screenshot==
-[[Image:ConvexHull.gif]]
+[[Image:Consumer.gif]]
 ==Download Model==
-Download the .scn, .sel and .lse files by clicking on the following link: [[Image:ConvexHull.zip]]
+Download the .scn, .sel and .lse files by clicking on the following link: [[Image:ProducerConsumer.zip]]
 ==Model Code Exploration==
@@ Line 14: / Line 14: @@
 In the following sections we will examine all of the model files for this model. Note that instead of downloading the zip file above, you could just copy the text in the boxes below into a text editor and save it with the appropriate name (Section title). Opening the resulting .scn file in the SELES simulator would run this model.
-===MLM.scn===
+===Consume.scn===
+  Seles Scenario
   Model.sel
   SimPriority Low Priority
 ===Model.sel===
   Seles Model
-  Time Units: Step kiloStep 1000 100000
+  Time Units: Step kiloStep 1000 1000
   Landscape Events:
- ConvexHull.lse
+  Producer.lse DEBUG
- Spatial Constants:
+  Consumer.lse DEBUG
- StudyArea
+  Update.lse DEBUG
   Spatial Variables:
- ConvexHullPoints[1] <= 0
+  LSState[3]
- ConvexHull[1] <= 0
+  PrevLSState[3]
   Global Variables:
-  TRLoc = 0
+  pInitialProducer = 0.01
   Output Frequency: 1
-===ConvexHull.lse===
+===consumer.lse===
- // Sweep-line approach to identifying convex hull using a worker agent
+===update.lse===
 ==Suggested Experiments==
 To explore this cellular automata model further, try the following:

Fei at 22:57, 8 November 2006

2006-11-08T22:57:11Z

New page

==Summary==

==Screenshot==

[[Image:ConvexHull.gif]]

==Download Model==

Download the .scn, .sel and .lse files by clicking on the following link: [[Image:ConvexHull.zip]]

==Model Code Exploration==

In the following sections we will examine all of the model files for this model. Note that instead of downloading the zip file above, you could just copy the text in the boxes below into a text editor and save it with the appropriate name (Section title). Opening the resulting .scn file in the SELES simulator would run this model.

===MLM.scn===

SELES Scenario
$gisData$ = ..\..\CaseStudy_Morice\gisData\cell
StudyArea = $gisData$\BEC
Model Dimensions: StudyArea
Model.sel
SimPriority Low Priority
//SimStart 100000 1

===Model.sel===

Seles Model
Time Units: Step kiloStep 1000 100000
Landscape Events:
ConvexHull.lse
Spatial Constants:
StudyArea
Spatial Variables:
ConvexHullPoints[1] <= 0
ConvexHull[1] <= 0
Global Variables:
TRLoc = 0
Output Frequency: 1

===ConvexHull.lse===

// Sweep-line approach to identifying convex hull using a worker agent
LSAGENT: ConvexHull
DEFINITIONS
LAYER: StudyArea, ConvexHull, ConvexHullPoints
GLOBAL VARIABLE: TRLoc
LOCAL VARIABLE: EndingPivotLoc
AGENT VARIABLE: PivotLocation, EdgeLocation, EdgeRow, EdgeCol
ENDDEF
INITIALSTATE
INITIALSTATE = 1
// Find the top-right-most point
TRLoc = 0
OVER REGION WHOLE MAP
DECISION StudyArea > 0
TRLoc = MAX(TRLoc, Location)
ENDFN
ConvexHull = 1 // Initialize convex hull layer to all ones (will erase areas outside hull)
ENDIS
NUMAGENTS = 1 // Only need a single worker
// Start on the first pivot vertex on convex hull
AGENTLOCATION
REGION LOCATION(TRLoc)
ENDAL
PROBINIT
PROBINIT = 1
// Set up first pivot vertex, and first raster boundary location
PivotLocation = Location
EdgeRow = NUMROWS-1
EdgeCol = COL(Location)
EdgeLocation = LOCATION(EdgeRow, EdgeCol)
EndingPivotLoc = -1
ENDPI
TRANSITIONS
// Continue until we reach the ending pivot location (which is actually second time we hit the second pivot vertex,
// since the first pivot vertex gets only partially processed on first pass)
TRANSITIONS = (PivotLocation NEQ EndingPivotLoc)
ConvexHullPoints = 1
// Walk line from boundary to vertex and erase until a new pivot vertex is hit (then move there)
oldPivot = PivotLocation
hitPivot = FALSE
OVER REGION VECTOR(EdgeLocation, PivotLocation)
DECISION !hitPivot
IF (StudyArea > 0) // hit new vertex
hitPivot = TRUE
IF (EndingPivotLoc EQ -1) AND (Location NEQ PivotLocation)
EndingPivotLoc = PivotLocation
ENDFN
PivotLocation = Location
ELSE
ConvexHull = 0 // erase
ENDFN
ENDFN
// Move boundary location clockwise around outside edge of raster
IF (EdgeRow EQ NUMROWS-1) AND (EdgeCol < NUMCOLS-1)
EdgeCol = EdgeCol + 1
ELSE IF (EdgeCol EQ NUMCOLS-1) AND (EdgeRow > 0)
EdgeRow = EdgeRow - 1
ELSE IF (EdgeRow EQ 0) AND (EdgeCol > 0)
EdgeCol = EdgeCol - 1
ELSE
EdgeRow = EdgeRow + 1
ENDFN
EdgeLocation = LOCATION(EdgeRow, EdgeCol)
ENDTR
POPULATIONTIME = 1
// Move to current/new pivot location
MOVELOCATION
REGION LOCATION(PivotLocation)
ENDML

==Suggested Experiments==

To explore this cellular automata model further, try the following: