common.hpp

Go to the documentation of this file.

//
// Copyright (c) 2022 INRIA
//
#ifndef PROXSUITE_HELPERS_COMMON_HPP
#define PROXSUITE_HELPERS_COMMON_HPP
 
#include "proxsuite/config.hpp"
#include <limits>
 
namespace proxsuite {
namespace helpers {
 
template<typename Scalar>
struct infinite_bound
{
  static Scalar value()
  {
    using namespace std;
    return sqrt(std::numeric_limits<Scalar>::max());
  }
};
 
#define PROXSUITE_DEDUCE_RET(...)                                              \
  noexcept(noexcept(__VA_ARGS__))                                              \
    ->typename std::remove_const<decltype(__VA_ARGS__)>::type                  \
  {                                                                            \
    return __VA_ARGS__;                                                        \
  }                                                                            \
  static_assert(true, ".")
 
template<typename T, typename Scalar>
auto
at_most(T const& expr, const Scalar value) PROXSUITE_DEDUCE_RET(
  (expr.array() < value).select(expr, T::Constant(expr.rows(), value)));
 
template<typename T, typename Scalar>
auto
at_least(T const& expr, const Scalar value) PROXSUITE_DEDUCE_RET(
  (expr.array() > value).select(expr, T::Constant(expr.rows(), value)));
 
template<typename T>
auto
positive_part(T const& expr)
  PROXSUITE_DEDUCE_RET((expr.array() > 0).select(expr, T::Zero(expr.rows())));
 
template<typename T>
auto
negative_part(T const& expr)
  PROXSUITE_DEDUCE_RET((expr.array() < 0).select(expr, T::Zero(expr.rows())));
 
template<typename Condition, typename T, typename Scalar>
auto
select(Condition const& condition, T const& expr, const Scalar value)
  PROXSUITE_DEDUCE_RET((condition).select(expr,
                                          T::Constant(expr.rows(), value)));
 
} // helpers
} // proxsuite
 
#endif // ifndef PROXSUITE_HELPERS_COMMON_HPP