33 for (
unsigned int i=0;i<DIM;i++) {
37 for (
unsigned int i=0;i<DIM;i++) {
60 if (!c->
energy) c->
energy=
new Genfun::EnergyFunction(*
this);
65 double defStartingValue,
66 double startingValueMin,
67 double startingValueMax)
const {
const Component & momenta() const
const Component & coordinates() const
double startValue(const Genfun::Variable &component) const
Clockwork(Genfun::GENFUNCTION gH, const PhaseSpace &mphaseSpace)
Genfun::RKIntegrator * integrator
std::vector< Genfun::Parameter * > startingQ
const Classical::PhaseSpace & phaseSpace
Genfun::EnergyFunction * energy
std::vector< Genfun::Parameter * > startingP
virtual Genfun::GENFUNCTION hamiltonian() const
virtual Genfun::Parameter * createControlParameter(const std::string &variableName="anon", double defStartingValue=0.0, double startingValueMin=0.0, double startingValueMax=0.0) const
virtual Genfun::GENFUNCTION equationOf(const Genfun::Variable &v) const
virtual Genfun::GENFUNCTION energy() const
virtual Genfun::Parameter * takeP0(unsigned int index)
virtual Genfun::Parameter * takeQ0(unsigned int index)
RungeKuttaSolver(Genfun::GENFUNCTION H, const PhaseSpace &phaseSpace, const Genfun::RKIntegrator::RKStepper *stepper=NULL)
virtual const PhaseSpace & phaseSpace() const
virtual Derivative partial(unsigned int) const
Parameter * addDiffEquation(const AbsFunction *diffEquation, const std::string &variableName="anon", double defStartingValue=0.0, double startingValueMin=0.0, double startingValueMax=0.0)
Parameter * createControlParameter(const std::string &variableName="anon", double defStartingValue=0.0, double startingValueMin=0.0, double startingValueMax=0.0)
const RKFunction * getFunction(unsigned int i) const
unsigned int index() const