44 #ifndef ROL_TYPEB_PRIMALINTERIORPOINTALGORITHM_DEF_HPP
45 #define ROL_TYPEB_PRIMALINTERIORPOINTALGORITHM_DEF_HPP
52 template<
typename Real>
55 list_(list), subproblemIter_(0), print_(false) {
61 ParameterList& steplist = list.sublist(
"Step").sublist(
"Interior Point");
62 state_->searchSize = steplist.get(
"Initial Barrier Parameter", 1.0);
63 mumin_ = steplist.get(
"Minimum Barrier Parameter", 1e-4);
64 mumax_ = steplist.get(
"Maximum Barrier Parameter", 1e8);
65 rho_ = steplist.get(
"Barrier Penalty Reduction Factor", 0.5);
67 kappaD_ = steplist.get(
"Linear Damping Coefficient", 1.e-4);
68 print_ = steplist.sublist(
"Subproblem").get(
"Print History",
false);
70 gtol_ = steplist.sublist(
"Subproblem").get(
"Initial Optimality Tolerance", 1e-2);
72 int maxit = steplist.sublist(
"Subproblem").get(
"Iteration Limit", 1000);
73 list_.sublist(
"Status Test").set(
"Iteration Limit", maxit);
75 gtolrate_ = steplist.sublist(
"Subproblem").get(
"Optimality Tolerance Reduction Factor", 0.1);
76 mingtol_ =
static_cast<Real
>(1e-2)*list.sublist(
"Status Test").get(
"Gradient Tolerance", 1e-8);
78 stepname_ = steplist.sublist(
"Subproblem").get(
"Step Type",
"Augmented Lagrangian");
81 verbosity_ = list.sublist(
"General").get(
"Output Level", 0);
87 template<
typename Real>
93 std::ostream &outStream) {
95 if (proj_ == nullPtr) {
96 proj_ = makePtr<PolyhedralProjection<Real>>(makePtrFromRef(bnd));
99 proj_->project(x,outStream);
106 updateState(x,ipobj,bnd,pwa);
110 template<
typename Real>
115 std::ostream &outStream) {
117 Real zerotol = std::sqrt(ROL_EPSILON<Real>());
119 if (state_-> iter == 0) {
128 ipobj.
gradient(*state_->gradientVec, x, zerotol);
131 pwa.
axpy(-one,state_->gradientVec->dual());
132 proj_->project(pwa,outStream);
134 state_->gnorm = pwa.
norm();
140 template<
typename Real>
145 std::ostream &outStream ) {
147 Ptr<Vector<Real>> pwa = x.
clone();
150 x,g,useLinearDamping_,kappaD_,
152 initialize(x,g,ipobj,bnd,*pwa,outStream);
153 Ptr<TypeU::Algorithm<Real>> algo;
156 if (verbosity_ > 0) writeOutput(outStream,
true);
158 while (status_->check(*state_)) {
160 list_.sublist(
"Status Test").set(
"Gradient Tolerance", gtol_);
161 list_.sublist(
"Status Test").set(
"Step Tolerance", stol_);
162 algo = TypeU::AlgorithmFactory<Real>(list_,secant_);
163 if (hasPolyProj_) algo->run(x,g,ipobj,
164 *proj_->getLinearConstraint(),
165 *proj_->getMultiplier(),
166 *proj_->getResidual(),outStream);
167 else algo->run(x,g,ipobj,outStream);
168 subproblemIter_ = algo->getState()->iter;
169 state_->nfval += algo->getState()->nfval;
170 state_->ngrad += algo->getState()->ngrad;
173 state_->stepVec->set(x);
174 state_->stepVec->axpy(-one,*state_->iterateVec);
175 state_->snorm = state_->stepVec->norm();
178 state_->iterateVec->set(x);
185 if( (rho_< one && state_->searchSize > mumin_) || (rho_ > one && state_->searchSize < mumax_) ) {
186 state_->searchSize *= rho_;
189 gtol_ *= gtolrate_; gtol_ = std::max(gtol_,mingtol_);
190 stol_ =
static_cast<Real
>(1e-6)*gtol_;
194 updateState(x,ipobj,bnd,*pwa,outStream);
197 if (verbosity_ > 0) writeOutput(outStream,writeHeader_);
202 template<
typename Real>
204 std::ios_base::fmtflags osFlags(os.flags());
205 if (verbosity_ > 1) {
206 os << std::string(109,
'-') << std::endl;
207 os <<
"Interior Point Solver";
208 os <<
" status output definitions" << std::endl << std::endl;
209 os <<
" iter - Number of iterates (steps taken)" << std::endl;
210 os <<
" fval - Objective function value" << std::endl;
211 os <<
" gnorm - Norm of the gradient" << std::endl;
212 os <<
" snorm - Norm of the step (update to optimization vector)" << std::endl;
213 os <<
" penalty - Penalty parameter for bound constraints" << std::endl;
214 os <<
" #fval - Cumulative number of times the objective function was evaluated" << std::endl;
215 os <<
" #grad - Cumulative number of times the gradient was computed" << std::endl;
216 os <<
" optTol - Subproblem optimality tolerance" << std::endl;
217 os <<
" subiter - Number of subproblem iterations" << std::endl;
218 os << std::string(109,
'-') << std::endl;
222 os << std::setw(6) << std::left <<
"iter";
223 os << std::setw(15) << std::left <<
"fval";
224 os << std::setw(15) << std::left <<
"gnorm";
225 os << std::setw(15) << std::left <<
"snorm";
226 os << std::setw(10) << std::left <<
"penalty";
227 os << std::setw(8) << std::left <<
"#fval";
228 os << std::setw(8) << std::left <<
"#grad";
229 os << std::setw(10) << std::left <<
"optTol";
230 os << std::setw(8) << std::left <<
"subIter";
235 template<
typename Real>
237 std::ios_base::fmtflags osFlags(os.flags());
238 os << std::endl <<
"Interior Point Solver (Type B, Bound Constraints)";
240 os <<
"Subproblem Solver: " << stepname_ << std::endl;
244 template<
typename Real>
246 std::ios_base::fmtflags osFlags(os.flags());
247 os << std::scientific << std::setprecision(6);
248 if ( state_->iter == 0 ) writeName(os);
249 if ( write_header ) writeHeader(os);
250 if ( state_->iter == 0 ) {
252 os << std::setw(6) << std::left << state_->iter;
253 os << std::setw(15) << std::left << state_->value;
254 os << std::setw(15) << std::left << state_->gnorm;
255 os << std::setw(15) << std::left <<
"---";
256 os << std::scientific << std::setprecision(2);
257 os << std::setw(10) << std::left << state_->searchSize;
258 os << std::setw(8) << std::left << state_->nfval;
259 os << std::setw(8) << std::left << state_->ngrad;
260 os << std::setw(10) << std::left <<
"---";
261 os << std::setw(8) << std::left <<
"---";
266 os << std::setw(6) << std::left << state_->iter;
267 os << std::setw(15) << std::left << state_->value;
268 os << std::setw(15) << std::left << state_->gnorm;
269 os << std::setw(15) << std::left << state_->snorm;
270 os << std::scientific << std::setprecision(2);
271 os << std::setw(10) << std::left << state_->searchSize;
272 os << std::scientific << std::setprecision(6);
273 os << std::setw(8) << std::left << state_->nfval;
274 os << std::setw(8) << std::left << state_->ngrad;
275 os << std::scientific << std::setprecision(2);
276 os << std::setw(10) << std::left << gtol_;
277 os << std::scientific << std::setprecision(6);
278 os << std::setw(8) << std::left << subproblemIter_;
Provides the interface to evaluate objective functions.
virtual ROL::Ptr< Vector > clone() const =0
Clone to make a new (uninitialized) vector.
virtual void projectInterior(Vector< Real > &x)
Project optimization variables into the interior of the feasible set.
virtual void axpy(const Real alpha, const Vector &x)
Compute where .
InteriorPointAlgorithm(ParameterList &list, const Ptr< Secant< Real >> &secant=nullPtr)
void run(Vector< Real > &x, const Vector< Real > &g, Objective< Real > &obj, BoundConstraint< Real > &bnd, std::ostream &outStream=std::cout) override
Run algorithm on bound constrained problems (Type-B). This general interface supports the use of dual...
void writeName(std::ostream &os) const override
Print step name.
virtual void writeExitStatus(std::ostream &os) const
Defines the linear algebra or vector space interface.
void updatePenalty(const Real mu)
void initialize(Vector< Real > &x, const Vector< Real > &g, InteriorPointObjective< Real > &ipobj, BoundConstraint< Real > &bnd, Vector< Real > &pwa, std::ostream &outStream=std::cout)
void writeOutput(std::ostream &os, const bool write_header=false) const override
Print iterate status.
Provides an interface to run bound constrained optimization algorithms.
const Ptr< AlgorithmState< Real > > state_
Provides interface for and implements limited-memory secant operators.
void writeHeader(std::ostream &os) const override
Print iterate header.
Provides an interface to check status of optimization algorithms.
void update(const Vector< Real > &x, UpdateType type, int iter=-1)
Update objective function.
Provides the interface to apply upper and lower bound constraints.
void initialize(const Vector< Real > &x, const Vector< Real > &g)
void gradient(Vector< Real > &g, const Vector< Real > &x, Real &tol)
Compute gradient.
virtual void set(const Vector &x)
Set where .
virtual Real norm() const =0
Returns where .
void updateState(const Vector< Real > &x, InteriorPointObjective< Real > &ipobj, BoundConstraint< Real > &bnd, Vector< Real > &pwa, std::ostream &outStream=std::cout)
Real getObjectiveValue(const Vector< Real > &x, Real &tol)
const Ptr< CombinedStatusTest< Real > > status_