CAPD DynSys Library
5.2.0
|
#include <capd/dynsys/CnOdeSolver.h>
Public Types | |
typedef MapT | VectorFieldType |
typedef StepControlT | StepControlPolicy |
typedef EnclosurePolicyT | EnclosurePolicy |
typedef BasicCnOdeSolver< MapT, StepControlPolicy, CurveT > | BaseTaylor |
typedef VectorFieldType::ScalarType | ScalarType |
typedef VectorFieldType::VectorType | VectorType |
typedef VectorFieldType::MatrixType | MatrixType |
typedef VectorFieldType::HessianType | HessianType |
typedef VectorFieldType::FunctionType | FunctionType |
typedef BaseTaylor::JetType | JetType |
typedef MatrixType::size_type | size_type |
typedef diffAlgebra::C1TimeJet< MatrixType > | C1TimeJetType |
typedef diffAlgebra::C2TimeJet< MatrixType > | C2TimeJetType |
typedef diffAlgebra::CnTimeJet< MatrixType, 0 > | CnTimeJetType |
typedef CurveT | SolutionCurve |
typedef capd::vectalg::Norm< VectorType, MatrixType > | NormType |
typedef capd::dynsys::IEncFoundStepControl | StepControlType |
typedef JetType::Multipointer | Multipointer |
typedef JetType::Multiindex | Multiindex |
typedef TypeTraits< ScalarType >::Real | Real |
typedef __difference_type | difference_type |
Public Member Functions | |
template<typename SetType > | |
void | operator() (SetType &set) |
This operator computes image of the set (in given representation) using set.move function, see capd/dynsys/Move.h for details This template together with SetTraits prevent usage of various types of jets with incompatible solvers. The user will get an exception at runtime with clear message instead of unreadable compiler error. In this case a specialization CnSetMove is used meaning that this solver can integrate C^0, C^1 and C^2 sets only. Moreover, it cannot integrate nonrigorous jets (for user safety). More... | |
template<typename SetType > | |
void | operator() (SetType &set, SetType &result) |
CnOdeSolver (VectorFieldType &vectorField, size_type order, const StepControlT &stepControl=StepControlT()) | |
VectorType | Phi (const ScalarType &t, const VectorType &iv) |
MatrixType | JacPhi (const ScalarType &t, const VectorType &iv) |
VectorType | enclosure (const ScalarType &t, const VectorType &x) |
VectorType | Remainder (const ScalarType &t, const VectorType &iv, VectorType &o_enc) |
void | encloseC0Map (const ScalarType &t, const VectorType &x0, const VectorType &x, VectorType &o_phi, VectorType &o_rem, VectorType &o_enc, MatrixType &o_jacPhi) |
void | JacRemainder (const ScalarType &t, const VectorType &vecEnclosure, const MatrixType &jacEnclosure, VectorType &Remainder, MatrixType &jacRemainder) |
void | encloseC1Map (const ScalarType &t, const VectorType &x0, const VectorType &x, VectorType &o_phi, VectorType &o_rem, VectorType &o_enc, MatrixType &o_jacPhi, MatrixType &o_jacRem, MatrixType &o_jacEnc) |
void | encloseC2Map (const ScalarType &t, const VectorType &x, const VectorType &xx, VectorType &o_phi, VectorType &o_rem, VectorType &o_enc, MatrixType &o_jacPhi, MatrixType &o_jacRem, MatrixType &o_jacEnc, HessianType &o_hessianPhi, HessianType &o_hessianRem, HessianType &o_hessianEnc) |
void | c2Remainder (const VectorType &Enc, const MatrixType &jacEnc, const HessianType &hessianEnc, VectorType &o_Rem, MatrixType &o_jacRem, HessianType &o_hessianRem) |
template<class JetT > | |
VectorType | encloseCnMap (const ScalarType &t, const VectorType &x, const VectorType &xx, JetT &phi, JetT &rem, JetT &enc) |
template<class JetT > | |
VectorType | cnEnclosure (const ScalarType &t, const VectorType &x, JetT &result) |
template<class JetT > | |
void | cnRemainder (const JetT &enc, JetT &result) |
ScalarType | getStep () const |
ScalarType | getCoeffNorm (size_type i, size_type degree) const |
void | computeRemainderCoefficients (const VectorType &x) |
void | computeRemainderCoefficients (const VectorType &x, const MatrixType &M) |
void | computeRemainderCoefficients (ScalarType t, const VectorType &x) |
void | computeRemainderCoefficients (ScalarType t, const VectorType &x, const MatrixType &M) |
virtual void | computeRemainder (ScalarType t, const VectorType &xx, VectorType &o_enc, VectorType &o_rem) |
virtual void | computeRemainder (ScalarType t, const VectorType &xx, C1TimeJetType &o_enc, C1TimeJetType &o_rem) |
virtual void | computeRemainder (ScalarType t, const VectorType &xx, C2TimeJetType &o_enc, C2TimeJetType &o_rem) |
template<class JetT > | |
void | computeRemainder (ScalarType t, const VectorType &xx, JetT &o_enc, JetT &o_rem) |
void | sumTaylorSeries (VectorType &o_phi, MatrixType &o_jacPhi) |
void | sumTaylorSeries (VectorType &o_phi, MatrixType &o_jacPhi, HessianType &o_hessianPhi) |
template<class JetT > | |
void | sumTaylorSeries (JetT &o_phi) |
void | computeTimeStep (const ScalarType &t, const VectorType &x) |
const VectorFieldType & | getVectorField () const |
Returns vector field. More... | |
VectorFieldType & | getVectorField () |
void | setOrder (size_type order) |
Sets the order of the Taylor method. More... | |
void | setStep (const ScalarType &newStep) |
Sets time step and turns off step control. More... | |
void | setCurrentTime (const ScalarType &a_time) const |
const ScalarType & | getCurrentTime () const |
const SolutionCurve & | getCurve () |
VectorType | encloseCnMap (const ScalarType &, const VectorType &, const VectorType &, JetT &, JetT &, JetT &) |
virtual void | encloseC2Map (const ScalarType &t, const VectorType &x, const VectorType &xx, VectorType &o_phi, VectorType &o_rem, VectorType &o_enc, MatrixType &o_jacPhi, MatrixType &o_jacRem, MatrixType &o_jacEnc, HessianType &o_hessianPhi, HessianType &o_hessianRem, HessianType &o_hessianEnc)=0 |
virtual void | encloseC1Map (const ScalarType &t, const VectorType &x, const VectorType &xx, VectorType &o_phi, VectorType &o_rem, VectorType &o_enc, MapT::MatrixType &o_jacPhi, MapT::MatrixType &o_jacRem, MapT::MatrixType &o_jacEnc)=0 |
virtual VectorType | Phi (const ScalarType &t, const VectorType &iv)=0 |
Computes value of function (usually numerical scheme for an ODE) at time t and point iv. More... | |
virtual MatrixType | JacPhi (const ScalarType &t, const VectorType &iv)=0 |
Computes derivative of function (usually numerical scheme for an ODE) at time t and point iv. More... | |
virtual VectorType | Remainder (const ScalarType &t, const VectorType &iv, VectorType &out_enc)=0 |
Computes and returns bound for local error of a function (for instance if Phi is finite Taylor series of a map then this method computes bound for Lagrange remainder). If DynSys is an ODE, then out_enc contains enclosure of trajectories over the time step. If the function cannot validate existence of solutions to ODE over the time step, out_enc might be in an inconsistent state. More... | |
virtual ScalarType | Lipschitz (const ScalarType &t, const VectorType &iv, NormType &n) |
this function returns Lipschitz contants for maps. It should be overriden in classes implementing flows. More... | |
virtual VectorType | enclosure (const ScalarType &t, const VectorType &x)=0 |
Used for ODEs. It verifies the existence of solutions to IVP at time t and set of initial conditions x over the time step. If the function succeeds, a rigorous bound for the trajectories is returned. Otherwise, an exception is thrown. More... | |
virtual void | encloseC0Map (const ScalarType &t, const VectorType &x, const VectorType &xx, VectorType &o_phi, VectorType &o_rem, VectorType &o_enc, MatrixType &o_jacPhi)=0 |
For given set xx, time t and a point x from the set xx It simultaneously computes and returns enclosures for: More... | |
virtual ScalarType | getStep () const |
Returns time step of the dynamical system. By default it returns one - time step for discrete DS (maps). Shall be overridden in classes that implement numerical schemes for ODEs. More... | |
VectorType | operator() (VectorType) |
Computes image of vector v after one time step. More... | |
VectorType | operator() (ScalarType &t, const VectorType &) |
Computes image of vector v after one time step. The argument t is updated in this procedure. More... | |
VectorType | operator() (VectorType, MatrixType &o_resultDerivative) |
Computes image of vector v and derivatives of the flow with respect to init condition (v,identity). Version for autonomous systems. More... | |
VectorType | operator() (ScalarType &t, const VectorType &, MatrixType &o_resultDerivative) |
Computes image of vector v and derivatives of the flow with respect to init condition (v,identity). Version for nonautonomous systems. The argument t is updated in this procedure. More... | |
VectorType | operator() (VectorType, const MatrixType &derivative, MatrixType &o_resultDerivative) |
Computes image of vector v and derivatives of a flow with respect to init condition (v, derivative) More... | |
VectorType | operator() (ScalarType &t, const VectorType &v, const MatrixType &derivative, MatrixType &o_resultDerivative) |
Computes image of vector v and derivatives of a flow with respect to init condition (v, derivative). The argument t is updated in this procedure. More... | |
void | operator() (JetT &jet) |
This operator computes image of the set (in given representation) using set.move function, see capd/dynsys/Move.h for details This template together with SetTraits prevent usage of various types of jets with incompatible solvers. The user will get an exception at runtime with clear message instead of unreadable compiler error. In this case a specialization CnJetMove is used meaning that this solver can integrate any type of jets. More... | |
VectorType | operator() (VectorType, MatrixType &, HessianType &) |
Computes next point on the trajectory, first and second order derivatives with respect to initial conditions. Initial conditions for variational equations are Id and zero, respectively. More... | |
VectorType | operator() (ScalarType &t, const VectorType &, MatrixType &, HessianType &) |
Computes next point on the trajectory of a nonautonomous system, first and second order derivatives with respect to initial conditions. Initial conditions for variational equations are Id and zero, respectively. More... | |
VectorType | operator() (VectorType, const MatrixType &V, const HessianType &H, MatrixType &, HessianType &) |
The routine computes next point, derivatives and second order derivatives of a flow. Initial conditions for variational equations are V and H, respectively. More... | |
VectorType | operator() (ScalarType &t, const VectorType &x, const MatrixType &V, const HessianType &H, MatrixType &, HessianType &) |
The routine computes next point, derivatives and second order derivatives of a nonautonomous flow. Initial conditions for variational equations are V and H, respectively. More... | |
VectorType | operator() (ScalarType &t, JetType &) |
Initial conditions for the trajectory and for variational equations up to given degree are given in argument in_out_coeffs. The full result is stored in in_out_coeffs and also the new point on the trajectory is returned. Note: CnCoeffType is a data structure that stores current time. More... | |
VectorType | operator() (const ScalarType &h) const |
void | computeCoefficientsAtCenter (const VectorType &x, size_type order) |
void | computeCoefficientsAtCenter (ScalarType t, const VectorType &x, size_type order) |
void | computeCoefficients (const VectorType &x, size_type order) |
void | computeCoefficients (ScalarType t, const VectorType &x, size_type order) |
void | computeCoefficients (const VectorType &x, const MatrixType &M, size_type order) |
void | computeCoefficients (ScalarType t, const VectorType &x, const MatrixType &M, size_type order) |
void | setOrder (size_type order) |
Sets the order of the Taylor method. More... | |
ScalarType | getStep () const |
Returns current time step. More... | |
void | setStep (const ScalarType &newStep) |
Sets time step and turns off step control. More... | |
size_type | degree () const |
ScalarType | getCoeffNorm (size_type i, size_type degree) const |
VectorType | enclosure (const ScalarType &, const VectorType &) |
TODO. More... | |
void | adjustTimeStep (const ScalarType &newStep) |
sets time step but does not change step control settings (compare setStep) More... | |
void | setMask (Iterator b, Iterator e) |
The iterator range [b,e) should contain a range of Multiinideces the user requires to compute. More... | |
void | addMultiindexToMask (const capd::vectalg::Multiindex &mi) |
Adds new multiindex (along with dependencies) to the existing mask. More... | |
void | resetMask () |
Resets the mask of derivatives. More... | |
void | turnOnStepControl () |
void | turnOffStepControl () |
void | onOffStepControl (bool _onOffStepControl) |
const StepControlType & | getStepControl () const |
void | setStepControl (const StepControlType &stepControl) |
bool | isStepChangeAllowed () const |
void | setAbsoluteTolerance (TolScalarType tol) |
void | setRelativeTolerance (TolScalarType tol) |
TolScalarType | getAbsoluteTolerance () const |
TolScalarType | getRelativeTolerance () const |
MapT::ScalarType | getMaxStep () const |
void | setMaxStep (MapT::ScalarType maxStep) |
HessianType | hessian (const ScalarType &h) const |
JetType | jet (const ScalarType &h) const |
void | eval (ScalarType h, JetType &v) const |
virtual void | eval (ScalarType, JetType &) const |
VectorType | timeDerivative (const ScalarType &h) const |
MatrixType | derivative (const ScalarType &h) const |
MatrixType | operator[] (const ScalarType &h) const |
size_type | getOrder () const |
Returns the order of Taylor interpolation. More... | |
size_type | getAllocatedOrder () const |
Returns maximal allocated order - used to avoid memory reallocation. More... | |
size_type | dimension () const |
Returns the dimension in which the parametric curve is embedded. More... | |
void | setDegree (size_type degree) |
Sets new maximal degree of the jet and reallocates memory. More... | |
void | clearCoefficients () |
sets all coefficients to zero More... | |
const ScalarType & | centerCoefficient (size_type i, size_type k) const |
ScalarType & | centerCoefficient (size_type i, size_type k) |
const ScalarType & | coefficient (size_type i, size_type k) const |
const ScalarType & | coefficient (size_type i, size_type j, size_type k) const |
const ScalarType & | coefficient (size_type i, size_type j, size_type c, size_type k) const |
const ScalarType & | coefficient (size_type i, size_type j, size_type c, size_type s, size_type k) const |
ScalarType & | coefficient (size_type i, size_type k) |
ScalarType & | coefficient (size_type i, size_type j, size_type k) |
ScalarType & | coefficient (size_type i, size_type j, size_type c, size_type k) |
ScalarType & | coefficient (size_type i, size_type j, size_type c, size_type s, size_type k) |
const ScalarType & | remainderCoefficient (size_type i, size_type j) const |
const ScalarType & | remainderCoefficient (size_type i, size_type j, size_type k) const |
const ScalarType & | remainderCoefficient (size_type i, size_type j, size_type c, size_type k) const |
const ScalarType & | remainderCoefficient (size_type i, size_type j, size_type c, size_type s, size_type k) const |
ScalarType & | remainderCoefficient (size_type i, size_type j) |
ScalarType & | remainderCoefficient (size_type i, size_type j, size_type k) |
ScalarType & | remainderCoefficient (size_type i, size_type j, size_type c, size_type k) |
ScalarType & | remainderCoefficient (size_type i, size_type j, size_type c, size_type s, size_type k) |
void | allocate (size_type dimension, size_type degree) |
void | deallocate () |
void | copyData (const BasicCnCurve &c) |
const VectorType * | getCoefficientsAtCenter () const |
VectorType * | getCoefficientsAtCenter () |
const JetType * | getCoefficients () const |
JetType * | getCoefficients () |
const JetType * | getRemainderCoefficients () const |
JetType * | getRemainderCoefficients () |
void | setInitMatrix (const MatrixType &) |
void | setInitHessian (const HessianType &) |
void | setInitJet (const JetType &) |
std::runtime_error | domainErrorMessage (std::string msg, ScalarType h, Real left, Real right) const |
virtual void | setDomain (Real left, Real right) |
virtual Real | getLeftDomain () const |
virtual Real | getRightDomain () const |
Static Public Member Functions | |
static TolScalarType | getEffectiveTolerance (Solver &solver, const SetType &s) |
Public Attributes | |
VectorType * | m_coefficientsAtCenter |
JetType * | m_coefficients |
JetType * | m_remainderCoefficients |
size_type | m_order |
size_type | m_allocatedOrder |
typedef BasicCnOdeSolver<MapT,StepControlPolicy,CurveT> capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::BaseTaylor |
typedef diffAlgebra::C1TimeJet<MatrixType> capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::C1TimeJetType |
typedef diffAlgebra::C2TimeJet<MatrixType> capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::C2TimeJetType |
typedef diffAlgebra::CnTimeJet<MatrixType,0> capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::CnTimeJetType |
|
inherited |
typedef EnclosurePolicyT capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::EnclosurePolicy |
typedef VectorFieldType::FunctionType capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::FunctionType |
typedef VectorFieldType::HessianType capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::HessianType |
typedef BaseTaylor::JetType capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::JetType |
typedef VectorFieldType::MatrixType capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::MatrixType |
|
inherited |
|
inherited |
|
inherited |
|
inherited |
typedef VectorFieldType::ScalarType capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::ScalarType |
typedef MatrixType::size_type capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::size_type |
typedef CurveT capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::SolutionCurve |
typedef StepControlT capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::StepControlPolicy |
|
inherited |
typedef MapT capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::VectorFieldType |
typedef VectorFieldType::VectorType capd::dynsys::CnOdeSolver< MapT, StepControlT, EnclosurePolicyT, CurveT >::VectorType |
|
inlineprotected |
|
inlineinherited |
Adds new multiindex (along with dependencies) to the existing mask.
mi | multiindex to be added to the mask |
|
inlineinherited |
sets time step but does not change step control settings (compare setStep)
|
inherited |
|
inlineinherited |
|
inlineinherited |
|
inherited |
sets all coefficients to zero
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inline |
|
inlineprotectedinherited |
|
inherited |
|
inherited |
|
inlineinherited |
|
virtualinherited |
|
inlineinherited |
Returns the dimension in which the parametric curve is embedded.
|
inlineinherited |
|
pure virtualinherited |
For given set xx, time t and a point x from the set xx It simultaneously computes and returns enclosures for:
|
pure virtualinherited |
|
pure virtualinherited |
|
inlineinherited |
t | - time in ODE | |
x | - a centre of current set | |
xx | - the set | |
[out] | out_phi | - should contain coefficients of solution jet at (t,xx) |
[out] | out_rem | - should contain remainder coefficients computed at enclosure |
[out] | out_enc | - should contain enclosure for the main part and variational equations |
|
inlineinherited |
TODO.
|
inline |
|
pure virtualinherited |
Used for ODEs. It verifies the existence of solutions to IVP at time t and set of initial conditions x over the time step. If the function succeeds, a rigorous bound for the trajectories is returned. Otherwise, an exception is thrown.
Implemented in capd::dynsys::OdeSolver< MapT, StepControlPolicyT, EnclosurePolicy, CurveT >, and capd::dynsys::OdeSolver< MapT, capd::dynsys::ILastTermsStepControl, HighOrderEnclosure, capd::diffAlgebra::C2Curve< capd::diffAlgebra::BasicC2Curve< typename MapT::MatrixType > > >.
|
inherited |
|
inlinevirtualinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
inlineinherited |
|
inlineinherited |
Returns maximal allocated order - used to avoid memory reallocation.
|
inherited |
|
inlineinherited |
|
inherited |
|
inlineinherited |
|
inherited |
|
inline |
const BasicCnOdeSolver< MapType, StepControlType, CurveT >::SolutionCurve & capd::dynsys::BasicCnOdeSolver< MapType, StepControlType, CurveT >::getCurve |
|
inlinestaticinherited |
|
inlinevirtualinherited |
|
inlineprotectedinherited |
|
inlineprotectedinherited |
|
inlineprotectedinherited |
|
inlineinherited |
|
inlineinherited |
Returns the order of Taylor interpolation.
|
inlineinherited |
|
inherited |
|
inlineinherited |
|
inlinevirtualinherited |
|
inlinevirtualinherited |
Returns time step of the dynamical system. By default it returns one - time step for discrete DS (maps). Shall be overridden in classes that implement numerical schemes for ODEs.
|
inlineinherited |
Returns current time step.
|
inline |
|
inlineinherited |
VectorFieldType& capd::dynsys::BasicCnOdeSolver< MapT, StepControlT, CurveT >::getVectorField |
|
inline |
Returns vector field.
|
virtualinherited |
|
inlineinherited |
|
pure virtualinherited |
Computes derivative of function (usually numerical scheme for an ODE) at time t and point iv.
Implemented in capd::dynsys::OdeSolver< MapT, StepControlPolicyT, EnclosurePolicy, CurveT >, and capd::dynsys::OdeSolver< MapT, capd::dynsys::ILastTermsStepControl, HighOrderEnclosure, capd::diffAlgebra::C2Curve< capd::diffAlgebra::BasicC2Curve< typename MapT::MatrixType > > >.
|
virtualinherited |
|
virtualinherited |
this function returns Lipschitz contants for maps. It should be overriden in classes implementing flows.
|
inlineinherited |
|
virtualinherited |
|
inlineinherited |
This operator computes image of the set (in given representation) using set.move function, see capd/dynsys/Move.h for details This template together with SetTraits prevent usage of various types of jets with incompatible solvers. The user will get an exception at runtime with clear message instead of unreadable compiler error. In this case a specialization CnJetMove is used meaning that this solver can integrate any type of jets.
|
inlineinherited |
Computes image of vector v after one time step. The argument t is updated in this procedure.
|
inlineinherited |
Computes next point on the trajectory of a nonautonomous system, first and second order derivatives with respect to initial conditions. Initial conditions for variational equations are Id and zero, respectively.
|
inlineinherited |
Computes image of vector v and derivatives of the flow with respect to init condition (v,identity). Version for nonautonomous systems. The argument t is updated in this procedure.
|
inlineinherited |
Computes image of vector v and derivatives of a flow with respect to init condition (v, derivative). The argument t is updated in this procedure.
|
inlineinherited |
The routine computes next point, derivatives and second order derivatives of a nonautonomous flow. Initial conditions for variational equations are V and H, respectively.
|
inherited |
Initial conditions for the trajectory and for variational equations up to given degree are given in argument in_out_coeffs. The full result is stored in in_out_coeffs and also the new point on the trajectory is returned. Note: CnCoeffType is a data structure that stores current time.
|
inline |
This operator computes image of the set (in given representation) using set.move function, see capd/dynsys/Move.h for details This template together with SetTraits prevent usage of various types of jets with incompatible solvers. The user will get an exception at runtime with clear message instead of unreadable compiler error. In this case a specialization CnSetMove is used meaning that this solver can integrate C^0, C^1 and C^2 sets only. Moreover, it cannot integrate nonrigorous jets (for user safety).
|
inline |
|
inherited |
Computes image of vector v after one time step.
|
inherited |
Computes image of vector v and derivatives of a flow with respect to init condition (v, derivative)
|
inherited |
The routine computes next point, derivatives and second order derivatives of a flow. Initial conditions for variational equations are V and H, respectively.
|
inherited |
Computes next point on the trajectory, first and second order derivatives with respect to initial conditions. Initial conditions for variational equations are Id and zero, respectively.
|
inherited |
Computes image of vector v and derivatives of the flow with respect to init condition (v,identity). Version for autonomous systems.
|
inlineprotected |
|
inlinevirtualinherited |
|
pure virtualinherited |
Computes value of function (usually numerical scheme for an ODE) at time t and point iv.
Implemented in capd::dynsys::OdeSolver< MapT, StepControlPolicyT, EnclosurePolicy, CurveT >, and capd::dynsys::OdeSolver< MapT, capd::dynsys::ILastTermsStepControl, HighOrderEnclosure, capd::diffAlgebra::C2Curve< capd::diffAlgebra::BasicC2Curve< typename MapT::MatrixType > > >.
|
pure virtualinherited |
Computes and returns bound for local error of a function (for instance if Phi is finite Taylor series of a map then this method computes bound for Lagrange remainder). If DynSys is an ODE, then out_enc contains enclosure of trajectories over the time step. If the function cannot validate existence of solutions to ODE over the time step, out_enc might be in an inconsistent state.
Implemented in capd::dynsys::OdeSolver< MapT, StepControlPolicyT, EnclosurePolicy, CurveT >, and capd::dynsys::OdeSolver< MapT, capd::dynsys::ILastTermsStepControl, HighOrderEnclosure, capd::diffAlgebra::C2Curve< capd::diffAlgebra::BasicC2Curve< typename MapT::MatrixType > > >.
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
Resets the mask of derivatives.
In consequence, full jet of derivatives will be computed after call to any method that computes derivative, hessian or jet.
|
inlineprotected |
|
inlineprotected |
|
inlineprotected |
|
inlineprotected |
|
inlineinherited |
|
inlineprotectedinherited |
|
inline |
|
inherited |
Sets new maximal degree of the jet and reallocates memory.
|
inlinevirtualinherited |
|
inlineinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
inlineinherited |
|
inlineinherited |
|
inlineinherited |
The iterator range [b,e) should contain a range of Multiinideces the user requires to compute.
The method automatically adds all the depending partial derivatives to this collection and defines a mask for computation of partial derivtives.
[b,e) | iterator range which contains collection of multiindices |
|
inlineinherited |
|
virtualinherited |
Sets the order of the Taylor method.
Reimplemented from capd::diffAlgebra::BasicCnCurve< typename MapT::MatrixType >.
void capd::dynsys::BasicCnOdeSolver< MapType, StepControlType, CurveT >::setOrder |
Sets the order of the Taylor method.
|
inlineinherited |
|
inlineinherited |
Sets time step and turns off step control.
|
inline |
Sets time step and turns off step control.
|
inlineinherited |
|
protectedinherited |
|
inherited |
|
inlineinherited |
|
inlineinherited |
|
protectedinherited |
|
inherited |
|
inherited |
|
inherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
protectedinherited |
|
inherited |
|
protectedinherited |
|
inherited |
|
protectedinherited |
domain
|
protected |
|
protectedinherited |
|
protected |