ROL_UpperBoundToConstraint_Def.hpp

Go to the documentation of this file.

// @HEADER
// *****************************************************************************
//               Rapid Optimization Library (ROL) Package
//
// Copyright 2014 NTESS and the ROL contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER
 
#ifndef ROL_UPPER_BOUND_TO_CONSTRAINT_DEF_H
#define ROL_UPPER_BOUND_TO_CONSTRAINT_DEF_H
 
namespace ROL {
 
template<typename Real>
UpperBoundToConstraint<Real>::UpperBoundToConstraint(BoundConstraint<Real> &bnd) {
  up_ = bnd.getUpperBound()->clone();
  up_->set(*bnd.getUpperBound());
}
 
template<typename Real>
UpperBoundToConstraint<Real>::UpperBoundToConstraint(const Vector<Real> &up) {
  up_ = up.clone();
  up_->set(up);
}
 
template<typename Real>
void UpperBoundToConstraint<Real>::value(Vector<Real> &c, const Vector<Real> &x, Real &tol) {
  const Real one(1);
  c.set(*up_);
  c.axpy(-one,x);
}
 
template<typename Real>
void UpperBoundToConstraint<Real>::applyJacobian(Vector<Real> &jv,
                                                 const Vector<Real> &v,
                                                 const Vector<Real> &x,
                                                 Real &tol) {
  const Real one(1);
  jv.set(v);
  jv.scale(-one);
}
 
template<typename Real>
void UpperBoundToConstraint<Real>::applyAdjointJacobian(Vector<Real> &ajv,
                                                        const Vector<Real> &v,
                                                        const Vector<Real> &x,
                                                        Real &tol) {
  const Real one(1);
  ajv.set(v);
  ajv.scale(-one);
}
 
template<typename Real>
void UpperBoundToConstraint<Real>::applyAdjointHessian(Vector<Real> &ahuv,
                                                       const Vector<Real> &u,
                                                       const Vector<Real> &v,
                                                       const Vector<Real> &x,
                                                       Real &tol) {
  ahuv.zero();
}
 
}
 
#endif