low-order leapfrog method provided for treecode, which has discontinuous potential/forces and is poorly handled by Runge-Kutta integrator. More...

#include <odeint.h>

Inheritance diagram for smile::COdeIntegratorLeapfrog:

Collaboration diagram for smile::COdeIntegratorLeapfrog:

[legend]

Public Member Functions
	COdeIntegratorLeapfrog (COdeSystem *_odeSystem, double _accTreecode, bool _symmetrizeTimestep)

virtual void	integrateToTime (double timeEnd)
	perform integration

virtual double	getInterpolatedSolution (unsigned int c, double t) const
	return interpolated value of c-th variable at time t, where t must lie within current timestep interval; called from various runtime functions that need to obtain x(t) at arbitrary times

Public Member Functions inherited from smile::CBasicOdeIntegrator
	CBasicOdeIntegrator (COdeSystem *_odeSystem)

double	getTime () const
	return the time to which the integration has proceeded so far

Static Public Member Functions
static const char *	myName ()

Private Attributes
double	accTreecode
	accuracy parameter (timestep is accTreecode*timeToCollision)

bool	symmetrizeTimestep
	whether to use Hut et al.(1995) timestep symmetrization algorithm

OdeStateType	stateCurr

OdeStateType	statePrev
	ODE system state at the beginning and the end of timestep.

Additional Inherited Members
Public Types inherited from smile::CBasicOdeIntegrator
enum	STEPPERKIND { SK_DEFAULT, SK_LEAPFROG_NB, SK_DOP853, SK_IAS15, SK_HERMITE, SK_ODEINT_CK5, SK_ODEINT_DP5, SK_ODEINT_BS3, SK_ODEINT_BS, SK_ODEINT_RK4, SK_ODEINT_SYMPL4 }
	list of various ODE stepper types

Protected Attributes inherited from smile::CBasicOdeIntegrator
COdeSystem *	odeSystem

double	timePrev

double	timeCurr

Detailed Description

low-order leapfrog method provided for treecode, which has discontinuous potential/forces and is poorly handled by Runge-Kutta integrator.

It does not compute Lyapunov exponent either.

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

Public Member Functions