Loading...
Searching...
No Matches
BGSDObjectiveFunction Class Reference
Inheritance diagram for BGSDObjectiveFunction:

Public Member Functions

 BGSDObjectiveFunction (Matter &matterRef, double reactantEnergyPassed, double bgsdAlphaPassed, const Parameters &parametersPassed)
 ~BGSDObjectiveFunction ()=default
double getEnergy ()
VectorXd getGradient (bool fdstep=false)
double getGradientnorm ()
void setPositions (const VectorXd &x)
VectorXd getPositions ()
int degreesOfFreedom ()
bool isConverged ()
bool isConvergedH ()
bool isConvergedV ()
bool isConvergedIP ()
double getConvergence ()
double getConvergenceH ()
double getConvergenceV ()
VectorXd difference (const VectorXd &a, const VectorXd &b)
Public Member Functions inherited from eonc::ObjectiveFunction
 ObjectiveFunction (const Parameters &paramsPassed)
virtual ~ObjectiveFunction ()

Private Attributes

Mattermatter
double reactantEnergy
double bgsdAlpha

Additional Inherited Members

Protected Attributes inherited from eonc::ObjectiveFunction
const Parametersparams

Detailed Description

Definition at line 26 of file BiasedGradientSquaredDescent.cpp.

Constructor & Destructor Documentation

◆ BGSDObjectiveFunction()

BGSDObjectiveFunction::BGSDObjectiveFunction ( Matter & matterRef,
double reactantEnergyPassed,
double bgsdAlphaPassed,
const Parameters & parametersPassed )
inline

Definition at line 30 of file BiasedGradientSquaredDescent.cpp.

33 : ObjectiveFunction(parametersPassed),
34 matter{matterRef} {
35 bgsdAlpha = bgsdAlphaPassed;
36 reactantEnergy = reactantEnergyPassed;
37 }
ObjectiveFunction(const Parameters &paramsPassed)

◆ ~BGSDObjectiveFunction()

BGSDObjectiveFunction::~BGSDObjectiveFunction ( )
default

Member Function Documentation

◆ degreesOfFreedom()

int BGSDObjectiveFunction::degreesOfFreedom ( )
inlinevirtual

Implements eonc::ObjectiveFunction.

Definition at line 79 of file BiasedGradientSquaredDescent.cpp.

79{ return 3 * matter.numberOfFreeAtoms(); }

◆ difference()

VectorXd BGSDObjectiveFunction::difference ( const VectorXd & a,
const VectorXd & b )
inlinevirtual

Implements eonc::ObjectiveFunction.

Definition at line 94 of file BiasedGradientSquaredDescent.cpp.

94 {
95 return matter.pbcV(a - b);
96 }

◆ getConvergence()

double BGSDObjectiveFunction::getConvergence ( )
inlinevirtual

Implements eonc::ObjectiveFunction.

Definition at line 91 of file BiasedGradientSquaredDescent.cpp.

91{ return getEnergy() && getGradient().norm(); }
VectorXd getGradient(bool fdstep=false)

◆ getConvergenceH()

double BGSDObjectiveFunction::getConvergenceH ( )
inline

Definition at line 92 of file BiasedGradientSquaredDescent.cpp.

92{ return getGradient().norm(); }

◆ getConvergenceV()

double BGSDObjectiveFunction::getConvergenceV ( )
inline

Definition at line 93 of file BiasedGradientSquaredDescent.cpp.

93{ return getEnergy(); }

◆ getEnergy()

double BGSDObjectiveFunction::getEnergy ( )
inlinevirtual

Implements eonc::ObjectiveFunction.

Definition at line 41 of file BiasedGradientSquaredDescent.cpp.

41 {
42 VectorXd Vforce = matter.getForcesFreeV();
43 double Henergy = 0.5 * Vforce.dot(Vforce) +
44 0.5 * bgsdAlpha *
45 (matter.getPotentialEnergy() -
46 (reactantEnergy + params.bgsd_options.beta)) *
47 (matter.getPotentialEnergy() -
48 (reactantEnergy + params.bgsd_options.beta));
49 return Henergy;
50 }
const Parameters & params

◆ getGradient()

VectorXd BGSDObjectiveFunction::getGradient ( bool fdstep = false)
inlinevirtual

Implements eonc::ObjectiveFunction.

Definition at line 52 of file BiasedGradientSquaredDescent.cpp.

52 {
53 VectorXd Vforce = matter.getForcesFreeV();
54 double magVforce = Vforce.norm();
55 VectorXd normVforce = Vforce / magVforce;
56 VectorXd Vpositions = matter.getPositionsFreeV();
57 matter.setPositionsFreeV(
58 matter.getPositionsFreeV() -
59 normVforce * params.bgsd_options.gradient_finite_difference);
60 VectorXd Vforcenew = matter.getForcesFreeV();
61 matter.setPositionsFreeV(Vpositions);
62 VectorXd Hforce = magVforce * (Vforcenew - Vforce) /
63 params.bgsd_options.gradient_finite_difference +
64 bgsdAlpha *
65 (matter.getPotentialEnergy() -
66 (reactantEnergy + params.bgsd_options.beta)) *
67 Vforce;
68 return -Hforce;
69 }

◆ getGradientnorm()

double BGSDObjectiveFunction::getGradientnorm ( )
inline

Definition at line 71 of file BiasedGradientSquaredDescent.cpp.

71 {
72 VectorXd Hforce = getGradient();
73 double Hnorm = Hforce.norm();
74 return Hnorm;
75 }

◆ getPositions()

VectorXd BGSDObjectiveFunction::getPositions ( )
inlinevirtual

Implements eonc::ObjectiveFunction.

Definition at line 78 of file BiasedGradientSquaredDescent.cpp.

78{ return matter.getPositionsFreeV(); }

◆ isConverged()

bool BGSDObjectiveFunction::isConverged ( )
inlinevirtual

◆ isConvergedH()

bool BGSDObjectiveFunction::isConvergedH ( )
inline

Definition at line 81 of file BiasedGradientSquaredDescent.cpp.

81 {
82 return getConvergenceH() < params.bgsd_options.h_force_convergence;
83 }

◆ isConvergedIP()

bool BGSDObjectiveFunction::isConvergedIP ( )
inline

Definition at line 87 of file BiasedGradientSquaredDescent.cpp.

87 {
88 return getConvergenceH() < params.bgsd_options.grad2force_convergence;
89 }

◆ isConvergedV()

bool BGSDObjectiveFunction::isConvergedV ( )
inline

Definition at line 84 of file BiasedGradientSquaredDescent.cpp.

84 {
85 return getConvergenceV() < params.bgsd_options.grad2energy_convergence;
86 }

◆ setPositions()

void BGSDObjectiveFunction::setPositions ( const VectorXd & x)
inlinevirtual

Implements eonc::ObjectiveFunction.

Definition at line 77 of file BiasedGradientSquaredDescent.cpp.

77{ matter.setPositionsFreeV(x); }

Member Data Documentation

◆ bgsdAlpha

double BGSDObjectiveFunction::bgsdAlpha
private

Definition at line 100 of file BiasedGradientSquaredDescent.cpp.

◆ matter

Matter& BGSDObjectiveFunction::matter
private

Definition at line 27 of file BiasedGradientSquaredDescent.cpp.

◆ reactantEnergy

double BGSDObjectiveFunction::reactantEnergy
private

Definition at line 99 of file BiasedGradientSquaredDescent.cpp.


The documentation for this class was generated from the following file: