ROL_LinearOperatorFromConstraint.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_LINEAROPERATOR_FROM_EQUALITYCONSTRAINT_H
#define ROL_LINEAROPERATOR_FROM_EQUALITYCONSTRAINT_H
 
#include "ROL_Vector.hpp"
#include "ROL_Constraint.hpp"
#include "ROL_LinearOperator.hpp"
 
namespace ROL {
 
template <class Real>
class LinearOperatorFromConstraint : public LinearOperator<Real> {
private:
  const ROL::Ptr<const Vector<Real> > x_;
  ROL::Ptr<Constraint<Real> > con_;
 
 
public:
 
  LinearOperatorFromConstraint( const ROL::Ptr<const Vector<Real> > &x, 
                                        const ROL::Ptr<Constraint<Real> > &con ) : 
                                        x_(x), con_(con) {
  }
 
  virtual ~LinearOperatorFromConstraint() {}
 
  virtual void apply( Vector<Real> &Hv, const Vector<Real> &v, Real &tol ) const {
    con_->applyJacobian(Hv,v,*x_,tol);
  }
 
  // Not implemented for generic equality constraint
  virtual void applyInverse( Vector<Real> &Hv, const Vector<Real> &v, Real &tol ) const {
    Hv.set(v);
  }
 
}; // class LinearOperator
 
} // namespace ROL
 
#endif