#include <ROL_Reduced_Objective_SimOpt.hpp>

Inheritance diagram for ROL::Reduced_Objective_SimOpt< Real >:

Public Member Functions
	Reduced_Objective_SimOpt (const ROL::Ptr< Objective_SimOpt< Real > > &obj, const ROL::Ptr< Constraint_SimOpt< Real > > &con, const ROL::Ptr< Vector< Real > > &state, const ROL::Ptr< Vector< Real > > &control, const ROL::Ptr< Vector< Real > > &adjoint, const bool storage=true, const bool useFDhessVec=false)
	Constructor. More...

	Reduced_Objective_SimOpt (const ROL::Ptr< Objective_SimOpt< Real > > &obj, const ROL::Ptr< Constraint_SimOpt< Real > > &con, const ROL::Ptr< Vector< Real > > &state, const ROL::Ptr< Vector< Real > > &control, const ROL::Ptr< Vector< Real > > &adjoint, const ROL::Ptr< Vector< Real > > &dualstate, const ROL::Ptr< Vector< Real > > &dualcontrol, const ROL::Ptr< Vector< Real > > &dualadjoint, const bool storage=true, const bool useFDhessVec=false)
	Secondary, general constructor for use with dual optimization vector spaces where the user does not define the dual() method. More...

	Reduced_Objective_SimOpt (const ROL::Ptr< Objective_SimOpt< Real > > &obj, const ROL::Ptr< Constraint_SimOpt< Real > > &con, const ROL::Ptr< SimController< Real > > &stateStore, const ROL::Ptr< Vector< Real > > &state, const ROL::Ptr< Vector< Real > > &control, const ROL::Ptr< Vector< Real > > &adjoint, const bool storage=true, const bool useFDhessVec=false)
	Constructor. More...

	Reduced_Objective_SimOpt (const ROL::Ptr< Objective_SimOpt< Real > > &obj, const ROL::Ptr< Constraint_SimOpt< Real > > &con, const ROL::Ptr< SimController< Real > > &stateStore, const ROL::Ptr< Vector< Real > > &state, const ROL::Ptr< Vector< Real > > &control, const ROL::Ptr< Vector< Real > > &adjoint, const ROL::Ptr< Vector< Real > > &dualstate, const ROL::Ptr< Vector< Real > > &dualcontrol, const ROL::Ptr< Vector< Real > > &dualadjoint, const bool storage=true, const bool useFDhessVec=false)
	Secondary, general constructor for use with dual optimization vector spaces where the user does not define the dual() method. More...

void	update (const Vector< Real > &z, bool flag=true, int iter=-1)
	Update the SimOpt objective function and equality constraint. More...

Real	value (const Vector< Real > &z, Real &tol)
	Given \(z\in\mathcal{Z}\), evaluate the objective function \(\widehat{J}(z) = J(u(z),z)\) where \(u=u(z)\in\mathcal{U}\) solves \(e(u,z) = 0\). More...

void	gradient (Vector< Real > &g, const Vector< Real > &z, Real &tol)
	Given \(z\in\mathcal{Z}\), evaluate the gradient of the objective function \(\nabla\widehat{J}(z) = J_z(z) + c_z(u,z)^\lambda\) where \(\lambda=\lambda(u,z)\in\mathcal{C}^\) solves \(e_u(u,z)^*\lambda+J_u(u,z) = 0\). More...

void	hessVec (Vector< Real > &hv, const Vector< Real > &v, const Vector< Real > &z, Real &tol)
	Given \(z\in\mathcal{Z}\), evaluate the Hessian of the objective function \(\nabla^2\widehat{J}(z)\) in the direction \(v\in\mathcal{Z}\). More...

virtual void	precond (Vector< Real > &Pv, const Vector< Real > &v, const Vector< Real > &z, Real &tol)
	Apply a reduced Hessian preconditioner. More...

void	setParameter (const std::vector< Real > &param)

Public Member Functions inherited from ROL::Objective< Real >
virtual	~Objective ()

virtual Real	dirDeriv (const Vector< Real > &x, const Vector< Real > &d, Real &tol)
	Compute directional derivative. More...

virtual void	invHessVec (Vector< Real > &hv, const Vector< Real > &v, const Vector< Real > &x, Real &tol)
	Apply inverse Hessian approximation to vector. More...

virtual std::vector < std::vector< Real > >	checkGradient (const Vector< Real > &x, const Vector< Real > &d, const bool printToStream=true, std::ostream &outStream=std::cout, const int numSteps=ROL_NUM_CHECKDERIV_STEPS, const int order=1)
	Finite-difference gradient check. More...

virtual std::vector < std::vector< Real > >	checkGradient (const Vector< Real > &x, const Vector< Real > &g, const Vector< Real > &d, const bool printToStream=true, std::ostream &outStream=std::cout, const int numSteps=ROL_NUM_CHECKDERIV_STEPS, const int order=1)
	Finite-difference gradient check. More...

virtual std::vector < std::vector< Real > >	checkGradient (const Vector< Real > &x, const Vector< Real > &d, const std::vector< Real > &steps, const bool printToStream=true, std::ostream &outStream=std::cout, const int order=1)
	Finite-difference gradient check with specified step sizes. More...

virtual std::vector < std::vector< Real > >	checkGradient (const Vector< Real > &x, const Vector< Real > &g, const Vector< Real > &d, const std::vector< Real > &steps, const bool printToStream=true, std::ostream &outStream=std::cout, const int order=1)
	Finite-difference gradient check with specified step sizes. More...

virtual std::vector < std::vector< Real > >	checkHessVec (const Vector< Real > &x, const Vector< Real > &v, const bool printToStream=true, std::ostream &outStream=std::cout, const int numSteps=ROL_NUM_CHECKDERIV_STEPS, const int order=1)
	Finite-difference Hessian-applied-to-vector check. More...

virtual std::vector < std::vector< Real > >	checkHessVec (const Vector< Real > &x, const Vector< Real > &hv, const Vector< Real > &v, const bool printToStream=true, std::ostream &outStream=std::cout, const int numSteps=ROL_NUM_CHECKDERIV_STEPS, const int order=1)
	Finite-difference Hessian-applied-to-vector check. More...

virtual std::vector < std::vector< Real > >	checkHessVec (const Vector< Real > &x, const Vector< Real > &v, const std::vector< Real > &steps, const bool printToStream=true, std::ostream &outStream=std::cout, const int order=1)
	Finite-difference Hessian-applied-to-vector check with specified step sizes. More...

virtual std::vector < std::vector< Real > >	checkHessVec (const Vector< Real > &x, const Vector< Real > &hv, const Vector< Real > &v, const std::vector< Real > &steps, const bool printToStream=true, std::ostream &outStream=std::cout, const int order=1)
	Finite-difference Hessian-applied-to-vector check with specified step sizes. More...

virtual std::vector< Real >	checkHessSym (const Vector< Real > &x, const Vector< Real > &v, const Vector< Real > &w, const bool printToStream=true, std::ostream &outStream=std::cout)
	Hessian symmetry check. More...

virtual std::vector< Real >	checkHessSym (const Vector< Real > &x, const Vector< Real > &hv, const Vector< Real > &v, const Vector< Real > &w, const bool printToStream=true, std::ostream &outStream=std::cout)
	Hessian symmetry check. More...

Private Member Functions
void	solve_state_equation (const Vector< Real > &z, Real &tol)

void	solve_adjoint_equation (const Vector< Real > &z, Real &tol)
	Given \((u,z)\in\mathcal{U}\times\mathcal{Z}\) which solves the state equation, solve the adjoint equation \(c_u(u,z)^\lambda + J_u(u,z) = 0\) for \(\lambda=\lambda(u,z)\in\mathcal{C}^\). More...

void	solve_state_sensitivity (const Vector< Real > &v, const Vector< Real > &z, Real &tol)
	Given \((u,z)\in\mathcal{U}\times\mathcal{Z}\) which solves the state equation and a direction \(v\in\mathcal{Z}\), solve the state senstivity equation \(c_u(u,z)s + c_z(u,z)v = 0\) for \(s=u_z(z)v\in\mathcal{U}\). More...

void	solve_adjoint_sensitivity (const Vector< Real > &v, const Vector< Real > &z, Real &tol)
	Given \((u,z)\in\mathcal{U}\times\mathcal{Z}\), the adjoint variable \(\lambda\in\mathcal{C}^\), and a direction \(v\in\mathcal{Z}\), solve the adjoint sensitvity equation \(c_u(u,z)^p + J_{uu}(u,z)s + J_{uz}(u,z)v + c_{uu}(u,z)(\cdot,s)^\lambda + c_{zu}(u,z)(\cdot,v)^\lambda = 0\) for \(p = \lambda_z(u(z),z)v\in\mathcal{C}^*\). More...

Private Attributes
const ROL::Ptr < Objective_SimOpt< Real > >	obj_

const ROL::Ptr < Constraint_SimOpt< Real > >	con_

ROL::Ptr< SimController< Real > >	stateStore_

ROL::Ptr< SimController< Real > >	adjointStore_

ROL::Ptr< Vector< Real > >	state_

ROL::Ptr< Vector< Real > >	adjoint_

ROL::Ptr< Vector< Real > >	state_sens_

ROL::Ptr< Vector< Real > >	adjoint_sens_

ROL::Ptr< Vector< Real > >	dualstate_

ROL::Ptr< Vector< Real > >	dualstate1_

ROL::Ptr< Vector< Real > >	dualadjoint_

ROL::Ptr< Vector< Real > >	dualcontrol_

const bool	storage_

const bool	useFDhessVec_

bool	updateFlag_

int	updateIter_

Additional Inherited Members
Protected Member Functions inherited from ROL::Objective< Real >
const std::vector< Real >	getParameter (void) const

Detailed Description

template<class Real>
class ROL::Reduced_Objective_SimOpt< Real >

Definition at line 55 of file ROL_Reduced_Objective_SimOpt.hpp.

Constructor & Destructor Documentation

template<class Real>

ROL::Reduced_Objective_SimOpt< Real >::Reduced_Objective_SimOpt	(	const ROL::Ptr< Objective_SimOpt< Real > > &	obj,
		const ROL::Ptr< Constraint_SimOpt< Real > > &	con,
		const ROL::Ptr< Vector< Real > > &	state,
		const ROL::Ptr< Vector< Real > > &	control,
		const ROL::Ptr< Vector< Real > > &	adjoint,
		const bool	storage = `true`,
		const bool	useFDhessVec = `false`
	)

inline

Constructor.

Parameters

[in]	obj	is a pointer to a SimOpt objective function.
[in]	con	is a pointer to a SimOpt equality constraint.
[in]	state	is a pointer to a state space vector, \(\mathcal{U}\).
[in]	control	is a pointer to a optimization space vector, \(\mathcal{Z}\).
[in]	adjoint	is a pointer to a dual constraint space vector, \(\mathcal{C}^*\).
[in]	storage	is a flag whether or not to store computed states and adjoints.
[in]	useFDhessVec	is a flag whether or not to use a finite-difference Hessian approximation.

Definition at line 171 of file ROL_Reduced_Objective_SimOpt.hpp.

template<class Real>

ROL::Reduced_Objective_SimOpt< Real >::Reduced_Objective_SimOpt	(	const ROL::Ptr< Objective_SimOpt< Real > > &	obj,
		const ROL::Ptr< Constraint_SimOpt< Real > > &	con,
		const ROL::Ptr< Vector< Real > > &	state,
		const ROL::Ptr< Vector< Real > > &	control,
		const ROL::Ptr< Vector< Real > > &	adjoint,
		const ROL::Ptr< Vector< Real > > &	dualstate,
		const ROL::Ptr< Vector< Real > > &	dualcontrol,
		const ROL::Ptr< Vector< Real > > &	dualadjoint,
		const bool	storage = `true`,
		const bool	useFDhessVec = `false`
	)

inline

Secondary, general constructor for use with dual optimization vector spaces where the user does not define the dual() method.

Parameters

[in]	obj	is a pointer to a SimOpt objective function.
[in]	con	is a pointer to a SimOpt equality constraint.
[in]	state	is a pointer to a state space vector, \(\mathcal{U}\).
[in]	control	is a pointer to a optimization space vector, \(\mathcal{Z}\).
[in]	adjoint	is a pointer to a dual constraint space vector, \(\mathcal{C}^*\).
[in]	dualstate	is a pointer to a dual state space vector, \(\mathcal{U}^*\).
[in]	dualadjoint	is a pointer to a constraint space vector, \(\mathcal{C}\).
[in]	storage	is a flag whether or not to store computed states and adjoints.
[in]	useFDhessVec	is a flag whether or not to use a finite-difference Hessian approximation.

Definition at line 207 of file ROL_Reduced_Objective_SimOpt.hpp.

template<class Real>

ROL::Reduced_Objective_SimOpt< Real >::Reduced_Objective_SimOpt	(	const ROL::Ptr< Objective_SimOpt< Real > > &	obj,
		const ROL::Ptr< Constraint_SimOpt< Real > > &	con,
		const ROL::Ptr< SimController< Real > > &	stateStore,
		const ROL::Ptr< Vector< Real > > &	state,
		const ROL::Ptr< Vector< Real > > &	control,
		const ROL::Ptr< Vector< Real > > &	adjoint,
		const bool	storage = `true`,
		const bool	useFDhessVec = `false`
	)

inline

Constructor.

Parameters

[in]	obj	is a pointer to a SimOpt objective function.
[in]	con	is a pointer to a SimOpt equality constraint.
[in]	stateStore	is a pointer to a SimController object.
[in]	state	is a pointer to a state space vector, \(\mathcal{U}\).
[in]	control	is a pointer to a optimization space vector, \(\mathcal{Z}\).
[in]	adjoint	is a pointer to a dual constraint space vector, \(\mathcal{C}^*\).
[in]	storage	is a flag whether or not to store computed states and adjoints.
[in]	useFDhessVec	is a flag whether or not to use a finite-difference Hessian approximation.

Definition at line 244 of file ROL_Reduced_Objective_SimOpt.hpp.

References ROL::Reduced_Objective_SimOpt< Real >::adjoint_, ROL::Reduced_Objective_SimOpt< Real >::adjoint_sens_, ROL::Reduced_Objective_SimOpt< Real >::adjointStore_, ROL::Reduced_Objective_SimOpt< Real >::dualadjoint_, ROL::Reduced_Objective_SimOpt< Real >::dualcontrol_, ROL::Reduced_Objective_SimOpt< Real >::dualstate1_, ROL::Reduced_Objective_SimOpt< Real >::dualstate_, ROL::Reduced_Objective_SimOpt< Real >::state_, and ROL::Reduced_Objective_SimOpt< Real >::state_sens_.

template<class Real>

ROL::Reduced_Objective_SimOpt< Real >::Reduced_Objective_SimOpt	(	const ROL::Ptr< Objective_SimOpt< Real > > &	obj,
		const ROL::Ptr< Constraint_SimOpt< Real > > &	con,
		const ROL::Ptr< SimController< Real > > &	stateStore,
		const ROL::Ptr< Vector< Real > > &	state,
		const ROL::Ptr< Vector< Real > > &	control,
		const ROL::Ptr< Vector< Real > > &	adjoint,
		const ROL::Ptr< Vector< Real > > &	dualstate,
		const ROL::Ptr< Vector< Real > > &	dualcontrol,
		const ROL::Ptr< Vector< Real > > &	dualadjoint,
		const bool	storage = `true`,
		const bool	useFDhessVec = `false`
	)

inline

Secondary, general constructor for use with dual optimization vector spaces where the user does not define the dual() method.

Parameters

[in]	obj	is a pointer to a SimOpt objective function.
[in]	con	is a pointer to a SimOpt equality constraint.
[in]	stateStore	is a pointer to a SimController object.
[in]	state	is a pointer to a state space vector, \(\mathcal{U}\).
[in]	control	is a pointer to a optimization space vector, \(\mathcal{Z}\).
[in]	adjoint	is a pointer to a dual constraint space vector, \(\mathcal{C}^*\).
[in]	dualstate	is a pointer to a dual state space vector, \(\mathcal{U}^*\).
[in]	dualadjoint	is a pointer to a constraint space vector, \(\mathcal{C}\).
[in]	storage	is a flag whether or not to store computed states and adjoints.
[in]	useFDhessVec	is a flag whether or not to use a finite-difference Hessian approximation.

Definition at line 281 of file ROL_Reduced_Objective_SimOpt.hpp.

References ROL::Reduced_Objective_SimOpt< Real >::adjoint_, ROL::Reduced_Objective_SimOpt< Real >::adjoint_sens_, ROL::Reduced_Objective_SimOpt< Real >::adjointStore_, ROL::Reduced_Objective_SimOpt< Real >::dualadjoint_, ROL::Reduced_Objective_SimOpt< Real >::dualcontrol_, ROL::Reduced_Objective_SimOpt< Real >::dualstate1_, ROL::Reduced_Objective_SimOpt< Real >::dualstate_, ROL::Reduced_Objective_SimOpt< Real >::state_, and ROL::Reduced_Objective_SimOpt< Real >::state_sens_.

Member Function Documentation

template<class Real>

void ROL::Reduced_Objective_SimOpt< Real >::solve_state_equation	(	const Vector< Real > &	z,
		Real &	tol
	)

inlineprivate

Definition at line 80 of file ROL_Reduced_Objective_SimOpt.hpp.

References ROL::Reduced_Objective_SimOpt< Real >::con_, ROL::Reduced_Objective_SimOpt< Real >::dualadjoint_, ROL::Reduced_Objective_SimOpt< Real >::obj_, ROL::Reduced_Objective_SimOpt< Real >::state_, ROL::Reduced_Objective_SimOpt< Real >::stateStore_, ROL::Reduced_Objective_SimOpt< Real >::storage_, ROL::Reduced_Objective_SimOpt< Real >::updateFlag_, and ROL::Reduced_Objective_SimOpt< Real >::updateIter_.

Referenced by ROL::Reduced_Objective_SimOpt< Real >::gradient(), ROL::Reduced_Objective_SimOpt< Real >::hessVec(), and ROL::Reduced_Objective_SimOpt< Real >::value().

template<class Real>

void ROL::Reduced_Objective_SimOpt< Real >::solve_adjoint_equation	(	const Vector< Real > &	z,
		Real &	tol
	)

inlineprivate

Given \((u,z)\in\mathcal{U}\times\mathcal{Z}\) which solves the state equation, solve the adjoint equation \(c_u(u,z)^*\lambda + J_u(u,z) = 0\) for \(\lambda=\lambda(u,z)\in\mathcal{C}^*\).

Definition at line 107 of file ROL_Reduced_Objective_SimOpt.hpp.

References ROL::Reduced_Objective_SimOpt< Real >::adjoint_, ROL::Reduced_Objective_SimOpt< Real >::adjointStore_, ROL::Reduced_Objective_SimOpt< Real >::con_, ROL::Reduced_Objective_SimOpt< Real >::dualstate_, ROL::Reduced_Objective_SimOpt< Real >::obj_, ROL::Reduced_Objective_SimOpt< Real >::state_, and ROL::Reduced_Objective_SimOpt< Real >::storage_.

Referenced by ROL::Reduced_Objective_SimOpt< Real >::gradient(), and ROL::Reduced_Objective_SimOpt< Real >::hessVec().

template<class Real>

void ROL::Reduced_Objective_SimOpt< Real >::solve_state_sensitivity	(	const Vector< Real > &	v,
		const Vector< Real > &	z,
		Real &	tol
	)

inlineprivate

Given \((u,z)\in\mathcal{U}\times\mathcal{Z}\) which solves the state equation and a direction \(v\in\mathcal{Z}\), solve the state senstivity equation \(c_u(u,z)s + c_z(u,z)v = 0\) for \(s=u_z(z)v\in\mathcal{U}\).

Definition at line 131 of file ROL_Reduced_Objective_SimOpt.hpp.

References ROL::Reduced_Objective_SimOpt< Real >::con_, ROL::Reduced_Objective_SimOpt< Real >::dualadjoint_, ROL::Reduced_Objective_SimOpt< Real >::state_, and ROL::Reduced_Objective_SimOpt< Real >::state_sens_.

Referenced by ROL::Reduced_Objective_SimOpt< Real >::hessVec().

template<class Real>

void ROL::Reduced_Objective_SimOpt< Real >::solve_adjoint_sensitivity	(	const Vector< Real > &	v,
		const Vector< Real > &	z,
		Real &	tol
	)

inlineprivate

Given \((u,z)\in\mathcal{U}\times\mathcal{Z}\), the adjoint variable \(\lambda\in\mathcal{C}^*\), and a direction \(v\in\mathcal{Z}\), solve the adjoint sensitvity equation \(c_u(u,z)^*p + J_{uu}(u,z)s + J_{uz}(u,z)v + c_{uu}(u,z)(\cdot,s)^*\lambda + c_{zu}(u,z)(\cdot,v)^*\lambda = 0\) for \(p = \lambda_z(u(z),z)v\in\mathcal{C}^*\).

Definition at line 145 of file ROL_Reduced_Objective_SimOpt.hpp.

References ROL::Reduced_Objective_SimOpt< Real >::adjoint_, ROL::Reduced_Objective_SimOpt< Real >::adjoint_sens_, ROL::Reduced_Objective_SimOpt< Real >::con_, ROL::Reduced_Objective_SimOpt< Real >::dualstate1_, ROL::Reduced_Objective_SimOpt< Real >::dualstate_, ROL::Reduced_Objective_SimOpt< Real >::obj_, ROL::Reduced_Objective_SimOpt< Real >::state_, and ROL::Reduced_Objective_SimOpt< Real >::state_sens_.

Referenced by ROL::Reduced_Objective_SimOpt< Real >::hessVec().

template<class Real>

void ROL::Reduced_Objective_SimOpt< Real >::update	(	const Vector< Real > &	z,
		bool	flag = `true`,
		int	iter = `-1`
	)

inlinevirtual

Update the SimOpt objective function and equality constraint.

Reimplemented from ROL::Objective< Real >.

Definition at line 309 of file ROL_Reduced_Objective_SimOpt.hpp.

References ROL::Reduced_Objective_SimOpt< Real >::adjointStore_, ROL::Reduced_Objective_SimOpt< Real >::stateStore_, ROL::Reduced_Objective_SimOpt< Real >::updateFlag_, and ROL::Reduced_Objective_SimOpt< Real >::updateIter_.

template<class Real>

Real ROL::Reduced_Objective_SimOpt< Real >::value	(	const Vector< Real > &	z,
		Real &	tol
	)

inlinevirtual

Given \(z\in\mathcal{Z}\), evaluate the objective function \(\widehat{J}(z) = J(u(z),z)\) where \(u=u(z)\in\mathcal{U}\) solves \(e(u,z) = 0\).

Implements ROL::Objective< Real >.

Definition at line 320 of file ROL_Reduced_Objective_SimOpt.hpp.

References ROL::Reduced_Objective_SimOpt< Real >::obj_, ROL::Reduced_Objective_SimOpt< Real >::solve_state_equation(), and ROL::Reduced_Objective_SimOpt< Real >::state_.

template<class Real>

void ROL::Reduced_Objective_SimOpt< Real >::gradient	(	Vector< Real > &	g,
		const Vector< Real > &	z,
		Real &	tol
	)

inlinevirtual

Given \(z\in\mathcal{Z}\), evaluate the gradient of the objective function \(\nabla\widehat{J}(z) = J_z(z) + c_z(u,z)^*\lambda\) where \(\lambda=\lambda(u,z)\in\mathcal{C}^*\) solves \(e_u(u,z)^*\lambda+J_u(u,z) = 0\).

Reimplemented from ROL::Objective< Real >.

Definition at line 332 of file ROL_Reduced_Objective_SimOpt.hpp.

References ROL::Reduced_Objective_SimOpt< Real >::adjoint_, ROL::Reduced_Objective_SimOpt< Real >::con_, ROL::Reduced_Objective_SimOpt< Real >::dualcontrol_, ROL::Reduced_Objective_SimOpt< Real >::obj_, ROL::Vector< Real >::plus(), ROL::Reduced_Objective_SimOpt< Real >::solve_adjoint_equation(), ROL::Reduced_Objective_SimOpt< Real >::solve_state_equation(), and ROL::Reduced_Objective_SimOpt< Real >::state_.

template<class Real>