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 Common Public License.
//
// $Id: IpAugRestoSystemSolver.hpp 759 2006-07-07 03:07:08Z andreasw $
//
// 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,
      double W_factor,
      const Vector* D_x,
      double delta_x,
      const Vector* D_s,
      double delta_s,
      const Matrix* J_c,
      const Vector* D_c,
      double delta_c,
      const Matrix* J_d,
      const Vector* D_d,
      double 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