constraint-set.hxx

Go to the documentation of this file.

#pragma once
 
#include "proxsuite-nlp/constraint-set.hpp"
 
namespace proxsuite {
namespace nlp {
 
template <typename Scalar>
void ConstraintSetTpl<Scalar>::applyProjectionJacobian(const ConstVectorRef &z,
                                                       MatrixRef Jout) const {
  const int nr = (int)z.size();
  assert(nr == Jout.rows());
  ActiveType active_set(nr);
  computeActiveSet(z, active_set);
  for (int i = 0; i < nr; i++) {
    if (active_set(i)) {
      Jout.row(i).setZero();
    }
  }
}
 
template <typename Scalar>
void ConstraintSetTpl<Scalar>::applyNormalConeProjectionJacobian(
    const ConstVectorRef &z, MatrixRef Jout) const {
  const int nr = (int)z.size();
  assert(nr == Jout.rows());
  ActiveType active_set(nr);
  computeActiveSet(z, active_set);
  for (int i = 0; i < nr; i++) {
    if (!active_set(i)) {
      Jout.row(i).setZero();
    }
  }
}
 
} // namespace nlp
} // namespace proxsuite