Provides an interface to run the projected secant algorithm. More...

#include <ROL_TypeB_NewtonKrylovAlgorithm.hpp>

Inheritance diagram for ROL::TypeB::NewtonKrylovAlgorithm< Real >:

Classes
class	HessianPNK

class	PrecondPNK

Public Member Functions
	NewtonKrylovAlgorithm (ParameterList &list, const Ptr< Secant< Real >> &secant=nullPtr)

	NewtonKrylovAlgorithm (ParameterList &list, const Ptr< Krylov< Real >> &krylov, 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	run (Problem< Real > &problem, std::ostream &outStream=std::cout) override
	Run algorithm on bound constrained problems (Type-B). This is the primary Type-B interface. More...

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) override
	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...

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 (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, 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)

void	parseParameterList (ParameterList &list)

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	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 line search steps (default: 20) More...

Real	alpha0_
	Initial line search step size (default: 1.0) More...

Real	rhodec_
	Backtracking rate (default: 0.5) More...

Real	c1_
	Sufficient Decrease Parameter (default: 1e-4) More...

bool	useralpha_
	Flag to use user-defined initial step size (default: false) More...

bool	usePrevAlpha_
	Flag to use previous step size as next initial step size (default: false) More...

int	ls_nfval_

int	verbosity_

bool	writeHeader_

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_

Detailed Description

template<typename Real>
class ROL::TypeB::NewtonKrylovAlgorithm< Real >

Provides an interface to run the projected secant algorithm.

Definition at line 59 of file ROL_TypeB_NewtonKrylovAlgorithm.hpp.

Constructor & Destructor Documentation

template<typename Real >

ROL::TypeB::NewtonKrylovAlgorithm< Real >::NewtonKrylovAlgorithm	(	ParameterList &	list,
		const Ptr< Secant< Real >> &	secant = `nullPtr`
	)

Definition at line 51 of file ROL_TypeB_NewtonKrylovAlgorithm_Def.hpp.

References ROL::TypeB::NewtonKrylovAlgorithm< Real >::ekv_, ROL::TypeB::NewtonKrylovAlgorithm< Real >::esec_, ROL::TypeB::NewtonKrylovAlgorithm< Real >::krylov_, ROL::TypeB::NewtonKrylovAlgorithm< Real >::krylovName_, ROL::TypeB::NewtonKrylovAlgorithm< Real >::parseParameterList(), ROL::TypeB::NewtonKrylovAlgorithm< Real >::secant_, ROL::TypeB::NewtonKrylovAlgorithm< Real >::secantName_, ROL::StringToEKrylov(), and ROL::StringToESecant().

template<typename Real >

ROL::TypeB::NewtonKrylovAlgorithm< Real >::NewtonKrylovAlgorithm	(	ParameterList &	list,
		const Ptr< Krylov< Real >> &	krylov,
		const Ptr< Secant< Real >> &	secant = `nullPtr`
	)

Definition at line 72 of file ROL_TypeB_NewtonKrylovAlgorithm_Def.hpp.

References ROL::TypeB::NewtonKrylovAlgorithm< Real >::esec_, ROL::TypeB::NewtonKrylovAlgorithm< Real >::krylovName_, ROL::TypeB::NewtonKrylovAlgorithm< Real >::parseParameterList(), ROL::TypeB::NewtonKrylovAlgorithm< Real >::secant_, ROL::TypeB::NewtonKrylovAlgorithm< Real >::secantName_, and ROL::StringToESecant().

Member Function Documentation

template<typename Real >

void ROL::TypeB::NewtonKrylovAlgorithm< Real >::initialize	(	Vector< Real > &	x,
		const Vector< Real > &	g,
		Objective< Real > &	obj,
		BoundConstraint< Real > &	bnd,
		std::ostream &	outStream = `std::cout`
	)

private

Definition at line 115 of file ROL_TypeB_NewtonKrylovAlgorithm_Def.hpp.

References ROL::Objective< Real >::gradient(), ROL::Initial, ROL::TypeB::Algorithm< Real >::initialize(), ROL::Objective< Real >::update(), and ROL::Objective< Real >::value().

template<typename Real >

void ROL::TypeB::NewtonKrylovAlgorithm< Real >::parseParameterList ( ParameterList & list )

private

Definition at line 93 of file ROL_TypeB_NewtonKrylovAlgorithm_Def.hpp.

Referenced by ROL::TypeB::NewtonKrylovAlgorithm< Real >::NewtonKrylovAlgorithm().

template<typename Real >

void ROL::TypeB::NewtonKrylovAlgorithm< Real >::run	(	Vector< Real > &	x,
		const Vector< Real > &	g,
		Objective< Real > &	obj,
		BoundConstraint< Real > &	bnd,
		std::ostream &	outStream = `std::cout`
	)

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 144 of file ROL_TypeB_NewtonKrylovAlgorithm_Def.hpp.

References ROL::Accept, ROL::Vector< Real >::axpy(), ROL::Vector< Real >::clone(), ROL::Objective< Real >::gradient(), ROL::Vector< Real >::set(), ROL::Trial, ROL::Objective< Real >::update(), ROL::Objective< Real >::value(), and ROL::TypeB::Algorithm< Real >::writeExitStatus().

template<typename Real >

void ROL::TypeB::NewtonKrylovAlgorithm< Real >::run	(	Problem< Real > &	problem,
		std::ostream &	outStream = `std::cout`
	)

overridevirtual

Run algorithm on bound constrained problems (Type-B). This is the primary Type-B interface.

Reimplemented from ROL::TypeB::Algorithm< Real >.

Definition at line 247 of file ROL_TypeB_NewtonKrylovAlgorithm_Def.hpp.

References ROL::Problem< Real >::getPolyhedralProjection(), and ROL::TypeB::Algorithm< Real >::run().

template<typename Real >

void ROL::TypeB::NewtonKrylovAlgorithm< Real >::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`
	)

overridevirtual

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.

Reimplemented from ROL::TypeB::Algorithm< Real >.

Definition at line 258 of file ROL_TypeB_NewtonKrylovAlgorithm_Def.hpp.

template<typename Real >

void ROL::TypeB::NewtonKrylovAlgorithm< Real >::writeHeader ( std::ostream & os ) const

overridevirtual

Print iterate header.

Reimplemented from ROL::TypeB::Algorithm< Real >.

Definition at line 270 of file ROL_TypeB_NewtonKrylovAlgorithm_Def.hpp.

References ROL::CG_FLAG_SUCCESS, ROL::CG_FLAG_UNDEFINED, ROL::ECGFlagToString(), and ROL::NumberToString().

template<typename Real >

void ROL::TypeB::NewtonKrylovAlgorithm< Real >::writeName ( std::ostream & os ) const

overridevirtual