|
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) |
| Computes value of function (usually numerical scheme for an ODE) at time t and point iv. More...
|
|
MatrixType | JacPhi (const ScalarType &t, const VectorType &iv) |
| Computes derivative of function (usually numerical scheme for an ODE) at time t and point iv. More...
|
|
VectorType | enclosure (const ScalarType &t, const VectorType &x) |
| 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...
|
|
VectorType | Remainder (const ScalarType &t, const VectorType &iv, VectorType &o_enc) |
| 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...
|
|
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 |
| 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...
|
|
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 () |
|
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 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 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...
|
|
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) |
|
const VectorFieldType & | getVectorField () const |
| Returns vector field. More...
|
|
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...
|
|
size_type | degree () const |
|
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 |
|
Scalar | getMaxStep () const |
|
void | setMaxStep (Scalar maxStep) |
|
HessianType | hessian (const ScalarType &h) const |
|
JetType | jet (const ScalarType &h) const |
|
void | eval (ScalarType h, JetType &v) const |
|
VectorType | timeDerivative (const ScalarType &h) const |
|
MatrixType | derivative (const ScalarType &h) const |
|
MatrixType | operator[] (const ScalarType &h) const |
|
virtual void | setDomain (Real left, Real right) |
|
virtual Real | getLeftDomain () const |
|
virtual Real | getRightDomain () const |
|