class PerturbDistanceAngle : public feasst::PerturbDistance
Put first site in selection, i, in a truncated conical shell that is formed by the intersection of a spherical shell about the first site in the anchor, and a conical shell with axis of symmetry along the vector connecting the first anchor site to the second anchor site. For angle potentials, the equilibrium angle and spring constant are as described in Random::bond_angle Currently implemented for harmonic bonds (exponent: 2), but could add an optional exponent model parameter to generalize this in the future.
void precompute(TrialSelect *select, System *system)
Same as PerturbDistance, but also obtain the equilibrium angle and spring constant. If 2D, angles are positive when clockwise. Thus, when 2D and reverse (e.g., kji instead of ijk), angle = 2pi - angle.
double angle_type() const
Return the angle.
double old_angle_energy(const System &system, const TrialSelect *select)
Return the current energy of the existing, old angle.
double random_angle_radians(const System &system, const TrialSelect *select, Random *random, double *bond_energy)
Return the randomly selected angle from the potential.
- void precompute(TrialSelect *select, System *system)