Jump to content

The IBM model: Difference between revisions

From SELESwiki
Newer edit →
Fei (talk | contribs)
Administrators
860 edits
VisualWikitext
No edit summary
 
No edit summary
Line 25: Line 25:
   Producer.lse DEBUG
   Producer.lse DEBUG
   Consumer.lse DEBUG
   Consumer.lse DEBUG
  Update.lse DEBUG
 
  Spatial Variables:
  Spatial Variables:
   LSState[3]
   LSState[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 39:
   pDeathProducer = 0.01
   pDeathProducer = 0.01
   pDeathConsumer = 0.05
   pDeathConsumer = 0.05
//  SpreadRateProducer = 1
//  SpreadRateConsumer = 1
  Output Frequency: 1
  Output Frequency: 1


===consumer.lse===
===consumer.lse===
   // Process based definition for consumber behaviour in producer/consumer model
   // Process based definition for consumber behaviour in producer/consumer model
  LSEVENT: Consumder
  LSAGENT: Consumer
  DEFINITIONS
  DEFINITIONS
   LAYER: LSState, PrevLSState
   LAYER: LSState
   GLOBAL CONSTANT: Uncolonised, Producer, Consumer
   GLOBAL CONSTANT: Uncolonised, Producer, Consumer
   GLOBAL VARIABLE: pInitialConsumer,   pSpreadConsumer, pDeathConsumer //, SpreadRateConsumer
   GLOBAL VARIABLE: pInitialConsumer, pSpreadConsumer, pDeathConsumer //, SpreadRateConsumer
  AGENT VARIABLE: isSedentary
  ENDDEF
  ENDDEF
  // Just start once
  // Just immigrate once
  RETURNTIME = 0
  IMMIGRATIONTIME = 0
NUMAGENTS = NUMCELLS * pInitialConsumer
  // Initialize in an empty cell
  // Initialize in an empty cell
  PROBINIT = IF LSState EQ Uncolonised THEN pInitialConsumer ELSE 0
  PROBINIT
  PROBINIT = LSState EQ Uncolonised
  LSState = Consumer
  isSedentary = TRUE
  ENDPI
  TRANSITIONS
  TRANSITIONS
   // Decide if this cell survives
   // Decide if this cell survives
   Dies = UNIFORM(0,1) < pDeathConsumer
   Dies = UNIFORM(0,1) < pDeathConsumer
   // If the consumer is dying
   // If the consumer is dying
   IF Dies
   IF Dies
Line 62: Line 69:
   ENDFN
   ENDFN
   TRANSITIONS = !Dies
   TRANSITIONS = !Dies
  LSState = Consumer
  ENDTR
  ENDTR
  SPREADTIMESTEP = 1 //SpreadRateConsumer
  POPULATIONTIME = 1
  SPREADLOCATION
  MOVELOCATION
   REGION CENTRED(0,1, WRAPPED)
  d = IF isSedentary THEN 0 ELSE 1
  ENDSL
   REGION CENTRED(d, d, WRAPPED)
  SPREADPROB
  ENDML
   // Always spread to the same cell  (self), and spread probabilistical to neighbouring producer cells
  MOVEPROB
   SPREADPROB = IF Location EQ (SOURCE  Location) THEN 1
   // Parent stays in place and offspring must move to producer neighbour
              ELSE IF PrevLSState EQ  Producer THEN pSpreadConsumer
   MOVEPROB = IF isSedentary THEN 1
              ELSE 0
            ELSE LSState EQ Producer
  ENDSP
   isSedentary = TRUE
  LSState = Consumer
ENDMP
  NUMOFFSPRING
  n = 0
  OVER REGION CENTRED(1,1,WRAPPED)
      DECISION LSState EQ Uncolonised
      n = n + (UNIFORM(0,1) < pSpreadConsumer)
  ENDFN
 
  NUMOFFSPRING = n
  isSedentary = FALSE
  ENDNO


===producer.lse===
===producer.lse===
  // Process based definiton of producer behaviour in producer/consumer model
  // Agent based definiton of producer behaviour in producer/consumer model
  LSEVENT: Producer
  LSAGENT: Producer
  DEFINITIONS
  DEFINITIONS
   LAYER: LSState, PrevLSState
   LAYER: LSState
   GLOBAL CONSTANT: Uncolonised, Producer, Consumer
   GLOBAL CONSTANT: Uncolonised, Producer, Consumer
   GLOBAL VARIABLE: pInitialProducer, pSpreadProducer, pDeathProducer //,  SpreadRateProducer
   GLOBAL VARIABLE: pInitialProducer, pSpreadProducer, pDeathProducer
  AGENT VARIABLE: isSedentary
  ENDDEF
  ENDDEF
  // Just start once
  // Just immigrate once
  RETURNTIME = 0
  IMMIGRATIONTIME = 0
NUMAGENTS = NUMCELLS * pInitialProducer
  // Initialize in an empty cell
  // Initialize in an empty cell
  PROBINIT = IF LSState EQ Uncolonised THEN pInitialProducer ELSE 0
  PROBINIT
  PROBINIT = LSState EQ Uncolonised
  LSState = Producer
  isSedentary = TRUE
  ENDPI
  TRANSITIONS
  TRANSITIONS
   // Decide if this cell survives
   // Decide if this cell survives
   Dies = UNIFORM(0,1) < pDeathProducer
   Dies = (UNIFORM(0,1) < pDeathProducer) OR (LSState NEQ Producer)
   // If the producer is dying
   // If the producer is dying
   IF Dies
   IF Dies
Line 95: Line 119:
   ENDFN
   ENDFN
   TRANSITIONS = (!Dies)
   TRANSITIONS = (!Dies)
  LSState = Producer
  ENDTR
  ENDTR
  SPREADTIMESTEP = 1 //SpreadRateProducer
  POPULATIONTIME = 1
  SPREADLOCATION
  MOVELOCATION
   notEaten = PrevLSState EQ Producer //  check that the producer wasn't eaten  during the last step
   d = IF isSedentary THEN 0 ELSE 1
   REGION CENTRED(0, 1, WRAPPED)
   REGION CENTRED(d, d, WRAPPED)
      DECISION notEaten // don't spread if  eaten (i.e. stop process)
  ENDML
  ENDSL
  MOVEPROB
  SPREADPROB
   // Parent stays in place and offspring must move to unoccupied neighbour
   // Always spread to the same cell  (self), and spread probabilistically to empty neighbouring cells
   MOVEPROB = IF isSedentary THEN 1
   SPREADPROB = IF Location EQ (SOURCE  Location) THEN 1 // spread to same cell
            ELSE LSState EQ Uncolonised
              ELSE IF PrevLSState EQ Uncolonised THEN pSpreadProducer
  isSedentary = TRUE
              ELSE 0
  LSState = Producer
ENDSP
  ENDMP
 
  NUMOFFSPRING
===update.lse===
  n = 0
  // Event to update time step
   OVER REGION CENTRED(1,1,WRAPPED)
  LSEVENT: Update
      DECISION LSState EQ Uncolonised
DEFINITIONS
      n = n + (UNIFORM(0,1) < pSpreadProducer)
   LAYER: LSState, PrevLSState
  ENDFN
ENDDEF
   NUMOFFSPRING = n
RETURNTIME
   isSedentary = FALSE
   RETURNTIME = 1
  ENDNO
   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:

Revision as of 21:13, 15 November 2006

Summary

Screenshot

Download Model

Download the .scn, .sel and .lse files by clicking on the following link: File: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

Spatial Variables:
 LSState[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
Output Frequency: 1


consumer.lse

 // Process based definition for consumber behaviour in producer/consumer model
LSAGENT: Consumer
DEFINITIONS
  LAYER: LSState
  GLOBAL CONSTANT: Uncolonised, Producer, Consumer
  GLOBAL VARIABLE: pInitialConsumer, pSpreadConsumer, pDeathConsumer //, SpreadRateConsumer
  AGENT VARIABLE: isSedentary
ENDDEF
// Just immigrate once
IMMIGRATIONTIME = 0
NUMAGENTS = NUMCELLS * pInitialConsumer
// Initialize in an empty cell
PROBINIT
  PROBINIT = LSState EQ Uncolonised
  LSState = Consumer
  isSedentary = TRUE
ENDPI
TRANSITIONS
  // Decide if this cell survives
  Dies = UNIFORM(0,1) < pDeathConsumer

  // If the consumer is dying
  IF Dies
      LSState = Uncolonised
  ENDFN
  TRANSITIONS = !Dies
ENDTR
POPULATIONTIME = 1
MOVELOCATION
  d = IF isSedentary THEN 0 ELSE 1
  REGION CENTRED(d, d, WRAPPED)
ENDML
MOVEPROB
  // Parent stays in place and offspring must move to producer neighbour
 MOVEPROB = IF isSedentary THEN 1
            ELSE LSState EQ Producer
 isSedentary = TRUE
 LSState = Consumer
ENDMP
NUMOFFSPRING
  n = 0
  OVER REGION CENTRED(1,1,WRAPPED)
     DECISION LSState EQ Uncolonised
     n = n + (UNIFORM(0,1) < pSpreadConsumer)
  ENDFN

  NUMOFFSPRING = n
  isSedentary = FALSE
ENDNO

producer.lse

// Agent based definiton of producer behaviour in producer/consumer model
LSAGENT: Producer
DEFINITIONS
  LAYER: LSState
  GLOBAL CONSTANT: Uncolonised, Producer, Consumer
  GLOBAL VARIABLE: pInitialProducer, pSpreadProducer, pDeathProducer
  AGENT VARIABLE: isSedentary
ENDDEF
// Just immigrate once
IMMIGRATIONTIME = 0
NUMAGENTS = NUMCELLS * pInitialProducer 
// Initialize in an empty cell
PROBINIT
  PROBINIT = LSState EQ Uncolonised
  LSState = Producer
  isSedentary = TRUE
ENDPI
TRANSITIONS
  // Decide if this cell survives
  Dies = (UNIFORM(0,1) < pDeathProducer) OR (LSState NEQ Producer)
  // If the producer is dying
  IF Dies
     LSState = Uncolonised
  ENDFN
  TRANSITIONS = (!Dies)
ENDTR
POPULATIONTIME = 1
MOVELOCATION
  d = IF isSedentary THEN 0 ELSE 1
  REGION CENTRED(d, d, WRAPPED)
ENDML
MOVEPROB
  // Parent stays in place and offspring must move to unoccupied neighbour
 MOVEPROB = IF isSedentary THEN 1
            ELSE LSState EQ Uncolonised
 isSedentary = TRUE
 LSState = Producer
ENDMP
NUMOFFSPRING
  n = 0
  OVER REGION CENTRED(1,1,WRAPPED)
     DECISION LSState EQ Uncolonised
     n = n + (UNIFORM(0,1) < pSpreadProducer)
  ENDFN
  NUMOFFSPRING = n
  isSedentary = FALSE
ENDNO

Suggested Experiments

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

Retrieved from "https://seles.info/index.php?title=The_IBM_model&oldid=3184"