Ipopt  3.11.9
IpOrigIpoptNLP.hpp
Go to the documentation of this file.
1 // Copyright (C) 2004, 2010 International Business Machines and others.
2 // All Rights Reserved.
3 // This code is published under the Eclipse Public License.
4 //
5 // $Id: IpOrigIpoptNLP.hpp 2160 2012-12-26 19:14:42Z stefan $
6 //
7 // Authors: Carl Laird, Andreas Waechter IBM 2004-08-13
8 
9 #ifndef __IPORIGIPOPTNLP_HPP__
10 #define __IPORIGIPOPTNLP_HPP__
11 
12 #include "IpIpoptNLP.hpp"
13 #include "IpException.hpp"
14 #include "IpTimingStatistics.hpp"
15 
16 namespace Ipopt
17 {
18 
21  EXACT=0,
23  };
24 
29  };
30 
37  class OrigIpoptNLP : public IpoptNLP
38  {
39  public:
43  const SmartPtr<NLP>& nlp,
44  const SmartPtr<NLPScalingObject>& nlp_scaling);
45 
47  virtual ~OrigIpoptNLP();
49 
51  virtual bool Initialize(const Journalist& jnlst,
52  const OptionsList& options,
53  const std::string& prefix);
54 
58  bool init_x,
59  SmartPtr<Vector>& y_c,
60  bool init_y_c,
61  SmartPtr<Vector>& y_d,
62  bool init_y_d,
63  SmartPtr<Vector>& z_L,
64  bool init_z_L,
65  SmartPtr<Vector>& z_U,
66  bool init_z_U,
67  SmartPtr<Vector>& v_L,
68  SmartPtr<Vector>& v_U
69  );
70 
72  virtual bool GetWarmStartIterate(IteratesVector& warm_start_iterate)
73  {
74  return nlp_->GetWarmStartIterate(warm_start_iterate);
75  }
79  virtual Number f(const Vector& x);
80 
83  virtual Number f(const Vector& x, Number mu);
84 
86  virtual SmartPtr<const Vector> grad_f(const Vector& x);
87 
90  virtual SmartPtr<const Vector> grad_f(const Vector& x, Number mu);
91 
93  virtual SmartPtr<const Vector> c(const Vector& x);
94 
96  virtual SmartPtr<const Matrix> jac_c(const Vector& x);
97 
100  virtual SmartPtr<const Vector> d(const Vector& x);
101 
103  virtual SmartPtr<const Matrix> jac_d(const Vector& x);
104 
106  virtual SmartPtr<const SymMatrix> h(const Vector& x,
107  Number obj_factor,
108  const Vector& yc,
109  const Vector& yd
110  );
111 
114  virtual SmartPtr<const SymMatrix> h(const Vector& x,
115  Number obj_factor,
116  const Vector& yc,
117  const Vector& yd,
118  Number mu);
119 
124 
126  virtual SmartPtr<const Vector> x_L() const
127  {
128  return x_L_;
129  }
130 
132  virtual SmartPtr<const Matrix> Px_L() const
133  {
134  return Px_L_;
135  }
136 
138  virtual SmartPtr<const Vector> x_U() const
139  {
140  return x_U_;
141  }
142 
144  virtual SmartPtr<const Matrix> Px_U() const
145  {
146  return Px_U_;
147  }
148 
150  virtual SmartPtr<const Vector> d_L() const
151  {
152  return d_L_;
153  }
154 
156  virtual SmartPtr<const Matrix> Pd_L() const
157  {
158  return Pd_L_;
159  }
160 
162  virtual SmartPtr<const Vector> d_U() const
163  {
164  return d_U_;
165  }
166 
168  virtual SmartPtr<const Matrix> Pd_U() const
169  {
170  return Pd_U_;
171  }
172 
174  {
175  return h_space_;
176  }
178 
180  virtual void GetSpaces(SmartPtr<const VectorSpace>& x_space,
183  SmartPtr<const VectorSpace>& x_l_space,
184  SmartPtr<const MatrixSpace>& px_l_space,
185  SmartPtr<const VectorSpace>& x_u_space,
186  SmartPtr<const MatrixSpace>& px_u_space,
187  SmartPtr<const VectorSpace>& d_l_space,
188  SmartPtr<const MatrixSpace>& pd_l_space,
189  SmartPtr<const VectorSpace>& d_u_space,
190  SmartPtr<const MatrixSpace>& pd_u_space,
191  SmartPtr<const MatrixSpace>& Jac_c_space,
192  SmartPtr<const MatrixSpace>& Jac_d_space,
193  SmartPtr<const SymMatrixSpace>& Hess_lagrangian_space);
194 
197  virtual void AdjustVariableBounds(const Vector& new_x_L,
198  const Vector& new_x_U,
199  const Vector& new_d_L,
200  const Vector& new_d_U);
201 
204  virtual Index f_evals() const
205  {
206  return f_evals_;
207  }
208  virtual Index grad_f_evals() const
209  {
210  return grad_f_evals_;
211  }
212  virtual Index c_evals() const
213  {
214  return c_evals_;
215  }
216  virtual Index jac_c_evals() const
217  {
218  return jac_c_evals_;
219  }
220  virtual Index d_evals() const
221  {
222  return d_evals_;
223  }
224  virtual Index jac_d_evals() const
225  {
226  return jac_d_evals_;
227  }
228  virtual Index h_evals() const
229  {
230  return h_evals_;
231  }
233 
236  void FinalizeSolution(SolverReturn status,
237  const Vector& x, const Vector& z_L, const Vector& z_U,
238  const Vector& c, const Vector& d,
239  const Vector& y_c, const Vector& y_d,
240  Number obj_value,
241  const IpoptData* ip_data,
244  Index iter, Number obj_value,
245  Number inf_pr, Number inf_du,
246  Number mu, Number d_norm,
247  Number regularization_size,
248  Number alpha_du, Number alpha_pr,
249  Index ls_trials,
253 
257  static void RegisterOptions(SmartPtr<RegisteredOptions> roptions);
259 
262  {
263  return nlp_;
264  }
265 
268 
270  void ResetTimes();
271 
272  void PrintTimingStatistics(Journalist& jnlst,
273  EJournalLevel level,
274  EJournalCategory category) const;
275 
276  const TimedTask& f_eval_time() const
277  {
278  return f_eval_time_;
279  }
281  {
282  return grad_f_eval_time_;
283  }
284  const TimedTask& c_eval_time() const
285  {
286  return c_eval_time_;
287  }
288  const TimedTask& jac_c_eval_time() const
289  {
290  return jac_c_eval_time_;
291  }
292  const TimedTask& d_eval_time() const
293  {
294  return d_eval_time_;
295  }
296  const TimedTask& jac_d_eval_time() const
297  {
298  return jac_d_eval_time_;
299  }
300  const TimedTask& h_eval_time() const
301  {
302  return h_eval_time_;
303  }
304 
309 
310  private:
313 
316 
333 
338 
342 
345 
348 
352 
356 
360 
364 
367 
370 
373 
376 
379 
382 
385 
388 
391 
394 
398 
408  OrigIpoptNLP();
409 
411  OrigIpoptNLP(const OrigIpoptNLP&);
412 
414  void operator=(const OrigIpoptNLP&);
416 
423  void relax_bounds(Number bound_relax_factor, Vector& bounds);
427 
454 
465 
468 
479  };
480 
481 } // namespace Ipopt
482 
483 #endif
Number * x
Input: Starting point Output: Optimal solution.
CachedResults< SmartPtr< const Vector > > c_cache_
Equality constraint residuals.
SmartPtr< const Vector > orig_x_U_
Original unmodified upper bounds on x.
virtual SmartPtr< const Vector > c(const Vector &x)
Equality constraint residual.
bool jac_d_constant_
Flag indicating if we need to ask for inequality constraint Jacobians only once.
const TimedTask & f_eval_time() const
virtual SmartPtr< const Matrix > jac_d(const Vector &x)
Jacobian Matrix for inequality constraints.
Specialized CompoundVector class specifically for the algorithm iterates.
Class for all IPOPT specific calculated quantities.
CachedResults< SmartPtr< const Vector > > d_cache_
Inequality constraint residual (reformulated as equalities with slacks.
virtual SmartPtr< const Vector > d_L() const
Lower bounds on d.
void ResetTimes()
Reset the timing statistics.
SmartPtr< const Matrix > Px_L_
Permutation matrix (x_L_ -> x)
const TimedTask & jac_c_eval_time() const
bool jac_c_constant_
Flag indicating if we need to ask for equality constraint Jacobians only once.
virtual SmartPtr< const Vector > x_L() const
Lower bounds on x.
SmartPtr< const MatrixSpace > jac_d_space_
This class is used to collect timing information for a particular task.
Definition: IpTimedTask.hpp:18
SmartPtr< const MatrixSpace > px_u_space_
virtual bool GetWarmStartIterate(IteratesVector &warm_start_iterate)
Method accessing the GetWarmStartIterate of the NLP.
virtual Index f_evals() const
const TimedTask & h_eval_time() const
AlgorithmMode
enum to indicate the mode in which the algorithm is
const TimedTask & jac_d_eval_time() const
const TimedTask & c_eval_time() const
double Number
Type of all numbers.
Definition: IpTypes.hpp:17
Vector Base Class.
Definition: IpVector.hpp:47
SmartPtr< const VectorSpace > d_l_space_
SmartPtr< const SymMatrixSpace > scaled_h_space_
SmartPtr< const VectorSpace > x_u_space_
SmartPtr< const MatrixSpace > scaled_jac_c_space_
bool initialized_
Flag indicating if initialization method has been called.
EJournalLevel
Print Level Enum.
SmartPtr< const VectorSpace > x_l_space_
virtual SmartPtr< const Vector > d(const Vector &x)
Inequality constraint residual (reformulated as equalities with slacks.
CachedResults< SmartPtr< const Vector > > unscaled_x_cache_
Unscaled version of x vector.
HessianApproximationType hessian_approximation_
Flag indicating what Hessian information is to be used.
virtual ~OrigIpoptNLP()
Default destructor.
virtual void AdjustVariableBounds(const Vector &new_x_L, const Vector &new_x_U, const Vector &new_d_L, const Vector &new_d_U)
Method for adapting the variable bounds.
SmartPtr< const MatrixSpace > scaled_jac_d_space_
CachedResults< SmartPtr< const Matrix > > jac_d_cache_
Jacobian Matrix for inequality constraints (current iteration)
bool honor_original_bounds_
Flag indicating whether the primal variables should be projected back into original bounds are optimi...
virtual bool Initialize(const Journalist &jnlst, const OptionsList &options, const std::string &prefix)
Initialize - overloaded from IpoptNLP.
virtual SmartPtr< const Matrix > Pd_L() const
Permutation matrix (d_L_ -> d)
Number TotalFunctionEvaluationWallclockTime() const
Template class for Smart Pointers.
Definition: IpSmartPtr.hpp:182
This class stores a list of user set options.
virtual SmartPtr< const SymMatrix > uninitialized_h()
Provides a Hessian matrix from the correct matrix space with uninitialized values.
SolverReturn
enum for the return from the optimize algorithm (obviously we need to add more)
Definition: IpAlgTypes.hpp:22
bool IntermediateCallBack(AlgorithmMode mode, Index iter, Number obj_value, Number inf_pr, Number inf_du, Number mu, Number d_norm, Number regularization_size, Number alpha_du, Number alpha_pr, Index ls_trials, SmartPtr< const IpoptData > ip_data, SmartPtr< IpoptCalculatedQuantities > ip_cq)
SmartPtr< const Vector > x_L_
Lower bounds on x.
virtual SmartPtr< const Vector > x_U() const
Upper bounds on x.
SmartPtr< const Vector > d_L_
Lower bounds on d.
SmartPtr< NLP > nlp_
Pointer to the NLP.
virtual bool InitializeStructures(SmartPtr< Vector > &x, bool init_x, SmartPtr< Vector > &y_c, bool init_y_c, SmartPtr< Vector > &y_d, bool init_y_d, SmartPtr< Vector > &z_L, bool init_z_L, SmartPtr< Vector > &z_U, bool init_z_U, SmartPtr< Vector > &v_L, SmartPtr< Vector > &v_U)
Initialize (create) structures for the iteration data.
virtual SmartPtr< const SymMatrixSpace > HessianMatrixSpace() const
Accessor method to obtain the MatrixSpace for the Hessian matrix (or it&#39;s approximation) ...
HessianApproximationType
enumeration for the Hessian information type.
SmartPtr< const Matrix > Pd_L_
Permutation matrix (d_L_ -> d)
SmartPtr< const Journalist > jnlst_
journalist
SmartPtr< const Vector > d_U_
Upper bounds on d.
virtual Index h_evals() const
const TimedTask & grad_f_eval_time() const
virtual SmartPtr< const Matrix > Px_L() const
Permutation matrix (x_L_ -> x)
HessianApproximationSpace
enumeration for the Hessian approximation space.
SmartPtr< const MatrixSpace > pd_u_space_
Class to organize all the data required by the algorithm.
Definition: IpIpoptData.hpp:83
SmartPtr< const SymMatrixSpace > h_space_
SmartPtr< NLP > nlp()
Accessor method to the underlying NLP.
void PrintTimingStatistics(Journalist &jnlst, EJournalLevel level, EJournalCategory category) const
virtual Index d_evals() const
virtual Number f(const Vector &x)
Accessor methods for model data.
virtual SmartPtr< const Matrix > Px_U() const
Permutation matrix (x_U_ -> x.
SmartPtr< const Vector > get_unscaled_x(const Vector &x)
Method for getting the unscaled version of the x vector.
void operator=(const OrigIpoptNLP &)
Overloaded Equals Operator.
int Index
Type of all indices of vectors, matrices etc.
Definition: IpTypes.hpp:19
virtual SmartPtr< const SymMatrix > h(const Vector &x, Number obj_factor, const Vector &yc, const Vector &yd)
Hessian of the Lagrangian.
SmartPtr< const VectorSpace > d_space_
CachedResults< SmartPtr< const Matrix > > jac_c_cache_
Jacobian Matrix for equality constraints (current iteration)
SmartPtr< const Matrix > Px_U_
Permutation matrix (x_U_ -> x.
virtual Index jac_d_evals() const
void FinalizeSolution(SolverReturn status, const Vector &x, const Vector &z_L, const Vector &z_U, const Vector &c, const Vector &d, const Vector &y_c, const Vector &y_d, Number obj_value, const IpoptData *ip_data, IpoptCalculatedQuantities *ip_cq)
Solution Routines - overloaded from IpoptNLP.
virtual SmartPtr< const Vector > d_U() const
Upper bounds on d.
Number TotalFunctionEvaluationCpuTime() const
CachedResults< SmartPtr< const Vector > > grad_f_cache_
Gradient of the objective function.
CachedResults< Number > f_cache_
Objective function.
virtual void GetSpaces(SmartPtr< const VectorSpace > &x_space, SmartPtr< const VectorSpace > &c_space, SmartPtr< const VectorSpace > &d_space, SmartPtr< const VectorSpace > &x_l_space, SmartPtr< const MatrixSpace > &px_l_space, SmartPtr< const VectorSpace > &x_u_space, SmartPtr< const MatrixSpace > &px_u_space, SmartPtr< const VectorSpace > &d_l_space, SmartPtr< const MatrixSpace > &pd_l_space, SmartPtr< const VectorSpace > &d_u_space, SmartPtr< const MatrixSpace > &pd_u_space, SmartPtr< const MatrixSpace > &Jac_c_space, SmartPtr< const MatrixSpace > &Jac_d_space, SmartPtr< const SymMatrixSpace > &Hess_lagrangian_space)
Accessor method for vector/matrix spaces pointers.
This class maps the traditional NLP into something that is more useful by Ipopt.
Number TotalFunctionEvaluationSysTime() const
Class responsible for all message output.
OrigIpoptNLP()
Default Constructor.
bool hessian_constant_
Flag indicating if we need to ask for Hessian only once.
SmartPtr< const Vector > orig_x_L_
Original unmodified lower bounds on x.
static void RegisterOptions(SmartPtr< RegisteredOptions > roptions)
Called by IpoptType to register the options.
virtual Index grad_f_evals() const
virtual Index jac_c_evals() const
SmartPtr< const VectorSpace > d_u_space_
void relax_bounds(Number bound_relax_factor, Vector &bounds)
relax the bounds by a relative move of relax_bound_factor.
bool check_derivatives_for_naninf_
Flag indicating whether it is desired to check if there are Nan or Inf entries in first and second de...
bool warm_start_same_structure_
Flag indicating whether the TNLP with identical structure has already been solved before...
CachedResults< SmartPtr< const SymMatrix > > h_cache_
Hessian of the lagrangian (current iteration)
SmartPtr< const MatrixSpace > px_l_space_
SmartPtr< const VectorSpace > x_space_
Necessary Vector/Matrix spaces.
virtual SmartPtr< const Matrix > Pd_U() const
Permutation matrix (d_U_ -> d.
SmartPtr< const MatrixSpace > pd_l_space_
virtual SmartPtr< const Vector > grad_f(const Vector &x)
Gradient of the objective.
const TimedTask & d_eval_time() const
SmartPtr< const MatrixSpace > jac_c_space_
virtual Index c_evals() const
This is the abstract base class for classes that map the traditional NLP into something that is more ...
Definition: IpIpoptNLP.hpp:28
Number bound_relax_factor_
relaxation factor for the bounds
SmartPtr< const VectorSpace > c_space_
HessianApproximationSpace hessian_approximation_space_
Flag indicating in which space Hessian is to be approximated.
EJournalCategory
Category Selection Enum.
virtual SmartPtr< const Matrix > jac_c(const Vector &x)
Jacobian Matrix for equality constraints.
SmartPtr< const Vector > x_U_
Upper bounds on x.
SmartPtr< const Matrix > Pd_U_
Permutation matrix (d_U_ -> d.