/home/coin/SVN-release/CoinAll-1.1.0/Ipopt/src/contrib/CGPenalty/IpCGPerturbationHandler.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 Common Public License.
//
// $Id: IpPDPerturbationHandler.hpp 778 2006-08-16 18:11:13Z andreasw $
//
// Authors:  Carl Laird, Andreas Waechter              IBM    2005-08-04

#ifndef __IPCGPERTURBATIONHANDLER_HPP__
#define __IPCGPERTURBATIONHANDLER_HPP__

#include "IpPDPerturbationHandler.hpp"

namespace Ipopt
{

  class CGPerturbationHandler: public PDPerturbationHandler
  {
  public:
    CGPerturbationHandler();
    virtual ~CGPerturbationHandler()
    {}

    /* 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:
    CGPerturbationHandler(const CGPerturbationHandler&);

    void operator=(const CGPerturbationHandler&);

    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_;

    Number penalty_max_;
    Number mult_diverg_feasibility_tol_;

    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

---