15 EnergyFunction::~EnergyFunction() {
18 EnergyFunction::EnergyFunction(
const EnergyFunction & right)
19 :AbsFunction(right), solver(right.solver)
23 double EnergyFunction::operator() (
double t)
const {
25 const unsigned int DIM=phaseSpace.
dim();
27 for (
unsigned int i=0;i<phaseSpace.
dim();i++) {
28 arg[i] =solver.equationOf(phaseSpace.
coordinates()[i])(t);
29 arg[i+DIM]=solver.equationOf(phaseSpace.
momenta()[i])(t);
#define FUNCTION_OBJECT_IMP(classname)
const Component & momenta() const
const Component & coordinates() const