Standard RK integrator. More...
#include <rk.hpp>
Public Types | |
typedef ActionDataAbstractTpl< Scalar > | ActionDataAbstract |
typedef IntegratedActionModelAbstractTpl< Scalar > | Base |
typedef ControlParametrizationModelAbstractTpl< Scalar > | ControlParametrizationModelAbstract |
typedef IntegratedActionDataRKTpl< Scalar > | Data |
typedef DifferentialActionModelAbstractTpl< Scalar > | DifferentialActionModelAbstract |
typedef MathBaseTpl< Scalar > | MathBase |
typedef MathBase::MatrixXs | MatrixXs |
typedef MathBase::VectorXs | VectorXs |
Public Types inherited from IntegratedActionModelAbstractTpl< _Scalar > | |
typedef ActionDataAbstractTpl< Scalar > | ActionDataAbstract |
typedef ActionModelAbstractTpl< Scalar > | Base |
typedef ControlParametrizationModelAbstractTpl< Scalar > | ControlParametrizationModelAbstract |
typedef IntegratedActionDataAbstractTpl< Scalar > | Data |
typedef DifferentialActionModelAbstractTpl< Scalar > | DifferentialActionModelAbstract |
typedef MathBaseTpl< Scalar > | MathBase |
typedef MathBase::MatrixXs | MatrixXs |
typedef MathBase::VectorXs | VectorXs |
Public Types inherited from ActionModelAbstractTpl< _Scalar > | |
typedef ActionDataAbstractTpl< Scalar > | ActionDataAbstract |
typedef MathBaseTpl< Scalar > | MathBase |
typedef StateAbstractTpl< Scalar > | StateAbstract |
typedef MathBase::VectorXs | VectorXs |
Public Member Functions | |
IntegratedActionModelRKTpl (boost::shared_ptr< DifferentialActionModelAbstract > model, boost::shared_ptr< ControlParametrizationModelAbstract > control, const RKType rktype, const Scalar time_step=Scalar(1e-3), const bool with_cost_residual=true) | |
Initialize the RK integrator. More... | |
IntegratedActionModelRKTpl (boost::shared_ptr< DifferentialActionModelAbstract > model, const RKType rktype, const Scalar time_step=Scalar(1e-3), const bool with_cost_residual=true) | |
Initialize the RK integrator. More... | |
virtual void | calc (const boost::shared_ptr< ActionDataAbstract > &data, const Eigen::Ref< const VectorXs > &x) |
Integrate the total cost value for nodes that depends only on the state using RK scheme. More... | |
virtual void | calc (const boost::shared_ptr< ActionDataAbstract > &data, const Eigen::Ref< const VectorXs > &x, const Eigen::Ref< const VectorXs > &u) |
Integrate the differential action model using RK scheme. More... | |
virtual void | calcDiff (const boost::shared_ptr< ActionDataAbstract > &data, const Eigen::Ref< const VectorXs > &x) |
Compute the partial derivatives of the cost. More... | |
virtual void | calcDiff (const boost::shared_ptr< ActionDataAbstract > &data, const Eigen::Ref< const VectorXs > &x, const Eigen::Ref< const VectorXs > &u) |
Compute the partial derivatives of the RK integrator. More... | |
virtual bool | checkData (const boost::shared_ptr< ActionDataAbstract > &data) |
Checks that a specific data belongs to this model. | |
virtual boost::shared_ptr< ActionDataAbstract > | createData () |
Create the RK integrator data. More... | |
std::size_t | get_ni () const |
Return the number of nodes of the integrator. | |
virtual void | print (std::ostream &os) const |
Print relevant information of the RK integrator model. More... | |
virtual void | quasiStatic (const boost::shared_ptr< ActionDataAbstract > &data, Eigen::Ref< VectorXs > u, const Eigen::Ref< const VectorXs > &x, const std::size_t maxiter=100, const Scalar tol=Scalar(1e-9)) |
Computes the quasic static commands. More... | |
Public Member Functions inherited from IntegratedActionModelAbstractTpl< _Scalar > | |
IntegratedActionModelAbstractTpl (boost::shared_ptr< DifferentialActionModelAbstract > model, boost::shared_ptr< ControlParametrizationModelAbstract > control, const Scalar time_step=Scalar(1e-3), const bool with_cost_residual=true) | |
Initialize the integrator. More... | |
IntegratedActionModelAbstractTpl (boost::shared_ptr< DifferentialActionModelAbstract > model, const Scalar time_step=Scalar(1e-3), const bool with_cost_residual=true) | |
Initialize the integrator. More... | |
DEPRECATED ("The DifferentialActionModel should be set at construction time", void set_differential(boost::shared_ptr< DifferentialActionModelAbstract > model)) | |
const boost::shared_ptr< ControlParametrizationModelAbstract > & | get_control () const |
Return the control parametrization model associated to this integrated action model. | |
const boost::shared_ptr< DifferentialActionModelAbstract > & | get_differential () const |
Return the differential action model associated to this integrated action model. | |
const Scalar | get_dt () const |
Return the time step used for the integration. | |
virtual const VectorXs & | get_g_lb () const |
Return the lower bound of the inequality constraints. | |
virtual const VectorXs & | get_g_ub () const |
Return the upper bound of the inequality constraints. | |
virtual std::size_t | get_ng () const |
Return the number of inequality constraints. | |
virtual std::size_t | get_nh () const |
Return the number of equality constraints. | |
void | set_dt (const Scalar dt) |
Set the time step for the integration. | |
Public Member Functions inherited from ActionModelAbstractTpl< _Scalar > | |
ActionModelAbstractTpl (boost::shared_ptr< StateAbstract > state, const std::size_t nu, const std::size_t nr=0, const std::size_t ng=0, const std::size_t nh=0) | |
Initialize the action model. More... | |
bool | get_has_control_limits () const |
Indicates if there are defined control limits. | |
std::size_t | get_nr () const |
Return the dimension of the cost-residual vector. | |
std::size_t | get_nu () const |
Return the dimension of the control input. | |
const boost::shared_ptr< StateAbstract > & | get_state () const |
Return the state. | |
const VectorXs & | get_u_lb () const |
Return the control lower bound. | |
const VectorXs & | get_u_ub () const |
Return the control upper bound. | |
VectorXs | quasiStatic_x (const boost::shared_ptr< ActionDataAbstract > &data, const VectorXs &x, const std::size_t maxiter=100, const Scalar tol=Scalar(1e-9)) |
void | set_g_lb (const VectorXs &g_lb) |
Modify the lower bound of the inequality constraints. | |
void | set_g_ub (const VectorXs &g_ub) |
Modify the upper bound of the inequality constraints. | |
void | set_u_lb (const VectorXs &u_lb) |
Modify the control lower bounds. | |
void | set_u_ub (const VectorXs &u_ub) |
Modify the control upper bounds. | |
Public Attributes | |
EIGEN_MAKE_ALIGNED_OPERATOR_NEW typedef _Scalar | Scalar |
Public Attributes inherited from IntegratedActionModelAbstractTpl< _Scalar > | |
EIGEN_MAKE_ALIGNED_OPERATOR_NEW typedef _Scalar | Scalar |
Public Attributes inherited from ActionModelAbstractTpl< _Scalar > | |
EIGEN_MAKE_ALIGNED_OPERATOR_NEW typedef _Scalar | Scalar |
Protected Attributes | |
std::size_t | ng_ |
< Differential action model | |
std::size_t | nh_ |
< Number of inequality constraints | |
std::size_t | nu_ |
< Number of equality constraints | |
boost::shared_ptr< StateAbstract > | state_ |
< Dimension of the control | |
Protected Attributes inherited from IntegratedActionModelAbstractTpl< _Scalar > | |
boost::shared_ptr< ControlParametrizationModelAbstract > | control_ |
Model of the control parametrization. | |
boost::shared_ptr< DifferentialActionModelAbstract > | differential_ |
Differential action model that is integrated. | |
bool | has_control_limits_ |
std::size_t | nr_ |
std::size_t | nu_ |
< Dimension of the cost residual | |
boost::shared_ptr< StateAbstract > | state_ |
< Dimension of the control | |
Scalar | time_step2_ |
Square of the time step used for integration. | |
Scalar | time_step_ |
Time step used for integration. | |
VectorXs | u_lb_ |
< Model of the state | |
VectorXs | u_ub_ |
< Lower control limits | |
bool | with_cost_residual_ |
Flag indicating whether a cost residual is used. | |
Protected Attributes inherited from ActionModelAbstractTpl< _Scalar > | |
VectorXs | g_lb_ |
Lower bound of the inequality constraints. | |
VectorXs | g_ub_ |
Lower bound of the inequality constraints. | |
bool | has_control_limits_ |
std::size_t | ng_ |
Number of inequality constraints. | |
std::size_t | nh_ |
Number of equality constraints. | |
std::size_t | nr_ |
Dimension of the cost residual. | |
std::size_t | nu_ |
Control dimension. | |
boost::shared_ptr< StateAbstract > | state_ |
Model of the state. | |
VectorXs | u_lb_ |
Lower control limits. | |
VectorXs | u_ub_ |
Upper control limits. | |
VectorXs | unone_ |
Neutral state. | |
Additional Inherited Members | |
Protected Member Functions inherited from IntegratedActionModelAbstractTpl< _Scalar > | |
void | init () |
< Upper control limits | |
Protected Member Functions inherited from ActionModelAbstractTpl< _Scalar > | |
void | update_has_control_limits () |
Update the status of the control limits (i.e. if there are defined limits) | |
Standard RK integrator.
It applies a standard RK integration schemes to a differential (i.e., continuous time) action model. The available integrators are: RK2, RK3, and RK4.
This standard RK scheme introduces also the possibility to parametrize the control trajectory inside an integration step, for instance using polynomials. This requires introducing some notation to clarify the difference between the control inputs of the differential model and the control inputs to the integrated model. We have decided to use \(\mathbf{w}\) to refer to the control inputs of the differential model and \(\mathbf{u}\) for the control inputs of the integrated action model.
IntegratedActionModelRKTpl | ( | boost::shared_ptr< DifferentialActionModelAbstract > | model, |
boost::shared_ptr< ControlParametrizationModelAbstract > | control, | ||
const RKType | rktype, | ||
const Scalar | time_step = Scalar(1e-3) , |
||
const bool | with_cost_residual = true |
||
) |
Initialize the RK integrator.
[in] | model | Differential action model |
[in] | control | Control parametrization |
[in] | rktype | Type of RK integrator |
[in] | time_step | Step time (default 1e-3) |
[in] | with_cost_residual | Compute cost residual (default true) |
IntegratedActionModelRKTpl | ( | boost::shared_ptr< DifferentialActionModelAbstract > | model, |
const RKType | rktype, | ||
const Scalar | time_step = Scalar(1e-3) , |
||
const bool | with_cost_residual = true |
||
) |
Initialize the RK integrator.
This initialization uses ControlParametrizationPolyZeroTpl
for the control parametrization.
[in] | model | Differential action model |
[in] | rktype | Type of RK integrator |
[in] | time_step | Step time (default 1e-3) |
[in] | with_cost_residual | Compute cost residual (default true) |
|
virtual |
Integrate the differential action model using RK scheme.
[in] | data | RK integrator data |
[in] | x | State point \(\mathbf{x}\in\mathbb{R}^{ndx}\) |
[in] | u | Control input \(\mathbf{u}\in\mathbb{R}^{nu}\) |
Implements ActionModelAbstractTpl< _Scalar >.
|
virtual |
Integrate the total cost value for nodes that depends only on the state using RK scheme.
It computes the total cost and defines the next state as the current one. This function is used in the terminal nodes of an optimal control problem.
[in] | data | RK integrator data |
[in] | x | State point \(\mathbf{x}\in\mathbb{R}^{ndx}\) |
Reimplemented from ActionModelAbstractTpl< _Scalar >.
|
virtual |
Compute the partial derivatives of the RK integrator.
[in] | data | RK integrator data |
[in] | x | State point \(\mathbf{x}\in\mathbb{R}^{ndx}\) |
[in] | u | Control input \(\mathbf{u}\in\mathbb{R}^{nu}\) |
Implements ActionModelAbstractTpl< _Scalar >.
|
virtual |
Compute the partial derivatives of the cost.
It updates the derivatives of the cost function with respect to the state only. This function is used in the terminal nodes of an optimal control problem.
[in] | data | RK integrator data |
[in] | x | State point \(\mathbf{x}\in\mathbb{R}^{ndx}\) |
Reimplemented from ActionModelAbstractTpl< _Scalar >.
|
virtual |
Create the RK integrator data.
Reimplemented from IntegratedActionModelAbstractTpl< _Scalar >.
|
virtual |
Computes the quasic static commands.
The quasic static commands are the ones produced for a the reference posture as an equilibrium point, i.e. for \(\mathbf{f^q_x}\delta\mathbf{q}+\mathbf{f_u}\delta\mathbf{u}=\mathbf{0}\)
[in] | data | RK integrator data |
[out] | u | Quasic static commands |
[in] | x | State point (velocity has to be zero) |
[in] | maxiter | Maximum allowed number of iterations |
[in] | tol | Tolerance |
Reimplemented from ActionModelAbstractTpl< _Scalar >.
|
virtual |
Print relevant information of the RK integrator model.
[out] | os | Output stream object |
Reimplemented from ActionModelAbstractTpl< _Scalar >.