|
Simpatico
v1.10
|
Implementation template for a AnglePotential.
Definition at line 33 of file ddMd/potentials/angle/AnglePotentialImpl.h.
#include <AnglePotentialImpl.h>
Public Member Functions | |
| AnglePotentialImpl (Simulation &simulation) | |
| Constructor. More... | |
| AnglePotentialImpl () | |
| Default constructor. More... | |
| virtual | ~AnglePotentialImpl () |
| Destructor. More... | |
| virtual void | setNAngleType (int nAtomType) |
| Set the maximum number of atom types. More... | |
| virtual void | readParameters (std::istream &in) |
| Read potential energy parameters. More... | |
| virtual void | loadParameters (Serializable::IArchive &ar) |
| Load internal state from an archive. More... | |
| virtual void | save (Serializable::OArchive &ar) |
| Save internal state to an archive. More... | |
Interaction interface | |
| double | angleEnergy (double cosTheta, int type) const |
| Returns potential energy for one angle. More... | |
| void | angleForce (const Vector &R1, const Vector &R2, Vector &F1, Vector &F2, int type) const |
| Returns forces along two bonds at the angle, for use in MD and stress calculation. More... | |
| void | set (std::string name, int type, double value) |
| Modify an angle parameter, identified by a string. More... | |
| double | get (std::string name, int type) const |
| Get a parameter value, identified by a string. More... | |
| virtual std::string | interactionClassName () const |
| Return pair interaction class name (e.g., "CosineAngle"). More... | |
| const Interaction & | interaction () const |
| Return angle interaction by const reference. More... | |
| Interaction & | interaction () |
| Return angle interaction by reference. More... | |
Total Energy, Force and Stress | |
| virtual void | computeForces () |
| Compute angle forces for atoms. More... | |
| virtual void | computeEnergy (MPI::Intracomm &communicator) |
| Compute the total angle energy for all processors. More... | |
| virtual void | computeStress (MPI::Intracomm &communicator) |
| Compute the covalent bond stress. More... | |
Public Member Functions inherited from DdMd::AnglePotential | |
| AnglePotential (Simulation &simulation) | |
| Constructor. More... | |
| AnglePotential () | |
| Default constructor (for unit testing). More... | |
| ~AnglePotential () | |
| Destructor. More... | |
| void | associate (Boundary &boundary, GroupStorage< 3 > &storage) |
| Associate with related objects. More... | |
| double | angleEnergy (double cosTheta, int type) const |
| Returns potential energy for one angle. More... | |
| void | angleForce (const Vector &R1, const Vector &R2, Vector &F1, Vector &F2, int type) const |
| Computes forces along two bonds within the angle. More... | |
Public Member Functions inherited from DdMd::Potential | |
| Potential () | |
| Constructor. More... | |
| virtual | ~Potential () |
| Destructor. More... | |
| void | setReverseUpdateFlag (bool reverseUpdateFlag) |
| Set flag to identify if reverse communication is enabled. More... | |
| bool | reverseUpdateFlag () const |
| Get flag to identify if reverse communication is enabled. More... | |
| double | energy () const |
| Return the total potential, from all processors. More... | |
| void | unsetEnergy () |
| Mark the energy as unknown (nullify). More... | |
| bool | isEnergySet () const |
| Is the energy set (known)? More... | |
| virtual void | computeForcesAndStress (MPI::Intracomm &communicator) |
| Compute forces and stress for all processors. More... | |
| Tensor | stress () const |
| Return the stress tensor. More... | |
| double | pressure () const |
| Return the pressure. More... | |
| void | unsetStress () |
| Mark the stress as unknown (nullify). More... | |
| bool | isStressSet () const |
| Is the stress set (known)? More... | |
| virtual bool | isValid (MPI::Intracomm &communicator) const |
| Is the potential in a valid internal state? More... | |
Public Member Functions inherited from Util::ParamComposite | |
| 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... | |
Public Member Functions inherited from Util::ParamComponent | |
| 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... | |
Public Member Functions inherited from Util::Serializable | |
| virtual | ~Serializable () |
| Destructor. More... | |
Public Member Functions inherited from Util::MpiFileIo | |
| 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... | |
Additional Inherited Members | |
Public Types inherited from Util::Serializable | |
| 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 Public Member Functions inherited from Util::ParamComponent | |
| 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... | |
Protected Member Functions inherited from DdMd::AnglePotential | |
| Boundary & | boundary () const |
| Return boundary by reference. More... | |
| GroupStorage< 3 > & | storage () const |
| Return bond storage by reference. More... | |
Protected Member Functions inherited from DdMd::Potential | |
| void | setEnergy (double energy) |
| Set a value for the total energy. More... | |
| void | setStress (const Tensor &stress) |
| Set a value for the total stress. More... | |
| void | incrementPairStress (const Vector &f, const Vector &dr, Tensor &stress) const |
| Add a pair contribution to the stress tensor. More... | |
| void | reduceEnergy (double localEnergy, MPI::Intracomm &communicator) |
| Add local energies from all processors, set energy on master. More... | |
| void | reduceStress (Tensor &localStress, MPI::Intracomm &communicator) |
| Add local stresses from all processors, set total on master. More... | |
Protected Member Functions inherited from Util::ParamComposite | |
| 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... | |
Protected Member Functions inherited from Util::ParamComponent | |
| ParamComponent () | |
| Constructor. More... | |
| ParamComponent (const ParamComponent &other) | |
| Copy constructor. More... | |
| DdMd::AnglePotentialImpl< Interaction >::AnglePotentialImpl | ( | Simulation & | simulation | ) |
Constructor.
Definition at line 215 of file ddMd/potentials/angle/AnglePotentialImpl.h.
References DdMd::Simulation::nAngleType(), and DdMd::AnglePotentialImpl< Interaction >::setNAngleType().
| DdMd::AnglePotentialImpl< Interaction >::AnglePotentialImpl | ( | ) |
Default constructor.
Definition at line 228 of file ddMd/potentials/angle/AnglePotentialImpl.h.
|
virtual |
Destructor.
Definition at line 238 of file ddMd/potentials/angle/AnglePotentialImpl.h.
|
virtual |
Set the maximum number of atom types.
Implements DdMd::AnglePotential.
Definition at line 249 of file ddMd/potentials/angle/AnglePotentialImpl.h.
References DdMd::AnglePotentialImpl< Interaction >::interaction().
Referenced by DdMd::AnglePotentialImpl< Interaction >::AnglePotentialImpl().
|
virtual |
Read potential energy parameters.
This method reads the angle potential Interaction parameter block. Before calling Interaction::readParameters(), it passes simulation().nAngleType() to Interaction::setNAtomType().
Reimplemented from Util::ParamComposite.
Definition at line 256 of file ddMd/potentials/angle/AnglePotentialImpl.h.
References Util::ParamComposite::addParamComposite(), DdMd::AnglePotentialImpl< Interaction >::interaction(), and UTIL_CHECK.
|
virtual |
Load internal state from an archive.
Precondition: setNAngleType must have been called before this.
| ar | input/loading archive |
Reimplemented from Util::ParamComposite.
Definition at line 270 of file ddMd/potentials/angle/AnglePotentialImpl.h.
References Util::ParamComposite::addParamComposite(), DdMd::AnglePotentialImpl< Interaction >::interaction(), and UTIL_CHECK.
|
virtual |
Save internal state to an archive.
| ar | output/saving archive |
Reimplemented from Util::ParamComposite.
Definition at line 283 of file ddMd/potentials/angle/AnglePotentialImpl.h.
References DdMd::AnglePotentialImpl< Interaction >::interaction().
| double DdMd::AnglePotentialImpl< Interaction >::angleEnergy | ( | double | cosTheta, |
| int | type | ||
| ) | const |
Returns potential energy for one angle.
| cosTheta | cosine of the bend angle |
| type | type index of angle group |
Definition at line 291 of file ddMd/potentials/angle/AnglePotentialImpl.h.
References DdMd::AnglePotentialImpl< Interaction >::interaction().
Referenced by DdMd::AnglePotentialImpl< Interaction >::computeEnergy().
| void DdMd::AnglePotentialImpl< Interaction >::angleForce | ( | const Vector & | R1, |
| const Vector & | R2, | ||
| Vector & | F1, | ||
| Vector & | F2, | ||
| int | type | ||
| ) | const |
Returns forces along two bonds at the angle, for use in MD and stress calculation.
| R1 | bond vector from atom 1 to 2 |
| R2 | bond vector from atom 2 to 3 |
| F1 | return force along R1 direction |
| F2 | return force along R2 direction |
| type | type index for angle group |
Definition at line 299 of file ddMd/potentials/angle/AnglePotentialImpl.h.
References DdMd::AnglePotentialImpl< Interaction >::interaction().
|
inlinevirtual |
Modify an angle parameter, identified by a string.
| name | parameter variable name |
| type | type index for angle group |
| value | new value of parameter |
Implements DdMd::AnglePotential.
Definition at line 114 of file ddMd/potentials/angle/AnglePotentialImpl.h.
|
inlinevirtual |
Get a parameter value, identified by a string.
| name | parameter variable name |
| type | type index of angle group |
Implements DdMd::AnglePotential.
Definition at line 123 of file ddMd/potentials/angle/AnglePotentialImpl.h.
|
virtual |
Return pair interaction class name (e.g., "CosineAngle").
Implements DdMd::AnglePotential.
Definition at line 309 of file ddMd/potentials/angle/AnglePotentialImpl.h.
References DdMd::AnglePotentialImpl< Interaction >::interaction().
|
inline |
Return angle interaction by const reference.
Get Interaction by const reference.
Definition at line 323 of file ddMd/potentials/angle/AnglePotentialImpl.h.
Referenced by DdMd::AnglePotentialImpl< Interaction >::angleEnergy(), DdMd::AnglePotentialImpl< Interaction >::angleForce(), DdMd::AnglePotentialImpl< Interaction >::computeEnergy(), DdMd::AnglePotentialImpl< Interaction >::computeForces(), DdMd::AnglePotentialImpl< Interaction >::computeStress(), DdMd::AnglePotentialImpl< Interaction >::interactionClassName(), DdMd::AnglePotentialImpl< Interaction >::loadParameters(), DdMd::AnglePotentialImpl< Interaction >::readParameters(), DdMd::AnglePotentialImpl< Interaction >::save(), and DdMd::AnglePotentialImpl< Interaction >::setNAngleType().
|
inline |
Return angle interaction by reference.
Get Interaction by reference.
Definition at line 316 of file ddMd/potentials/angle/AnglePotentialImpl.h.
|
virtual |
Compute angle forces for atoms.
Adds angle forces to existing Atom forces.
Implements DdMd::Potential.
Definition at line 330 of file ddMd/potentials/angle/AnglePotentialImpl.h.
References DdMd::GroupStorage< N >::begin(), DdMd::AnglePotential::boundary(), Simp::OrthorhombicBoundary::distanceSq(), DdMd::Atom::force(), DdMd::AnglePotentialImpl< Interaction >::interaction(), DdMd::Atom::isGhost(), Util::PArrayIterator< Group< N > >::notEnd(), DdMd::Atom::position(), and DdMd::AnglePotential::storage().
|
virtual |
Compute the total angle energy for all processors.
Call on all processors (MPI reduce operation).
Implements DdMd::Potential.
Definition at line 373 of file ddMd/potentials/angle/AnglePotentialImpl.h.
References DdMd::AnglePotentialImpl< Interaction >::angleEnergy(), DdMd::GroupStorage< N >::begin(), DdMd::AnglePotential::boundary(), Simp::OrthorhombicBoundary::distanceSq(), Util::Vector::dot(), DdMd::AnglePotentialImpl< Interaction >::interaction(), DdMd::Potential::isEnergySet(), DdMd::Atom::isGhost(), Util::PArrayIterator< Group< N > >::notEnd(), DdMd::Atom::position(), DdMd::Potential::reduceEnergy(), and DdMd::AnglePotential::storage().
|
virtual |
Compute the covalent bond stress.
Call on all processors.
Reimplemented from DdMd::Potential.
Definition at line 424 of file ddMd/potentials/angle/AnglePotentialImpl.h.
References DdMd::GroupStorage< N >::begin(), DdMd::AnglePotential::boundary(), Simp::OrthorhombicBoundary::distanceSq(), DdMd::Potential::incrementPairStress(), DdMd::AnglePotentialImpl< Interaction >::interaction(), DdMd::Atom::isGhost(), DdMd::Potential::isStressSet(), Util::PArrayIterator< Group< N > >::notEnd(), DdMd::Atom::position(), DdMd::Potential::reduceStress(), DdMd::AnglePotential::storage(), Simp::OrthorhombicBoundary::volume(), and Util::Tensor::zero().
1.8.11