cost-wrap.hpp

Go to the documentation of this file.

#pragma once
 
#include "aligator/compat/crocoddyl/fwd.hpp"
#include "aligator/core/cost-abstract.hpp"
#include <crocoddyl/core/cost-base.hpp>
#include <crocoddyl/core/action-base.hpp>
 
namespace aligator::compat::croc {
 
template <typename _Scalar>
struct CrocCostModelWrapperTpl : CostAbstractTpl<_Scalar> {
  using Scalar = _Scalar;
  ALIGATOR_DYNAMIC_TYPEDEFS(Scalar);
  using CrocCostModel = crocoddyl::CostModelAbstractTpl<Scalar>;
  using CrocActionModel = crocoddyl::ActionModelAbstractTpl<Scalar>;
  using Base = CostAbstractTpl<Scalar>;
  using BaseData = CostDataAbstractTpl<Scalar>;
  using StateWrap = StateWrapperTpl<Scalar>;
 
  shared_ptr<CrocCostModel> croc_cost_;
  shared_ptr<CrocActionModel> action_model_;

  explicit CrocCostModelWrapperTpl(shared_ptr<CrocCostModel> cost)
      : Base(StateWrap(cost->get_state()), (int)cost->get_nu())
      , croc_cost_(cost) {}

  explicit CrocCostModelWrapperTpl(shared_ptr<CrocActionModel> action_model)
      : Base(StateWrap(action_model->get_state()), (int)action_model->get_nu())
      , action_model_(action_model) {}
 
  void evaluate(const ConstVectorRef &x, const ConstVectorRef &u,
                BaseData &data) const {
    using Data = CrocCostDataWrapperTpl<Scalar>;
    Data &d = static_cast<Data &>(data);
    if (croc_cost_ != 0) {
      croc_cost_->calc(d.croc_cost_data_, x, u);
      d.value_ = d.croc_cost_data_->cost;
    } else {
      action_model_->calc(d.croc_act_data_, x);
      d.value_ = d.croc_act_data_->cost;
    }
  }
 
  void computeGradients(const ConstVectorRef &x, const ConstVectorRef &u,
                        BaseData &data) const {
    using Data = CrocCostDataWrapperTpl<Scalar>;
    Data &d = static_cast<Data &>(data);
    if (croc_cost_ != 0) {
      croc_cost_->calcDiff(d.croc_cost_data_, x, u);
      d.Lx_ = d.croc_cost_data_->Lx;
      d.Lu_ = d.croc_cost_data_->Lu;
    } else {
      action_model_->calcDiff(d.croc_act_data_, x);
      d.Lx_ = d.croc_act_data_->Lx;
      d.Lu_ = d.croc_act_data_->Lu;
    }
  }
 
  void computeHessians(const ConstVectorRef &x, const ConstVectorRef &u,
                       BaseData &data) const {
    using Data = CrocCostDataWrapperTpl<Scalar>;
    Data &d = static_cast<Data &>(data);
    if (croc_cost_ != 0) {
      croc_cost_->calcDiff(d.croc_cost_data_, x, u);
      d.Lxx_ = d.croc_cost_data_->Lxx;
      d.Lxu_ = d.croc_cost_data_->Lxu;
      d.Luu_ = d.croc_cost_data_->Luu;
    } else {
      action_model_->calcDiff(d.croc_act_data_, x);
      d.Lxx_ = d.croc_act_data_->Lxx;
      d.Lxu_ = d.croc_act_data_->Lxu;
      d.Luu_ = d.croc_act_data_->Luu;
    }
  }
 
  shared_ptr<BaseData> createData() const {
    if (action_model_ != 0) {
      shared_ptr<crocoddyl::ActionDataAbstractTpl<Scalar>> am_data =
          action_model_->createData();
      return std::make_shared<CrocCostDataWrapperTpl<Scalar>>(am_data);
    } else {
      ALIGATOR_DOMAIN_ERROR("Invalid call. Cannot build Data from"
                            "crocoddyl cost model only.");
    }
  }
};
 
template <typename Scalar>
struct CrocCostDataWrapperTpl : CostDataAbstractTpl<Scalar> {
  using CostData = ::crocoddyl::CostDataAbstractTpl<Scalar>;
  using ActionData = ::crocoddyl::ActionDataAbstractTpl<Scalar>;
  using Base = CostDataAbstractTpl<Scalar>;
  shared_ptr<CostData> croc_cost_data_;
  shared_ptr<ActionData> croc_act_data_;
 
  explicit CrocCostDataWrapperTpl(const shared_ptr<CostData> &crocdata)
      : Base((int)crocdata->Lx.rows(), (int)crocdata->Lu.rows())
      , croc_cost_data_(crocdata) {}
 
  explicit CrocCostDataWrapperTpl(const shared_ptr<ActionData> &actdata)
      : Base((int)actdata->Lx.rows(), (int)actdata->Lu.rows())
      , croc_act_data_(actdata) {}
};
 
#ifdef ALIGATOR_ENABLE_TEMPLATE_INSTANTIATION
extern template struct CrocCostModelWrapperTpl<context::Scalar>;
extern template struct CrocCostDataWrapperTpl<context::Scalar>;
#endif
 
} // namespace aligator::compat::croc