Classical::RungeKuttaSolver Class Reference

#include <RungeKuttaClassicalSolver.hh>

Inheritance diagram for Classical::RungeKuttaSolver:

Classes
class	Clockwork

Public Member Functions
	RungeKuttaSolver (Genfun::GENFUNCTION H, const PhaseSpace &phaseSpace, const Genfun::RKIntegrator::RKStepper *stepper=NULL)
	~RungeKuttaSolver ()
virtual Genfun::GENFUNCTION	equationOf (const Genfun::Variable &v) const
virtual const PhaseSpace &	phaseSpace () const
virtual Genfun::GENFUNCTION	hamiltonian () const
virtual Genfun::GENFUNCTION	energy () const
virtual Genfun::Parameter *	takeQ0 (unsigned int index)
virtual Genfun::Parameter *	takeP0 (unsigned int index)
virtual Genfun::Parameter *	createControlParameter (const std::string &variableName="anon", double defStartingValue=0.0, double startingValueMin=0.0, double startingValueMax=0.0) const
Public Member Functions inherited from Classical::Solver
	Solver ()
virtual	~Solver ()
virtual Genfun::GENFUNCTION	equationOf (const Genfun::Variable &v) const =0
virtual const PhaseSpace &	phaseSpace () const =0
virtual Genfun::GENFUNCTION	hamiltonian () const =0
virtual Genfun::GENFUNCTION	energy () const =0
virtual Genfun::Parameter *	takeQ0 (unsigned int index)=0
virtual Genfun::Parameter *	takeP0 (unsigned int index)=0
virtual Genfun::Parameter *	createControlParameter (const std::string &variableName="anon", double defStartingValue=0.0, double startingValueMin=0.0, double startingValueMax=0.0) const =0

Detailed Description

Definition at line 21 of file RungeKuttaClassicalSolver.hh.

Constructor & Destructor Documentation

RungeKuttaSolver()

Classical::RungeKuttaSolver::RungeKuttaSolver	(	Genfun::GENFUNCTION	H,
		const PhaseSpace &	phaseSpace,
		const Genfun::RKIntegrator::RKStepper *	stepper = NULL
	)

Definition at line 25 of file RungeKuttaClassicalSolver.cc.

                                                                                                                                    :c(new Clockwork(gH,mphaseSpace)){
    c->integrator=new Genfun::RKIntegrator(stepper);
    //
    // Dimension (of coords, or phase space)
    //
    const unsigned int DIM=c->phaseSpace.dim();
    //
    // Equations of motion via hamilton's equations: 
    //
    const Classical::PhaseSpace::Component & X=c->phaseSpace.coordinates();
    const Classical::PhaseSpace::Component & P=c->phaseSpace.momenta();
    
    for (unsigned int i=0;i<DIM;i++) {
      Genfun::GENFUNCTION DXDT  =  c->H.partial(P[i].index());
      c->startingQ.push_back(c->integrator->addDiffEquation(&DXDT,"X",c->phaseSpace.startValue(X[i])));
    }
    for (unsigned int i=0;i<DIM;i++) {
      Genfun::GENFUNCTION DPDT  = -c->H.partial(X[i].index());
      c->startingP.push_back(c->integrator->addDiffEquation(&DPDT,"P",c->phaseSpace.startValue(P[i])));
    }
    c->energy=NULL;
    
  }

~RungeKuttaSolver()

Classical::RungeKuttaSolver::~RungeKuttaSolver

(

)

Definition at line 48 of file RungeKuttaClassicalSolver.cc.

                                     {
    delete c->integrator;
    delete c->energy;
    delete c;
  }

Member Function Documentation

createControlParameter()

Genfun::Parameter * Classical::RungeKuttaSolver::createControlParameter ( const std::string &  variableName = "anon", double  defStartingValue = 0.0, double  startingValueMin = 0.0, double  startingValueMax = 0.0  ) const

virtual

Implements Classical::Solver.

Definition at line 68 of file RungeKuttaClassicalSolver.cc.

                                                                     {
    return c->integrator->createControlParameter(variableName, defStartingValue, startingValueMin, startingValueMax) ;
  }

energy()

Genfun::GENFUNCTION Classical::RungeKuttaSolver::energy ( ) const

virtual

Implements Classical::Solver.

Definition at line 63 of file RungeKuttaClassicalSolver.cc.

                                                   {
    if (!c->energy) c->energy=new Genfun::EnergyFunction(*this);
    return *c->energy;
  }

equationOf()

Genfun::GENFUNCTION Classical::RungeKuttaSolver::equationOf ( const Genfun::Variable &  v ) const

virtual

Implements Classical::Solver.

Definition at line 54 of file RungeKuttaClassicalSolver.cc.

                                                                               {
    return *c->integrator->getFunction(v.index());
  }

hamiltonian()

Genfun::GENFUNCTION Classical::RungeKuttaSolver::hamiltonian ( ) const

virtual

Implements Classical::Solver.

Definition at line 57 of file RungeKuttaClassicalSolver.cc.

                                                        {
    return c->H;
  }

phaseSpace()

const Classical::PhaseSpace & Classical::RungeKuttaSolver::phaseSpace ( ) const

virtual

Implements Classical::Solver.

Definition at line 60 of file RungeKuttaClassicalSolver.cc.

                                                                 {
    return c->phaseSpace;
  }

takeP0()

Genfun::Parameter * Classical::RungeKuttaSolver::takeP0 ( unsigned int  index )

virtual

Implements Classical::Solver.

Definition at line 78 of file RungeKuttaClassicalSolver.cc.

                                                              {
    return c->startingP[index];
  }

takeQ0()

Genfun::Parameter * Classical::RungeKuttaSolver::takeQ0 ( unsigned int  index )

virtual

Implements Classical::Solver.

Definition at line 75 of file RungeKuttaClassicalSolver.cc.

                                                              {
    return c->startingQ[index];
  }

---

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

• RungeKuttaClassicalSolver.hh
• RungeKuttaClassicalSolver.cc