AMDIntegrator¶

class
simtk.openmm.amd.
AMDIntegrator
(dt, alpha, E)¶ AMDIntegrator implements the aMD integration algorithm.
The system is integrated based on a modified potential. Whenever the energy V(r) is less than a cutoff value E, the following effective potential is used:
V*(r) = V(r) + (EV(r))^2 / (alpha+EV(r))
For details, see Hamelberg et al., J. Chem. Phys. 127, 155102 (2007).

__init__
(dt, alpha, E)¶ Create an AMDIntegrator.
Parameters:  dt (time) – The integration time step to use
 alpha (energy) – The alpha parameter to use
 E (energy) – The energy cutoff to use
Methods
__init__
(dt, alpha, E)Create an AMDIntegrator. addComputeGlobal
(self, variable, expression)Add a step to the integration algorithm that computes a global value. addComputePerDof
(self, variable, expression)Add a step to the integration algorithm that computes a perDOF value. addComputeSum
(self, variable, expression)Add a step to the integration algorithm that computes a sum over degrees of freedom. addConstrainPositions
(self)Add a step to the integration algorithm that updates particle positions so all constraints are satisfied. addConstrainVelocities
(self)Add a step to the integration algorithm that updates particle velocities so the net velocity along all constraints is 0. addGlobalVariable
(self, name, initialValue)Define a new global variable. addPerDofVariable
(self, name, initialValue)Define a new perDOF variable. addTabulatedFunction
(self, name, function)Add a tabulated function that may appear in expressions. addUpdateContextState
(self)Add a step to the integration algorithm that allows Forces to update the context state. beginIfBlock
(self, condition)Add a step which begins a new “if” block. beginWhileBlock
(self, condition)Add a step which begins a new “while” block. endBlock
(self)Add a step which marks the end of the most recently begun “if” or “while” block. getAlpha
()Get the value of alpha for the integrator. getComputationStep
(self, index)Get the details of a computation step that has been added to the integration algorithm. getConstraintTolerance
(self)Get the distance tolerance within which constraints are maintained, as a fraction of the constrained distance. getE
()Get the energy threshold E for the integrator. getEffectiveEnergy
(energy)Given the actual potential energy of the system, return the value of the effective potential. getGlobalVariable
(self, index)Get the current value of a global variable. getGlobalVariableByName
(self, name)Get the current value of a global variable, specified by name. getGlobalVariableName
(self, index)Get the name of a global variable. getKineticEnergyExpression
(self)Get the expression to use for computing the kinetic energy. getNumComputations
(self)Get the number of computation steps that have been added. getNumGlobalVariables
(self)Get the number of global variables that have been defined. getNumPerDofVariables
(self)Get the number of perDOF variables that have been defined. getNumTabulatedFunctions
(self)Get the number of tabulated functions that have been defined. getPerDofVariable
(self, index)getPerDofVariable(self, index) > PyObject * getPerDofVariableByName
(self, name)Get the value of a perDOF variable, specified by name. getPerDofVariableName
(self, index)Get the name of a perDOF variable. getRandomNumberSeed
(self)Get the random number seed. getStepSize
(self)Get the size of each time step, in picoseconds. getTabulatedFunction
(self, index)getTabulatedFunction(self, index) > TabulatedFunction getTabulatedFunctionName
(self, index)Get the name of a tabulated function that may appear in expressions. setAlpha
(alpha)Set the value of alpha for the integrator. setConstraintTolerance
(self, tol)Set the distance tolerance within which constraints are maintained, as a fraction of the constrained distance. setE
(E)Set the energy threshold E for the integrator. setGlobalVariable
(self, index, value)Set the value of a global variable. setGlobalVariableByName
(self, name, value)Set the value of a global variable, specified by name. setKineticEnergyExpression
(self, expression)Set the expression to use for computing the kinetic energy. setPerDofVariable
(self, index, values)Set the value of a perDOF variable. setPerDofVariableByName
(self, name, values)Set the value of a perDOF variable, specified by name. setRandomNumberSeed
(self, seed)Set the random number seed. setStepSize
(self, size)Set the size of each time step, in picoseconds. step
(self, steps)Advance a simulation through time by taking a series of time steps. Attributes
BlockEnd
ComputeGlobal
ComputePerDof
ComputeSum
ConstrainPositions
ConstrainVelocities
IfBlockStart
UpdateContextState
WhileBlockStart

getAlpha
()¶ Get the value of alpha for the integrator.

setAlpha
(alpha)¶ Set the value of alpha for the integrator.

getE
()¶ Get the energy threshold E for the integrator.

setE
(E)¶ Set the energy threshold E for the integrator.

getEffectiveEnergy
(energy)¶ Given the actual potential energy of the system, return the value of the effective potential.

__copy__
(self) → Integrator¶

addComputeGlobal
(self, variable, expression) → int¶ Add a step to the integration algorithm that computes a global value.
Parameters:  variable (string) – the global variable to store the computed value into
 expression (string) – a mathematical expression involving only global variables. In each integration step, its value is computed and stored into the specified variable.
Returns: the index of the step that was added
Return type: int

addComputePerDof
(self, variable, expression) → int¶ Add a step to the integration algorithm that computes a perDOF value.
Parameters:  variable (string) – the perDOF variable to store the computed value into
 expression (string) – a mathematical expression involving both global and perDOF variables. In each integration step, its value is computed for every degree of freedom and stored into the specified variable.
Returns: the index of the step that was added
Return type: int

addComputeSum
(self, variable, expression) → int¶ Add a step to the integration algorithm that computes a sum over degrees of freedom.
Parameters:  variable (string) – the global variable to store the computed value into
 expression (string) – a mathematical expression involving both global and perDOF variables. In each integration step, its value is computed for every degree of freedom. Those values are then added together, and the sum is stored in the specified variable.
Returns: the index of the step that was added
Return type: int

addConstrainPositions
(self) → int¶ Add a step to the integration algorithm that updates particle positions so all constraints are satisfied.
Returns: the index of the step that was added Return type: int

addConstrainVelocities
(self) → int¶ Add a step to the integration algorithm that updates particle velocities so the net velocity along all constraints is 0.
Returns: the index of the step that was added Return type: int

addGlobalVariable
(self, name, initialValue) → int¶ Define a new global variable.
Parameters:  name (string) – the name of the variable
 initialValue (double) – the variable will initially be set to this value
Returns: the index of the variable that was added
Return type: int

addPerDofVariable
(self, name, initialValue) → int¶ Define a new perDOF variable.
Parameters:  name (string) – the name of the variable
 initialValue (double) – the variable will initially be set to this value for all degrees of freedom
Returns: the index of the variable that was added
Return type: int

addTabulatedFunction
(self, name, function) → int¶ Add a tabulated function that may appear in expressions.
Parameters:  name (string) – the name of the function as it appears in expressions
 function (TabulatedFunction *) – a TabulatedFunction object defining the function. The TabulatedFunction should have been created on the heap with the “new” operator. The integrator takes over ownership of it, and deletes it when the integrator itself is deleted.
Returns: the index of the function that was added
Return type: int

addUpdateContextState
(self) → int¶ Add a step to the integration algorithm that allows Forces to update the context state.
Returns: the index of the step that was added Return type: int

beginIfBlock
(self, condition) → int¶ Add a step which begins a new “if” block.
Parameters: condition (string) – a mathematical expression involving a comparison operator and global variables. All steps between this one and the end of the block are executed only if the condition is true. Returns: the index of the step that was added Return type: int

beginWhileBlock
(self, condition) → int¶ Add a step which begins a new “while” block.
Parameters: condition (string) – a mathematical expression involving a comparison operator and global variables. All steps between this one and the end of the block are executed repeatedly as long as the condition remains true. Returns: the index of the step that was added Return type: int

endBlock
(self) → int¶ Add a step which marks the end of the most recently begun “if” or “while” block.
Returns: the index of the step that was added Return type: int

getComputationStep
(self, index)¶ Get the details of a computation step that has been added to the integration algorithm.
Parameters: index (int) – the index of the computation step to get Returns:  type (ComputationType) – the type of computation this step performs
 variable (string) – the variable into which this step stores its result. If this step does not store a result in a variable, this will be an empty string.
 expression (string) – the expression this step evaluates. If this step does not evaluate an expression, this will be an empty string.

getConstraintTolerance
(self) → double¶ Get the distance tolerance within which constraints are maintained, as a fraction of the constrained distance.

getGlobalVariable
(self, index) → double¶ Get the current value of a global variable.
Parameters: index (int) – the index of the variable to get Returns: the current value of the variable Return type: double

getGlobalVariableByName
(self, name) → double¶ Get the current value of a global variable, specified by name.
Parameters: name (string) – the name of the variable to get Returns: the current value of the parameter Return type: double

getGlobalVariableName
(self, index) → std::string const &¶ Get the name of a global variable.
Parameters: index (int) – the index of the variable to get Returns: the name of the variable Return type: string

getKineticEnergyExpression
(self) → std::string const &¶ Get the expression to use for computing the kinetic energy. The expression is evaluated for every degree of freedom. Those values are then added together, and the sum is reported as the current kinetic energy.

getNumComputations
(self) → int¶ Get the number of computation steps that have been added.

getNumGlobalVariables
(self) → int¶ Get the number of global variables that have been defined.

getNumPerDofVariables
(self) → int¶ Get the number of perDOF variables that have been defined.

getNumTabulatedFunctions
(self) → int¶ Get the number of tabulated functions that have been defined.

getPerDofVariable
(self, index)¶ getPerDofVariable(self, index) > PyObject *

getPerDofVariableByName
(self, name)¶ Get the value of a perDOF variable, specified by name.
Parameters: name (string) – the name of the variable to get Returns: values – the values of the variable for all degrees of freedom are stored into this Return type: vector< Vec3 >

getPerDofVariableName
(self, index) → std::string const &¶ Get the name of a perDOF variable.
Parameters: index (int) – the index of the variable to get Returns: the name of the variable Return type: string

getRandomNumberSeed
(self) → int¶ Get the random number seed. See setRandomNumberSeed() for details.

getStepSize
(self) → double¶ Get the size of each time step, in picoseconds. If this integrator uses variable time steps, the size of the most recent step is returned.
Returns: the step size, measured in ps Return type: double

getTabulatedFunction
(self, index) → TabulatedFunction¶ getTabulatedFunction(self, index) > TabulatedFunction
Get a reference to a tabulated function that may appear in expressions.
Parameters: index (int) – the index of the function to get Returns: the TabulatedFunction object defining the function Return type: TabulatedFunction

getTabulatedFunctionName
(self, index) → std::string const &¶ Get the name of a tabulated function that may appear in expressions.
Parameters: index (int) – the index of the function to get Returns: the name of the function as it appears in expressions Return type: string

setConstraintTolerance
(self, tol)¶ Set the distance tolerance within which constraints are maintained, as a fraction of the constrained distance.

setGlobalVariable
(self, index, value)¶ Set the value of a global variable.
Parameters:  index (int) – the index of the variable to set
 value (double) – the new value of the variable

setGlobalVariableByName
(self, name, value)¶ Set the value of a global variable, specified by name.
Parameters:  name (string) – the name of the variable to set
 value (double) – the new value of the variable

setKineticEnergyExpression
(self, expression)¶ Set the expression to use for computing the kinetic energy. The expression is evaluated for every degree of freedom. Those values are then added together, and the sum is reported as the current kinetic energy.

setPerDofVariable
(self, index, values)¶ Set the value of a perDOF variable.
Parameters:  index (int) – the index of the variable to set
 values (vector< Vec3 >) – the new values of the variable for all degrees of freedom

setPerDofVariableByName
(self, name, values)¶ Set the value of a perDOF variable, specified by name.
Parameters:  name (string) – the name of the variable to set
 values (vector< Vec3 >) – the new values of the variable for all degrees of freedom

setRandomNumberSeed
(self, seed)¶ Set the random number seed. The precise meaning of this parameter is undefined, and is left up to each Platform to interpret in an appropriate way. It is guaranteed that if two simulations are run with different random number seeds, the sequence of random numbers will be different. On the other hand, no guarantees are made about the behavior of simulations that use the same seed. In particular, Platforms are permitted to use nondeterministic algorithms which produce different results on successive runs, even if those runs were initialized identically.
If seed is set to 0 (which is the default value assigned), a unique seed is chosen when a Context is created from this Force. This is done to ensure that each Context receives unique random seeds without you needing to set them explicitly.

setStepSize
(self, size)¶ Set the size of each time step, in picoseconds. If this integrator uses variable time steps, the effect of calling this method is undefined, and it may simply be ignored.
Parameters: size (double) – the step size, measured in ps

step
(self, steps)¶ Advance a simulation through time by taking a series of time steps.
Parameters: steps (int) – the number of time steps to take
