Ipopt  3.12.12
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros
IpLimMemQuasiNewtonUpdater.hpp
Go to the documentation of this file.
1 // Copyright (C) 2005, 2010 International Business Machines and others.
2 // All Rights Reserved.
3 // This code is published under the Eclipse Public License.
4 //
5 // $Id: IpLimMemQuasiNewtonUpdater.hpp 2269 2013-05-05 11:32:40Z stefan $
6 //
7 // Authors: Andreas Waechter IBM 2005-12-26
8 
9 #ifndef __IPLIMMEMQUASINEWTONUPDATER_HPP__
10 #define __IPLIMMEMQUASINEWTONUPDATER_HPP__
11 
12 #include "IpHessianUpdater.hpp"
14 #include "IpMultiVectorMatrix.hpp"
15 #include "IpDenseVector.hpp"
16 #include "IpDenseGenMatrix.hpp"
17 #include "IpDenseSymMatrix.hpp"
18 
19 namespace Ipopt
20 {
21 
26  {
27  public:
31  LimMemQuasiNewtonUpdater(bool update_for_resto);
32 
35  {}
37 
39  virtual bool InitializeImpl(const OptionsList& options,
40  const std::string& prefix);
41 
44  virtual void UpdateHessian();
45 
48  static void RegisterOptions(SmartPtr<RegisteredOptions> roptions);
50 
51  private:
62 
66 
69 
76  {
77  BFGS=0,
79  };
84  {
90  };
107 
113  const bool update_for_resto_;
131 
135 
221 
228  bool CheckSkippingBFGS(Vector& s_new, Vector& y_new);
233  bool UpdateInternalData(const Vector& s_new, const Vector& y_new,
234  SmartPtr<Vector> ypart_new);
241  const Vector& v_new);
248  Number v_new);
257  const MultiVectorMatrix& S,
258  const MultiVectorMatrix& Y);
266  const MultiVectorMatrix& S);
275  const MultiVectorMatrix& S,
276  const MultiVectorMatrix& DRS);
277 
282  void ShiftMultiVector(SmartPtr<MultiVectorMatrix>& V, const Vector& v_new);
296  const MultiVectorMatrix& S,
297  const MultiVectorMatrix& Y);
304  const MultiVectorMatrix& S);
311  const MultiVectorMatrix& S,
312  const MultiVectorMatrix& DRS);
315  void RecalcY(Number eta, const Vector& DR_x,
317  MultiVectorMatrix& Ypart,
320  void RecalcD(MultiVectorMatrix& S,
324  void RecalcL(MultiVectorMatrix& S,
336  bool SplitEigenvalues(DenseGenMatrix& Q, const DenseVector& E,
337  SmartPtr<DenseGenMatrix>& Qminus,
338  SmartPtr<DenseGenMatrix>& Qplus);
351  void SetW();
353 
354  };
355 
356 } // namespace Ipopt
357 
358 #endif
SmartPtr< DenseGenMatrix > L_old_
Matrix L_k for compact formulation from last update (backup).
SmartPtr< DenseSymMatrix > STDRS_
For efficient implementation, we store the S^T S DR * S.
SmartPtr< Vector > B0_
First term (starting matrix) for the approximation.
SmartPtr< DenseSymMatrix > STDRS_old_
For efficient implementation, we store the S^T S DR * S.
Number sigma_old_
First term (starting matrix) for the approximation.
virtual bool InitializeImpl(const OptionsList &options, const std::string &prefix)
overloaded from AlgorithmStrategyObject
SmartPtr< MultiVectorMatrix > V_old_
V in LowRankUpdateMatrix from last update (backup)
Index curr_lm_memory_old_
current size of limited memory
SmartPtr< const Matrix > last_jac_c_
Jacobian for equality constraints w.r.t x at x_last.
SmartPtr< DenseVector > D_old_
Diagonal elements D_k for compact formulation from last update (backup).
SmartPtr< Vector > B0_old_
First term (starting matrix) for the approximation (backup).
LMInitialization limited_memory_initialization_
How to choose B0 in the low-rank update.
Index limited_memory_max_skipping_
Number of successive iterations of skipped updates after which the approximation is reset...
Number limited_memory_init_val_
Value of B0 (as this multiple of the identity in certain situations.)
LMUpdateType
enumeration for the Hessian update type.
Number sigma_safe_min_
Minimal safeguard value for sigma.
SmartPtr< const Vector > curr_red_DR_x_
Current DR_x scaling factors in the restoration phase objective function in the smaller space for the...
double Number
Type of all numbers.
Definition: IpTypes.hpp:17
Vector Base Class.
Definition: IpVector.hpp:47
Index limited_memory_max_history_
Size of memory for limited memory update.
void ShiftMultiVector(SmartPtr< MultiVectorMatrix > &V, const Vector &v_new)
Given a MutliVector V, get rid of the first column, shift all other columns to the left...
virtual void UpdateHessian()
Update the Hessian based on the current information in IpData.
Dense Vector Implementation.
void RecalcY(Number eta, const Vector &DR_x, MultiVectorMatrix &S, MultiVectorMatrix &Ypart, SmartPtr< MultiVectorMatrix > &Y)
Method for recomputing Y from scratch, using Ypart (only for restoration phase)
LMUpdateType limited_memory_update_type_
Type of Hessian update.
SmartPtr< const LowRankUpdateSymMatrixSpace > h_space_
Matrix space for the low-rank Hessian approximation.
void ReleaseInternalDataBackup()
Release anything that we allocated for StoreInternalDataBackup and is no longer needed.
SmartPtr< const Vector > last_x_
Primal variables x from most recent update.
TaggedObject::Tag curr_DR_x_tag_
Tag for curr_DR_x_.
Number sigma_safe_max_
Maximal safeguard value for sigma.
SmartPtr< MultiVectorMatrix > S_old_
s pairs for the recent iterations (backup)
void RecalcD(MultiVectorMatrix &S, MultiVectorMatrix &Y, SmartPtr< DenseVector > &D)
Method for recomputing D from S and Y.
bool SplitEigenvalues(DenseGenMatrix &Q, const DenseVector &E, SmartPtr< DenseGenMatrix > &Qminus, SmartPtr< DenseGenMatrix > &Qplus)
Split the eigenvectors into negative and positive ones.
bool limited_memory_special_for_resto_
Flag indicating if Hessian approximation should be done in a special manner for the restoration phase...
bool CheckSkippingBFGS(Vector &s_new, Vector &y_new)
Method deciding whether the BFGS update should be skipped.
Template class for Smart Pointers.
Definition: IpSmartPtr.hpp:172
This class stores a list of user set options.
bool SdotS_uptodate_
Flag indicating whether SdotS_ is update to date from most recent update.
void AugmentSdotSMatrix(SmartPtr< DenseSymMatrix > &V, const MultiVectorMatrix &S)
Given a DenseSymMatrix V, create a new DenseGenMatrixSpace with one more dimension, and return V as a member of that space, consisting of all previous elements, and in addition elements s_i^Ts_j for the new entries, where s are the vectors in the MultiVector S.
void RecalcL(MultiVectorMatrix &S, MultiVectorMatrix &Y, SmartPtr< DenseGenMatrix > &L)
Method for recomputing L from S and Y.
bool UpdateInternalData(const Vector &s_new, const Vector &y_new, SmartPtr< Vector > ypart_new)
Update the internal data, such as the S, Y, L, D etc matrices and vectors that are required for compu...
bool SdotS_uptodate_old_
Flag indicating whether SdotS_ is update to date from most recent update (backup).
void StoreInternalDataBackup()
Store a copy of the pointers to the internal data (S, Y, D, L, SdotS, curr_lm_memory) This is called ...
void AugmentMultiVector(SmartPtr< MultiVectorMatrix > &V, const Vector &v_new)
Given a MutliVector V, create a new MultiVectorSpace with one more column, and return V as a member o...
virtual ~LimMemQuasiNewtonUpdater()
Default destructor.
Index curr_lm_memory_
current size of limited memory
void ShiftDenseVector(SmartPtr< DenseVector > &V, Number v_new)
Given a DenseVector V, get rid of the first element, shift all other elements one position to the top...
Class for Matrices with few columns that consists of Vectors.
SmartPtr< MultiVectorMatrix > U_
U in LowRankUpdateMatrix from last update.
Implementation of the HessianUpdater for limit-memory quasi-Newton approximation of the Lagrangian He...
SmartPtr< MultiVectorMatrix > V_
V in LowRankUpdateMatrix from last update.
int Index
Type of all indices of vectors, matrices etc.
Definition: IpTypes.hpp:19
SmartPtr< const Vector > last_grad_f_
Gradient of objective function w.r.t.
SmartPtr< MultiVectorMatrix > S_
s pairs for the recent iterations
SmartPtr< DenseVector > D_
Diagonal elements D_k for compact formulation from last update.
void AugmentDenseVector(SmartPtr< DenseVector > &V, Number v_new)
Given a DenseVector V, create a new DenseVectorSpace with one more row, and return V as a member of t...
SmartPtr< const Matrix > last_jac_d_
Jacobian for inequality constraints w.r.t x at x_last.
void ShiftLMatrix(SmartPtr< DenseGenMatrix > &V, const MultiVectorMatrix &S, const MultiVectorMatrix &Y)
Given a strictly-lower triangular square DenseGenMatrix V, shift everything one row and column up...
void AugmentLMatrix(SmartPtr< DenseGenMatrix > &V, const MultiVectorMatrix &S, const MultiVectorMatrix &Y)
Given a strictly-lower triangular square DenseGenMatrix V, create a new DenseGenMatrixSpace with one ...
SmartPtr< DenseSymMatrix > SdotS_
For efficient implementation, we store the pairwise products for s&#39;s.
LMInitialization
enumeration for the Hessian initialization.
SmartPtr< MultiVectorMatrix > Ypart_
For restoration phase update: Y without the quadratic objective function part.
SmartPtr< MultiVectorMatrix > U_old_
U in LowRankUpdateMatrix from last update (backup)
SmartPtr< MultiVectorMatrix > DRS_
DR * S (only for restoration phase)
Index lm_skipped_iter_
Counter for successive iterations in which the update was skipped.
unsigned int Tag
Type for the Tag values.
Number sigma_
First term (starting matrix) for the approximation.
void SetW()
Set the W field in IpData based on the current values of B0_, V_, and U_.
Number last_eta_
Most recent value for eta in the restoration phase objective function (only for update_for_resto_ = t...
Class for dense general matrices.
SmartPtr< DenseSymMatrix > SdotS_old_
For efficient implementation, we store the pairwise products for s&#39;s (backup).
void operator=(const LimMemQuasiNewtonUpdater &)
Overloaded Equals Operator.
Number curr_eta_
Current value of weighing factor eta in the restoration phase objective function (only for update_for...
SmartPtr< MultiVectorMatrix > Ypart_old_
For restoration phase update: Y without the quadratic objective function part (backup) ...
static void RegisterOptions(SmartPtr< RegisteredOptions > roptions)
Methods for OptionsList.
void RestoreInternalDataBackup()
Restore the copy of the pointers to the internal data most recently stored with StoreInternalDataBack...
SmartPtr< MultiVectorMatrix > Y_old_
y pairs for the recent iterations.
void ShiftSdotSMatrix(SmartPtr< DenseSymMatrix > &V, const MultiVectorMatrix &S)
Given a DenseSymMatrix V, shift everything up one row and column, and fill the new entries as s_i^Ts_...
SmartPtr< DenseGenMatrix > L_
Matrix L_k for compact formulation from last update.
Abstract base class for objects responsible for updating the Hessian information. ...
SmartPtr< const Vector > curr_DR_x_
Current DR_x scaling factors in the restoration phase objective function (only for update_for_resto_ ...
LimMemQuasiNewtonUpdater(bool update_for_resto)
Default Constructor.
void AugmentSTDRSMatrix(SmartPtr< DenseSymMatrix > &V, const MultiVectorMatrix &S, const MultiVectorMatrix &DRS)
Given a DenseSymMatrix V, create a new DenseGenMatrixSpace with one more dimension, and return V as a member of that space, consisting of all previous elements, and in addition elements s_i^TDRs_j for the new entries, where s are the vectors in the MultiVector S, and DRs are the vectors in DRS.
const bool update_for_resto_
Flag indicating if the update is to be done for the original NLP or for the restoration phase NLP...
SmartPtr< MultiVectorMatrix > DRS_old_
DR * S (only for restoration phase) (backup)
void ShiftSTDRSMatrix(SmartPtr< DenseSymMatrix > &V, const MultiVectorMatrix &S, const MultiVectorMatrix &DRS)
Given a DenseSymMatrix V, shift everything up one row and column, and fill the new entries as s_i^TDR...
SmartPtr< MultiVectorMatrix > Y_
y pairs for the recent iterations.