positive_quadrant.h

Go to the documentation of this file.

#ifndef SOPT_PROJECTED_ALGORITHM_H
#define SOPT_PROJECTED_ALGORITHM_H
 
#include "sopt/linear_transform.h"
#include "sopt/types.h"
 
#include <utility> // for std::forward
 
namespace sopt::algorithm {
 
template <typename ALGORITHM>
class PositiveQuadrant {
 public:
  using Algorithm = ALGORITHM;
  using Scalar = typename Algorithm::Scalar;
  using t_Vector = typename Algorithm::t_Vector;
  using t_IsConverged = typename Algorithm::t_IsConverged;
  using Diagnostic = typename ALGORITHM::Diagnostic;
  using DiagnosticAndResult = typename ALGORITHM::DiagnosticAndResult;
 
  PositiveQuadrant(Algorithm const &algo) : algorithm_(algo) {}
  PositiveQuadrant(Algorithm &&algo) : algorithm_(std::move(algo)) {}
 
  Algorithm &algorithm() { return algorithm_; }
  Algorithm const &algorithm() const { return algorithm_; }
 
  template <typename... T>
  Diagnostic operator()(t_Vector &out, T const &... args) const {
    auto const diagnostic = algorithm()(out, std::forward<T const &>(args)...);
    out = positive_quadrant(out);
    return diagnostic;
  };
 
  template <typename... T>
  DiagnosticAndResult operator()(T const &... args) const {
    auto result = algorithm()(std::forward<T const &>(args)...);
    result.x = positive_quadrant(result.x);
    return result;
  };
 
 protected:
  Algorithm algorithm_;
};
 
template <typename ALGORITHM>
PositiveQuadrant<ALGORITHM> positive_quadrant(ALGORITHM const &algo) {
  return {algo};
}
} // namespace sopt::algorithm
 
#endif