IpPDPerturbationHandler.hpp

Go to the documentation of this file.

// Copyright (C) 2005, 2007 International Business Machines and others.
// All Rights Reserved.
// This code is published under the Eclipse Public License.
//
// $Id: IpPDPerturbationHandler.hpp 1861 2010-12-21 21:34:47Z andreasw $
//
// Authors:  Carl Laird, Andreas Waechter              IBM    2005-08-04
 
#ifndef __IPPDPERTURBATIONHANDLER_HPP__
#define __IPPDPERTURBATIONHANDLER_HPP__
 
#include "IpAlgStrategy.hpp"
 
namespace Ipopt
{
 
  class PDPerturbationHandler: public AlgorithmStrategyObject
  {
  public:
    PDPerturbationHandler();
    virtual ~PDPerturbationHandler()
    {}
 
    /* overloaded from AlgorithmStrategyObject */
    virtual bool InitializeImpl(const OptionsList& options,
                                const std::string& prefix);
 
    virtual bool ConsiderNewSystem(Number& delta_x, Number& delta_s,
                                   Number& delta_c, Number& delta_d);
 
    virtual bool PerturbForSingularity(Number& delta_x, Number& delta_s,
                                       Number& delta_c, Number& delta_d);
 
    virtual bool PerturbForWrongInertia(Number& delta_x, Number& delta_s,
                                        Number& delta_c, Number& delta_d);
 
    virtual void CurrentPerturbation(Number& delta_x, Number& delta_s,
                                     Number& delta_c, Number& delta_d);
 
    static void RegisterOptions(SmartPtr<RegisteredOptions> roptions);
 
  protected:
    PDPerturbationHandler(const PDPerturbationHandler&);
 
    void operator=(const PDPerturbationHandler&);
 
    Number delta_x_last_;
    Number delta_s_last_;
    Number delta_c_last_;
    Number delta_d_last_;
 
    Number delta_x_curr_;
    Number delta_s_curr_;
    Number delta_c_curr_;
    Number delta_d_curr_;
 
    bool get_deltas_for_wrong_inertia_called_;
 
    enum DegenType
    {
      NOT_YET_DETERMINED,
      NOT_DEGENERATE,
      DEGENERATE
    };
 
    DegenType hess_degenerate_;
 
    DegenType jac_degenerate_;
 
    Index degen_iters_;
 
    enum TrialStatus
    {
      NO_TEST,
      TEST_DELTA_C_EQ_0_DELTA_X_EQ_0,
      TEST_DELTA_C_GT_0_DELTA_X_EQ_0,
      TEST_DELTA_C_EQ_0_DELTA_X_GT_0,
      TEST_DELTA_C_GT_0_DELTA_X_GT_0
    };
 
    TrialStatus test_status_;
 
    Number delta_xs_max_;
    Number delta_xs_min_;
    Number delta_xs_first_inc_fact_;
    Number delta_xs_inc_fact_;
    Number delta_xs_dec_fact_;
    Number delta_xs_init_;
    Number delta_cd_val_;
    Number delta_cd_exp_;
    bool reset_last_;
    Index degen_iters_max_;
    bool perturb_always_cd_;
 
    bool get_deltas_for_wrong_inertia(Number& delta_x, Number& delta_s,
                                      Number& delta_c, Number& delta_d);
 
    void finalize_test();
    Number delta_cd();
 
  };
 
} // namespace Ipopt
 
#endif