Simpatico
v1.10
|
A two-step velocity-Verlet style integrator.
Definition at line 24 of file TwoStepIntegrator.h.
#include <TwoStepIntegrator.h>
Public Member Functions | |
TwoStepIntegrator (Simulation &simulation) | |
Constructor. More... | |
~TwoStepIntegrator () | |
Destructor. More... | |
void | run (int nStep) |
Run a simulation. More... | |
![]() | |
Integrator (Simulation &simulation) | |
Constructor. More... | |
~Integrator () | |
Destructor. More... | |
void | readParameters (std::istream &in) |
Read saveInterval and saveFileName. More... | |
void | loadParameters (Serializable::IArchive &ar) |
Load saveInterval and saveFileName from restart archive. More... | |
void | save (Serializable::OArchive &ar) |
Save saveInterval and saveFileName from restart archive. More... | |
virtual void | clear () |
Set integrator to initial state and clears all statistics. More... | |
virtual void | outputStatistics (std::ostream &out) |
Output timing statistics from a run. More... | |
double | time () const |
Get average time per processor of previous run. More... | |
int | iStep () const |
Get current time step index. More... | |
![]() | |
ParamComposite () | |
Constructor. More... | |
ParamComposite (const ParamComposite &other) | |
Copy constructor. More... | |
ParamComposite (int capacity) | |
Constructor. More... | |
virtual | ~ParamComposite () |
Virtual destructor. More... | |
void | resetParam () |
Resets ParamComposite to its empty state. More... | |
virtual void | readParam (std::istream &in) |
Read the parameter file block. More... | |
virtual void | readParamOptional (std::istream &in) |
Read optional parameter file block. More... | |
virtual void | writeParam (std::ostream &out) |
Write all parameters to an output stream. More... | |
virtual void | load (Serializable::IArchive &ar) |
Load all parameters from an input archive. More... | |
virtual void | loadOptional (Serializable::IArchive &ar) |
Load an optional ParamComposite. More... | |
void | saveOptional (Serializable::OArchive &ar) |
Saves isActive flag, and then calls save() iff isActive is true. More... | |
void | readParamComposite (std::istream &in, ParamComposite &child, bool next=true) |
Add and read a required child ParamComposite. More... | |
void | readParamCompositeOptional (std::istream &in, ParamComposite &child, bool next=true) |
Add and attempt to read an optional child ParamComposite. More... | |
template<typename Type > | |
ScalarParam< Type > & | read (std::istream &in, const char *label, Type &value) |
Add and read a new required ScalarParam < Type > object. More... | |
template<typename Type > | |
ScalarParam< Type > & | readOptional (std::istream &in, const char *label, Type &value) |
Add and read a new optional ScalarParam < Type > object. More... | |
template<typename Type > | |
CArrayParam< Type > & | readCArray (std::istream &in, const char *label, Type *value, int n) |
Add and read a required C array parameter. More... | |
template<typename Type > | |
CArrayParam< Type > & | readOptionalCArray (std::istream &in, const char *label, Type *value, int n) |
Add and read an optional C array parameter. More... | |
template<typename Type > | |
DArrayParam< Type > & | readDArray (std::istream &in, const char *label, DArray< Type > &array, int n) |
Add and read a required DArray < Type > parameter. More... | |
template<typename Type > | |
DArrayParam< Type > & | readOptionalDArray (std::istream &in, const char *label, DArray< Type > &array, int n) |
Add and read an optional DArray < Type > parameter. More... | |
template<typename Type , int N> | |
FArrayParam< Type, N > & | readFArray (std::istream &in, const char *label, FArray< Type, N > &array) |
Add and read a required FArray < Type, N > array parameter. More... | |
template<typename Type , int N> | |
FArrayParam< Type, N > & | readOptionalFArray (std::istream &in, const char *label, FArray< Type, N > &array) |
Add and read an optional FArray < Type, N > array parameter. More... | |
template<typename Type > | |
CArray2DParam< Type > & | readCArray2D (std::istream &in, const char *label, Type *value, int m, int n, int np) |
Add and read a required CArray2DParam < Type > 2D C-array. More... | |
template<typename Type > | |
CArray2DParam< Type > & | readOptionalCArray2D (std::istream &in, const char *label, Type *value, int m, int n, int np) |
Add and read an optional CArray2DParam < Type > 2D C-array parameter. More... | |
template<typename Type > | |
DMatrixParam< Type > & | readDMatrix (std::istream &in, const char *label, DMatrix< Type > &matrix, int m, int n) |
Add and read a required DMatrix < Type > matrix parameter. More... | |
template<typename Type > | |
DMatrixParam< Type > & | readOptionalDMatrix (std::istream &in, const char *label, DMatrix< Type > &matrix, int m, int n) |
Add and read an optional DMatrix < Type > matrix parameter. More... | |
template<typename Type > | |
DSymmMatrixParam< Type > & | readDSymmMatrix (std::istream &in, const char *label, DMatrix< Type > &matrix, int n) |
Add and read a required symmetrix DMatrix. More... | |
template<typename Type > | |
DSymmMatrixParam< Type > & | readOptionalDSymmMatrix (std::istream &in, const char *label, DMatrix< Type > &matrix, int n) |
Add and read an optional DMatrix matrix parameter. More... | |
Begin & | readBegin (std::istream &in, const char *label, bool isRequired=true) |
Add and read a class label and opening bracket. More... | |
End & | readEnd (std::istream &in) |
Add and read the closing bracket. More... | |
Blank & | readBlank (std::istream &in) |
Add and read a new Blank object, representing a blank line. More... | |
void | loadParamComposite (Serializable::IArchive &ar, ParamComposite &child, bool next=true) |
Add and load a required child ParamComposite. More... | |
void | loadParamCompositeOptional (Serializable::IArchive &ar, ParamComposite &child, bool next=true) |
Add and load an optional child ParamComposite if isActive. More... | |
template<typename Type > | |
ScalarParam< Type > & | loadParameter (Serializable::IArchive &ar, const char *label, Type &value, bool isRequired) |
Add and load a new ScalarParam < Type > object. More... | |
template<typename Type > | |
ScalarParam< Type > & | loadParameter (Serializable::IArchive &ar, const char *label, Type &value) |
Add and load new required ScalarParam < Type > object. More... | |
template<typename Type > | |
CArrayParam< Type > & | loadCArray (Serializable::IArchive &ar, const char *label, Type *value, int n, bool isRequired) |
Add a C array parameter and load its elements. More... | |
template<typename Type > | |
CArrayParam< Type > & | loadCArray (Serializable::IArchive &ar, const char *label, Type *value, int n) |
Add and load a required CArrayParam< Type > array parameter. More... | |
template<typename Type > | |
DArrayParam< Type > & | loadDArray (Serializable::IArchive &ar, const char *label, DArray< Type > &array, int n, bool isRequired) |
Add an load a DArray < Type > array parameter. More... | |
template<typename Type > | |
DArrayParam< Type > & | loadDArray (Serializable::IArchive &ar, const char *label, DArray< Type > &array, int n) |
Add and load a required DArray< Type > array parameter. More... | |
template<typename Type , int N> | |
FArrayParam< Type, N > & | loadFArray (Serializable::IArchive &ar, const char *label, FArray< Type, N > &array, bool isRequired) |
Add and load an FArray < Type, N > fixed-size array parameter. More... | |
template<typename Type , int N> | |
FArrayParam< Type, N > & | loadFArray (Serializable::IArchive &ar, const char *label, FArray< Type, N > &array) |
Add and load a required FArray < Type > array parameter. More... | |
template<typename Type > | |
CArray2DParam< Type > & | loadCArray2D (Serializable::IArchive &ar, const char *label, Type *value, int m, int n, int np, bool isRequired) |
Add and load a CArray2DParam < Type > C 2D array parameter. More... | |
template<typename Type > | |
CArray2DParam< Type > & | loadCArray2D (Serializable::IArchive &ar, const char *label, Type *value, int m, int n, int np) |
Add and load a required < Type > matrix parameter. More... | |
template<typename Type > | |
DMatrixParam< Type > & | loadDMatrix (Serializable::IArchive &ar, const char *label, DMatrix< Type > &matrix, int m, int n, bool isRequired) |
Add and load a DMatrixParam < Type > matrix parameter. More... | |
template<typename Type > | |
DMatrixParam< Type > & | loadDMatrix (Serializable::IArchive &ar, const char *label, DMatrix< Type > &matrix, int m, int n) |
Add and load a required DMatrixParam < Type > matrix parameter. More... | |
template<typename Type > | |
DSymmMatrixParam< Type > & | loadDSymmMatrix (Serializable::IArchive &ar, const char *label, DMatrix< Type > &matrix, int n, bool isRequired) |
Add and load a symmetric DSymmMatrixParam < Type > matrix parameter. More... | |
template<typename Type > | |
DSymmMatrixParam< Type > & | loadDSymmMatrix (Serializable::IArchive &ar, const char *label, DMatrix< Type > &matrix, int n) |
Add and load a required DSymmMatrixParam < Type > matrix parameter. More... | |
void | addParamComposite (ParamComposite &child, bool next=true) |
Add a child ParamComposite object to the format array. More... | |
Begin & | addBegin (const char *label) |
Add a Begin object representing a class name and bracket. More... | |
End & | addEnd () |
Add a closing bracket. More... | |
Blank & | addBlank () |
Create and add a new Blank object, representing a blank line. More... | |
std::string | className () const |
Get class name string. More... | |
bool | isRequired () const |
Is this ParamComposite required in the input file? More... | |
bool | isActive () const |
Is this parameter active? More... | |
![]() | |
virtual | ~ParamComponent () |
Destructor. More... | |
void | setIndent (const ParamComponent &parent, bool next=true) |
Set indent level. More... | |
std::string | indent () const |
Return indent string for this object (string of spaces). More... | |
template<class Archive > | |
void | serialize (Archive &ar, const unsigned int version) |
Serialize this ParamComponent as a string. More... | |
![]() | |
virtual | ~Serializable () |
Destructor. More... | |
![]() | |
MpiFileIo () | |
Constructor. More... | |
MpiFileIo (const MpiFileIo &other) | |
Copy constructor. More... | |
bool | isIoProcessor () const |
Can this processor do file I/O ? More... | |
void | setIoCommunicator (MPI::Intracomm &communicator) |
Set the communicator. More... | |
void | clearCommunicator () |
Clear (nullify) the communicator. More... | |
bool | hasIoCommunicator () const |
Does this object have an associated MPI communicator? More... | |
MPI::Intracomm & | ioCommunicator () const |
Get the MPI communicator by reference. More... | |
![]() | |
SimulationAccess (Simulation &simulation) | |
Constructor. More... | |
~SimulationAccess () | |
Destructor. More... | |
Simulation & | simulation () |
Get the parent simulation. More... | |
Boundary & | boundary () |
Get the Boundary. More... | |
AtomStorage & | atomStorage () |
Get the AtomStorage. More... | |
GroupStorage< 2 > & | bondStorage () |
Get the BondStorage. More... | |
GroupStorage< 3 > & | angleStorage () |
Get the AngleStorage. More... | |
GroupStorage< 4 > & | dihedralStorage () |
Get the angleStorage. More... | |
PairPotential & | pairPotential () |
Get the PairPotential. More... | |
BondPotential & | bondPotential () |
Get the PairPotential. More... | |
AnglePotential & | anglePotential () |
Get the AnglePotential. More... | |
DihedralPotential & | dihedralPotential () |
Get the DihedralPotential. More... | |
ExternalPotential & | externalPotential () |
Get the ExternalPotential. More... | |
EnergyEnsemble & | energyEnsemble () |
Get the EnergyEnsemble. More... | |
BoundaryEnsemble & | boundaryEnsemble () |
Get the BoundaryEnsemble. More... | |
Random & | random () |
Get the Random number generator. More... | |
Domain & | domain () |
Get the Domain. More... | |
Exchanger & | exchanger () |
Get the Exchanger. More... | |
FileMaster & | fileMaster () |
Get the FileMaster. More... | |
int | nAtomType () |
Get maximum number of atom types. More... | |
int | nBondType () |
Get maximum number of bond types. More... | |
int | nAngleType () |
Get maximum number of angle types. More... | |
int | nDihedralType () |
Get maximum number of dihedral types. More... | |
bool | hasExternal () |
Does this simulation have an external potential? More... | |
AtomType & | atomType (int i) |
Get an AtomType descriptor for atomtype i. More... | |
MaskPolicy | maskedPairPolicy () const |
Return the value of the mask policy (MaskNone or MaskBonded). More... | |
bool | reverseUpdateFlag () const |
Are forces evaluated by reverse communication (true) or not (false)? More... | |
Protected Member Functions | |
virtual void | integrateStep1 ()=0 |
Execute first step of two-step integrator. More... | |
virtual void | integrateStep2 ()=0 |
Execute secodn step of two-step integrator. More... | |
![]() | |
virtual void | initDynamicalState () |
Set any internal dynamical to default initial values. More... | |
virtual void | setup ()=0 |
Setup integrator just before integration. More... | |
void | setIsSetup () |
Mark the integrator as having been setup at least once. More... | |
bool | isSetup () const |
Has the setup() method been called at least once previously? More... | |
void | setupAtoms () |
Setup state of atoms just before integration. More... | |
void | computeForces () |
Compute forces for all local atoms, with timing. More... | |
void | computeForcesAndVirial () |
Compute forces for all local atoms and virial, with timing. More... | |
bool | isExchangeNeeded (double skin) |
Determine whether an atom exchange and reneighboring is needed. More... | |
const std::string & | saveFileName () const |
Get restart file base name. More... | |
int | saveInterval () const |
Get interval for writing a restart file. More... | |
![]() | |
void | setClassName (const char *className) |
Set class name string. More... | |
void | setIsRequired (bool isRequired) |
Set or unset the isActive flag. More... | |
void | setIsActive (bool isActive) |
Set or unset the isActive flag. More... | |
void | setParent (ParamComponent ¶m, bool next=true) |
Set this to the parent of a child component. More... | |
void | addComponent (ParamComponent ¶m, bool isLeaf=true) |
Add a new ParamComponent object to the format array. More... | |
template<typename Type > | |
ScalarParam< Type > & | add (std::istream &in, const char *label, Type &value, bool isRequired=true) |
Add a new required ScalarParam < Type > object. More... | |
template<typename Type > | |
CArrayParam< Type > & | addCArray (std::istream &in, const char *label, Type *value, int n, bool isRequired=true) |
Add (but do not read) a required C array parameter. More... | |
template<typename Type > | |
DArrayParam< Type > & | addDArray (std::istream &in, const char *label, DArray< Type > &array, int n, bool isRequired=true) |
Add (but do not read) a DArray < Type > parameter. More... | |
template<typename Type , int N> | |
FArrayParam< Type, N > & | addFArray (std::istream &in, const char *label, FArray< Type, N > &array, bool isRequired=true) |
Add (but do not read) a FArray < Type, N > array parameter. More... | |
template<typename Type > | |
CArray2DParam< Type > & | addCArray2D (std::istream &in, const char *label, Type *value, int m, int n, int np, bool isRequired=true) |
Add (but do not read) a CArray2DParam < Type > 2D C-array. More... | |
template<typename Type > | |
DMatrixParam< Type > & | addDMatrix (std::istream &in, const char *label, DMatrix< Type > &matrix, int m, int n, bool isRequired=true) |
Add and read a required DMatrix < Type > matrix parameter. More... | |
![]() | |
ParamComponent () | |
Constructor. More... | |
ParamComponent (const ParamComponent &other) | |
Copy constructor. More... | |
Additional Inherited Members | |
![]() | |
typedef BinaryFileOArchive | OArchive |
Type of output archive used by save method. More... | |
typedef BinaryFileIArchive | IArchive |
Type of input archive used by load method. More... | |
![]() | |
static void | initStatic () |
Initialize static echo member to false. More... | |
static void | setEcho (bool echo=true) |
Enable or disable echoing for all subclasses of ParamComponent. More... | |
static bool | echo () |
Get echo parameter. More... | |
![]() | |
![]() | |
int | iStep_ |
Current step number. More... | |
DdMd::TwoStepIntegrator::TwoStepIntegrator | ( | Simulation & | simulation | ) |
Constructor.
Definition at line 33 of file TwoStepIntegrator.cpp.
DdMd::TwoStepIntegrator::~TwoStepIntegrator | ( | ) |
Destructor.
Definition at line 40 of file TwoStepIntegrator.cpp.
|
virtual |
Run a simulation.
nStep | number of steps |
Implements DdMd::Integrator.
Definition at line 46 of file TwoStepIntegrator.cpp.
References DdMd::Simulation::analyzerManager(), DdMd::SimulationAccess::atomStorage(), DdMd::SimulationAccess::boundary(), DdMd::SimulationAccess::boundaryEnsemble(), DdMd::PairPotential::buildCellList(), DdMd::PairPotential::buildPairList(), DdMd::AtomStorage::clearSnapshot(), DdMd::Integrator::computeForces(), DdMd::Integrator::computeForcesAndVirial(), DdMd::AtomStorage::computeNAtomTotal(), DdMd::SimulationAccess::domain(), DdMd::ModifierManager::endOfStep(), DdMd::Exchanger::exchange(), DdMd::SimulationAccess::exchanger(), DdMd::Simulation::exchangeSignal(), Util::Log::file(), integrateStep1(), integrateStep2(), DdMd::Integrator::isExchangeNeeded(), DdMd::Integrator::iStep_, DdMd::Simulation::isValid(), DdMd::AtomStorage::makeSnapshot(), DdMd::Simulation::modifierManager(), DdMd::Simulation::modifySignal(), Util::Signal< T >::notify(), DdMd::SimulationAccess::pairPotential(), DdMd::ModifierManager::postExchange(), DdMd::ModifierManager::postForce(), DdMd::ModifierManager::postIntegrate1(), DdMd::ModifierManager::postNeighbor(), DdMd::ModifierManager::postUpdate(), DdMd::ModifierManager::preExchange(), DdMd::ModifierManager::preForce(), DdMd::ModifierManager::preIntegrate1(), DdMd::ModifierManager::preTransform(), DdMd::ModifierManager::preUpdate(), DdMd::AnalyzerManager::sample(), DdMd::Simulation::save(), DdMd::Integrator::saveFileName(), DdMd::Integrator::saveInterval(), DdMd::AnalyzerManager::setup(), DdMd::ModifierManager::setup(), DdMd::Integrator::setup(), DdMd::SimulationAccess::simulation(), DdMd::DdTimer::stamp(), DdMd::DdTimer::start(), DdMd::DdTimer::stop(), DdMd::Exchanger::timer(), DdMd::AtomStorage::transformCartToGen(), DdMd::AtomStorage::transformGenToCart(), DdMd::AtomStorage::unsetNAtomTotal(), DdMd::Exchanger::update(), and UTIL_THROW.
|
protectedpure virtual |
Execute first step of two-step integrator.
The function simulation().modifySignal().notify() is called immediately after this function in the run method, so this function need not send this signal.
Implemented in DdMd::NvtLangevinIntegrator, DdMd::NptIntegrator, DdMd::NveIntegrator, DdMd::NphIntegrator, and DdMd::NvtIntegrator.
Referenced by run().
|
protectedpure virtual |
Execute secodn step of two-step integrator.
This method must call simulation().velocitySignal().notify() if it modifies the velocities.
Implemented in DdMd::NvtLangevinIntegrator, DdMd::NptIntegrator, DdMd::NveIntegrator, DdMd::NphIntegrator, and DdMd::NvtIntegrator.
Referenced by run().