IpPDPerturbationHandler.hpp

Go to the documentation of this file.

// Copyright (C) 2005, 2006 International Business Machines and others.
// All Rights Reserved.
// This code is published under the Common Public License.
//
// $Id: IpPDPerturbationHandler.hpp 759 2006-07-07 03:07:08Z 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);

    bool ConsiderNewSystem(Number& delta_x, Number& delta_s,
                           Number& delta_c, Number& delta_d);

    bool PerturbForSingularity(Number& delta_x, Number& delta_s,
                               Number& delta_c, Number& delta_d);

    bool PerturbForWrongInertia(Number& delta_x, Number& delta_s,
                                Number& delta_c, Number& delta_d);

    void CurrentPerturbation(Number& delta_x, Number& delta_s,
                             Number& delta_c, Number& delta_d);

    static void RegisterOptions(SmartPtr<RegisteredOptions> roptions);

  private:
    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 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