IpGenAugSystemSolver.hpp

Go to the documentation of this file.

// Copyright (C) 2007 International Business Machines and others.
// All Rights Reserved.
// This code is published under the Eclipse Public License.
//
// $Id: IpGenAugSystemSolver.hpp 2269 2013-05-05 11:32:40Z stefan $
//
// Authors:  Andreas Waechter     IBM    2007-03-01

#ifndef __IP_STDAUGSYSTEMSOLVER_HPP__
#define __IP_STDAUGSYSTEMSOLVER_HPP__

#include "IpAugSystemSolver.hpp"
#include "IpGenKKTSolverInterface.hpp"

namespace Ipopt
{
  class GenAugSystemSolver : public AugSystemSolver
  {
  public:
    GenAugSystemSolver(GenKKTSolverInterface& SolverInterface);

    virtual ~GenAugSystemSolver();

    bool InitializeImpl(const OptionsList& options,
                        const std::string& prefix);

    virtual ESymSolverStatus MultiSolve(
      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,
      std::vector<SmartPtr<const Vector> >& rhs_xV,
      std::vector<SmartPtr<const Vector> >& rhs_sV,
      std::vector<SmartPtr<const Vector> >& rhs_cV,
      std::vector<SmartPtr<const Vector> >& rhs_dV,
      std::vector<SmartPtr<Vector> >& sol_xV,
      std::vector<SmartPtr<Vector> >& sol_sV,
      std::vector<SmartPtr<Vector> >& sol_cV,
      std::vector<SmartPtr<Vector> >& sol_dV,
      bool check_NegEVals,
      Index numberOfNegEVals);

    virtual Index NumberOfNegEVals() const;

    virtual bool ProvidesInertia() const;

    virtual bool IncreaseQuality();

  private:
    GenAugSystemSolver();
    GenAugSystemSolver(const GenAugSystemSolver&);

    void operator=(const GenAugSystemSolver&);

    bool AugmentedSystemChanged(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);

    void UpdateTags(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);

    SmartPtr<GenKKTSolverInterface> solver_interface_;

    TaggedObject::Tag w_tag_;
    double w_factor_;
    TaggedObject::Tag d_x_tag_;
    double delta_x_;
    TaggedObject::Tag d_s_tag_;
    double delta_s_;
    TaggedObject::Tag j_c_tag_;
    TaggedObject::Tag d_c_tag_;
    double delta_c_;
    TaggedObject::Tag j_d_tag_;
    TaggedObject::Tag d_d_tag_;
    double delta_d_;

    Number* dx_vals_copy_;
    Number* ds_vals_copy_;
    Number* dc_vals_copy_;
    Number* dd_vals_copy_;

    bool warm_start_same_structure_;
  };

} // namespace Ipopt

#endif