IpAugRestoSystemSolver.hpp

Go to the documentation of this file.

// Copyright (C) 2004, 2006 International Business Machines and others.
// All Rights Reserved.
// This code is published under the Eclipse Public License.
//
// Authors:  Carl Laird, Andreas Waechter     IBM    2004-08-13
 
#ifndef __IP_AUGRESTOSYSTEMSOLVER_HPP__
#define __IP_AUGRESTOSYSTEMSOLVER_HPP__
 
#include "IpAugSystemSolver.hpp"
 
namespace Ipopt
{
 
class AugRestoSystemSolver: public AugSystemSolver
{
public:
 
   AugRestoSystemSolver(
      AugSystemSolver& orig_aug_solver,
      bool             skip_orig_aug_solver_init = true
   );
 
   virtual ~AugRestoSystemSolver();
 
   bool InitializeImpl(
      const OptionsList& options,
      const std::string& prefix
   );
 
   virtual ESymSolverStatus Solve(
      const SymMatrix* W,
      Number           W_factor,
      const Vector*    D_x,
      Number           delta_x,
      const Vector*    D_s,
      Number           delta_s,
      const Matrix*    J_c,
      const Vector*    D_c,
      Number           delta_c,
      const Matrix*    J_d,
      const Vector*    D_d,
      Number           delta_d,
      const Vector&    rhs_x,
      const Vector&    rhs_s,
      const Vector&    rhs_c,
      const Vector&    rhs_d,
      Vector&          sol_x,
      Vector&          sol_s,
      Vector&          sol_c,
      Vector&          sol_d,
      bool             check_NegEVals,
      Index            numberOfNegEVals
   );
 
   virtual Index NumberOfNegEVals() const
   {
      return orig_aug_solver_->NumberOfNegEVals();
   }
 
   virtual bool ProvidesInertia() const
   {
      return orig_aug_solver_->ProvidesInertia();
   }
 
   virtual bool IncreaseQuality()
   {
      return orig_aug_solver_->IncreaseQuality();
   }
 
private:
 
   AugRestoSystemSolver();
 
   AugRestoSystemSolver(
      const AugRestoSystemSolver&
   );
 
   void operator=(
      const AugRestoSystemSolver&
   );
 
   CachedResults<SmartPtr<Vector> > neg_omega_c_plus_D_c_cache_;
   CachedResults<SmartPtr<Vector> > neg_omega_d_plus_D_d_cache_;
   CachedResults<SmartPtr<Vector> > sigma_tilde_n_c_inv_cache_;
   CachedResults<SmartPtr<Vector> > sigma_tilde_p_c_inv_cache_;
   CachedResults<SmartPtr<Vector> > sigma_tilde_n_d_inv_cache_;
   CachedResults<SmartPtr<Vector> > sigma_tilde_p_d_inv_cache_;
   CachedResults<SmartPtr<Vector> > d_x_plus_wr_d_cache_;
   CachedResults<SmartPtr<Vector> > rhs_cR_cache_;
   CachedResults<SmartPtr<Vector> > rhs_dR_cache_;
 
   SmartPtr<const Vector> Neg_Omega_c_plus_D_c(
      const SmartPtr<const Vector>& sigma_tilde_n_c_inv,
      const SmartPtr<const Vector>& sigma_tilde_p_c_inv,
      const Vector*                 D_c,
      const Vector&                 any_vec_in_c
   );
 
   SmartPtr<const Vector> Neg_Omega_d_plus_D_d(
      const Matrix&                 Pd_L,
      const SmartPtr<const Vector>& sigma_tilde_n_d_inv,
      const Matrix&                 neg_Pd_U,
      const SmartPtr<const Vector>& sigma_tilde_p_d_inv,
      const Vector*                 D_d,
      const Vector&                 any_vec_in_d
   );
 
   SmartPtr<const Vector> Sigma_tilde_n_c_inv(
      const SmartPtr<const Vector>& sigma_tilde_n_c,
      Number                        delta_x,
      const Vector&                 any_vec_in_n_c
   );
 
   SmartPtr<const Vector> Sigma_tilde_p_c_inv(
      const SmartPtr<const Vector>& sigma_tilde_p_c,
      Number                        delta_x,
      const Vector&                 any_vec_in_p_c
   );
 
   SmartPtr<const Vector> Sigma_tilde_n_d_inv(
      const SmartPtr<const Vector>& sigma_tilde_n_d,
      Number                        delta_x,
      const Vector&                 any_vec_in_n_d
   );
 
   SmartPtr<const Vector> Sigma_tilde_p_d_inv(
      const SmartPtr<const Vector>& sigma_tilde_p_d,
      Number                        delta_x,
      const Vector&                 any_vec_in_p_d
   );
 
   SmartPtr<const Vector> D_x_plus_wr_d(
      const SmartPtr<const Vector>& CD_x0,
      Number                        factor,
      const Vector&                 wr_d
   );
 
   SmartPtr<const Vector> Rhs_cR(
      const Vector&                 rhs_c,
      const SmartPtr<const Vector>& sigma_tilde_n_c_inv,
      const Vector&                 rhs_n_c,
      const SmartPtr<const Vector>& sigma_tilde_p_c_inv,
      const Vector&                 rhs_p_c
   );
 
   SmartPtr<const Vector> Rhs_dR(
      const Vector&                 rhs_d,
      const SmartPtr<const Vector>& sigma_tilde_n_d_inv,
      const Vector&                 rhs_n_d,
      const Matrix&                 pd_L,
      const SmartPtr<const Vector>& sigma_tilde_p_d_inv,
      const Vector&                 rhs_p_d,
      const Matrix&                 pd_U
   );
 
   SmartPtr<AugSystemSolver> orig_aug_solver_;
   bool skip_orig_aug_solver_init_;
};
 
} // namespace Ipopt
 
#endif