negative-orthant.hpp

#pragma once
 
#include "proxsuite-nlp/constraint-set.hpp"
 
namespace proxsuite {
namespace nlp {

template <typename _Scalar>
struct NegativeOrthantTpl : ConstraintSetTpl<_Scalar> {
  using Scalar = _Scalar;
  PROXSUITE_NLP_DYNAMIC_TYPEDEFS(Scalar);
 
  NegativeOrthantTpl() = default;
  NegativeOrthantTpl(const NegativeOrthantTpl &) = default;
  NegativeOrthantTpl &operator=(const NegativeOrthantTpl &) = default;
  NegativeOrthantTpl(NegativeOrthantTpl &&) = default;
  NegativeOrthantTpl &operator=(NegativeOrthantTpl &&) = default;
 
  using Base = ConstraintSetTpl<Scalar>;
  using ActiveType = typename Base::ActiveType;
 
  void projection(const ConstVectorRef &z, VectorRef zout) const {
    zout = z.cwiseMin(static_cast<Scalar>(0.));
  }
 
  void normalConeProjection(const ConstVectorRef &z, VectorRef zout) const {
    zout = z.cwiseMax(static_cast<Scalar>(0.));
  }

  void computeActiveSet(const ConstVectorRef &z,
                        Eigen::Ref<ActiveType> out) const {
    out.array() = (z.array() > static_cast<Scalar>(0.));
  }
};
 
template <typename Scalar>
using NegativeOrthant PROXSUITE_NLP_DEPRECATED_MESSAGE(
    "Use NegativeOrthantTpl<T> instead") = NegativeOrthantTpl<Scalar>;
 
#ifdef PROXSUITE_NLP_ENABLE_TEMPLATE_INSTANTIATION
extern template struct PROXSUITE_NLP_EXPLICIT_INSTANTIATION_DECLARATION_DLLAPI
    NegativeOrthantTpl<context::Scalar>;
#endif
 
} // namespace nlp
} // namespace proxsuite
 
#include "proxsuite-nlp/modelling/constraints.hpp"