ROL
|
Provides an interface to run the projected secant algorithm. More...
#include <ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp>
Classes | |
class | HessianPDAS |
class | HessianPDAS_Poly |
class | PrecondPDAS |
class | PrecondPDAS_Poly |
Public Member Functions | |
PrimalDualActiveSetAlgorithm (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 optimization vector spaces, where the user does not define the dual() method. More... | |
void | writeHeader (std::ostream &os) const override |
Print iterate header. More... | |
void | writeName (std::ostream &os) const override |
Print step name. More... | |
void | writeOutput (std::ostream &os, const bool write_header=false) const override |
Print iterate status. More... | |
Public Member Functions inherited from ROL::TypeB::Algorithm< Real > | |
virtual | ~Algorithm () |
Algorithm () | |
Constructor, given a step and a status test. More... | |
void | setStatusTest (const Ptr< StatusTest< Real >> &status, const bool combineStatus=false) |
virtual void | run (Problem< Real > &problem, std::ostream &outStream=std::cout) |
Run algorithm on bound constrained problems (Type-B). This is the primary Type-B interface. More... | |
virtual void | run (Vector< Real > &x, Objective< Real > &obj, BoundConstraint< Real > &bnd, std::ostream &outStream=std::cout) |
Run algorithm on bound constrained problems (Type-B). This is the primary Type-B interface. More... | |
virtual void | run (Vector< Real > &x, Objective< Real > &obj, BoundConstraint< Real > &bnd, Constraint< Real > &linear_econ, Vector< Real > &linear_emul, std::ostream &outStream=std::cout) |
Run algorithm on bound constrained problems with explicit linear constraints (Type-B). This is the primary Type-B with explicit linear constraints interface. More... | |
virtual void | run (Vector< Real > &x, const Vector< Real > &g, Objective< Real > &obj, BoundConstraint< Real > &bnd, Constraint< Real > &linear_econ, Vector< Real > &linear_emul, const Vector< Real > &linear_eres, std::ostream &outStream=std::cout) |
Run algorithm on bound constrained problems with explicit linear constraints (Type-B). This general interface supports the use of dual optimization vector spaces, where the user does not define the dual() method. More... | |
virtual void | run (Vector< Real > &x, Objective< Real > &obj, BoundConstraint< Real > &bnd, Constraint< Real > &linear_icon, Vector< Real > &linear_imul, BoundConstraint< Real > &linear_ibnd, std::ostream &outStream=std::cout) |
Run algorithm on bound constrained problems with explicit linear constraints (Type-B). This is the primary Type-B with explicit linear constraints interface. More... | |
virtual void | run (Vector< Real > &x, const Vector< Real > &g, Objective< Real > &obj, BoundConstraint< Real > &bnd, Constraint< Real > &linear_icon, Vector< Real > &linear_imul, BoundConstraint< Real > &linear_ibnd, const Vector< Real > &linear_ires, std::ostream &outStream=std::cout) |
Run algorithm on bound constrained problems with explicit linear constraints (Type-B). This general interface supports the use of dual optimization vector spaces, where the user does not define the dual() method. More... | |
virtual void | run (Vector< Real > &x, Objective< Real > &obj, BoundConstraint< Real > &bnd, Constraint< Real > &linear_econ, Vector< Real > &linear_emul, Constraint< Real > &linear_icon, Vector< Real > &linear_imul, BoundConstraint< Real > &linear_ibnd, std::ostream &outStream=std::cout) |
Run algorithm on bound constrained problems with explicit linear constraints (Type-B). This is the primary Type-B with explicit linear constraints interface. More... | |
virtual void | run (Vector< Real > &x, const Vector< Real > &g, Objective< Real > &obj, BoundConstraint< Real > &bnd, Constraint< Real > &linear_econ, Vector< Real > &linear_emul, const Vector< Real > &linear_eres, Constraint< Real > &linear_icon, Vector< Real > &linear_imul, BoundConstraint< Real > &linear_ibnd, const Vector< Real > &linear_ires, std::ostream &outStream=std::cout) |
Run algorithm on bound constrained problems with explicit linear constraints (Type-B). This general interface supports the use of dual optimization vector spaces, where the user does not define the dual() method. More... | |
virtual void | writeExitStatus (std::ostream &os) const |
Ptr< const AlgorithmState< Real > > | getState () const |
void | reset () |
Private Member Functions | |
void | initialize (Vector< Real > &x, const Vector< Real > &g, Objective< Real > &obj, BoundConstraint< Real > &bnd, std::ostream &outStream=std::cout) |
Private Attributes | |
Ptr< Secant< Real > > | secant_ |
Secant object (used for quasi-Newton) More... | |
ESecant | esec_ |
Secant type. More... | |
std::string | secantName_ |
Secant name. More... | |
Ptr< Krylov< Real > > | krylov_ |
Krylov solver object (used for inexact Newton) More... | |
EKrylov | ekv_ |
Krylov type. More... | |
std::string | krylovName_ |
Krylov name. More... | |
int | totalKrylov_ |
Total number of Krylov iterations per PDAS iteration. More... | |
int | iterKrylov_ |
Number of Krylov iterations (used for inexact Newton) More... | |
int | flagKrylov_ |
Termination flag for Krylov method (used for inexact Newton) More... | |
bool | useSecantHessVec_ |
Whether or not to use to a secant approximation as the Hessian. More... | |
bool | useSecantPrecond_ |
Whether or not to use a secant approximation to precondition inexact Newton. More... | |
int | maxit_ |
Maximum number of PDAS steps (default: 10) More... | |
int | iter_ |
PDAS iteration counter. More... | |
int | flag_ |
PDAS termination flag. More... | |
Real | stol_ |
PDAS minimum step size stopping tolerance (default: 1e-8) More... | |
Real | gtol_ |
PDAS gradient stopping tolerance (default: 1e-6) More... | |
Real | scale_ |
Scale for dual variables in the active set, \(c\) (default: 1) More... | |
Real | neps_ |
\(\epsilon\)-active set parameter More... | |
Real | itol_ |
Inexact Krylov tolerance. More... | |
Real | atolKrylov_ |
Absolute tolerance for Krylov solve (default: 1e-4) More... | |
Real | rtolKrylov_ |
Relative tolerance for Krylov solve (default: 1e-2) More... | |
int | maxitKrylov_ |
Maximum number of Krylov iterations (default: 100) More... | |
bool | feasible_ |
Flag whether the current iterate is feasible or not. More... | |
int | verbosity_ |
bool | writeHeader_ |
bool | hasPoly_ |
Additional Inherited Members | |
Protected Member Functions inherited from ROL::TypeB::Algorithm< Real > | |
void | initialize (const Vector< Real > &x, const Vector< Real > &g) |
Real | optimalityCriterion (const Vector< Real > &x, const Vector< Real > &g, Vector< Real > &primal, std::ostream &outStream=std::cout) const |
Protected Attributes inherited from ROL::TypeB::Algorithm< Real > | |
const Ptr< CombinedStatusTest < Real > > | status_ |
const Ptr< AlgorithmState< Real > > | state_ |
Ptr< PolyhedralProjection< Real > > | proj_ |
Provides an interface to run the projected secant algorithm.
Definition at line 59 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm | ( | ParameterList & | list, |
const Ptr< Secant< Real >> & | secant = nullPtr |
||
) |
Definition at line 51 of file ROL_TypeB_PrimalDualActiveSetAlgorithm_Def.hpp.
References ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::atolKrylov_, ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::esec_, ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::gtol_, ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::maxit_, ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::maxitKrylov_, ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::rtolKrylov_, ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::scale_, ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::secant_, ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::secantName_, ROL::TypeB::Algorithm< Real >::status_, ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::stol_, ROL::StringToESecant(), ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::useSecantHessVec_, ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::useSecantPrecond_, ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::verbosity_, and ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::writeHeader_.
|
private |
Definition at line 87 of file ROL_TypeB_PrimalDualActiveSetAlgorithm_Def.hpp.
References ROL::Objective< Real >::gradient(), ROL::Initial, ROL::TypeB::Algorithm< Real >::initialize(), ROL::StringToEKrylov(), ROL::Objective< Real >::update(), and ROL::Objective< Real >::value().
|
overridevirtual |
Run algorithm on bound constrained problems (Type-B). This general interface supports the use of dual optimization vector spaces, where the user does not define the dual() method.
Implements ROL::TypeB::Algorithm< Real >.
Definition at line 130 of file ROL_TypeB_PrimalDualActiveSetAlgorithm_Def.hpp.
References ROL::Accept, ROL::Vector< Real >::clone(), ROL::BoundConstraint< Real >::getLowerBound(), ROL::BoundConstraint< Real >::getUpperBound(), ROL::Objective< Real >::gradient(), ROL::Objective< Real >::hessVec(), ROL::BoundConstraint< Real >::isFeasible(), ROL::Vector< Real >::norm(), ROL::Vector< Real >::plus(), ROL::BoundConstraint< Real >::project(), ROL::BoundConstraint< Real >::pruneActive(), ROL::BoundConstraint< Real >::pruneInactive(), ROL::BoundConstraint< Real >::pruneLowerInactive(), ROL::BoundConstraint< Real >::pruneUpperInactive(), ROL::Vector< Real >::set(), ROL::Objective< Real >::update(), ROL::Objective< Real >::value(), ROL::TypeB::Algorithm< Real >::writeExitStatus(), and zero.
Referenced by main().
|
overridevirtual |
Print iterate header.
Reimplemented from ROL::TypeB::Algorithm< Real >.
Definition at line 304 of file ROL_TypeB_PrimalDualActiveSetAlgorithm_Def.hpp.
References ROL::CG_FLAG_SUCCESS, ROL::CG_FLAG_UNDEFINED, ROL::ECGFlagToString(), and ROL::NumberToString().
|
overridevirtual |
Print step name.
Reimplemented from ROL::TypeB::Algorithm< Real >.
Definition at line 360 of file ROL_TypeB_PrimalDualActiveSetAlgorithm_Def.hpp.
|
overridevirtual |
Print iterate status.
Reimplemented from ROL::TypeB::Algorithm< Real >.
Definition at line 373 of file ROL_TypeB_PrimalDualActiveSetAlgorithm_Def.hpp.
|
private |
Secant object (used for quasi-Newton)
Definition at line 61 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
Secant type.
Definition at line 62 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
Secant name.
Definition at line 63 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
Krylov solver object (used for inexact Newton)
Definition at line 65 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
|
private |
Krylov type.
Definition at line 66 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
|
private |
Krylov name.
Definition at line 67 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
|
private |
Total number of Krylov iterations per PDAS iteration.
Definition at line 69 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
|
private |
Number of Krylov iterations (used for inexact Newton)
Definition at line 70 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
|
private |
Termination flag for Krylov method (used for inexact Newton)
Definition at line 71 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
|
private |
Whether or not to use to a secant approximation as the Hessian.
Definition at line 73 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
Whether or not to use a secant approximation to precondition inexact Newton.
Definition at line 74 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
Maximum number of PDAS steps (default: 10)
Definition at line 76 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
PDAS iteration counter.
Definition at line 77 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
|
private |
PDAS termination flag.
Definition at line 78 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
|
private |
PDAS minimum step size stopping tolerance (default: 1e-8)
Definition at line 79 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
PDAS gradient stopping tolerance (default: 1e-6)
Definition at line 80 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
Scale for dual variables in the active set, \(c\) (default: 1)
Definition at line 81 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
\(\epsilon\)-active set parameter
Definition at line 82 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
|
private |
Inexact Krylov tolerance.
Definition at line 83 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
|
private |
Absolute tolerance for Krylov solve (default: 1e-4)
Definition at line 84 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
Relative tolerance for Krylov solve (default: 1e-2)
Definition at line 85 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
Maximum number of Krylov iterations (default: 100)
Definition at line 86 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
Flag whether the current iterate is feasible or not.
Definition at line 87 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
|
private |
Definition at line 89 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
Definition at line 90 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.
Referenced by ROL::TypeB::PrimalDualActiveSetAlgorithm< Real >::PrimalDualActiveSetAlgorithm().
|
private |
Definition at line 91 of file ROL_TypeB_PrimalDualActiveSetAlgorithm.hpp.