AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM > Class Template Reference

#include <AbstractCellBasedSimulation.hpp>

Inherits Identifiable.

Inherited by OffLatticeSimulation< ELEMENT_DIM, SPACE_DIM >.

Collaboration diagram for AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >:
Collaboration graph
[legend]

List of all members.

Public Member Functions

 AbstractCellBasedSimulation (AbstractCellPopulation< ELEMENT_DIM, SPACE_DIM > &rCellPopulation, bool deleteCellPopulationInDestructor=false, bool initialiseCells=true)
virtual ~AbstractCellBasedSimulation ()
void SetCellBasedPdeHandler (CellBasedPdeHandler< SPACE_DIM > *pCellBasedPdeHandler)
CellBasedPdeHandler< SPACE_DIM > * GetCellBasedPdeHandler ()
std::vector< doubleGetNodeLocation (const unsigned &rNodeIndex)
double GetDt ()
unsigned GetNumBirths ()
unsigned GetNumDeaths ()
std::string GetOutputDirectory ()
void SetDt (double dt)
void SetEndTime (double endTime)
void SetOutputDirectory (std::string outputDirectory)
void SetSamplingTimestepMultiple (unsigned samplingTimestepMultiple)
void SetNoBirth (bool noBirth)
void SetUpdateCellPopulationRule (bool updateCellPopulation)
bool GetUpdateCellPopulationRule ()
void AddCellKiller (boost::shared_ptr< AbstractCellKiller< SPACE_DIM > > pCellKiller)
void RemoveAllCellKillers ()
void AddSimulationModifier (boost::shared_ptr< AbstractCellBasedSimulationModifier< ELEMENT_DIM, SPACE_DIM > > pSimulationModifier)
std::vector< boost::shared_ptr
< AbstractCellBasedSimulationModifier
< ELEMENT_DIM, SPACE_DIM > > > * 
GetSimulationModifiers ()
void Solve ()
AbstractCellPopulation
< ELEMENT_DIM, SPACE_DIM > & 
rGetCellPopulation ()
const AbstractCellPopulation
< ELEMENT_DIM, SPACE_DIM > & 
rGetCellPopulation () const
bool GetOutputDivisionLocations ()
void SetOutputDivisionLocations (bool outputDivisionLocations)
bool GetOutputCellVelocities ()
void SetOutputCellVelocities (bool outputCellVelocities)
virtual void OutputSimulationParameters (out_stream &rParamsFile)=0

Protected Member Functions

virtual void WriteVisualizerSetupFile ()
virtual unsigned DoCellBirth ()
virtual c_vector< double,
SPACE_DIM > 
CalculateCellDivisionVector (CellPtr pParentCell)=0
unsigned DoCellRemoval ()
virtual void SetupSolve ()
virtual bool StoppingEventHasOccurred ()
virtual void UpdateCellPopulation ()
virtual void UpdateCellLocationsAndTopology ()=0
void OutputSimulationSetup ()
virtual void OutputAdditionalSimulationSetup (out_stream &rParamsFile)

Protected Attributes

double mDt
double mEndTime
AbstractCellPopulation
< ELEMENT_DIM, SPACE_DIM > & 
mrCellPopulation
bool mDeleteCellPopulationInDestructor
bool mInitialiseCells
bool mNoBirth
bool mUpdateCellPopulation
std::string mOutputDirectory
std::string mSimulationOutputDirectory
out_stream mpVizSetupFile
unsigned mNumBirths
unsigned mNumDeaths
bool mOutputDivisionLocations
out_stream mpDivisionLocationFile
bool mOutputCellVelocities
out_stream mpCellVelocitiesFile
std::vector< boost::shared_ptr
< AbstractCellKiller
< SPACE_DIM > > > 
mCellKillers
std::vector< boost::shared_ptr
< AbstractCellBasedSimulationModifier
< ELEMENT_DIM, SPACE_DIM > > > 
mSimulationModifiers
unsigned mSamplingTimestepMultiple
CellBasedPdeHandler< SPACE_DIM > * mpCellBasedPdeHandler

Private Member Functions

template<class Archive >
void serialize (Archive &archive, const unsigned int version)

Friends

class TestCryptSimulation2d
class TestOffLatticeSimulation3d
class TestOffLatticeSimulation
class boost::serialization::access

Detailed Description

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
class AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >

An abstract cell-based simulation class. This class contains common functionality from off lattice and on lattice simulations.

The AbstractCellBasedSimulation is constructed with a CellPopulation, which updates the correspondence between each Cell and its spatial representation and handles cell division (governed by the CellCycleModel associated with each cell). Once constructed, one or more CellKillers may be passed to the AbstractCellBasedSimulation object to specify conditions in which Cells may die,

Subclasses use one or more Force laws or update rules (Which are passed to the child class object) to define the mechanical properties of the CellPopulation.

Definition at line 66 of file AbstractCellBasedSimulation.hpp.


Constructor & Destructor Documentation

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::AbstractCellBasedSimulation ( AbstractCellPopulation< ELEMENT_DIM, SPACE_DIM > &  rCellPopulation,
bool  deleteCellPopulationInDestructor = false,
bool  initialiseCells = true 
) [inline]

Constructor.

Parameters:
rCellPopulation A cell population object
deleteCellPopulationInDestructor Whether to delete the cell population on destruction to free up memory (defaults to false)
initialiseCells Whether to initialise cells (defaults to true; set to false when loading from an archive)

Definition at line 50 of file AbstractCellBasedSimulation.cpp.

References RandomNumberGenerator::Instance(), AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mInitialiseCells, and AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mrCellPopulation.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::~AbstractCellBasedSimulation (  )  [inline, virtual]

Member Function Documentation

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::AddCellKiller ( boost::shared_ptr< AbstractCellKiller< SPACE_DIM > >  pCellKiller  )  [inline]

Add a cell killer to be used in this simulation.

Parameters:
pCellKiller pointer to a cell killer

Definition at line 275 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mCellKillers.

Referenced by OffLatticeSimulation< ELEMENT_DIM, SPACE_DIM >::OffLatticeSimulation().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::AddSimulationModifier ( boost::shared_ptr< AbstractCellBasedSimulationModifier< ELEMENT_DIM, SPACE_DIM > >  pSimulationModifier  )  [inline]

Add a SimulationModifier to be used in this simulation.

Parameters:
pSimulationModifier pointer to a SimulationModifier

Definition at line 287 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mSimulationModifiers.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
virtual c_vector<double, SPACE_DIM> AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::CalculateCellDivisionVector ( CellPtr  pParentCell  )  [protected, pure virtual]

Method for determining how cell division occurs. This method returns a vector which is then passed into the CellPopulation method AddCell().

As this method is pure virtual, it must be overridden in subclasses.

Parameters:
pParentCell the parent cell
Returns:
a vector containing information on cell division.

Implemented in OffLatticeSimulation< ELEMENT_DIM, SPACE_DIM >, OnLatticeSimulation< DIM >, CryptSimulation1d, CryptSimulation2d, OffLatticeSimulation< 1 >, and OffLatticeSimulation< 2 >.

Referenced by AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::DoCellBirth().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
unsigned AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::DoCellBirth (  )  [inline, protected, virtual]

During a simulation time step, process any cell divisions that need to occur. If the simulation includes cell birth, causes (almost) all cells that are ready to divide to produce daughter cells.

If mOutputDivisionLocations is set to true, then this method also writes the location of each cell division at the present time to mpDivisionLocationFile. This outputs a line of tab-separated values of the form: [time] [div 0 x-pos] [div 0 y-pos] [div 0 z-pos] [div 0 age] [div 1 x-pos] [div 1 y-pos] [div 1 z-pos] [div 1 age] ...

with [y-pos] and [z-pos] included for 2 and 3 dimensional simulations, respectively, and [...age] denoting the age of the dividing cell.

Returns:
the number of births that occurred.

Todo:
(#2578)

For consistency with the rest of the output code, consider removing the AbstractCellBasedSimulation member mOutputDivisionLocations, adding a new member mAgesAndLocationsOfDividingCells to AbstractCellPopulation, adding a new class CellDivisionLocationsWriter to the CellPopulationWriter hierarchy to output the content of mAgesAndLocationsOfDividingCells to file (remembering to clear mAgesAndLocationsOfDividingCells at each timestep), and replacing the following conditional statement with something like

if (mrCellPopulation.HasWriter<CellDivisionLocationsWriter>()) { mCellDivisionLocations.push_back(new_location); }

Definition at line 101 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::CalculateCellDivisionVector(), SimulationTime::GetTime(), SimulationTime::Instance(), AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mNoBirth, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mOutputDivisionLocations, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mpDivisionLocationFile, and AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mrCellPopulation.

Referenced by AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::UpdateCellPopulation().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
unsigned AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::DoCellRemoval (  )  [inline, protected]

During a simulation time step, process any cell sloughing or death

This uses the cell killers to remove cells and associated nodes from the facade class.

Returns:
the number of deaths that occurred.

Definition at line 172 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mCellKillers, and AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mrCellPopulation.

Referenced by AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::UpdateCellPopulation().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
CellBasedPdeHandler< SPACE_DIM > * AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::GetCellBasedPdeHandler (  )  [inline]
Returns:
mpCellBasedPdeHandler

Definition at line 95 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mpCellBasedPdeHandler.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
double AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::GetDt (  )  [inline]
Returns:
the timestep of the simulation

Definition at line 200 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mDt.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
std::vector< double > AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::GetNodeLocation ( const unsigned rNodeIndex  )  [inline]

Get a node's location (ONLY FOR TESTING).

Parameters:
rNodeIndex the node index
Returns:
the co-ordinates of this node.

Definition at line 299 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mrCellPopulation.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
unsigned AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::GetNumBirths (  )  [inline]
Returns:
the number of births that have occurred in the entire simulation (since t=0)

Definition at line 206 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mNumBirths.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
unsigned AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::GetNumDeaths (  )  [inline]
Returns:
the number of deaths that have occurred in the entire simulation (since t=0).

Definition at line 212 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mNumDeaths.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
bool AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::GetOutputCellVelocities (  )  [inline]
Returns:
mOutputCellVelocities

Definition at line 649 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mOutputCellVelocities.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
std::string AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::GetOutputDirectory (  )  [inline]
Returns:
the output directory of the simulation.

Definition at line 232 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mOutputDirectory.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
bool AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::GetOutputDivisionLocations (  )  [inline]
Returns:
mOutputDivisionLocations

Definition at line 637 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mOutputDivisionLocations.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
std::vector< boost::shared_ptr< AbstractCellBasedSimulationModifier< ELEMENT_DIM, SPACE_DIM > > > * AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::GetSimulationModifiers (  )  [inline]
Returns:
a pointer to the vector of SimulationModifiers used in this simulation.

Definition at line 293 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mSimulationModifiers.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
bool AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::GetUpdateCellPopulationRule (  )  [inline]

Return whether to update the topology of the cell population at each time step.

Returns:
whether to update the cell population each time step

Definition at line 263 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mUpdateCellPopulation.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
virtual void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::OutputAdditionalSimulationSetup ( out_stream &  rParamsFile  )  [inline, protected, virtual]

Helper method to output additional simulations parameters and information defined in subclasses to file.

Parameters:
rParamsFile the file stream to which the parameters are output

Reimplemented in OffLatticeSimulation< ELEMENT_DIM, SPACE_DIM >, OnLatticeSimulation< DIM >, OffLatticeSimulation< 1 >, and OffLatticeSimulation< 2 >.

Definition at line 270 of file AbstractCellBasedSimulation.hpp.

Referenced by AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::OutputSimulationSetup().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::OutputSimulationParameters ( out_stream &  rParamsFile  )  [inline, pure virtual]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::OutputSimulationSetup (  )  [inline, protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::RemoveAllCellKillers (  )  [inline]

Method to remove all the cell killers.

Definition at line 281 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mCellKillers.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
const AbstractCellPopulation< ELEMENT_DIM, SPACE_DIM > & AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::rGetCellPopulation (  )  const [inline]
Returns:
const reference to the cell population (used in archiving).

Definition at line 251 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mrCellPopulation.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
AbstractCellPopulation< ELEMENT_DIM, SPACE_DIM > & AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::rGetCellPopulation (  )  [inline]
Returns:
reference to the cell population.

Definition at line 245 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mrCellPopulation.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
template<class Archive >
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::serialize ( Archive &  archive,
const unsigned int  version 
) [inline, private]

Save or restore the simulation.

Parameters:
archive the archive
version the current version of this class

Reimplemented in OffLatticeSimulation< ELEMENT_DIM, SPACE_DIM >, OnLatticeSimulation< DIM >, CryptSimulation1d, CryptSimulation2d, OffLatticeSimulation< 1 >, and OffLatticeSimulation< 2 >.

Definition at line 85 of file AbstractCellBasedSimulation.hpp.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::SetCellBasedPdeHandler ( CellBasedPdeHandler< SPACE_DIM > *  pCellBasedPdeHandler  )  [inline]

Set mpCellBasedPdeHandler

Parameters:
pCellBasedPdeHandler pointer to a CellBasedPdeHandler object

Definition at line 89 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mpCellBasedPdeHandler.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::SetDt ( double  dt  )  [inline]

Set the timestep of the simulation.

Parameters:
dt the timestep to use

Definition at line 193 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mDt.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::SetEndTime ( double  endTime  )  [inline]

Set the end time and resets the timestep to be endtime/100.

Parameters:
endTime the end time to use

Definition at line 218 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mEndTime.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::SetNoBirth ( bool  noBirth  )  [inline]

Set the simulation to run with no birth.

Parameters:
noBirth whether to run with no birth

Definition at line 269 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mNoBirth.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::SetOutputCellVelocities ( bool  outputCellVelocities  )  [inline]

Set mOutputCellVelocities.

Parameters:
outputCellVelocities the new value of mOutputCellVelocities

Definition at line 655 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mOutputCellVelocities.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::SetOutputDirectory ( std::string  outputDirectory  )  [inline]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::SetOutputDivisionLocations ( bool  outputDivisionLocations  )  [inline]

Set mOutputDivisionLocations.

Parameters:
outputDivisionLocations the new value of mOutputDivisionLocations

Definition at line 643 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mOutputDivisionLocations.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::SetSamplingTimestepMultiple ( unsigned  samplingTimestepMultiple  )  [inline]

Set the ratio of the number of actual timesteps to the number of timesteps at which results are written to file. Default value is set to 1 by the constructor.

Parameters:
samplingTimestepMultiple the ratio to use

Definition at line 238 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mSamplingTimestepMultiple.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::SetUpdateCellPopulationRule ( bool  updateCellPopulation  )  [inline]

Set whether to update the topology of the cell population at each time step.

Parameters:
updateCellPopulation whether to update the cell population each time step

Definition at line 257 of file AbstractCellBasedSimulation.cpp.

References AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mUpdateCellPopulation.

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
virtual void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::SetupSolve (  )  [inline, protected, virtual]

A method for subclasses to do something at before the start of the time loop.

Reimplemented in OffLatticeSimulation< ELEMENT_DIM, SPACE_DIM >, CryptSimulation2d, OffLatticeSimulation< 1 >, and OffLatticeSimulation< 2 >.

Definition at line 228 of file AbstractCellBasedSimulation.hpp.

Referenced by AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::Solve().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::Solve (  )  [inline]

Main Solve() method, used to evolve the cell population. Note that prior to calling Solve() we must have called SetEndTime(). We may also have optionally called SetDt(); if not, then a default time step is used.

The Solve() method proceeds as follows.

Setting up:

First, we set up SimulationTime, which (i) provides a globally consistent time, accessible to all other classes in the cell_based code and (ii) handles any rounding issues when the time step does not exactly divide the end time.

Next, we create output files. We then call SetupSolve(), which is empty in the parent class but may be overridden, e.g. to open additional output files. We then call SetupSolve() on any member objects inheriting from AbstractCellBasedSimulationModifier. This class hierarchy allows the user to introduce new ways of updating the cell population within the simulation.

Next, we set up each cell by calling ReadyToDivide() on it, which updates the cell's age and cell cycle model. Finally, we call WriteVisualizerSetupFile() and OutputSimulationSetup(), as well as WriteResultsToFiles() on the cell population, to record the initial configration. This completes the set up process.

The main time loop:

At each time step, we begin by calling UpdateCellPopulation(), which implements any cell deaths and cell divisions through DoCellRemoval() and DoCellBirth() respectively. We then update the correspondence between cells and the mesh by calling Update() on the cell population.

If mOutputCellVelocities is set to true and we are at a printing time, then we als write the velocity of each cell at the present time to mpCellVelocitiesFile. This outputs a line of tab-separated values of the form: [time] [cell 0 x-pos] [cell 0 y-pos] [cell 0 z-pos] [cell 0 x-vel] [cell 0 y-vel] [cell 0 z-vel] ...

with [y-pos] and [z-pos] included for 2 and 3 dimensional simulations, respectively, and data for cells being ordered as given by the cell population Iterator.

Next, we call UpdateCellLocationsAndTopology(), which is pure virtual in the parent class so must be overridden. As the cell population has been updated, we then increment SimulationTime by one time step. We then call UpdateAtEndOfTimeStep() on any AbstractCellBasedSimulationModifiers present, e.g. to write additional output. In an analogous manner to the calls to SetupSolve() prior to entering the main time loop.

The last step within the main time loop is to output the present results to file.

Finishing up:

After exiting the main time loop, we call UpdateCellPopulation() in order to carry out a final update of the cell population. We also call UpdateAtEndOfSolve()} on any member objects inheriting from AbstractCellBasedSimulationModifier in an analogous manner to the aforementioned calls to SetupSolve() UpdateAtEndOfTimeStep(). Finally, we close output files. This completes the Solve() method.

Definition at line 310 of file AbstractCellBasedSimulation.cpp.

References PetscTools::AmMaster(), GenericEventHandler< 11, CellBasedEventHandler >::BeginEvent(), CellBasedPdeHandler< DIM >::CloseResultsFiles(), DOUBLE_UNSET, GenericEventHandler< 11, CellBasedEventHandler >::EndEvent(), EXCEPTION, SimulationTime::GetTime(), SimulationTime::GetTimeStepsElapsed(), SimulationTime::IncrementTimeOneStep(), SimulationTime::Instance(), SimulationTime::IsEndTimeAndNumberOfTimeStepsSetUp(), SimulationTime::IsFinished(), AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mDt, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mEndTime, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mOutputCellVelocities, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mOutputDirectory, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mOutputDivisionLocations, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mpCellBasedPdeHandler, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mpCellVelocitiesFile, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mpDivisionLocationFile, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mpVizSetupFile, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mrCellPopulation, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mSamplingTimestepMultiple, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mSimulationModifiers, AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mSimulationOutputDirectory, OutputFileHandler::OpenOutputFile(), CellBasedPdeHandler< DIM >::OpenResultsFiles(), AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::OutputSimulationSetup(), SimulationTime::ResetEndTimeAndNumberOfTimeSteps(), SimulationTime::SetEndTimeAndNumberOfTimeSteps(), AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::SetupSolve(), CellBasedPdeHandler< DIM >::SolvePdeAndWriteResultsToFile(), AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::StoppingEventHasOccurred(), AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::UpdateCellLocationsAndTopology(), AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::UpdateCellPopulation(), and AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::WriteVisualizerSetupFile().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
bool AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::StoppingEventHasOccurred (  )  [inline, protected, virtual]

A child class can overload this if they want the simulation to stop based on certain conditions before the specified end time (for example, run until a crypt becomes monoclonal).

Returns:
true if stopping event has occurred

Definition at line 596 of file AbstractCellBasedSimulation.cpp.

Referenced by AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::Solve().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
virtual void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::UpdateCellLocationsAndTopology (  )  [protected, pure virtual]

Update the cell locations and topology (connectivity) of the cell population. This method is called within the main time loop of Solve() and is pure virtual in the parent class, so must be overridden.

In the case of an OffLatticeSimulation, the method computes the force acting on each node (corresponding to a cell in centre-based models and to a vertex in vertex-based models) and integrates equations of motion to find the new position of each node.

In the case of an OnLatticeSimulation, the method performs Monte Carlo updating of the cell population, through a call to UpdateCellLocations() on the cell population object.

Implemented in OffLatticeSimulation< ELEMENT_DIM, SPACE_DIM >, OnLatticeSimulation< DIM >, OffLatticeSimulation< 1 >, and OffLatticeSimulation< 2 >.

Referenced by AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::Solve().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::UpdateCellPopulation (  )  [inline, protected, virtual]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
virtual void AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::WriteVisualizerSetupFile (  )  [inline, protected, virtual]

Friends And Related Function Documentation

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
friend class boost::serialization::access [friend]

Member Data Documentation

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
std::vector<boost::shared_ptr<AbstractCellKiller<SPACE_DIM> > > AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mCellKillers [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
bool AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mDeleteCellPopulationInDestructor [protected]

Whether to delete the cell population in the destructor.

Definition at line 120 of file AbstractCellBasedSimulation.hpp.

Referenced by AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::~AbstractCellBasedSimulation().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
double AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mDt [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
double AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mEndTime [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
bool AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mInitialiseCells [protected]

Whether to initialise the cells.

Definition at line 123 of file AbstractCellBasedSimulation.hpp.

Referenced by AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::AbstractCellBasedSimulation().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
bool AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mNoBirth [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
unsigned AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mNumBirths [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
unsigned AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mNumDeaths [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
bool AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mOutputCellVelocities [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
std::string AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mOutputDirectory [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
bool AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mOutputDivisionLocations [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
CellBasedPdeHandler<SPACE_DIM>* AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mpCellBasedPdeHandler [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
out_stream AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mpCellVelocitiesFile [protected]

Results file for cell velocities.

Definition at line 159 of file AbstractCellBasedSimulation.hpp.

Referenced by AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::Solve().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
out_stream AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mpDivisionLocationFile [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
out_stream AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mpVizSetupFile [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
AbstractCellPopulation<ELEMENT_DIM, SPACE_DIM>& AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mrCellPopulation [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
unsigned AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mSamplingTimestepMultiple [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
std::vector<boost::shared_ptr<AbstractCellBasedSimulationModifier<ELEMENT_DIM, SPACE_DIM> > > AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mSimulationModifiers [protected]
template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
std::string AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mSimulationOutputDirectory [protected]

Simulation Output directory either the same as mOutputDirectory or includes mOutputDirectory/results_from_time_[TIME].

Definition at line 135 of file AbstractCellBasedSimulation.hpp.

Referenced by AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::OutputSimulationSetup(), AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::SetOutputDirectory(), and AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::Solve().

template<unsigned ELEMENT_DIM, unsigned SPACE_DIM = ELEMENT_DIM>
bool AbstractCellBasedSimulation< ELEMENT_DIM, SPACE_DIM >::mUpdateCellPopulation [protected]

The documentation for this class was generated from the following files:

Generated by  doxygen 1.6.2