CellGPU
0.8.0
GPU-accelerated simulations of cells
|
A class that implements simple self-propelled particle dynamics in 2D. More...
#include <selfPropelledAligningParticleDynamics.h>
Public Member Functions | |
selfPropelledAligningParticleDynamics () | |
base constructor sets the default time step size | |
selfPropelledAligningParticleDynamics (int N) | |
additionally set the number of particles andinitialize things More... | |
virtual void | integrateEquationsOfMotion () |
the fundamental function that models will call, using vectors of different data structures More... | |
virtual void | integrateEquationsOfMotionCPU () |
call the CPU routine to integrate the e.o.m. More... | |
virtual void | integrateEquationsOfMotionGPU () |
call the GPU routine to integrate the e.o.m. More... | |
Dscalar | getMu () |
Get the inverse friction constant, mu. | |
void | setMu (Dscalar _mu) |
Set the number of degrees of freedom of the equation of motion. | |
void | setJ (Dscalar _J) |
virtual void | spatialSorting () |
call the Simple2DCell spatial vertex sorter, and re-index arrays of cell activity More... | |
virtual void | set2DModel (shared_ptr< Simple2DModel > _model) |
set the active model More... | |
int | getTimestep () |
get the number of timesteps run | |
Dscalar | getTime () |
get the current simulation time | |
virtual void | setDeltaT (Dscalar dt) |
Set the simulation time stepsize. | |
virtual void | performUpdate () |
performUpdate just maps to integrateEquationsOfMotion | |
virtual void | setReproducible (bool rep) |
Set whether the source of noise should always use the same random numbers. | |
void | reIndexRNG (GPUArray< curandState > &array) |
re-index the any RNGs associated with the e.o.m. | |
virtual void | Update (int timestep) |
The fundamental function that a controlling Simulation can call. | |
void | setPeriod (int _p) |
set the period | |
void | setPhase (int _p) |
set the phase | |
virtual void | setGPU () |
Enforce GPU-only operation. This is the default mode, so this method need not be called most of the time. | |
virtual void | setCPU () |
Enforce CPU-only operation. Derived classes might have to do more work when the CPU mode is invoked. | |
int | getNdof () |
Get the number of degrees of freedom of the equation of motion. | |
void | setNdof (int _n) |
Set the number of degrees of freedom of the equation of motion. | |
Public Attributes | |
shared_ptr< Simple2DModel > | model |
A pointer to a Simple2DModel that the updater acts on. | |
Protected Member Functions | |
void | reIndexArray (GPUArray< int > &array) |
Re-index cell arrays after a spatial sorting has occured. | |
void | reIndexArray (GPUArray< Dscalar > &array) |
why use templates when you can type more? | |
void | reIndexArray (GPUArray< Dscalar2 > &array) |
why use templates when you can type more? | |
Protected Attributes | |
shared_ptr< Simple2DActiveCell > | activeModel |
A shared pointer to a simple active model. | |
Dscalar | mu |
The value of the inverse friction constant. | |
Dscalar | J |
The value of the aligning couping. | |
int | Timestep |
Count the number of integration timesteps. | |
Dscalar | deltaT |
The time stepsize of the simulation. | |
GPUArray< Dscalar2 > | displacements |
an internal GPUArray for holding displacements | |
noiseSource | noise |
A source of noise for the equation of motion. | |
int | Period |
The period of the updater... the updater will work every Period timesteps. | |
int | Phase |
The phase of the updater... the updater will work every Period timesteps offset by a phase. | |
bool | GPUcompute |
whether the updater does its work on the GPU or not | |
int | Ndof |
The number of degrees of freedom the equations of motion need to know about. | |
vector< int > | reIndexing |
a vector of the re-indexing information | |
A class that implements simple self-propelled particle dynamics in 2D.
implements dr/dt = mu*F + v_0 {n}, where {n} = (cos(theta),sin(theta))
Additionally there is a tendancy for the directors to align with the particles current instantaneous velocity, theta/dt = -J sin(-)+(brownian noise), where = dr/dt
selfPropelledAligningParticleDynamics::selfPropelledAligningParticleDynamics | ( | int | _N | ) |
additionally set the number of particles andinitialize things
An extremely simple constructor that does nothing, but enforces default GPU operation
the | number of points in the system (cells or particles) |
References simpleEquationOfMotion::deltaT, simpleEquationOfMotion::displacements, updater::GPUcompute, noiseSource::initialize(), J, mu, updater::Ndof, updaterWithNoise::noise, GPUArray< T >::resize(), and simpleEquationOfMotion::Timestep.
|
virtual |
the fundamental function that models will call, using vectors of different data structures
Advances self-propelled dynamics with random noise in the director by one time step
Reimplemented from simpleEquationOfMotion.
References activeModel, simpleEquationOfMotion::displacements, updater::GPUcompute, noiseSource::initialize(), integrateEquationsOfMotionCPU(), integrateEquationsOfMotionGPU(), updater::Ndof, updaterWithNoise::noise, GPUArray< T >::resize(), and simpleEquationOfMotion::Timestep.
Referenced by selfPropelledAligningParticleDynamics().
|
virtual |
call the CPU routine to integrate the e.o.m.
The straightforward CPU implementation
References activeModel, simpleEquationOfMotion::displacements, access_location::host, updater::Ndof, access_mode::overwrite, and access_mode::read.
Referenced by integrateEquationsOfMotion(), and selfPropelledAligningParticleDynamics().
|
virtual |
call the GPU routine to integrate the e.o.m.
The straightforward GPU implementation
References activeModel, ArrayHandle< T >::data, simpleEquationOfMotion::deltaT, access_location::device, simpleEquationOfMotion::displacements, gpu_spp_aligning_eom_integration(), J, mu, updater::Ndof, updaterWithNoise::noise, access_mode::overwrite, access_mode::read, access_mode::readwrite, noiseSource::RNGs, and simpleEquationOfMotion::Timestep.
Referenced by integrateEquationsOfMotion(), and selfPropelledAligningParticleDynamics().
|
virtual |
call the Simple2DCell spatial vertex sorter, and re-index arrays of cell activity
When spatial sorting is performed, re-index the array of cuda RNGs... This function is currently commented out, for greater flexibility (i.e., to not require that the indexToTag (or Itt) be the re-indexing array), since that assumes cell and not particle-based dynamics
Reimplemented from updater.
|
virtual |
set the active model
Set the shared pointer of the base class to passed variable; cast it as an active cell model
Reimplemented from updater.
References activeModel, and updater::model.