44#ifndef ROL_AUGMENTEDLAGRANGIAN_SIMOPT_H
45#define ROL_AUGMENTEDLAGRANGIAN_SIMOPT_H
100 const ROL::Ptr<Objective_SimOpt<Real> >
obj_;
101 ROL::Ptr<QuadraticPenalty_SimOpt<Real> >
pen_;
129 const Real penaltyParameter,
133 ROL::ParameterList &parlist)
143 ROL::ParameterList& sublist = parlist.sublist(
"Step").sublist(
"Augmented Lagrangian");
145 int HessianApprox = sublist.get(
"Level of Hessian Approximation", 0);
147 pen_ = ROL::makePtr<QuadraticPenalty_SimOpt<Real>>(con,multiplier,penaltyParameter,simVec,optVec,conVec,
scaleLagrangian_,HessianApprox);
151 obj_->update(u,z,flag,iter);
152 pen_->update(u,z,flag,iter);
165 Real pval =
pen_->value(u,z,tol);
209 obj_->hessVec_11(hv,v,u,z,tol);
222 obj_->hessVec_12(hv,v,u,z,tol);
235 obj_->hessVec_21(hv,v,u,z,tol);
248 obj_->hessVec_22(hv,v,u,z,tol);
260 Real tol = std::sqrt(ROL_EPSILON<Real>());
271 pen_->getConstraintVec(c,u,z);
276 return pen_->getNumberConstraintEvaluations();
292 pen_->reset(multiplier,penaltyParameter);
Contains definitions of custom data types in ROL.
Provides the interface to evaluate the SimOpt augmented Lagrangian.
virtual Real getObjectiveValue(const Vector< Real > &u, const Vector< Real > &z)
virtual int getNumberFunctionEvaluations(void) const
AugmentedLagrangian_SimOpt(const ROL::Ptr< Objective_SimOpt< Real > > &obj, const ROL::Ptr< Constraint_SimOpt< Real > > &con, const Vector< Real > &multiplier, const Real penaltyParameter, const Vector< Real > &simVec, const Vector< Real > &optVec, const Vector< Real > &conVec, ROL::ParameterList &parlist)
virtual Real value(const Vector< Real > &u, const Vector< Real > &z, Real &tol)
Compute value.
virtual void getConstraintVec(Vector< Real > &c, const Vector< Real > &u, const Vector< Real > &z)
virtual int getNumberConstraintEvaluations(void) const
virtual void hessVec_22(Vector< Real > &hv, const Vector< Real > &v, const Vector< Real > &u, const Vector< Real > &z, Real &tol)
ROL::Ptr< Vector< Real > > gradient2_
ROL::Ptr< Vector< Real > > gradient1_
virtual void update(const Vector< Real > &u, const Vector< Real > &z, bool flag=true, int iter=-1)
Update objective function. u is an iterate, z is an iterate, flag = true if the iterate has changed...
virtual void gradient_2(Vector< Real > &g, const Vector< Real > &u, const Vector< Real > &z, Real &tol)
Compute gradient with respect to second component.
virtual void reset(const Vector< Real > &multiplier, const Real penaltyParameter)
bool isGradient2Computed_
ROL::Ptr< Vector< Real > > dualSimVector_
virtual int getNumberGradientEvaluations(void) const
ROL::Ptr< Vector< Real > > dualOptVector_
virtual void hessVec_12(Vector< Real > &hv, const Vector< Real > &v, const Vector< Real > &u, const Vector< Real > &z, Real &tol)
bool isGradient1Computed_
virtual void hessVec_11(Vector< Real > &hv, const Vector< Real > &v, const Vector< Real > &u, const Vector< Real > &z, Real &tol)
Apply Hessian approximation to vector.
virtual void hessVec_21(Vector< Real > &hv, const Vector< Real > &v, const Vector< Real > &u, const Vector< Real > &z, Real &tol)
const ROL::Ptr< Objective_SimOpt< Real > > obj_
ROL::Ptr< QuadraticPenalty_SimOpt< Real > > pen_
virtual void gradient_1(Vector< Real > &g, const Vector< Real > &u, const Vector< Real > &z, Real &tol)
Compute gradient with respect to first component.
Defines the constraint operator interface for simulation-based optimization.
Provides the interface to evaluate simulation-based objective functions.
Defines the linear algebra or vector space interface.
virtual void set(const Vector &x)
Set where .
virtual void scale(const Real alpha)=0
Compute where .
virtual const Vector & dual() const
Return dual representation of , for example, the result of applying a Riesz map, or change of basis,...
virtual void plus(const Vector &x)=0
Compute , where .
virtual ROL::Ptr< Vector > clone() const =0
Clone to make a new (uninitialized) vector.