The CA model - Revision history

Fei: /* consumerCA.lse */

2006-11-15T21:35:11Z

consumerCA.lse

← Older revision		Revision as of 21:35, 15 November 2006
Line 38:		Line 38:
	Output Frequency: 1		Output Frequency: 1

	===~~consumerCA~~.lse===		===consumeCA.lse===
	// Producer/Consumer conceptual model implemented as a cellular automata		// Producer/Consumer conceptual model implemented as a cellular automata
	// updating everything in a single event		// updating everything in a single event

Fei at 21:10, 15 November 2006

2006-11-15T21:10:25Z

← Older revision		Revision as of 21:10, 15 November 2006
Line 86:		Line 86:
	TRANSITIONS = FALSE		TRANSITIONS = FALSE
	ENDTR		ENDTR





	==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: /* update.lse */

2006-11-15T21:10:05Z

update.lse

← Older revision		Revision as of 21:10, 15 November 2006
Line 89:		Line 89:


	~~===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-15T21:09:48Z

producer.lse

@@ Line 87: / Line 87: @@
   ENDTR
-===producer.lse===
 ===update.lse===

Fei: /* consumer.lse */

2006-11-15T21:08:56Z

consumer.lse

← Older revision		Revision as of 21:08, 15 November 2006
Line 38:		Line 38:
	Output Frequency: 1		Output Frequency: 1

	===~~consumer~~.lse===		===consumerCA.lse===
	// Producer/Consumer conceptual model implemented as a cellular automata		// Producer/Consumer conceptual model implemented as a cellular automata
	// updating everything in a single event		// updating everything in a single event

Fei: /* consumer.lse */

2006-11-15T21:07:15Z

consumer.lse

← Older revision		Revision as of 21:07, 15 November 2006
Line 39:		Line 39:

	===consumer.lse===		===consumer.lse===
	// ~~Process based definition for~~ ~~consumber behaviour~~ in ~~producer/consumer model~~		// Producer/Consumer conceptual model implemented as a cellular automata
	LSEVENT: ~~Consumder~~		// updating everything in a single event
			LSEVENT: ConsumeCA
	DEFINITIONS		DEFINITIONS
	LAYER: LSState, PrevLSState		LAYER: LSState, PrevLSState
	GLOBAL CONSTANT: Uncolonised, Producer, Consumer		GLOBAL CONSTANT: Uncolonised, Producer, Consumer
	GLOBAL VARIABLE: pInitialConsumer, pSpreadConsumer, pDeathConsumer ~~//, SpreadRateConsumer~~		GLOBAL VARIABLE: pInitialProducer, pInitialConsumer
			GLOBAL VARIABLE: pSpreadProducer, pDeathProducer, pSpreadConsumer, pDeathConsumer
	ENDDEF		ENDDEF
	~~// Just start once~~		INITIALSTATE
	RETURNTIME = 0		LSState = IF (UNIFORM(0,1) < pInitialConsumer) THEN Consumer
	// ~~Initialize in an empty~~ cell		ELSE IF (UNIFORM(0,1) < pInitialProducer) THEN Producer
	~~PROBINIT = IF LSState EQ Uncolonised THEN pInitialConsumer ELSE 0~~		ELSE Uncolonised
			INITIALSTATE = 1
			ENDIS
			RETURNTIME
			RETURNTIME = 1
			PrevLSState = LSState
			ENDRT
			// Cellular automata: visit every cell and change state based on current state
			// and neighbourhood
	TRANSITIONS		TRANSITIONS
	// ~~Decide~~ if ~~this cell survives~~		// An empty cell can possibly be colonised by a producer
	~~Dies~~ = UNIFORM(0,1) < ~~pDeathConsumer~~		// if it is next to at least one producer
	// ~~If the~~ consumer ~~is dying~~		IF PrevLSState EQ Uncolonised
	IF ~~Dies~~		OVER REGION CENTRED(1,1, WRAPPED)
	LSState = Uncolonised		DECISION (PrevLSState EQ Producer) AND (UNIFORM(0,1) < pSpreadProducer)
	ENDFN		SOURCE LSState = Producer
	TRANSITIONS = ~~!Dies~~		ENDFN
	~~LSState = Consumer~~		// A producer can possibly die, or be eaten by a consumer
			ELSE IF PrevLSState EQ Producer
			IF (UNIFORM(0,1) < pDeathProducer) // mortality
			LSState = Uncolonised
			ELSE
			OVER REGION CENTRED(1,1, WRAPPED)
			DECISION (PrevLSState EQ Consumer) AND (UNIFORM(0,1) < pSpreadConsumer)
			SOURCE LSState = Consumer
			ENDFN
			ENDFN
			// A consumer may die
			ELSE // PrevLSState EQ Consumer
			IF UNIFORM(0,1) < pDeathConsumer
			LSState = Uncolonised
			ENDFN
			ENDFN
			TRANSITIONS = FALSE
	ENDTR		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: /* Model.sel */

2006-11-15T21:06:12Z

Model.sel

← Older revision		Revision as of 21:06, 15 November 2006
Line 23:		Line 23:
	Time Units: Step kiloStep 1000 1000		Time Units: Step kiloStep 1000 1000
	Landscape Events:		Landscape Events:
	~~Producer.lse DEBUG~~		ConsumeCA.lse DEBUG
	~~Consumer.lse DEBUG~~
	~~Update~~.lse DEBUG
	Spatial Variables:		Spatial Variables:
	LSState[3]		LSState[3]
	PrevLSState[3]		PrevLSState[3]
	Legends:		Legends:
	LSStateTypes = {0:Uncolonised, 1:Producer, 2:Consumer}		LSStateTypes = {0:Uncolonised, 1:Producer, 2:Consumer}
	Global Variables:		Global Variables:
	pInitialProducer = 0.01		pInitialProducer = 0.01
Line 38:		Line 36:
	pDeathProducer = 0.01		pDeathProducer = 0.01
	pDeathConsumer = 0.05		pDeathConsumer = 0.05
	~~// SpreadRateProducer = 1~~
	~~// SpreadRateConsumer = 1~~
	Output Frequency: 1		Output Frequency: 1

Fei at 21:03, 15 November 2006

2006-11-15T21:03:52Z

New page

==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: